Fixed comments for poisson_solve
glue code can now take an arbitrary RHS vector to pass to SiGMA