お世話になっております. 36サイトスピン系のLanczos法の計算をリスタートモードを使いながら計算すると La

<a class="user-mention notranslate" data-hovercard-type="user" data-hovercard-url="/us

対角化のリスタートについて2 about hphi HOT 4 CLOSED

issp-center-dev commented on September 25, 2024

対角化のリスタートについて2

from hphi.

Comments (4)

k-yoshimi commented on September 25, 2024

@ultimatile
ご報告ありがとうございます。
動作確認の上修正後ご連絡しますので、それまで今しばらくお待ちください。

from hphi.

k-yoshimi commented on September 25, 2024

@ultimatile
リスタート時にメモリ再確保する際の処理について、
三重対角化行列の固有ベクトルのメモリ再確保に漏れがありました。
テストでは偶然通っていたようです。ご指摘大変助かりました。

次のバージョンでは本修正を取り入れる予定ですが、取り急ぎは
Lanczos_EigenValue.cにあるReadTMComponents関数のメモリ再確保部分について、

 if(iFlg==0) {
    alpha = (double *) realloc(alpha, sizeof(double) * (i_max + X->Def.Lanczos_max + 1));
    beta = (double *) realloc(beta, sizeof(double) * (i_max + X->Def.Lanczos_max + 1));
    for (i = 0; i < X->Def.nvec + 1; i++) {
      vec[i] = (complex double *) realloc(vec[i], (i_max + X->Def.Lanczos_max + 1) * sizeof(complex double));
    }
  }
  else if(iFlg==1){
    alpha=(double*)realloc(alpha, sizeof(double)*(i_max + 1));
    beta=(double*)realloc(beta, sizeof(double)*(i_max + 1));
    for (i = 0; i < X->Def.nvec + 1; i++) {
      vec[i] = (complex double *) realloc(vec[i], (i_max + X->Def.Lanczos_max + 1) * sizeof(complex double));
    }
  }

と変更しvecについてメモリの再確保をするようにした上で、リビルドしてご使用いただけるでしょうか？
お手数ですがよろしくお願いします。

from hphi.