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