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

limiting exp()

This commit is contained in:
dwarning 2010-03-07 00:14:13 +00:00
parent 54046ea353
commit f73055d7be
1 changed files with 2 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -21,6 +21,7 @@ Author: 1987 Gary W. Ng
#include "ngspice.h" #include "ngspice.h"
#include "noisedef.h" #include "noisedef.h"
#define limexp(x) (x > 700 ? exp(700.0)*(1.0+x-700.0) : exp(x))
double double
Nintegrate (double noizDens, double lnNdens, double lnNlstDens, Ndata *data) Nintegrate (double noizDens, double lnNdens, double lnNlstDens, Ndata *data)
@ -32,7 +33,7 @@ Nintegrate (double noizDens, double lnNdens, double lnNlstDens, Ndata *data)
if ( fabs(exponent) < N_INTFTHRESH ) { if ( fabs(exponent) < N_INTFTHRESH ) {
return (noizDens * data->delFreq); return (noizDens * data->delFreq);
} else { } else {
a = exp(lnNdens - exponent*data->lnFreq); a = limexp(lnNdens - exponent*data->lnFreq);
exponent += 1.0; exponent += 1.0;
if (fabs(exponent) < N_INTUSELOG) { if (fabs(exponent) < N_INTUSELOG) {
return (a * (data->lnFreq - data->lnLastFreq)); return (a * (data->lnFreq - data->lnLastFreq));