Index: /issm/trunk-jpl/src/c/analyses/controltao_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16483)
+++ /issm/trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16484)
@@ -82,4 +82,6 @@
 	TaoSolve(tao);
 	TaoView(tao,PETSC_VIEWER_STDOUT_WORLD);
+
+	/*Save results*/
 	TaoGetSolutionVector(tao,&X->pvector->vector);
 	G=new Vector<IssmDouble>(0); VecFree(&G->pvector->vector);
@@ -87,7 +89,5 @@
 	SetControlInputsFromVectorx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,X);
 	ControlInputSetGradientx(femmodel->elements,femmodel->nodes,femmodel->vertices,femmodel->loads,femmodel->materials,femmodel->parameters,G);
-	for(int i=0;i<num_controls;i++){
-		InputToResultx(femmodel,control_list[i]);
-	}
+	femmodel->OutputControlsx(&femmodel->results);
 	femmodel->results->AddObject(new GenericExternalResult<double*>(femmodel->results->Size()+1,JEnum,user.J,maxiter+3,1,1,0));