Grids
calculate and LU factorise the matrix multiplying the electrostatic potential in quasineutrality this involves the factor 1-Gamma_0(kperp(alpha))
if using a modified Boltzmann response for the electrons obtain the response of phi_homogeneous to a unit perturbation in flux-surface-averaged phi