Carl R. Witty on Sun, 14 May 2000 13:33:26 -0700


[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

patch for previously reported polsturm bug


Here's a patch to fix the bug I just reported against polsturm.

After the patch, I get the correct answer for the previously buggy
case:

(13:28) gp > polsturm(x, 0, 1)
%1 = 0

--- polarit2.c~ Mon Feb  7 07:29:08 2000
+++ polarit2.c  Sun May 14 13:28:12 2000
@@ -2456,6 +2456,7 @@
   {
     s = b? gsigne(poleval(x,b)):  sl;
     t = a? gsigne(poleval(x,a)): -sl;
+    if (t == 0) return 0;
     avma = av; return (s == t)? 0: 1;
   }
   u=gdiv(x,content(x)); v=derivpol(x);