On Wed, Jun 22, 2011 at 9:36 AM, John Peterson <[email protected]> wrote: > 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.
I just checked in a trivial fix for this, TRI6 VTK output now seems to work for me. -- 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
