Exemplo n.º 1
0
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;
}