Improve convergence in Infineon Power MOS.

0**something is 0
1 year ago · 8c806912d0
1 changed files with 3 additions and 1 deletions
--- a/src/spicelib/parser/ptfuncs.c
+++ b/src/spicelib/parser/ptfuncs.c
@ -70,7 +70,9 @@ 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