Index: /issm/trunk-jpl/test/NightlyRun/test3001.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3001.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3001.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -16,4 +17,5 @@
 md.stressbalance.requested_outputs=['default','DeviatoricStressxx','DeviatoricStressyy','DeviatoricStressxy']
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3002.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3002.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3002.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -16,4 +17,5 @@
 md.cluster=generic('name',gethostname(),'np',1)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3003.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3003.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3003.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -17,4 +18,5 @@
 md.stressbalance.requested_outputs=['default','StressTensorxx','StressTensoryy','StressTensorzz','StressTensorxy','StressTensorxz','StressTensoryz']
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3004.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3004.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3004.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -16,4 +17,5 @@
 md.cluster=generic('name',gethostname(),'np',1)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3005.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3005.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3005.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -15,4 +16,5 @@
 md.cluster=generic('name',gethostname(),'np',1)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3006.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3006.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3006.py	(revision 23746)
@@ -8,4 +8,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -19,4 +20,5 @@
 md.masstransport.spcthickness=md.geometry.thickness
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3007.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3007.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3007.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -16,4 +17,5 @@
 md.cluster=generic('name',gethostname(),'np',1)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3008.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3008.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3008.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -17,4 +18,5 @@
 md.cluster=generic('name',gethostname(),'np',1)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Thermal')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3009.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3009.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3009.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -21,4 +22,5 @@
 md.transient.isgroundingline=False
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Transient')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3010.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3010.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3010.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -17,4 +18,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 md=solve(md,'Transient')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3101.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3101.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3101.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -17,4 +18,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3102.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3102.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3102.py	(revision 23746)
@@ -7,5 +7,5 @@
 from setflowequation import *
 from solve import *
-
+from issmmumpssolver import issmmumpssolver
 
 md=triangle(model(),'../Exp/Square.exp',180000.)
@@ -16,4 +16,5 @@
 md.cluster=generic('name',gethostname(),'np',3)
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 
 md=solve(md,'Stressbalance')
Index: /issm/trunk-jpl/test/NightlyRun/test3103.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3103.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3103.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -18,4 +19,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3104.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3104.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3104.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -17,4 +18,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Stressbalance')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3105.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3105.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3105.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -16,4 +17,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3106.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3106.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3106.py	(revision 23746)
@@ -8,4 +8,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -20,4 +21,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3107.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3107.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3107.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -17,4 +18,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Masstransport')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3108.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3108.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3108.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -18,4 +19,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Thermal')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3109.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3109.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3109.py	(revision 23746)
@@ -7,4 +7,5 @@
 from setflowequation import *
 from solve import *
+from issmmumpssolver import issmmumpssolver
 
 
@@ -22,4 +23,5 @@
 
 md.autodiff.isautodiff=True
+md.toolkits.DefaultAnalysis=issmmumpssolver()
 md=solve(md,'Transient')
 
Index: /issm/trunk-jpl/test/NightlyRun/test3119.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test3119.py	(revision 23745)
+++ /issm/trunk-jpl/test/NightlyRun/test3119.py	(revision 23746)
@@ -9,4 +9,5 @@
 from dependent import *
 from solve import *
+from issmgslsolver import issmgslsolver
 
 
@@ -19,7 +20,7 @@
 
 md.autodiff.isautodiff=True
-md.verbose.autodiff=True
+md.toolkits.DefaultAnalysis=issmgslsolver()
 
-#first run scalar reverse mode: 
+#first run scalar reverse mode:
 md.autodiff.independents=[independent('name','md.geometry.thickness','type','vertex','nods',md.mesh.numberofvertices)]
 md.autodiff.dependents=[dependent('name','MaxVel','type','scalar','fos_reverse_index',1)]
@@ -28,5 +29,5 @@
 md=solve(md,'Transient')
 
-#recover jacobian: 
+#recover jacobian:
 jac_reverse=md.results.TransientSolution[0].AutodiffJacobian