Commit f5faf502 authored by Andreas Marek's avatar Andreas Marek

Fix the setup of elpa_cholesky

parent de792054
......@@ -62,7 +62,6 @@
complex(kind=COMPLEX_DATATYPE) :: a(lda,matrixCols)
#endif
#endif
integer(kind=ik) :: nev
logical, intent(in) :: wantDebug
logical :: success
integer(kind=ik) :: successInternal
......@@ -86,7 +85,6 @@
e => elpa_allocate()
call e%set("na", na)
call e%set("nev", nev)
call e%set("local_nrows", lda)
call e%set("local_ncols", matrixCols)
call e%set("nblk", nblk)
......@@ -94,6 +92,13 @@
call e%set("mpi_comm_rows", mpi_comm_rows)
call e%set("mpi_comm_cols", mpi_comm_cols)
!! the elpa object needs nev to be set (in case the EVP-solver is
!! called later. Thus it is set by user, do nothing, otherwise,
!! set it to na as default
!if (e%is_set("nev")) .ne. 1) then
! call e%set("nev", na)
!endif
if (e%setup() .ne. ELPA_OK) then
print *, "Cannot setup ELPA instance"
success = .false.
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment