solves Delta * phi_hom = -delta_{ky,0} * ne/Te for phi_hom this is the vector describing the response of phi_hom to a unit impulse in phi_fsa it is the sum over ky and integral over kx of this that is needed, and this is stored in adiabatic_response_factor
Layouts Parameters Grids Calculations adiabatic_response_vector is initialised to be the rhs of the equation for the 'homogeneous' part of phi, with a unit impulse assumed for the flux-surface-averaged phi only the ky=0 component contributes to the flux-surface-averaged potential
pass in the rhs and overwrite with the solution for phi_homogeneous
obtain the flux surface average of the response vector