IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 2979


Ignore:
Timestamp:
Jan 13, 2005, 11:25:44 AM (22 years ago)
Author:
desonia
Message:

modified the checking to allow for float-point inaccuracies.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/psLib/test/dataManip/tst_psVectorFFT.c

    r2978 r2979  
    66*  @author Robert DeSonia, MHPCC
    77*
    8 *  @version $Revision: 1.11 $ $Name: not supported by cvs2svn $
    9 *  @date $Date: 2005-01-13 21:23:56 $
     8*  @version $Revision: 1.12 $ $Name: not supported by cvs2svn $
     9*  @date $Date: 2005-01-13 21:25:44 $
    1010*
    1111*  Copyright 2004 Maui High Performance Computing Center, University of Hawaii
     
    597597        val1 = (cabsf(vec3->data.C64[n]) * cabsf(vec3->data.C64[n]) +
    598598                cabsf(vec3->data.C64[100-n]) * cabsf(vec3->data.C64[100-n]))/100/100;
    599         if ( fabsf( val - vec2->data.F32[ n ] ) > FLT_EPSILON ) {
     599        if ( fabsf( val - vec2->data.F32[ n ] ) > 10*FLT_EPSILON ) {
    600600            psError(PS_ERR_UNKNOWN,true, "psVectorPowerSpectrum result is invalid (n=%d, %.2f %.2f)",
    601601                    n, vec2->data.F32[ n ], val );
    602602            return 2;
    603603        }
    604         if (fabsf(val1 - vec4->data.F64[n]) > 0.001) {
     604        if (fabsf(val1 - vec4->data.F64[n]) > 10*FLT_EPSILON) {
    605605            psError(PS_ERR_UNKNOWN,true,"psVectorPowerSpectrum result is invalid (n=%d, %.2f %.2f)",n,vec4->data.F64[n],val1);
    606606            return 22;
     
    609609
    610610    val = cabsf( vec->data.C32[ 50 ] ) * cabsf( vec->data.C32[ 50 ] ) / 100 / 100;
    611     if ( fabsf( vec2->data.F32[ 50 ] - val ) > FLT_EPSILON ) {
     611    if ( fabsf( vec2->data.F32[ 50 ] - val ) > 10*FLT_EPSILON ) {
    612612        psError(PS_ERR_UNKNOWN,true, "psVectorPowerSpectrum result is invalid (n=50, %.2f %.2f)",
    613613                vec2->data.F32[ 0 ], val );
Note: See TracChangeset for help on using the changeset viewer.