From 22869b6e3922940c845443e2c93b6e6fc03d4f84 Mon Sep 17 00:00:00 2001 From: dwarning Date: Thu, 4 Apr 2024 16:19:35 +0200 Subject: [PATCH] use limexp to prevent NaN with extreme exponents --- src/spicelib/analysis/ninteg.c | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/spicelib/analysis/ninteg.c b/src/spicelib/analysis/ninteg.c index dd347e2da..43c5b6732 100644 --- a/src/spicelib/analysis/ninteg.c +++ b/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)); } }