tests on symmetric laplacian
some mpi stuff
found a way to template constructors
added bessel.h