Browse Source

use limexp to prevent NaN with extreme exponents

pre-master-46
dwarning 2 years ago
committed by Holger Vogt
parent
commit
3684db0126
  1. 2
      src/spicelib/analysis/ninteg.c

2
src/spicelib/analysis/ninteg.c

@ -38,7 +38,7 @@ Nintegrate (double noizDens, double lnNdens, double lnNlstDens, Ndata *data)
if (fabs(exponent) < N_INTUSELOG) {
return (a * (data->lnFreq - data->lnLastFreq));
} else {
return (a * ((exp(exponent * data->lnFreq) - exp(exponent * data->lnLastFreq)) /
return (a * ((limexp(exponent * data->lnFreq) - limexp(exponent * data->lnLastFreq)) /
exponent));
}
}

Loading…
Cancel
Save