41 long ipHi, ipLo, typeOfRate;
48 for( ipHi=1; ipHi<=
iso_sp[ipISO][nelem].numLevels_max; ipHi++ )
51 for( ipLo=0; ipLo< ipHi; ipLo++ )
53 for( typeOfRate=0; typeOfRate<=1; typeOfRate++ )
55 if(
iso_sp[ipISO][nelem].
ex[ipHi][ipLo].Error[typeOfRate] >= 0. )
57 iso_sp[ipISO][nelem].ex[ipHi][ipLo].ErrorFactor[typeOfRate] =
59 ASSERT(
iso_sp[ipISO][nelem].
ex[ipHi][ipLo].ErrorFactor[typeOfRate] > 0. );
63 iso_sp[ipISO][nelem].ex[ipHi][ipLo].ErrorFactor[typeOfRate] = 1.0f;
70 iso_sp[ipISO][nelem].lgErrGenDone =
true;
void iso_put_error(long int ipISO, long int nelem, long int ipHi, long int ipLo, long int whichData, realnum errorOpt, realnum errorPess)