I guess, here the Poisson matrix is solved. Maybe it'll be more convenient to move it into a separate function?