gmres - restart and Gauss-Seidel

Fri Aug 15 21:41:36 CDT 2008

    I don't think this really solves the problem, I think it merely  
hides the problem.

    Barry

On Aug 15, 2008, at 6:45 AM, Lukasz Kaczmarczyk wrote:

> Thank for all, problem is solved, I add PCJACOBI,
> PCSetType(Composite_prec,PCCOMPOSITE);
> PCCompositeAddPC(Composite_prec,PCJACOBI);
> PCCompositeAddPC(Composite_prec,PCSHELL); // Gauss-Seidel COLOR  
> Forward
> PCCompositeAddPC(Composite_prec,PCSHELL); // MulltiGrid
> PCCompositeAddPC(Composite_prec,PCSHELL); // Gauss-Seidel COLOR  
> Forward
> PCCompositeAddPC(Composite_prec,PCJACOBI);
>
> Now it works :)
> -ksp_type fgmres -ksp_gmres_restart 10  -ksp_monitor_true_residual
>  0 KSP Residual norm 8.770931682456e+05
>  0 KSP preconditioned resid norm 8.770931682456e+05 true resid norm
> 8.770931682456e+05 ||Ae||/||Ax|| 2.025559907164e+06
>  1 KSP Residual norm 2.470886778565e-01
>  1 KSP preconditioned resid norm 2.470886778565e-01 true resid norm
> 2.470886778563e-01 ||Ae||/||Ax|| 5.706268586963e-01
>  2 KSP Residual norm 6.271486855102e-02
>  2 KSP preconditioned resid norm 6.271486855102e-02 true resid norm
> 6.271486855100e-02 ||Ae||/||Ax|| 1.448337849605e-01
>  3 KSP Residual norm 1.782615793912e-02
>  3 KSP preconditioned resid norm 1.782615793912e-02 true resid norm
> 1.782615793917e-02 ||Ae||/||Ax|| 4.116774833918e-02
>  4 KSP Residual norm 6.410478294417e-03
>  4 KSP preconditioned resid norm 6.410478294417e-03 true resid norm
> 6.410478294455e-03 ||Ae||/||Ax|| 1.480436547575e-02
>  5 KSP Residual norm 1.314224590977e-03
>  5 KSP preconditioned resid norm 1.314224590977e-03 true resid norm
> 1.314224590937e-03 ||Ae||/||Ax|| 3.035071685413e-03
>  6 KSP Residual norm 2.683633682600e-04
>  6 KSP preconditioned resid norm 2.683633682600e-04 true resid norm
> 2.683633681768e-04 ||Ae||/||Ax|| 6.197586514302e-04
>  7 KSP Residual norm 8.799918944041e-05
>  7 KSP preconditioned resid norm 8.799918944041e-05 true resid norm
> 8.799918940855e-05 ||Ae||/||Ax|| 2.032254227740e-04
>  8 KSP Residual norm 2.508995399544e-05
>  8 KSP preconditioned resid norm 2.508995399544e-05 true resid norm
> 2.508995405419e-05 ||Ae||/||Ax|| 5.794276690857e-05
>  9 KSP Residual norm 1.026378586306e-05
>  9 KSP preconditioned resid norm 1.026378586306e-05 true resid norm
> 1.026378592108e-05 ||Ae||/||Ax|| 2.370319825777e-05
> 10 KSP Residual norm 3.126671988424e-06
> 10 KSP preconditioned resid norm 3.126671988424e-06 true resid norm
> 3.126671988424e-06 ||Ae||/||Ax|| 7.220739656737e-06
> 11 KSP Residual norm 1.229059204228e-06
> 11 KSP preconditioned resid norm 1.229059204228e-06 true resid norm
> 1.229059216105e-06 ||Ae||/||Ax|| 2.838390677073e-06
> 12 KSP Residual norm 5.416643445915e-07
> 12 KSP preconditioned resid norm 5.416643445915e-07 true resid norm
> 5.416643499083e-07 ||Ae||/||Ax|| 1.250920232920e-06
> 13 KSP Residual norm 2.633538639496e-07
> 13 KSP preconditioned resid norm 2.633538639496e-07 true resid norm
> 2.633538571527e-07 ||Ae||/||Ax|| 6.081896812770e-07
> 14 KSP Residual norm 1.027652175716e-07
> 14 KSP preconditioned resid norm 1.027652175716e-07 true resid norm
> 1.027652231589e-07 ||Ae||/||Ax|| 2.373261170164e-07
> 15 KSP Residual norm 3.189221996600e-08
> 15 KSP preconditioned resid norm 3.189221996600e-08 true resid norm
> 3.189223368074e-08 ||Ae||/||Ax|| 7.365195880253e-08
> 16 KSP Residual norm 1.202828530862e-08
> 16 KSP preconditioned resid norm 1.202828530862e-08 true resid norm
> 1.202828930905e-08 ||Ae||/||Ax|| 2.777814428189e-08
> 17 KSP Residual norm 2.701944514090e-09
> 17 KSP preconditioned resid norm 2.701944514090e-09 true resid norm
> 2.701946176761e-09 ||Ae||/||Ax|| 6.239877409956e-09
> 18 KSP Residual norm 8.920433449729e-10
> 18 KSP preconditioned resid norm 8.920433449729e-10 true resid norm
> 8.920273046556e-10 ||Ae||/||Ax|| 2.060048817870e-09
> 19 KSP Residual norm 2.868707402558e-10
> 19 KSP preconditioned resid norm 2.868707402558e-10 true resid norm
> 2.868717856448e-10 ||Ae||/||Ax|| 6.625020106597e-10
> 20 KSP Residual norm 5.959506466900e-11
> 20 KSP preconditioned resid norm 5.959506466900e-11 true resid norm
> 5.959506466900e-11 ||Ae||/||Ax|| 1.376289065161e-10
> 21 KSP Residual norm 2.347311204935e-11
> 21 KSP preconditioned resid norm 2.347311204935e-11 true resid norm
> 2.346228685060e-11 ||Ae||/||Ax|| 5.418383051598e-11
> 22 KSP Residual norm 9.323415814482e-12
> 22 KSP preconditioned resid norm 9.323415814482e-12 true resid norm
> 9.323298967857e-12 ||Ae||/||Ax|| 2.153123667531e-11
> 23 KSP Residual norm 2.925536240498e-12
> 23 KSP preconditioned resid norm 2.925536240498e-12 true resid norm
> 2.974960765051e-12 ||Ae||/||Ax|| 6.870377594124e-12
> 24 KSP Residual norm 1.025371996314e-12
> 24 KSP preconditioned resid norm 1.025371996314e-12 true resid norm
> 1.141211075669e-12 ||Ae||/||Ax|| 2.635514086960e-12
> 25 KSP Residual norm 3.207279162194e-13
> 25 KSP preconditioned resid norm 3.207279162194e-13 true resid norm
> 6.238776238382e-13 ||Ae||/||Ax|| 1.440783656257e-12
> 26 KSP Residual norm 1.290295601110e-13
> 26 KSP preconditioned resid norm 1.290295601110e-13 true resid norm
> 4.880917397572e-13 ||Ae||/||Ax|| 1.127199589352e-12
> 27 KSP Residual norm 3.648895071464e-14
> 27 KSP preconditioned resid norm 3.648895071464e-14 true resid norm
> 4.915218264645e-13 ||Ae||/||Ax|| 1.135121035287e-12
> 28 KSP Residual norm 9.038317929266e-15
> 28 KSP preconditioned resid norm 9.038317929266e-15 true resid norm
> 4.525760014181e-13 ||Ae||/||Ax|| 1.045179504990e-12
> 29 KSP Residual norm 3.168044355090e-15
> 29 KSP preconditioned resid norm 3.168044355090e-15 true resid norm
> 4.800332481647e-13 ||Ae||/||Ax|| 1.108589300191e-12
> 30 KSP Residual norm 4.793973611579e-13
> 30 KSP preconditioned resid norm 4.793973611579e-13 true resid norm
> 4.793973611579e-13 ||Ae||/||Ax|| 1.107120782053e-12
> 31 KSP Residual norm 1.622163699810e-13
> 31 KSP preconditioned resid norm 1.622163699810e-13 true resid norm
> 5.177440837392e-13 ||Ae||/||Ax|| 1.195678744473e-12
> 32 KSP Residual norm 3.278455971610e-14
> 32 KSP preconditioned resid norm 3.278455971610e-14 true resid norm
> 5.095690679196e-13 ||Ae||/||Ax|| 1.176799354136e-12
> 33 KSP Residual norm 9.985582557377e-15
> 33 KSP preconditioned resid norm 9.985582557377e-15 true resid norm
> 5.216297041713e-13 ||Ae||/||Ax|| 1.204652200482e-12
> 34 KSP Residual norm 4.197972975568e-15
> 34 KSP preconditioned resid norm 4.197972975568e-15 true resid norm
> 5.581506198985e-13 ||Ae||/||Ax|| 1.288993642587e-12
> 35 KSP Residual norm 1.034379105713e-15
>