[cig-commits] commit: Add a triInds and tetInds mapping for Q2 elements. Seems to work completely for Q2P-1 elements.

[Mercurial]

changeset:   797:03c1c4e942ca
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Sep 21 23:30:19 2011 -0700
files:       Discretisation/src/C2Generator.cxx
description:
Add a triInds and tetInds mapping for Q2 elements.  Seems to work completely for Q2P-1 elements.


diff -r 52aeca915d0d -r 03c1c4e942ca Discretisation/src/C2Generator.cxx
--- a/Discretisation/src/C2Generator.cxx	Sun Sep 18 12:38:22 2011 -0700
+++ b/Discretisation/src/C2Generator.cxx	Wed Sep 21 23:30:19 2011 -0700
@@ -466,7 +466,6 @@ void C2Generator_GenElementTypes( void* 
 	C2Generator*	self = (C2Generator*)meshGenerator;
 	Stream*		stream;
 	unsigned	nDomainEls;
-	unsigned	vertMap[8] = {0, 2, 6, 8, 18, 20, 24, 26};
 	unsigned	e_i;
 
 	assert( self );
@@ -479,7 +478,7 @@ void C2Generator_GenElementTypes( void* 
 	mesh->elTypes = AllocArray( Mesh_ElementType*, mesh->nElTypes );
 	mesh->elTypes[0] = (Mesh_ElementType*)Mesh_HexType_New();
 	Mesh_ElementType_SetMesh( mesh->elTypes[0], mesh );
-	Mesh_HexType_SetVertexMap( mesh->elTypes[0], vertMap );
+	Mesh_HexType_SetQ2Inds( mesh->elTypes[0] );
 	nDomainEls = Mesh_GetDomainSize( mesh, Mesh_GetDimSize( mesh ) );
 	mesh->elTypeMap = AllocArray( unsigned, nDomainEls );
 	for( e_i = 0; e_i < nDomainEls; e_i++ )