PCNN preconditioner

[Andreas Grassl]

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