[cig-commits] r6105 - in long/3D/Gale/trunk/src/StgFEM: . SLE/LinearAlgebra/src

[walter at geodynamics.org]

Author: walter
Date: 2007-02-23 12:20:47 -0800 (Fri, 23 Feb 2007)
New Revision: 6105

Modified:
   long/3D/Gale/trunk/src/StgFEM/
   long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c
Log:
 r988 at earth (orig r739):  LukeHodkinson | 2007-02-12 14:19:03 -0800
 Modified the search algorithm so that it uses the
 faster optimised version for regular meshes.
 



Property changes on: long/3D/Gale/trunk/src/StgFEM
___________________________________________________________________
Name: svk:merge
   - 38867592-cf10-0410-9e16-a142ea72ac34:/cig:880
db209038-57f2-0310-97fa-b160e0ae9d04:/branches/decomp3d:738
db209038-57f2-0310-97fa-b160e0ae9d04:/trunk:669
   + 38867592-cf10-0410-9e16-a142ea72ac34:/cig:880
db209038-57f2-0310-97fa-b160e0ae9d04:/branches/decomp3d:739
db209038-57f2-0310-97fa-b160e0ae9d04:/trunk:669

Modified: long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c
===================================================================
--- long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c	2007-02-23 20:20:45 UTC (rev 6104)
+++ long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c	2007-02-23 20:20:47 UTC (rev 6105)
@@ -560,11 +560,8 @@
 			if( fEqNum == (unsigned)-1 )
 				continue;
 
-			if( !Mesh_Search( cMesh, Mesh_GetVertex( fMesh, n_i ), &dim, &ind ) ) {
-				printf( "*** Missed a node: %g, %g, %g\n", 
-					Mesh_GetVertex( fMesh, n_i )[0], Mesh_GetVertex( fMesh, n_i )[1], Mesh_GetVertex( fMesh, n_i )[2] );
-				assert( 0 );
-			}
+			insist( Mesh_SearchElements( cMesh, Mesh_GetVertex( fMesh, n_i ), &ind ) );
+			dim = nDims;
 			if( dim == MT_VERTEX ) {
 				cTopNode = self->topMaps[level - 1][ind];
 				for( dof_j = 0; dof_j < dofLayout->dofCounts[cTopNode]; dof_j++ ) {