Changeset 2033 for issm/trunk/src/c/parallel/convergence.cpp

Timestamp:

08/27/09 17:15:17 (16 years ago)

Author:

seroussi

Message:

fixed control stedy, but stokes not working yet

File:

• issm/trunk/src/c/parallel/convergence.cpp (modified) (2 diffs)

issm/trunk/src/c/parallel/convergence.cpp

@@ -67 +67 @@
         VecNorm(pf,NORM_2,&nF);
         res=nKUoldF/nF;
-        if (isnan(res)) throw ErrorException(__FUNCT__,exprintf("mechanical equilibrium convergence criterion is NaN! "));
+        if (isnan(res)){
+                PetscSynchronizedPrintf(MPI_COMM_WORLD,"norm nf %lf \n",nF);
+                PetscSynchronizedFlush(MPI_COMM_WORLD);
+
+                PetscSynchronizedPrintf(MPI_COMM_WORLD,"norm kuold %lf \n",nKUoldF);
+                PetscSynchronizedFlush(MPI_COMM_WORLD);
+                throw ErrorException(__FUNCT__,exprintf("mechanical equilibrium convergence criterion is NaN! "));
+        }
 
         //clean up
@@ -89 +96 @@
                 VecDuplicate(old_uf,&duf);VecCopy(old_uf,duf); VecAYPX(duf,-1.0,uf);
                 VecNorm(duf,NORM_2,&ndu); VecNorm(old_uf,NORM_2,&nu);
+
                 if (isnan(ndu) || isnan(nu)) throw ErrorException(__FUNCT__,exprintf("convergence criterion is NaN! "));