Changeset 14450

Timestamp:

03/13/13 20:30:28 (12 years ago)

Author:

habbalf

Message:

BUG: checking matrix is n=m, not m=m

File:

• issm/trunk-jpl/src/c/modules/Solverx/SolverxPetsc.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/modules/Solverx/SolverxPetsc.cpp

@@ -67 +67 @@
         /*First, check that f-set is not NULL, i.e. model is fully constrained:*/ 
         _assert_(Kff);
-        MatGetSize(Kff,&global_m,&global_n); _assert_(global_m==global_m);
+        MatGetSize(Kff,&global_m,&global_n); _assert_(global_m==global_n);
         if(!global_n){
                 *puf=NewVec(0,IssmComm::GetComm()); return;