Computational Thinking, Python, Linear Algebra, Linear Regression, Jupyter Notebook.
params, the function didn’t return what we thought, or sometimes the docs are even outdated) np.linalg.norm(R.T @ w - Atb) 1.1368683772161603e-13 coeffs_chol = scipy.linalg.solve_triangular(R, w, lower=False) np.linalg.norm(R @ coeffs_chol - w) 6.861429794408013e-14 def ls_chol(A, b): R =