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

[walter at geodynamics.org]

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

Modified:
   long/3D/Gale/trunk/src/StgFEM/
   long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c
Log:
 r987 at earth (orig r738):  LukeHodkinson | 2007-02-11 22:39:29 -0800
 Adding a debugging line.
 



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:731
db209038-57f2-0310-97fa-b160e0ae9d04:/trunk:669
   + 38867592-cf10-0410-9e16-a142ea72ac34:/cig:880
db209038-57f2-0310-97fa-b160e0ae9d04:/branches/decomp3d:738
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:19:55 UTC (rev 6103)
+++ long/3D/Gale/trunk/src/StgFEM/SLE/LinearAlgebra/src/SROpGenerator.c	2007-02-23 20:20:45 UTC (rev 6104)
@@ -560,7 +560,11 @@
 			if( fEqNum == (unsigned)-1 )
 				continue;
 
-			insist( Mesh_Search( cMesh, Mesh_GetVertex( fMesh, n_i ), &dim, &ind ) );
+			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 );
+			}
 			if( dim == MT_VERTEX ) {
 				cTopNode = self->topMaps[level - 1][ind];
 				for( dof_j = 0; dof_j < dofLayout->dofCounts[cTopNode]; dof_j++ ) {