t-sakashita / rokko

Integrated Interface for libraries of eigenvalue decomposition
Boost Software License 1.0
10 stars 2 forks source link

scalapack_pdsyevxがエラー #30

Open rigarash opened 11 years ago

rigarash commented 11 years ago

$ mpirun -np 3 benchmark/frank_matrix scalapack_pdsyevx 1000 として実行すると(VPS上で、3CPUのマシンです)、次のエラーで落ちました。 eigen_s, eigen_sx, scalapack, scalapack_pdsyevdは問題ないです。 Debianのscalapack(1.8.0), OpenMP(1.4.5), gcc(4.7.3)です。

これは、rokkoではなくて、scalapack自身のバグでしょうか?

... (長く表示されて)
Rank = 0  myrow=0 mycol=0
local_num_block_cols=1  local_rest_block_cols=1  rest_block_col=0
Rank = 1  myrow=0 mycol=1
local_num_block_cols=1  local_rest_block_cols=0  rest_block_col=0
Rank = 2  myrow=0 mycol=2
local_num_block_cols=1  local_rest_block_cols=0  rest_block_col=0
Rank = 1  myrow=0 mycol=1
local_num_block_rows=3  local_rest_block_rows=1  rest_block_row=0
Rank = 0  myrow=0 mycol=0
local_num_block_rows=3  local_rest_block_rows=1  rest_block_row=0
Rank = 2  myrow=0 mycol=2
local_num_block_rows=3  local_rest_block_rows=1  rest_block_row=0
proc=0
  mb=333  nb=333
  nprow=1  npcol=3
  m_local=1000 n_local=334
  myrow=0 mycol=0
  lld=1000
  length_array=334000
proc=1
  mb=333  nb=333
  nprow=1  npcol=3
  m_local=1000 n_local=333
  myrow=0 mycol=1
  lld=1000
  length_array=333000
proc=2
  mb=333  nb=333
  nprow=1  npcol=3
  m_local=1000 n_local=333
  myrow=0 mycol=2
  lld=1000
  length_array=333000
get_count:0
## error at pdsyevx function. info=2

mpirun has exited due to process rank 0 with PID 28307 on
node sakuravps2 exiting without calling "finalize". This may
have caused other processes in the application to be
## terminated by signals sent by mpirun (as reported here).
t-sakashita commented 11 years ago

abstol = 10-3しても,dim=80でエラー(info=2)

原因は未だ不明 返値info=2である.

以下のpdsyevxの引数infoの記述より. http://www.netlib.org/scalapack/double/pdsyevx.f

if (MOD(INFO,2).NE.0), then one or more eigenvectors

ifail等を詳細に調べてみる.