Starting to build cusp based sparse matrix solver for 3d poisson
Moved poisson build code into sperate function
Converted to template
Added test problem
Moved cusp build code into header. Cleaned up Makefile