get_phi_ffs Subroutine

private subroutine get_phi_ffs(rhs, phi)

change from rhs defined on grid with ky >=0 and kx from 0,...,kxmax,-kxmax,...,-dkx to rhs_swap defined on grid with ky = -kymax,...,kymax and kx >= 0 solve sum_s Z_s int d^3v = gam0*phi where sum_s Z_s int d^3v is initially passed in as rhs_swap and then rhs_swap is over-written with the solution to the linear system

swap back from the ordered grid in ky to the original (kx,ky) grid

Arguments

Type IntentOptional Attributes Name
complex, intent(in), dimension(:, :, -nzgrid:) :: rhs
complex, intent(out), dimension(:, :, -nzgrid:) :: phi