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

Fix a bug when evaluating -0.5^3. 

Use code which has been up to now only applied with LT compatibility.
This commit is contained in:
Holger Vogt 2026-04-11 14:12:52 +02:00
parent afd4eb64ed
commit b9221a07a3
1 changed files with 14 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

33
src/spicelib/parser/ptfuncs.c
View File

@ -69,24 +69,21 @@ double
PTpower(double arg1, double arg2)
{
double res;
if (newcompat.lt) {
if (arg1 == 0)
if (arg1 == 0)
res = 0;
else if(arg1 > 0)
res = pow(arg1, arg2);
else {
/* If arg2 is quasi an integer, round it to have pow not fail
when arg1 is negative. Takes into account the double
representation which sometimes differs in the last digit(s). */
if (AlmostEqualUlps(nearbyint(arg2), arg2, 10))
res = pow(arg1, round(arg2));
else
/* As per LTSPICE specification for ** */
res = 0;
else if(arg1 > 0)
res = pow(arg1, arg2);
else {
/* If arg2 is quasi an integer, round it to have pow not fail
when arg1 is negative. Takes into account the double
representation which sometimes differs in the last digit(s). */
if (AlmostEqualUlps(nearbyint(arg2), arg2, 10))
res = pow(arg1, round(arg2));
else
/* As per LTSPICE specification for ** */
res = 0;
}
}
else
res = pow(fabs(arg1), arg2);
return res;
}
@ -106,7 +103,7 @@ PTpowerH(double arg1, double arg2)
res = pow(arg1, arg2);
}
}
else if (newcompat.lt) {
else {
if (arg1 >= 0)
res = pow(arg1, arg2);
else {
@ -120,8 +117,6 @@ PTpowerH(double arg1, double arg2)
res = 0;
}
}
else
res = pow(fabs(arg1), arg2);
return res;
}