luke
/
ngspice
mirror of https://git.code.sf.net/p/ngspice/ngspice
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity

use limexp to prevent NaN with extreme exponents

This commit is contained in:
dwarning 2024-04-04 16:19:35 +02:00
parent f40285c2a2
commit 22869b6e39
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/spicelib/analysis/ninteg.c
View File

@ -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));
}
}