num::cholesky_update

Defined in header <linear/factorization/cholesky.hpp>

void cholesky_update(cholesky_result &factor, const vec<real> &update);

Replace A=LL^T by A+x*x^T in O(n^2).