Use the existing fn for Laplacian and check that matrix is symmetric.
Added some comments. Removed the Yen and Brand kernels which are only defined for symmetric matrices.