Use eigh for bmatrix

Assert that B-matrix is symmetrical and use eigh to make sure we get real values.