Changeset 21606

Timestamp:

03/11/17 16:50:19 (8 years ago)

Author:

Eric.Larour

Message:

CHG: speed up in case we have an empty vector.

File:

• issm/branches/trunk-larour-GRL2017/src/c/modules/Solverx/Solverx.cpp (modified) (1 diff)

issm/branches/trunk-larour-GRL2017/src/c/modules/Solverx/Solverx.cpp

@@ -22 +22 @@
         if(VerboseModule()) _printf0_("   Solving matrix system\n");
 
-        /*Initialize solver: */
-        solver=new Solver<IssmDouble>(Kff,pf,uf0,df,parameters);
+        if (pf->IsEmpty()){ 
+                /*no need to solve! return empty vector: */
+                uf=new Vector<IssmDouble>(0);
+        }
+        else{
+                /*Initialize solver: */
+                solver=new Solver<IssmDouble>(Kff,pf,uf0,df,parameters);
 
-        /*Solve:*/
-        uf=solver->Solve();
+                /*Solve:*/
+                uf=solver->Solve();
 
-        /*clean up and assign output pointers:*/
-        delete solver;
+                /*clean up: */
+                delete solver;
+        }
+                
+        /*assign output pointers:*/
         *puf=uf;
 }