Comment régler la tolérance de UMFPACK

Question

J'utilise umfpack en C++ pour résoudre une matrice clairsemée. J'ai du mal à trouver où le paramètre [UMFPACK PIVOT TOLERANCE] est défini. Quelqu'un sait-il comment définir ce paramètre de contrôle?

Mon code actuel est

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

Answer 1

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