SLEPcのLanczos使用時のエラー

[t-sakashita]

rokko/test/sparseで、以下を実行した。

mpirun -np 4 ./laplacian_mfree_mpi

routine=lanczos
[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------
[0]PETSC ERROR: Missing or incorrect user input 
[0]PETSC ERROR: You should explicitly provide the reorthogonalization type, e.g., -eps_lanczos_reorthog local
          ...   Note that the EPSLANCZOS solver is *NOT RECOMMENDED* for general use, because it uses
          ...   explicit restart which typically has slow convergence. The recommended solver is
          ...   EPSKRYLOVSCHUR (the default), which implements Lanczos with thick restart in the
          ...   case of symmetric/Hermitian problems
[0]PETSC ERROR: See https://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.15.1, Jun 17, 2021 

[t-sakashita]

SLEPcのver 3.15からは、Lanczos法の再直行化アルゴリズムのデフォルトを設定する必要がある。

EPSLANCZOS now requires that EPSLanczosSetReorthog() is called explicitly, or the corresponding command-line option is given. The intention is that inexperienced users realize that they should use EPSKRYLOVSCHUR instead of EPSLANCZOS.

https://gitlab.com/slepc/slepc/-/blob/main/CHANGELOG.md

Rokkoでは、デフォルトを"local"とした。 577bfded4fddacde416e1922abc2a19dab58dc67