Commit eaba0b63 authored by Gregor Michalicek's avatar Gregor Michalicek

Fixed a bug due optional parameters in franza call in diagonalization/eigen_diag.F90

Note that there are much more bugs of this type in diagonalization/eigen_diag.F90.
parent 2d82178b
......@@ -192,7 +192,13 @@ CONTAINS
jij,matind,nred,eig,ne_found,a_c,b_c,z_c)
endif
CASE (diag_lapack)
CALL franza(dimension%nbasfcn,ndim, lapw%nmat,(sym%l_zref.AND.(atoms%nlotot.EQ.0)), jij%l_j,matind,nred,input%gw,eig,ne_found,a_r,b_r,z_r,a_c,b_c,z_c)
IF (l_real) THEN
CALL franza(dimension%nbasfcn,ndim, lapw%nmat,(sym%l_zref.AND.(atoms%nlotot.EQ.0)),&
jij%l_j,matind,nred,input%gw,eig,ne_found,a_r=a_r,b_r=b_r,z_r=z_r)
ELSE
CALL franza(dimension%nbasfcn,ndim, lapw%nmat,(sym%l_zref.AND.(atoms%nlotot.EQ.0)),&
jij%l_j,matind,nred,input%gw,eig,ne_found,a_c=a_c,b_c=b_c,z_c=z_c)
END IF
CASE DEFAULT
!This should only happen if you select a solver by hand which was not compiled against
print*, "You selected a diagonalization scheme without compiling for it"
......
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