| 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);