first draft of code c1prealloc.c which will show Mat-directed preallocation; then c1poisson.c can emphasize KSP-level discussion