Changeset 950 for XIOS/trunk/extern

Timestamp:

09/23/16 14:23:16 (8 years ago)

Author:

ymipsl

Message:

remapper : fix rounding error which stop the remaping.

YM

File:

• XIOS/trunk/extern/remap/src/polyg.cpp (modified) (1 diff)

XIOS/trunk/extern/remap/src/polyg.cpp

@@ -173 +173 @@
                 t[1] = x[i];
                 t[2] = x[ii];
-                assert(scalarprod(crossprod(t[1] - t[0], t[2] - t[0]), t[0]) >= 0); // Error: tri a l'env (wrong orientation)
+                double sc=scalarprod(crossprod(t[1] - t[0], t[2] - t[0]), t[0]) ;
+                assert(sc >= -1e-10); // Error: tri a l'env (wrong orientation)
                 double area_gc = triarea(t[0], t[1], t[2]);
                 double area_sc_gc_moon = 0;