by markusm73 » Thu Jan 08, 2009 7:10 pm
The pivoted Cholesky recently introduced in LAPACK 3.2 (SPSTRF) does not seem to offer support for early stopping other than by providing a tolerance value. Wouldn't it seem useful to also offer a maximum rank as stopping criterion to essentially perform incomplete cholesky factorization with diagonal pivoting? I guess it's too late to change the API, but maybe this kind of function could be considered as a future addition.