[cig-commits] commit: Make _FiniteElementContext_DumpMeshHDF5 output quad and discontinuous elements in a non-jumbled order

Mercurial hg at geodynamics.org
Sun Oct 2 11:02:50 PDT 2011


changeset:   801:df953fda3eb6
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Sun Oct 02 11:01:01 2011 -0700
files:       SLE/SystemSetup/src/FiniteElementContext.cxx
description:
Make _FiniteElementContext_DumpMeshHDF5 output quad and discontinuous elements in a non-jumbled order


diff -r e24206efa24e -r df953fda3eb6 SLE/SystemSetup/src/FiniteElementContext.cxx
--- a/SLE/SystemSetup/src/FiniteElementContext.cxx	Sun Oct 02 10:58:42 2011 -0700
+++ b/SLE/SystemSetup/src/FiniteElementContext.cxx	Sun Oct 02 11:01:01 2011 -0700
@@ -798,13 +798,13 @@ void _FiniteElementContext_DumpMeshHDF5(
       
       for( node_I = 0 ; node_I < nodesPerEl ; node_I++ )
          buf_int[node_I] = Mesh_DomainToGlobal( mesh, MT_VERTEX, nodeList[node_I] );
-      if(nodesPerEl>3)
+      if((nDims == 2 && nodesPerEl==4) || (nDims == 3 && nodesPerEl==8))
         {
           /* some reordering is required to account for standard node ordering */ 
           buf_int[3] = Mesh_DomainToGlobal( mesh, MT_VERTEX, nodeList[2] );
           buf_int[2] = Mesh_DomainToGlobal( mesh, MT_VERTEX, nodeList[3] );
         }
-      if( nDims == 3 && nodesPerEl>7 ) {
+      if(nDims == 3 && nodesPerEl==8) {
          buf_int[7] = Mesh_DomainToGlobal( mesh, MT_VERTEX, nodeList[6] );
          buf_int[6] = Mesh_DomainToGlobal( mesh, MT_VERTEX, nodeList[7] );
       }



More information about the CIG-COMMITS mailing list