inline bool ls_solve(const mat &A, const vec &b, vec &result){ //result = A.jacobiSvd(ComputeThinU | ComputeThinV).solve(b); result = A.ldlt().solve(b); return true; }