[petsc-users] Computing determinant using MUMPS

[Hong Zhang]

Barbara :
I added mumps options ICNTL(30 - 33) to petsc-3.2.
With the latest petsc-3.2, you can get determinant with '-mat_mumps_icntl_33 1'
the result can be found using '-ksp_view'. For example,
petsc-3.2/src/ksp/ksp/examples/tutorials>./ex2 -pc_type lu
-pc_factor_mat_solver_package mumps -mat_mumps_icntl_33 1 -ksp_view
...
(RINFOG(12) RINFOG(13))*2^INFOG(34) (determinant): (0.900768,0)*(2^99)

Let us know if you need additional support - I guess you might want
MatGetMumpsRINFOG()?

Hong

>
> I am using version 3.2 of PETSC which uses MUMPS 4.10. When solving a system
> of linear equations with the option -mat_mumps_icntl_33 1 (= compute
> determinant) it is output: Option left: name:-mat_mumps_icntl_33 value: 1.
> Is this option not supported by PETSC?
>
> Thanks for your help,
> Barbara
>