61 INTEGER,
INTENT(IN) :: idir
62 REAL(
dp),
INTENT(IN) :: my_chi
64 INTEGER :: handle, ig, ng
68 CHARACTER(len=*),
PARAMETER :: routinen =
'mult_G_ov_G2_grid'
70 CALL timeset(routinen, handle)
72 CALL pw_pool%create_pw(influence_fn)
74 grid => influence_fn%pw_grid
75 DO ig = grid%first_gne0, grid%ngpts_cut_local
77 influence_fn%array(ig) =
gaussi*grid%g(idir, ig)/g2
79 IF (grid%have_g0) influence_fn%array(1) = 0.0_dp
84 funcg_times_rho%array(1:ng) = funcg_times_rho%array(1:ng)*influence_fn%array(1:ng)
85 IF (grid%have_g0) funcg_times_rho%array(1) = my_chi
87 CALL pw_pool%give_back_pw(influence_fn)
subroutine, public mult_g_ov_g2_grid(pw_pool, rho_gspace, funcg_times_rho, idir, my_chi)
Given the current density on the PW grid in reciprcal space (obtained by FFT), calculate the integral...