On 22/06/2011 18:19, John Peterson wrote:
> 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
>
Wonderful!
When compiling I get a ton of warnings like
In file included from src/mesh/vtk_io.C:73:
/usr/include/vtk-5.4/vtkCellArray.h:43:46: warning: invoking macro
vtkExportedTypeRevisionMacro argument 3: empty macro arguments are
undefined in ISO C90 and ISO C++98
but VTK output is definitely ok now (tested for TRI3, QUAD4, TRI6, QUAD8
and QUAD9).
Thanks a lot!
Matteo
--
Matteo Semplice Dip. di Fisica e Matematica
Phone: 031-2386132 Università dell'Insubria
Fax: 031-2386209 Via Valleggio, 11
22100 Como
------------------------------------------------------------------------------
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
