helmertのテストのエラー

[t-sakashita]

eigen_exaのeigen_sxを使用した場合にだけエラーが起こる。

rokko/test/diagonalize/helmert_mpi.cpp(57): error in "test_solver": difference{1.921%} between w.sum(){4952.9896907139155} and diag.sum(){5050} exceeds 0.0001

[t-sakashita]

考えられる原因：

• 行列・行列積で使用するBLACS周りが悪さをしている。（同サイズのFrank行列、ScaLAPACKに依存しないソルバでは起こらないから。）
  • BLACSに正しいlldが渡されているか。
• 行列の作り方は正しいか？　→表示して確認した。
• 逐次で作ったものをscatterルーチンで通信して確かめる。

[t-sakashita]

exampleの出力でも、微妙に許容できない誤差やnanが出ている。

mpirun -np 1 ./helmert_mpi eigen_exa:penta 9
largest eigenvalues: nan nan nan nan nan nan nan nan nan

上記は、-npが2の時もnan, -npが4だとうまくいく。

mpirun -np 1 ./helmert_mpi eigen_exa:penta 100
largest eigenvalues: 100 99 97.989691 96 95 94 93 92 91 90

mpirun -np 1 ./helmert_mpi eigen_exa:penta 20
largest eigenvalues: 20 19 17.94118 16 15 14 13 12 11 10 9 8 7 6 5 4 3 2 1 -7.041049e-15

[t-sakashita]

原因の切り分け： 生成した行列値を書き出し、生成関数を定義し、generateする。 この方法では、pdgemm, scatter等の並列処理を用いずに行列生成ができたことになる。 24602df56b282ce26a7ff836a46654fb138a6994 このような行列生成法でも、eigen_sxで上記の問題が起こる。

結論：

• RokkoやEigenExaのバグではない。
• この行列は、eigen_sxのアルゴリズムが苦手とする行列なのではないか？

[t-sakashita]

作者に報告する。

Rokkoを使わない場合、EigenExa単体でも発現することを確認。 ソルバのバグと判明。

[wistaria]

EigenExa バージョン 2.4b でもまだ治っていない test/diagonalize/helmert_mpi が fail する

[wistaria]

とりあえず eigenexa のテストをスキップ c26c082e38e67df11c312ea9e2e50ebfc2e08573