dpstf and dpstrf sometimes give different pivot results. Of course I can read the source code, but I don't get the idea from it. Since pivoting is for stability of the Cholesky decomposition, one of them is more stable. Which and why does the other exist?
What is the difference in the pivoting strategies between LAPACK's dpstrf and dpst2 and why?
linear algebra;lapack