如何设置UMFPACK的容差

时间:2011-07-04 00:31:14

标签: c++ linear-algebra numerical-methods umfpack

我在c ++中使用umfpack来解决稀疏矩阵。我无法找到[UMFPACK PIVOT TOLERANCE]参数的设置位置。有人知道如何设置这个控制参数吗?

我目前的代码是

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);

1 个答案:

答案 0 :(得分:0)

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