Changeset 6420 for issm/trunk/src/c/objects/Elements/Tria.cpp

Timestamp:

10/26/10 07:39:38 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added checks on solutions

File:

• issm/trunk/src/c/objects/Elements/Tria.cpp (modified) (5 diffs)

issm/trunk/src/c/objects/Elements/Tria.cpp

@@ -3722 +3722 @@
 
         /*Ok, we have vx and vy in values, fill in vx and vy arrays: */
-        for(i=0;i<numdof;i++) lambda[i]=values[i];
+        for(i=0;i<numdof;i++){
+                lambda[i]=values[i];
+                if(isnan(lambda[i])) _error_("NaN found in solution vector");
+        }
 
         /*Add vx and vy as inputs to the tria element: */
@@ -3752 +3755 @@
                 lambdax[i]=values[i*NDOF2+0];
                 lambday[i]=values[i*NDOF2+1];
+
+                /*Check solution*/
+                if(isnan(lambdax[i])) _error_("NaN found in solution vector");
+                if(isnan(lambday[i])) _error_("NaN found in solution vector");
         }
 
@@ -3789 +3796 @@
                 vx[i]=values[i*NDOF2+0];
                 vy[i]=values[i*NDOF2+1];
+
+                /*Check solution*/
+                if(isnan(vx[i])) _error_("NaN found in solution vector");
+                if(isnan(vy[i])) _error_("NaN found in solution vector");
         }
 
@@ -3848 +3859 @@
                 vx[i]=values[i*NDOF2+0];
                 vy[i]=values[i*NDOF2+1];
+
+                /*Check solution*/
+                if(isnan(vx[i])) _error_("NaN found in solution vector");
+                if(isnan(vy[i])) _error_("NaN found in solution vector");
         }
 
@@ -3901 +3916 @@
 
         /*Use the dof list to index into the solution vector: */
-        for(int i=0;i<numdof;i++) values[i]=solution[doflist[i]];
+        for(int i=0;i<numdof;i++){
+                values[i]=solution[doflist[i]];
+                if(isnan(values[i])) _error_("NaN found in solution vector");
+        }
 
         /*Add input to the element: */