Ignore:
Timestamp:
09/12/07 17:43:11 (17 years ago)
Author:
smasson
Message:

use of double precision to avoid rounding errors

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/SRC/Interpolation/compute_fromreg_imoms3_weigaddr.pro

    r262 r282  
    284284  weig[15, short] = wx3 * wy3 
    285285; 
    286   mi = min(total(weig[*, short], 1), max = ma) 
     286  mi = min(total(weig[*, short], 1, /double), max = ma) 
    287287  IF abs(mi-1) GE 1.e-6 THEN stop 
    288288  IF abs(ma-1) GE 1.e-6 THEN stop 
    … …  
    385385    weig[7, short] = wx3 * wy2 
    386386; 
    387     mi = min(total(weig[*, short], 1), max = ma) 
     387    mi = min(total(weig[*, short], 1, /double), max = ma) 
    388388    IF abs(mi-1) GE 1.e-6 THEN stop 
    389389    IF abs(ma-1) GE 1.e-6 THEN stop 
    … …  
    451451    weig[3, short] = wx3 
    452452; 
    453     mi = min(total(weig[*, short], 1), max = ma) 
     453    mi = min(total(weig[*, short], 1, /double), max = ma) 
    454454    IF abs(mi-1) GE 1.e-6 THEN stop 
    455455    IF abs(ma-1) GE 1.e-6 THEN stop 
    … …  
    517517    weig[3, short] = wx3 
    518518; 
    519     mi = min(total(weig[*, short], 1), max = ma) 
     519    mi = min(total(weig[*, short], 1, /double), max = ma) 
    520520    IF abs(mi-1) GE 1.e-6 THEN stop 
    521521    IF abs(ma-1) GE 1.e-6 THEN stop 
Note: See TracChangeset for help on using the changeset viewer.