Ignore:
Timestamp:
09/23/16 14:23:16 (5 years ago)
Author:
ymipsl
Message:

remapper : fix rounding error which stop the remaping.

YM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • XIOS/trunk/extern/remap/src/polyg.cpp

    r849 r950  
    173173                t[1] = x[i]; 
    174174                t[2] = x[ii]; 
    175                 assert(scalarprod(crossprod(t[1] - t[0], t[2] - t[0]), t[0]) >= 0); // Error: tri a l'env (wrong orientation) 
     175                double sc=scalarprod(crossprod(t[1] - t[0], t[2] - t[0]), t[0]) ; 
     176                assert(sc >= -1e-10); // Error: tri a l'env (wrong orientation) 
    176177                double area_gc = triarea(t[0], t[1], t[2]); 
    177178                double area_sc_gc_moon = 0; 
Note: See TracChangeset for help on using the changeset viewer.