Lorenzo, Thanks for the files. Although the output from ex14 looks correct in the vtu file (with your class), Paraview and Visit do not even plot the mesh correctly. So my guess is that this has something to do with how TRI6 connectivity is set in that element. Because the behavior seems to be independent of whether the VTKIO or VTKIO_NOVTK class is used to output. I will check this and try to get it working.
Roy, I will look in to vtk_io.C and see if that is the problem or if each element's conenctivity setting routines have something fishy going on. I will send a patch once I have fixed and tested it. Cheers, Vijay On Tue, Feb 24, 2009 at 1:47 PM, Lorenzo Botti <bottilore...@gmail.com> wrote: > Try this... > As John suggested some time ago I should have recasted this writer to make > it work if VTK is not linked to libmesh at compile time. > > Lorenzo > > 2009/2/24 Vijay S. Mahadevan <vija...@gmail.com> >> >> Hi, >> >> Do see my question below. >> >> Vijay >> >> PS: If you got the previous email already, then sorry for re-posting. >> I had not heard anything from anybody yet and so was not sure if this >> passed through the mailing list. >> >> ---------- Forwarded message ---------- >> From: Vijay S. Mahadevan <vija...@gmail.com> >> Date: Wed, Feb 18, 2009 at 7:10 PM >> Subject: VTK: Outputting higher order solution >> To: "libmesh-users@lists.sourceforge.net" >> <libmesh-users@lists.sourceforge.net> >> >> >> Hi, >> >> I have an Equation System that outputs the solution in VTK format with >> second order triangles (TRI6). For some reason, I think the >> connectivity is all messed up. Is this a known bug or do you think I >> am doing something weird here ? >> >> For example, the connectivity in the output vtu file : >> >> <DataArray type="Int32" Name="connectivity" format="ascii"> >> 0 3 5 0 -141730686 -1084317832 >> 2 4 8 0 -141730686 -1084317832 >> 6 10 8 0 -141730686 -1084317832 >> >> IMO, this looks like a corrupt array in memory and the output is >> garbage while accessing the connectivity array. Let me know if any of >> you have seen something like this. May be I am not setting something >> before outputting the data ?? Also if you guys have no problems >> outputting vtk with the current libmesh-dev version for higher order >> meshes, then do let me know. >> >> Thanks in advance, >> Vijay >> >> >> ------------------------------------------------------------------------------ >> Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, >> CA >> -OSBC tackles the biggest issue in open source: Open Sourcing the >> Enterprise >> -Strategies to boost innovation and cut costs with open source >> participation >> -Receive a $600 discount off the registration fee with the source code: >> SFAD >> http://p.sf.net/sfu/XcvMzF8H >> _______________________________________________ >> Libmesh-users mailing list >> Libmesh-users@lists.sourceforge.net >> https://lists.sourceforge.net/lists/listinfo/libmesh-users > > ------------------------------------------------------------------------------ Open Source Business Conference (OSBC), March 24-25, 2009, San Francisco, CA -OSBC tackles the biggest issue in open source: Open Sourcing the Enterprise -Strategies to boost innovation and cut costs with open source participation -Receive a $600 discount off the registration fee with the source code: SFAD http://p.sf.net/sfu/XcvMzF8H _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
