[cig-commits] commit: Rudimentary XDMF support for 3D discontinuous elements

Mercurial hg at geodynamics.org
Tue Oct 4 01:58:52 PDT 2011


changeset:   877:365bc6244be8
user:        Walter Landry <wlandry at caltech.edu>
date:        Mon Oct 03 20:41:47 2011 -0700
files:       Utils/src/XDMFGenerator.cxx
description:
Rudimentary XDMF support for 3D discontinuous elements


diff -r 63d1e39d385a -r 365bc6244be8 Utils/src/XDMFGenerator.cxx
--- a/Utils/src/XDMFGenerator.cxx	Mon Oct 03 20:29:44 2011 -0700
+++ b/Utils/src/XDMFGenerator.cxx	Mon Oct 03 20:41:47 2011 -0700
@@ -193,7 +193,7 @@ void _XDMFGenerator_WriteFieldSchema( Un
         }
       else if(element_family=="linear-inner")
         {
-          topologyType= (nDims==2) ? "Triangle" : "Hexahedron";
+          topologyType= (nDims==2) ? "Triangle" : "Tetrahedron";
         }
       /** first create the grid which is applicable to the
           checkpointed fevariables **/
@@ -402,7 +402,9 @@ void _XDMFGenerator_WriteFieldSchema( Un
             }
           else
             {
-              abort();
+              Journal_Printf( stream, "    <Topology Type=\"%s\" NumberOfElements=\"%u\"> \n", topologyType.c_str(), elementGlobalSize );
+              Journal_Printf( stream, "      <DataItem Format=\"HDF\" DataType=\"Int\"  Dimensions=\"%u %u\">Mesh.%s.%05d.h5:/connectivity</DataItem>\n", elementGlobalSize, maxNodes, feMesh->name, context->timeStep );
+              Journal_Printf( stream, "    </Topology>\n\n" );
             }
         }
       /* Continuous linear elements */



More information about the CIG-COMMITS mailing list