Working on simplified Laplace solver

"cyclic" solver using the CyclicReduction class for the tridiagonal solve.
Once tested should be able to replace serial "tri" and parallel "spt" implementations.

Main aim is to simplify the structure of the code, to enable extensions.