Commit ba22eb49 authored by Matteo Gerosa's avatar Matteo Gerosa
Browse files

Substituted variable fftdriver with 'Wave', only allowed value when using q-points.

parent f7b83c9a
......@@ -566,18 +566,18 @@ SUBROUTINE solve_gfreq_k(l_read_restart)
!
! Bring it to R-space
IF(noncolin) THEN
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,TRIM(fftdriver),igq_q(1,iq))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,'Wave',igq_q(1,iq))
DO ir=1,dffts%nnr
pertr(ir)=DCONJG(phase(ir))*psick_nc(ir,1)*DCONJG(pertr(ir))
ENDDO
CALL single_fwfft_k(dffts,npw,npwx,pertr,dvpsi(1,ip),'Wave',igk_k(1,current_k))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,TRIM(fftdriver),igq_q(1,iq))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,'Wave',igq_q(1,iq))
DO ir=1,dffts%nnr
pertr(ir)=DCONJG(phase(ir))*psick_nc(ir,2)*DCONJG(pertr(ir))
ENDDO
CALL single_fwfft_k(dffts,npw,npwx,pertr,dvpsi(1+npwx,ip),'Wave',igk_k(1,current_k))
ELSE
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,TRIM(fftdriver),igq_q(1,iq))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,'Wave',igq_q(1,iq))
DO ir=1,dffts%nnr
pertr(ir)=DCONJG(phase(ir))*psick(ir)*DCONJG(pertr(ir))
ENDDO
......
......@@ -950,7 +950,7 @@ SUBROUTINE solve_wfreq_k(l_read_restart,l_generate_plot)
pertr(ir)=phase(ir)*psick_nc(ir,1)*DCONJG(pertr(ir))
ENDDO
CALL single_fwfft_k(dffts,npw,npwx,pertr,dvpsi(1,ip),'Wave',igk_k(1,current_k))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,TRIM(fftdriver))
CALL single_invfft_k(dffts,npwq,npwqx,pertg(1),pertr,'Wave',igq_q(1,iq))
DO ir=1,dffts%nnr
pertr(ir)=phase(ir)*psick_nc(ir,2)*DCONJG(pertr(ir))
ENDDO
......
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