Index: ../trunk-jpl/src/c/analyses/controltao_core.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16483)
+++ ../trunk-jpl/src/c/analyses/controltao_core.cpp	(revision 16484)
@@ -81,14 +81,14 @@
 	if(VerboseControl()) _printf0_("   Starting optimization\n");
 	TaoSolve(tao);
 	TaoView(tao,PETSC_VIEWER_STDOUT_WORLD);
+
+	/*Save results*/
 	TaoGetSolutionVector(tao,&X->pvector->vector);
 	G=new Vector<IssmDouble>(0); VecFree(&G->pvector->vector);
 	TaoGetGradientVector(tao,&G->pvector->vector);
 	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));
 
 	/*Finalize*/