Closed lindonroberts closed 4 years ago
For interpolation linear solves, don't loop over m residuals (or n+1 Lagrange polynomials) - scipy.linalg.solve_triangular can handle multiple RHS (https://docs.scipy.org/doc/scipy/reference/generated/scipy.linalg.solve_triangular.html)
Should lead to a speed improvement.
For interpolation linear solves, don't loop over m residuals (or n+1 Lagrange polynomials) - scipy.linalg.solve_triangular can handle multiple RHS (https://docs.scipy.org/doc/scipy/reference/generated/scipy.linalg.solve_triangular.html)
Should lead to a speed improvement.