Commit 23919679 authored by Andreas Marek's avatar Andreas Marek

Pass wantDebug in legacy interface

parent 8d05d19c
......@@ -70,7 +70,7 @@
integer(kind=ik) :: nev
logical, intent(in) :: wantDebug
logical :: success
integer(kind=ik) :: successInternal, debug
integer(kind=ik) :: successInternal
class(elpa_t), pointer :: e
......
......@@ -64,12 +64,12 @@
#endif
implicit none
integer(kind=ik) :: na, lda, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols
integer(kind=ik) :: na, lda, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols
#if REALCASE == 1
#ifdef USE_ASSUMED_SIZE
real(kind=REAL_DATATYPE) :: a(lda,*)
real(kind=REAL_DATATYPE) :: a(lda,*)
#else
real(kind=REAL_DATATYPE) :: a(lda,matrixCols)
real(kind=REAL_DATATYPE) :: a(lda,matrixCols)
#endif
#endif
......@@ -81,9 +81,9 @@
#endif
#endif
! integer(kind=ik) :: my_prow, my_pcol, np_rows, np_cols, mpierr
! integer(kind=ik) :: l_cols, l_rows, l_col1, l_row1, l_colx, l_rowx
! integer(kind=ik) :: n, nc, i, info, ns, nb
! integer(kind=ik) :: my_prow, my_pcol, np_rows, np_cols, mpierr
! integer(kind=ik) :: l_cols, l_rows, l_col1, l_row1, l_colx, l_rowx
! integer(kind=ik) :: n, nc, i, info, ns, nb
!#if REALCASE == 1
! real(kind=REAL_DATATYPE), allocatable :: tmp1(:), tmp2(:,:), tmat1(:,:), tmat2(:,:)
!#endif
......@@ -92,8 +92,6 @@
!#endif
logical, intent(in) :: wantDebug
logical :: success
! integer(kind=ik) :: istat
! character(200) :: errorMessage
integer(kind=iK) :: successInternal
class(elpa_t), pointer :: e
......@@ -128,6 +126,10 @@
return
endif
if (wantDebug) then
call e%set("debug",1)
endif
call e%invert_triangular(a(1:lda,1:matrixCols), successInternal)
if (successInternal .ne. ELPA_OK) then
......
......@@ -64,15 +64,15 @@
use timings_dummy
#endif
implicit none
integer(kind=ik) :: na, nev, ldq, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols
real(kind=REAL_DATATYPE) :: d(na), e(na)
integer(kind=ik) :: na, nev, ldq, nblk, matrixCols, mpi_comm_rows, mpi_comm_cols
real(kind=REAL_DATATYPE) :: d(na), e(na)
#ifdef USE_ASSUMED_SIZE
real(kind=REAL_DATATYPE) :: q(ldq,*)
real(kind=REAL_DATATYPE) :: q(ldq,*)
#else
real(kind=REAL_DATATYPE) :: q(ldq,matrixCols)
real(kind=REAL_DATATYPE) :: q(ldq,matrixCols)
#endif
logical, intent(in) :: wantDebug
logical, intent(in) :: wantDebug
logical :: success
integer(kind=iK) :: successInternal
class(elpa_t), pointer :: obj
......@@ -105,6 +105,9 @@
return
endif
if (wantDebug) then
call obj%set("debug",1)
endif
call obj%solve_tridi(d(1:na), e(1:na), q(1:ldq,1:matrixCols), successInternal)
if (successInternal .ne. ELPA_OK) then
......
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