source: issm/oecreview/Archive/11681-11700/ISSM-11686-11687.diff@ 11991

proj/ice/larour/issm-uci-clean/trunk-jpl/src/c/modules/Solverx/Solverx.cpp

@@ -20 +20 @@
         uf=new Vector();
 
         #ifdef _HAVE_PETSC_
-        Vec uf0_vector=NULL;
-        if (uf0)uf0_vector=uf0->vector;
+        Vec uf0_vector = NULL;
+        Vec df_vector  = NULL;
+        if(uf0) uf0_vector=uf0->vector;
+        if(df)  df_vector = df->vector;
 
-        SolverxPetsc(&uf->vector, Kff->matrix, pf->vector, uf0_vector, df->vector, parameters);
+        SolverxPetsc(&uf->vector,Kff->matrix,pf->vector,uf0_vector,df_vector,parameters);
         VecGetSize(uf->vector,&uf->M);
         #else
         _error_("not supported yet!");