Adding Discrete Sine Transform option for Laplace solver

cyclic solver (type = cyclic) for Laplacian problems
can now use DSTs. This allows zero-value boundary conditions to be
set on the Z boundaries