c4poisson.c has important change from c2poisson: KSPSetOperators() --> KSPSetComputeOperators() so that we can use multigrid