using MUMPS with mg

Thomas Geenen geenen at gmail.com
Mon Jun 23 04:52:50 CDT 2008


dear Petsc users,

when using mg with Petsc (ml to be precise)
the default solver on the coarsest level is "redundant"
(I could not find much info in either the usermanual or de online
manuals about this solver)

it seems to be a sequential LU factorization running on all cpu's at
the same time?
doing the same thing?

i try to use MUMPS instead.
i try to invoke mumps by calling
ierr = PetscOptionsSetValue("-mg_coarse_ksp_type", "preonly");
ierr = PetscOptionsSetValue("-mg_coarse_mat_type", "aijmumps");
ierr = PetscOptionsSetValue("-mg_coarse_pc_type", "lu");

this gives me

[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: No support for this operation for this object type!
[0]PETSC ERROR: Matrix type mpiaij  symbolic

so apparently Petsc did not convert the matrix on the coarsest level.
if i remove the pc_type the conversion option is ignored and Petsc
uses redundant.



cheers
Thomas




More information about the petsc-users mailing list