first proof of concept - self-contained ripley system matrix with diagonal storage, CG solver, no preconditioner, no MPI, no blocks, quite a few quick'n'dirty hacks. Solves poisson faster than paso :-)