Recursive Formulation of Cholesky Algorithm in Fortran 90

Fortran 90 allows writing recursive procedures (see [7]). Recursion leads to automatic variable blocking for dense linear-algebra algorithms (see [5, 8]). The recursive way of programming algorithms eliminate the use of BLAS level 2 in the factorization steps. For this and other reasons recursion usually speed up the algorithms.