开发者

How to set UMFPACK's tolerance

开发者 https://www.devze.com 2023-03-16 21:21 出处:网络
I am using umfpack in c++ to solve a sparse matrix. I am having trouble finding wher开发者_运维问答e the [UMFPACK PIVOT TOLERANCE] parameter is set. Would anybody know how to set this control paramete

I am using umfpack in c++ to solve a sparse matrix. I am having trouble finding wher开发者_运维问答e the [UMFPACK PIVOT TOLERANCE] parameter is set. Would anybody know how to set this control parameter?

My current code is

umf::symbolic_type<double> Symbolic;
umf::numeric_type<double> Numeric;
umf::symbolic(J, Symbolic);
umf::numeric(J, Symbolic, Numeric);
umf::solve(J, unknowns, knowns, Numeric);


umf::control_type<ublas::compressed_matrix<double>::value_type > myControl;
myControl.ptr[UMFPACK_SYM_PIVOT_TOLERANCE]=0;
myControl.ptr[UMFPACK_PIVOT_TOLERANCE]=0;
0

精彩评论

暂无评论...
验证码 换一张
取 消