<josef.pktd at gmail.com> wrote: > It's in LINPACK but never made it into LAPACK. > Here is a cython version https://github.com/jcrudy/choldate The inner loop can be replaced with two calls to BLAS' *axpy subroutines. Indexing with casts to unsigned int just corrupts the readability. Sturla