Igor Schein on Sat, 22 May 1999 16:54:38 -0400


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

quadhilbert()


Hi,

when compiled with pgcc 1.1.3 under linux, quadhilbert(4364)
return "bug in quadhibert" error message.  The following patch hides
the problem.  Any idea how to fix the problem fundamentally? 

Thanks

Igor

--- src/modules/stark.c~	Thu May 20 13:05:27 1999
+++ src/modules/stark.c	Sat May 22 16:38:57 1999
@@ -265,11 +265,12 @@
 static GEN
 ComputeKernel0(GEN P, GEN DA, GEN DB, long nbA, long nbB)
 {
-  long rk, av = avma;
+  long rk, av = avma, j;
   GEN herm, mask1, mask2, U;
 
   herm  = hnfall(concat(P, DB));
-  rk = nbA + nbB + 1 - lg((GEN)herm[1]);
+  j = lg((GEN)herm[1]);
+  rk = nbA + nbB + 1 - j;
 
   mask1 = subis(shifti(gun, nbA), 1);
   mask2 = subis(shifti(gun, rk), 1);