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
<<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
