Hey John, thanks for your reply, here is the diff file.
Best regards Ping Am 07.09.2010 18:00, schrieb John Peterson:
Hi Ping, I'd be happy to take a look at your patch. Would you mind sending me the diffs in the 'svn diff' format? If you have an SVN checked-out copy of the source tree, just run 'svn diff', if not just send me your vtk_io.C file directly and I'll take a look. Thanks, John On Tue, Sep 7, 2010 at 10:53 AM, Ping Rong<[email protected]> wrote:Hello all, I have been recently trying to write my meshes into *.vtu files. In most cases it works pretty well, except the HEX27 elements. So I made a few modifications in the source code, and made it work. I guess maybe I just post it here, maybe someone can update it to the subversion, if nothing is wrong, of course. ------------------------------------ file: mesh\vtk_io.C line 215 ------------------------------------ ... case HEX27: celltype = VTK_TRIQUADRATIC_HEXAHEDRON; //Added line break; // 29 Added line .... ------------------------------------ file: mesh\cell_hex27.C line 460 -560 ------------------------------------ // according to the VTK documentation http://www.vtk.org/doc/nightly/html/classvtkTriQuadraticHexahedron.html // I change the connectivtiy mapping. Substitute the code from line 460 to 560 with the following: ... case 0: conn.resize(27); conn[0] = this->node(0); conn[1] = this->node(1); conn[2] = this->node(2); conn[3] = this->node(3); conn[4] = this->node(4); conn[5] = this->node(5); conn[6] = this->node(6); conn[7] = this->node(7); conn[8] = this->node(8); conn[9] = this->node(9); conn[10] = this->node(10); conn[11] = this->node(11); // conn[12] = this->node(16); conn[13] = this->node(17); conn[14] = this->node(18); conn[15] = this->node(19); conn[16] = this->node(12); conn[17] = this->node(13); // conn[18] = this->node(14); conn[19] = this->node(15); conn[20] = this->node(24); conn[21] = this->node(22); conn[22] = this->node(21); conn[23] = this->node(23); conn[24] = this->node(20); conn[25] = this->node(25); conn[26] = this->node(26); return; ... -- Mit freundlichen Grüßen / With best regards Ping Rong
-- Mit freundlichen Grüßen / With best regards Ping Rong
Index: src/geom/cell_hex27.C
===================================================================
--- src/geom/cell_hex27.C (revision 3940)
+++ src/geom/cell_hex27.C (working copy)
@@ -458,7 +458,35 @@
switch (sc)
{
case 0:
-
+ conn.resize(27);
+ conn[0] = this->node(0);
+ conn[1] = this->node(1);
+ conn[2] = this->node(2);
+ conn[3] = this->node(3);
+ conn[4] = this->node(4);
+ conn[5] = this->node(5);
+ conn[6] = this->node(6);
+ conn[7] = this->node(7);
+ conn[8] = this->node(8);
+ conn[9] = this->node(9);
+ conn[10] = this->node(10);
+ conn[11] = this->node(11); //
+ conn[12] = this->node(16);
+ conn[13] = this->node(17);
+ conn[14] = this->node(18);
+ conn[15] = this->node(19);
+ conn[16] = this->node(12);
+ conn[17] = this->node(13); //
+ conn[18] = this->node(14);
+ conn[19] = this->node(15);
+ conn[20] = this->node(24);
+ conn[21] = this->node(22);
+ conn[22] = this->node(21);
+ conn[23] = this->node(23);
+ conn[24] = this->node(20);
+ conn[25] = this->node(25);
+ conn[26] = this->node(26);
+ /*
conn[0] = this->node(0);
conn[1] = this->node(8);
conn[2] = this->node(20);
@@ -558,7 +586,7 @@
conn[5] = this->node(17);
conn[6] = this->node(6);
conn[7] = this->node(18);
-
+ */
return;
default:
Index: src/mesh/vtk_io.C
===================================================================
--- src/mesh/vtk_io.C (revision 3940)
+++ src/mesh/vtk_io.C (working copy)
@@ -213,6 +213,8 @@
#endif
case EDGE4:
case HEX27:
+ celltype = VTK_TRIQUADRATIC_HEXAHEDRON;
//Added line
+ break; // 29
case INFEDGE2:
case INFQUAD4:
case INFQUAD6:
<<attachment: ping_rong.vcf>>
------------------------------------------------------------------------------ This SF.net Dev2Dev email is sponsored by: Show off your parallel programming skills. Enter the Intel(R) Threading Challenge 2010. http://p.sf.net/sfu/intel-thread-sfd
_______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
