Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-08-04 Thread Wolfgang Bangerth
On 07/26/2017 01:18 PM, sml.imf...@gmail.com wrote: I must admit that I don't really remember ... I think it was GridOut, since I specifically wanted to export only the grid and not some solution data. What dou you mean with "GridOut should work"? I was wrong thinking that the GridOut::write_v

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-07-26 Thread sml . imfeld
I must admit that I don't really remember ... I think it was GridOut, since I specifically wanted to export only the grid and not some solution data. What dou you mean with "GridOut should work"? On Wednesday, July 26, 2017 at 9:04:16 PM UTC+2, Wolfgang Bangerth wrote: > > On 07/26/2017 12:55 PM

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-07-26 Thread Wolfgang Bangerth
On 07/26/2017 12:55 PM, sml.imf...@gmail.com wrote: Exactly, I generated it using GridOut. I see the problem now. I replaced all duplicate vertices in the vtk file, now it works as it should. Thanks for the clarification! Great. But just to clarify, you used GridOut or DataOut? The former sho

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-07-26 Thread sml . imfeld
Exactly, I generated it using GridOut. I see the problem now. I replaced all duplicate vertices in the vtk file, now it works as it should. Thanks for the clarification! Best, Samuel On Wednesday, July 26, 2017 at 8:38:39 PM UTC+2, Wolfgang Bangerth wrote: > > On 07/26/2017 09:18 AM, sml.i...@g

Re: [deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-07-26 Thread Wolfgang Bangerth
On 07/26/2017 09:18 AM, sml.imf...@gmail.com wrote: I am trying to read in a triangulation from a vtk file (generated by the dealii library). It works, but it seems like no information about neighboring cells is generated: iterating through the cells and calling cell->neighbor_index(face) alw

[deal.II] Triangulation neighbor information when using GridIn::read_vtk

2017-07-26 Thread sml . imfeld
Hello all I am trying to read in a triangulation from a vtk file (generated by the dealii library). It works, but it seems like no information about neighboring cells is generated: iterating through the cells and calling cell->neighbor_index(face) always gives -1, indicating that the face is at