[petsc-users] fieldsplit doesn't pass prefix to inner ksp
anton
popov at uni-mainz.de
Fri Sep 26 07:29:40 CDT 2014
Create preconditioner:
PCCreate(PETSC_COMM_WORLD, &pc);
PCSetOptionsPrefix(pc, "bf_");
PCSetFromOptions(pc);
Define fieldsplit options:
-bf_pc_type fieldsplit
-bf_pc_fieldsplit_type SCHUR
-bf_pc_fieldsplit_schur_factorization_type UPPER
Works OK.
Set options for the first field solver:
-bf_fieldsplit_0_ksp_type preonly
-bf_fieldsplit_0_pc_type lu
Doesn't work (ignored), because "bf_" prefix isn't pass to inner solver
ksp (checked in the debugger).
Indeed, the following works:
-fieldsplit_0_ksp_type preonly
-fieldsplit_0_pc_type lu
Observed with 3.5 but not with 3.4
Thanks.
Anton
More information about the petsc-users
mailing list