Index: /issm/trunk/src/c/PenaltyConstraintsx/PenaltyConstraintsx.cpp
===================================================================
--- /issm/trunk/src/c/PenaltyConstraintsx/PenaltyConstraintsx.cpp	(revision 3888)
+++ /issm/trunk/src/c/PenaltyConstraintsx/PenaltyConstraintsx.cpp	(revision 3889)
@@ -5,5 +5,4 @@
 #include "./PenaltyConstraintsx.h"
 #include "./RiftConstraints.h"
-
 #include "../shared/shared.h"
 #include "../include/include.h"
@@ -28,4 +27,5 @@
 
 	/*First, get nodes and loads configured: */
+	elements->Configure(elements, loads, nodes,vertices, materials,parameters);
 	nodes->Configure(elements, loads, nodes,vertices, materials,parameters);
 	loads->Configure(elements, loads, nodes,vertices, materials,parameters);
@@ -44,5 +44,4 @@
 		converged=1;
 	}
-
 		
 	/*Assign output pointers: */
Index: /issm/trunk/src/c/parallel/diagnostic_core.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3888)
+++ /issm/trunk/src/c/parallel/diagnostic_core.cpp	(revision 3889)
@@ -140,5 +140,4 @@
 	}
 	
-	
 	if (dim==3){
 
Index: /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp
===================================================================
--- /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 3888)
+++ /issm/trunk/src/c/parallel/diagnostic_core_linear.cpp	(revision 3889)
@@ -36,11 +36,10 @@
 
 	if (verbose) _printf_("   Generating penalty matrices\n");
-		//*Generate penalty system matrices
-		PenaltySystemMatricesx(Kgg, pg,NULL,fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,kflag,pflag,analysis_type,sub_analysis_type); 
+	//*Generate penalty system matrices
+	PenaltySystemMatricesx(Kgg, pg,NULL,fem->elements,fem->nodes,fem->vertices,fem->loads,fem->materials,fem->parameters,kflag,pflag,analysis_type,sub_analysis_type); 
 
 	/*!Reduce matrix from g to f size:*/
 	if (verbose) _printf_("   reducing matrix from g to f set\n");
 	Reducematrixfromgtofx(&Kff,&Kfs,Kgg,fem->Gmn,fem->nodesets); MatFree(&Kgg);
-
 
 	/*!Reduce load from g to f size: */