Use symmetric solver from SciPy