We should look into using the LinearSolve.jl API in the FEASolver instead of implementing our own API for A \ b.