Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18148)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18149)
@@ -7,5 +7,5 @@
 #include "../cores/cores.h"
 
-//#define FSANALYTICAL 21
+#define FSANALYTICAL 4
 
 /*Model processing*/
@@ -3064,4 +3064,5 @@
 
 	if(element->IsFloating() || !element->IsOnBase()) return NULL;
+	return NULL;
 
 	/*If on water or not FS, skip stiffness: */
@@ -3074,4 +3075,5 @@
 	IssmDouble  alpha2,Jdet;
 	IssmDouble  x_coord,y_coord,z_coord;
+	IssmDouble *xyz_list      = NULL;
 	IssmDouble *xyz_list_base = NULL;
 	Gauss*      gauss         = NULL;
@@ -3091,4 +3093,5 @@
 
 	/*Retrieve all inputs and parameters*/
+	element->GetVerticesCoordinates(&xyz_list);
 	element->GetVerticesCoordinatesBase(&xyz_list_base);
 	Input* vx_input         = element->GetInput(VxEnum);      _assert_(vx_input);
@@ -3123,4 +3126,5 @@
 	/*Clean up and return*/
 	delete gauss;
+	xDelete<IssmDouble>(xyz_list);
 	xDelete<IssmDouble>(xyz_list_base);
 	xDelete<IssmDouble>(B);
Index: /issm/trunk-jpl/src/c/shared/Matrix/MatrixUtils.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Matrix/MatrixUtils.cpp	(revision 18148)
+++ /issm/trunk-jpl/src/c/shared/Matrix/MatrixUtils.cpp	(revision 18149)
@@ -3,5 +3,4 @@
  */ 
 
-#include "../Elements/elements.h"
 /*Headers*/
 /*{{{*/
@@ -461,8 +460,5 @@
 	/*Compute determinant*/
 	Matrix4x4Determinant(&det,A);
-	if(fabs(det) < DBL_EPSILON){
-		printarray(A,4,4);
-		_error_("Determinant smaller than machine epsilon");
-	}
+	if(fabs(det) < DBL_EPSILON) _error_("Determinant smaller than machine epsilon");
 
 	/*Compute adjoint matrix*/