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/quadrilateral2square.pro

    r242 r282  
    3838; Can be scalar or array. 
    3939; 
    40 ; @keyword PERF 
     40; @keyword PERF {type=salar 0 or 1}{default=0} 
     41; activate to print the elapsed time spent within quadrilateral2square 
     42; 
     43; @keyword DOUBLE {type=salar 0 or 1}{default=0} 
     44; activate to perform double precision computation 
    4145; 
    4246; @returns 
    … …  
    7579;- 
    7680; 
    77 FUNCTION quadrilateral2square, x0in, y0in, x1in, y1in, x2in, y2in, x3in, y3in, xxin, yyin, PERF = perf 
     81FUNCTION quadrilateral2square, x0in, y0in, x1in, y1in, x2in, y2in, x3in, y3in, xxin, yyin, PERF = perf, DOUBLE = double 
    7882; 
    7983  compile_opt idl2, strictarrsubs 
    … …  
    111115; get the matrix A 
    112116; 
    113   a = square2quadrilateral(x0in, y0in, x1in, y1in, x2in, y2in, x3in, y3in) 
     117  a = square2quadrilateral(x0in, y0in, x1in, y1in, x2in, y2in, x3in, y3in, DOUBLE = double) 
    114118; 
    115119; compute the adjoint matrix 
Note: See TracChangeset for help on using the changeset viewer.