Prototype of (potentially faster) version of SAFT using SciPy sparse matrix multiplication