Rust, Design Patterns, Matrix Factorization, Numerical Methods, Software Engineering.
implements a trait PositiveDefinite with a method for computing the Cholesky decomposition. Along the way, multiple checks are performed to ensure that the computation goes smoothly: ensuring the matrix is square — we did not show it here in the article for brevity but in the code you can see it.