[cig-commits] commit:

[Mercurial]

changeset:   18:ee4b415d47d7
user:        LukeHodkinson
date:        Fri Jan 11 00:38:28 2008 +0000
files:       Mesh/src/Mesh_RegularAlgorithms.c
description:
The recent fix I made to the particle ownership tests
caused the use of quadratic element types to break.
I've modified the ownership tests to be based on
mesh elements instead of mesh vertices.


diff -r bc405ee8e151 -r ee4b415d47d7 Mesh/src/Mesh_RegularAlgorithms.c
--- a/Mesh/src/Mesh_RegularAlgorithms.c	Wed Jan 09 01:05:55 2008 +0000
+++ b/Mesh/src/Mesh_RegularAlgorithms.c	Fri Jan 11 00:38:28 2008 +0000
@@ -147,7 +147,7 @@ void Mesh_RegularAlgorithms_Update( void
 void Mesh_RegularAlgorithms_Update( void* algorithms ) {
 	Mesh_RegularAlgorithms*	self = (Mesh_RegularAlgorithms*)algorithms;
 	unsigned		nDims;
-	Grid*			vGrid;
+	Grid*			eGrid;
 	int			ii;
 
 	assert( self && Stg_CheckType( self, Mesh_RegularAlgorithms ) );
@@ -162,9 +162,9 @@ void Mesh_RegularAlgorithms_Update( void
 	Mesh_GetGlobalCoordRange( self->mesh, self->minCrd, self->maxCrd );
 
 	self->sep = AllocArray( double, nDims );
-	vGrid = *Mesh_GetExtension( self->mesh, Grid**, "vertexGrid" );
+	eGrid = *Mesh_GetExtension( self->mesh, Grid**, "elementGrid" );
 	for( ii = 0; ii < nDims; ii++ )
-		self->sep[ii] = (self->maxCrd[ii] - self->minCrd[ii]) / (vGrid->sizes[ii] - 1);
+		self->sep[ii] = (self->maxCrd[ii] - self->minCrd[ii]) / eGrid->sizes[ii];
 }
 
 Bool _Mesh_RegularAlgorithms_Search( void* algorithms, void* _mesh, double* point,