|Jeroen Demeyer on Sun, 11 Oct 2009 17:01:10 +0200|
[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]
|??gaussred is wrong|
Hello pari-dev,In the ??gaussred help, it is said that the result will be upper-triangular but this is not true. I suppose this happens when the given matrix is not (positive) definite. The formula which is written in ??gaussred is correct though, if one replaces j>i by j!=i. So I think the docs needs to clarified.
gp> qfgaussred([0,1;1,0]) %302 = [1/2 1] [-1 -1/2]