some comments in poisson.c