Index: /issm/trunk-jpl/src/jl/core/toolkits.jl
===================================================================
--- /issm/trunk-jpl/src/jl/core/toolkits.jl	(revision 26728)
+++ /issm/trunk-jpl/src/jl/core/toolkits.jl	(revision 26729)
@@ -127,6 +127,10 @@
 function Solverx(A::IssmMatrix,b::IssmVector,xold::IssmVector) #{{{
 
-	x = IssmVector(GetSize(xold))
-	x.vector = (A.matrix\b.vector)
+	#Initialize output
+	#x = IssmVector(GetSize(xold))
+	x = IssmVector(0)
+
+	#Solve linear system
+	x.vector = A.matrix\b.vector
 
 	return x
@@ -134,3 +138,2 @@
 
 end#}}}
-