Bill Allombert on Fri, 12 Nov 1999 13:41:44 +0100 (MET)


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

Re: matdet() bug 

A patch for matdet bug.

However, I do not see why nbco==1 was a special case.

(CVS version updated)
Index: src/basemath/alglin1.c
===================================================================
RCS file: /home/megrez/cvsroot/pari/src/basemath/alglin1.c,v
retrieving revision 1.2
diff -u -r1.2 alglin1.c
--- src/basemath/alglin1.c      1999/10/31 17:59:44     1.2
+++ src/basemath/alglin1.c      1999/11/12 12:39:44
@@ -1971,7 +1971,7 @@
     if (DEBUGLEVEL > 7) msgtimer("det, col %ld / %ld",i,nbco-1);
   }
   p = gcoeff(a,nbco,nbco);
-  if (s < 0) p = gneg(p); else if (nbco==1) p = gcopy(p);
+  if (s < 0) p = gneg(p); else p = gcopy(p);
   return gerepileupto(av, p);
 }