Index: /issm/trunk-jpl/src/m/solvers/asmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/asmoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/asmoptions.m	(revision 14610)
@@ -10,4 +10,5 @@
 
 %default asm options
+asm.toolkit='petsc';
 asm.mat_type=getfieldvalue(options,'mat_type','aij');
 asm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
Index: /issm/trunk-jpl/src/m/solvers/asmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/asmoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/asmoptions.py	(revision 14610)
@@ -12,5 +12,5 @@
 	arguments=pairoptions(*args) 
 	
-	options=[['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]];
+	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','lu'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-30]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/iluasmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/iluasmoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/iluasmoptions.m	(revision 14610)
@@ -10,4 +10,5 @@
 
 %default iluasm options
+iluasm.toolkit='petsc';
 iluasm.mat_type=getfieldvalue(options,'mat_type','aij');
 iluasm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
Index: /issm/trunk-jpl/src/m/solvers/iluasmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/iluasmoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/iluasmoptions.py	(revision 14610)
@@ -15,4 +15,5 @@
 
 	#default iluasm options
+	iluasm['toolkit']='petsc'
 	iluasm['mat_type']=options.getfieldvalue('mat_type','aij')
 	iluasm['ksp_type']=options.getfieldvalue('ksp_type','gmres')
Index: /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.m	(revision 14610)
@@ -10,4 +10,5 @@
 
 %default jacobiasm options
+jacobiasm.toolkit='petsc';
 jacobiasm.mat_type=getfieldvalue(options,'mat_type','aij');
 jacobiasm.ksp_type=getfieldvalue(options,'ksp_type','gmres');
Index: /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/jacobiasmoptions.py	(revision 14610)
@@ -12,5 +12,5 @@
 	arguments=pairoptions(*args) 
 	
-	options=[['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]];
+	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','gmres'],['pc_type','asm'],['sub_pc_type','jacobi'],['pc_asm_overlap',3],['ksp_max_it',100],['ksp_rtol',1e-15]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/jacobicgoptions.m	(revision 14610)
@@ -10,4 +10,5 @@
 
 %default jacobiasm options
+jacobicg.toolkit='petsc';
 jacobicg.mat_type=getfieldvalue(options,'mat_type','aij');
 jacobicg.ksp_type=getfieldvalue(options,'ksp_type','cg');
Index: /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/jacobicgoptions.py	(revision 14610)
@@ -12,5 +12,5 @@
 	arguments=pairoptions(*args) 
 	
-	options=[['mat_type','aij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]];
+	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['ksp_max_it',100],['ksp_rtol',1e-15]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/matlaboptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/matlaboptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/matlaboptions.m	(revision 14610)
@@ -10,3 +10,4 @@
 
 %default matlab options
+maltab.toolkit='petsc';
 maltab.ksp_type='matlab';
Index: /issm/trunk-jpl/src/m/solvers/matlaboptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/matlaboptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/matlaboptions.py	(revision 14610)
@@ -12,5 +12,5 @@
 	arguments=pairoptions(*args) 
 	
-	options=[['ksp_type','matlab']];
+	options=[['toolkit','petsc'],['ksp_type','matlab']];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/mumpsoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/mumpsoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/mumpsoptions.m	(revision 14610)
@@ -12,4 +12,5 @@
 PETSC_VERSION=petscversion();
 if PETSC_VERSION==2,
+	mumps.toolkit='petsc';
 	mumps.mat_type=getfieldvalue(options,'mat_type','aijmumps');
 	mumps.ksp_type=getfieldvalue(options,'ksp_type','preonly');
@@ -19,4 +20,5 @@
 end
 if PETSC_VERSION==3,
+	mumps.toolkit='petsc';
 	mumps.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 	mumps.ksp_type=getfieldvalue(options,'ksp_type','preonly');
Index: /issm/trunk-jpl/src/m/solvers/mumpsoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/mumpsoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/mumpsoptions.py	(revision 14610)
@@ -18,4 +18,5 @@
 	PETSC_VERSION=petscversion()
 	if PETSC_VERSION==2:
+		mumps['toolkit']='petsc'
 		mumps['mat_type']=options.getfieldvalue('mat_type','aijmumps')
 		mumps['ksp_type']=options.getfieldvalue('ksp_type','preonly')
@@ -24,4 +25,5 @@
 		mumps['pc_factor_shift_positive_definite']=options.getfieldvalue('pc_factor_shift_positive_definite','true')
 	if PETSC_VERSION==3:
+		mumps['toolkit']='petsc'
 		mumps['mat_type']=options.getfieldvalue('mat_type','mpiaij')
 		mumps['ksp_type']=options.getfieldvalue('ksp_type','preonly')
Index: /issm/trunk-jpl/src/m/solvers/soroptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/soroptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/soroptions.m	(revision 14610)
@@ -10,4 +10,5 @@
 
 %default sor options
+sor.toolkit='petsc';
 sor.mat_type=getfieldvalue(options,'mat_type','aij');
 sor.ksp_type=getfieldvalue(options,'ksp_type','cg');
Index: /issm/trunk-jpl/src/m/solvers/soroptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/soroptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/soroptions.py	(revision 14610)
@@ -12,5 +12,5 @@
 	arguments=pairoptions(*args) 
 	
-	options=[['mat_type','aij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]];
+	options=[['toolkit','petsc'],['mat_type','aij'],['ksp_type','cg'],['pc_type','sor'],['pc_sor_omega',1.1],['pc_sor_its',2]];
 
 	#now, go through our arguments, and write over default options.
Index: /issm/trunk-jpl/src/m/solvers/stokesoptions.m
===================================================================
--- /issm/trunk-jpl/src/m/solvers/stokesoptions.m	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/stokesoptions.m	(revision 14610)
@@ -15,4 +15,5 @@
 end
 if PETSC_VERSION==3,
+	stokes.toolkit='petsc';
 	stokes.mat_type=getfieldvalue(options,'mat_type','mpiaij');
 	stokes.ksp_max_it=getfieldvalue(options,'ksp_max_it',1000);
Index: /issm/trunk-jpl/src/m/solvers/stokesoptions.py
===================================================================
--- /issm/trunk-jpl/src/m/solvers/stokesoptions.py	(revision 14609)
+++ /issm/trunk-jpl/src/m/solvers/stokesoptions.py	(revision 14610)
@@ -19,5 +19,5 @@
 		raise RuntimeError('stokesoptions error message: multi-physics options not supported in Petsc 2')
 	if PETSC_VERSION==3:
-		options=[['mat_type','mpiaij'],['ksp_max_it',1000],['ksp_type','gmres'],['pc_type','fieldsplit'],['pc_field_split_type','schur'],\
+		options=[['toolkit','petsc'],['mat_type','mpiaij'],['ksp_max_it',1000],['ksp_type','gmres'],['pc_type','fieldsplit'],['pc_field_split_type','schur'],\
 	['fieldsplit_0_pc_type','hypre'],['fieldsplit_0_ksp_type','gmres'],['fieldsplit_0_pc_hypre_type','boomerang'],\
 	['fieldsplit_1_pc_type','jacobi'],['fieldsplit_1_ksp_type','preonly'],['issm_option_solver','stokes']]