Commit ab666d92 authored by Matthias Redies's avatar Matthias Redies

remove unused variables from read_z interface

parent 68f05c5e
......@@ -113,7 +113,7 @@ CONTAINS
CALL z%init(olap%l_real, nbasfcn, input%neig)
CALL read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nk, jsp, z)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nk, jsp, z)
! calculate exchange contribution of current k-point nk to total energy (te_hfex)
! in the case of a spin-unpolarized calculation the factor 2 is added in eigen.F90
......
......@@ -228,7 +228,7 @@
allocate(carr2(maxval(hybdat%nbands), (atoms%lmaxd + 1)**2))
allocate(carr3(maxval(hybdat%nbands), (atoms%lmaxd + 1)**2))
DO ikpt = 1, nkpti
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, ikpt, jsp, z(ikpt))
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, ikpt, jsp, z(ikpt))
END DO
iatom = 0
......
......@@ -186,7 +186,7 @@ CONTAINS
IF (input%neig < hybdat%nbands(nk)) call judft_error(' mhsfock: neigd < nbands(nk) ;trafo from wavefunctions to APW requires at least nbands(nk)')
call z%init(olap%l_real, nbasfcn, input%neig)
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nk, jsp, z)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nk, jsp, z)
z%matsize2 = hybdat%nbands(nk) ! reduce "visible matsize" for the following computations
call olap%multiply(z, trafo)
......
......@@ -100,7 +100,7 @@ CONTAINS
! read in z coefficient from direct access file z at k-point nk
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nk, jsp, z)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nk, jsp, z)
! construct local orbital consisting of radial function times spherical harmonic
! where the radial function vanishes on the MT sphere boundary
......@@ -837,7 +837,7 @@ CONTAINS
! read in z coefficients from direct access file z at kpoint nk
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nk, jsp, z)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nk, jsp, z)
!CALL intgrf_init(atoms%ntype,atoms%jmtd,atoms%jri,atoms%dx,atoms%rmsh,hybdat%gridf)
gpt(:, 1:lapw%nv(jsp)) = lapw%gvec(:, 1:lapw%nv(jsp), jsp)
......
......@@ -122,8 +122,8 @@ CONTAINS
call z_kqpt%alloc(.true., nbasfcn, input%neig)
! read in z at k-point nk and nkqpt
CALL read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nk, jsp, z_nk)
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nkqpt, jsp, z_kqpt)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nk, jsp, z_nk)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nkqpt, jsp, z_kqpt)
g = maxval(abs(lapw%gvec(:, :lapw%nv(jsp), jsp)), dim=2) &
+ maxval(abs(lapw_nkqpt%gvec(:, :lapw_nkqpt%nv(jsp), jsp)), dim=2)&
......
......@@ -125,8 +125,8 @@ CONTAINS
call z_kqpt%alloc(.false., nbasfcn, input%neig)
! read in z at k-point ik and nkqpt
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, ik, jsp, z_nk)
call read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, nkqpt, jsp, z_kqpt)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, ik, jsp, z_nk)
call read_z(atoms, cell, hybdat, kpts, sym, noco, input, nkqpt, jsp, z_kqpt)
g = maxval(abs(lapw%gvec(:,:lapw%nv(jsp), jsp)), dim=2) &
+ maxval(abs(lapw_nkqpt%gvec(:,:lapw_nkqpt%nv(jsp), jsp)), dim=2)&
......
......@@ -179,16 +179,14 @@ contains
CALL write_matrix(mat, rec, id_olap)
END subroutine write_olap
subroutine read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, ik, jsp, z_out)
subroutine read_z(atoms, cell, hybdat, kpts, sym, noco, input, ik, jsp, z_out)
USE m_eig66_io
use m_types
use m_trafo
implicit none
type(t_atoms), intent(in) :: atoms
type(t_cell), intent(in) :: cell
type(t_mpdata), intent(in) :: mpdata
type(t_hybdat), intent(in) :: hybdat
type(t_hybinp), intent(in) :: hybinp
type(t_kpts), intent(in) :: kpts
type(t_sym), intent(in) :: sym
type(t_noco), intent(in) :: noco
......
......@@ -564,7 +564,7 @@ SUBROUTINE rdmft(eig_id,mpi,input,kpts,banddos,sliceplot,cell,atoms,enpara,stars
CALL read_eig(eig_id,ikpt,jspin,list=[(i,i=1,hybdat%nbands(ikpt))],neig=nbands,zmat=zMat)
! CALL read_z(atoms, cell, mpdata, hybdat, hybinp, kpts, sym, noco, input, ikpt, jsp, zMat)
! call read_z(atoms, cell, hybdat, kpts, sym, noco, input, ikpt, jsp, zMat)
zMat%matsize2 = hybdat%nbands(ikpt) ! reduce "visible matsize" for the following computations
CALL olap%multiply(zMat,trafo)
......
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