On Wed, Jun 22, 2011 at 9:16 AM, Matteo Semplice <[email protected]> wrote: > > However, while I was testing, I found another problem. > The VTK output is correct only for first order elements (QUAD4 and TRI3). > In fact, if I use 2 TRI6 elements on a unit square domain, the output of > VTKIO::write_equation_systems contains 9 nodes (which includes the midpoints > of the edges), but the "connectivity" says > > <DataArray type="Int32" Name="connectivity" format="ascii" > RangeMin="-2147483648" RangeMax="8"> > 0 3 5 -2147483648 0 -2147483648 > 0 5 8 -2147483648 0 -2147483648 > </DataArray>
Definitely looks like a bug with quadratic elements in VTK...for one thing the Tri6::connectivity() function is designed to write out 4 linear child Tri3's, and it does not appear that VTK is expecting any kind of sub-triangle output, it only asks for the connectivity of child 0. > and this makes paraview crash. Since you mention paraview: if you are only interested in visualizing something in paraview, you might try writing the Exodus format instead. It's more heavily-used and better-tested than the VTK stuff... -- John ------------------------------------------------------------------------------ Simplify data backup and recovery for your virtual environment with vRanger. Installation's a snap, and flexible recovery options mean your data is safe, secure and there when you need it. Data protection magic? Nope - It's vRanger. Get your free trial download today. http://p.sf.net/sfu/quest-sfdev2dev _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
