Symmetric Positive Definite Skyline Matrix Factorization

Has anybody implemented one of these that they would like to share?

Somebody is going to implement (and then parallelize) an old Fortran code I have.I was looking for something in Chapel to give that person a second take on the problem. Even serial is good.

Thanks