solve a linear system with a cached LU decomposition
github.com/substack/ndarray-linear-solve
substack/ndarray-linear-solve