pari-gp - PARI/GP の matsolve メソッドの潜在的なバグ

Question

メソッドに不十分なサイズの行列を指定matsolve(M, B)すると、警告の代わりにセグメンテーション違反が発生しました。

  *** matsolve: bug in PARI/GP (Segmentation Fault), please report.

行列のサイズが正しいと、間違った結果が得られました。ここに私のテストコードがあります:

{
w = ffgen(w^7 + w^2 + Mod(2,3)*w + Mod(1,3));

info = [
  2*w^6 + 2*w^5 + 2*w^4 + 2*w^3 + 2*w^2 + w, w^6 + w^5 + w^4 + 2*w^2 + 2;
  w^4 + w^3 + 2*w^2 + w + 1, w^6 + 2*w^5 + 2*w^4 + w^3 + 2*w^2;
  2*w^6 + w^5 + w^4 + w^3 + 2, 2*w^5 + 2*w^4 + 2*w^3 + 2*w + 2;
  w^5 + w^3 + 2*w^2 + w, w^6 + 2*w^5 + 2*w^4 + 2*w^3 + w^2 + 1;
  w^3 + w^2 + 2*w + 1, 2*w^6 + w^5 + 2*w^3 + w;
  2*w^6 + 2*w^5 + w^4 + 2*w^3 + 2*w + 2, w^6 + 2*w^4 + 2*w^2 + 1;
  w^6 + 2*w^5 + 2*w^4 + 2*w^3 + 2*w^2 + 2*w + 1, w^6 + 2*w^4 + 2*w^3 + w^2 + w + 1
];

result = [
  w^6 + 2*w^5 + w^3 + 2*w;
  2*w^6 + w^2 + w;
  w^4 + 2*w^2 + 2*w + 2;
  2*w^6 + 2*w^5 + w^4 + w^3 + 2*w^2;
  2*w^5 + w^4 + w;
  w^5 + 2*w^2 + w + 2;
  2*w^6 + w^5 + 2*w^4 + 2*w^3 + 2*w^2 + w
];

\\ lamda * info = result
lamda = matsolve(info, result);

calculated_result = mattranspose(mattranspose(lamda) *
                                 mattranspose(info));
}

出力は次のとおりです。

? result
[                w^6 + 2*w^5 + w^3 + 2*w]
[                        2*w^6 + w^2 + w]
[                  w^4 + 2*w^2 + 2*w + 2]
[      2*w^6 + 2*w^5 + w^4 + w^3 + 2*w^2]
[                        2*w^5 + w^4 + w]
[                    w^5 + 2*w^2 + w + 2]
[2*w^6 + w^5 + 2*w^4 + 2*w^3 + 2*w^2 + w]

? calculated_result
[              w^6 + 2*w^5 + w^3 + 2*w]
[                      2*w^6 + w^2 + w]
[    2*w^5 + w^4 + w^3 + w^2 + 2*w + 2]
[w^6 + w^5 + 2*w^4 + 2*w^3 + 2*w^2 + w]
[      w^6 + w^4 + w^3 + w^2 + 2*w + 1]
[        2*w^4 + 2*w^3 + w^2 + 2*w + 1]
[                  2*w^4 + w^3 + w + 1]

結果の最初の 2 行のみが同じで、残りは異なります。特定の info(nxk) および result(kx 1) 行列について、GP は info 行列の最後の (n - k) 行を無視すると思います。