C/Visual C - Factorizacion Cholesky

Hola!
Necesito la factorizacion de Cholesky implementada en lenguaje C, agradeceria me ayudaseis a encontrarla.
De todas maneras, escrita en Matlab es como sigue:
n = length(A);
for k=1:n
A(k, k) = sqrt(A(k, k) - sum(A(k, 1:k-1).^2));
for i = k+1:n
A(i, k) = (A(i, k) - sum(A(i, 1:k-1) .* A(k, 1:k-1))) / A(k, k);
end
end

Alguien se atreve a pasarla a C??
Muchas gracias por adelantado.
Pedro

Mh, aunque no debiera ser demasiado complicado, no he programado en Matlab (bueno, muy poquito) y no entiendo muy bien qué es A, la sintaxis de los paréntesis que usas para A (supongo que será un vector o matriz), y otras cosillas. O también, la descripción de la factorización de Cholesky y si existe algún método numérico para calcularla (que probablemente sería más fácil de implementar en un lenguaje de programación).

Busca en Google, o directamente aquí:

http://www.saiensu.co.jp/book_support/4-7819-0855-1/cholesky.txt
http://www-b2.is.tokushima-u.ac.jp/~ikeda/num/cholesky.c
http://perso.wanadoo.fr/jean-pierre.moreau/Cplus/choles_cpp.txt

Saludos,
José Jorge (Geo).