symmetric reordering and incomplete factorization with tolerance?

Barry Smith bsmith at mcs.anl.gov
Tue Sep 16 10:37:00 CDT 2008


    You should use rcm+icc if you want to keep a symmetric  
preconditioner.

    Depending on your matrix you might want to use KSPCR or KSPMINRES  
or even KSPSYMMLQ
instead of bicgstab?

    We don't have a drop tolerance ICC and I do not recommend our drop  
tolerance ILU.

    Barry


On Sep 16, 2008, at 10:23 AM, zhifeng sheng wrote:

> Dear all
>
> I used the reordering scheme in petsc, and I would like to know  
> whether they are symmetric or not.
>
> I have a symmetric matrix (for 3D system) and I tried to solve it  
> with petsc and with matlab.
>
> In petsc I used rcm + ILU + bicgstab. while in matlab I used symrcm 
> +ICC+bicgstab.
>
> it seems that with symmetric reordering, the factorization is  
> somewhat faster... So I am wondering if it is possible to apply a  
> symmetric reordering?
>
> PS: in Matlab, an incomplete factorization can be done with a  
> tolerance (e.g. ICC(1e-4) ), can I do something like that with Petsc?
>
> Thanks a lot
> Best regards
> Zhifeng Sheng
>
>




More information about the petsc-users mailing list