portelli/Grid
1
0
Fork 0
mirror of https://github.com/paboyle/Grid.git synced 2024-11-15 02:05:37 +00:00
Code Releases Activity

Further fixes from multidimensional array

Browse Source
This commit is contained in:
Chulwoo Jung 2017-06-05 23:13:41 -04:00
parent 927c7ae3ed
commit b1b15f0b70
1 changed files with 10 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
lib/algorithms/iterative/ImplicitlyRestartedLanczos.h
View File

@ -297,11 +297,15 @@ public:
if (j==(i-1)) EE[j] = lme[j]; if (j==(i-1)) EE[j] = lme[j];
} }
LAPACK_INT evals_found; LAPACK_INT evals_found;
LAPACK_INT lwork = ( (18*NN) > (1+4*NN+NN*NN)? (18*NN):(1+4*NN+NN*NN)) ; // LAPACK_INT lwork = ( (18*NN) > (1+4*NN+NN*NN)? (18*NN):(1+4*NN+NN*NN)) ;
LAPACK_INT lwork = 1+(18*NN) ;
LAPACK_INT liwork = 3+NN*10 ; LAPACK_INT liwork = 3+NN*10 ;
LAPACK_INT iwork[liwork]; // LAPACK_INT iwork[liwork];
double work[lwork]; // double work[lwork];
LAPACK_INT isuppz[2*NN]; // LAPACK_INT isuppz[2*NN];
std::vector<LAPACK_INT> iwork(liwork);
std::vector<double> work(lwork);
std::vector<LAPACK_INT> isuppz(2*NN);
char jobz = 'V'; // calculate evals & evecs char jobz = 'V'; // calculate evals & evecs
char range = 'I'; // calculate all evals char range = 'I'; // calculate all evals
// char range = 'A'; // calculate all evals // char range = 'A'; // calculate all evals
@ -332,8 +336,8 @@ public:
&vl, &vu, &il, &iu, // these four are ignored if second parameteris 'A' &vl, &vu, &il, &iu, // these four are ignored if second parameteris 'A'
&tol, // tolerance &tol, // tolerance
&evals_found, evals_tmp.data(), evec_tmp.data(), &NN, &evals_found, evals_tmp.data(), evec_tmp.data(), &NN,
isuppz, isuppz.data(),
work, &lwork, iwork, &liwork, work.data(), &lwork, iwork.data(), &liwork,
&info); &info);
for (int i = iu-1; i>= il-1; i--){ for (int i = iu-1; i>= il-1; i--){
printf("node=%d evals_found=%d evals_tmp[%d] = %g\n",node,evals_found, i - (il-1),evals_tmp[i - (il-1)]); printf("node=%d evals_found=%d evals_tmp[%d] = %g\n",node,evals_found, i - (il-1),evals_tmp[i - (il-1)]);