Changeset 5010

Timestamp:

08/06/10 08:29:39 (15 years ago)

Author:

Mathieu Morlighem

Message:

better error message

File:

• issm/trunk/src/c/objects/Bamg/Triangles.cpp (modified) (3 diffs)

issm/trunk/src/c/objects/Bamg/Triangles.cpp

@@ -4571 +4571 @@
 
         // computation of the det 
-        int Nberr=0;
+        int number_of_errors=0;
         for (i=0;i<nbt;i++) {
                 MeshVertex & v0 = triangles[i][0];
@@ -4579 +4579 @@
                 //If this is not a boundary triangle
                 if ( &v0 && &v1 &&  &v2 ){
+
+                        /*Compute determinant*/
                         triangles[i].det= det(v0,v1,v2);
-                        if (triangles[i].det <=0 && Nberr++ <10){
-                                if(Nberr==1)
-                                 if (strfrom){
-                                         ISSMERROR("Fatal error %s (SetInCoor) : area of Triangle %i < 0",strfrom,i);
-                                 }
-                                 else{
-                                         ISSMERROR("Fatal error (SetInCoor) : area of Triangle %i < 0",i);
-                                 }
+
+                        /*Check that determinant is positive*/
+                        if (triangles[i].det <=0){
+
+                                /*increase number_of_errors and print error only for the first 20 triangles*/
+                                number_of_errors++;
+                                if (number_of_errors<20){
+                                        printf("Area of Triangle %i < 0 (det=%i)\n",i+1,triangles[i].det);
+                                }
                         }
                 }
@@ -4594 +4597 @@
                 else triangles[i].det=-1;
         }
+
+        if (number_of_errors) ISSMERROR("Fatal error: some triangles have negative areas, see above");
 }
 /*}}}1*/