Barry Smith schrieb: > > For the subdomain solves the prefix are is_localD_ and is_localN_ > so you should use the options > -is_localD_pc_factor_shift_positive_definite and > -is_localN_pc_factor_shift_positive_definite
With both options it is working now. > > There is currently no subroutine that "pulls out" the inner KSP's for > the > Neuman and Dirichlet problems for use in the code; though there should be > PCISGetDPC() and PCISGetNPC() that would get the pointer to ksp_N and ksp_D > objects inside the PC_IS data structured defined in > src/ksp/pc/impls/is/pcis.h > You can easily add these routines. Then use them to get the inner PC and > set > the shift option (and anything else you want to set). I'll try later, for now I'm happy with the options Thank you for helping, I'll give some feedback cheers ando -- /"\ Grassl Andreas \ / ASCII Ribbon Campaign Uni Innsbruck Institut f. Mathematik X against HTML email Technikerstr. 13 Zi 709 / \ +43 (0)512 507 6091
