Replace the existing discretization routines with their matrix-inverse-free counterparts, which are faster and more stable.