gmres - restart and Gauss-Seidel
Lukasz Kaczmarczyk
lukasz at civil.gla.ac.uk
Fri Aug 15 17:40:12 CDT 2008
Barry Smith wrote:
> 9 KSP Residual norm 3.942287795560e-03
> 10 KSP Residual norm 4.611215310185e+00
> ^^^^^^^^^^^^^^^^^^^^^^
>
> This should never happen. This means there is an error
> in your preconditioner, your smoother is not a linear operator
> on its input variables. With GMRES, at restart the residual
> norm should hardly change at all. The fact that it shoots
> up means that the norm indicated in iteration 9 is totally wrong.
>
> If you run the gmres version with -ksp_monitor_true_residual
> you will see that the true residual is not actually decreasing
> like you think it is.
>
Thank You for answer, result for which you asking are in my previous
email - look at end of the previous email. If You look there at true
residual, You will see that results are ok.
Regards,
Lukasz
2) -ksp_type fgmres -ksp_gmres_restart 100 -ksp_monitor_true_residual
0 KSP preconditioned resid norm 8.770931682456e+05 true resid norm
8.770931682456e+05 ||Ae||/||Ax|| 2.025559907164e+06
1 KSP preconditioned resid norm 1.466659768959e+00 true resid norm
1.466659768959e+00 ||Ae||/||Ax|| 3.387105649673e+00
2 KSP preconditioned resid norm 2.135783619677e-01 true resid norm
2.135783619678e-01 ||Ae||/||Ax|| 4.932380991006e-01
3 KSP preconditioned resid norm 1.342354228152e-01 true resid norm
1.342354228151e-01 ||Ae||/||Ax|| 3.100034299884e-01
4 KSP preconditioned resid norm 7.872906297489e-02 true resid norm
7.872906297491e-02 ||Ae||/||Ax|| 1.818169828064e-01
5 KSP preconditioned resid norm 3.414496482864e-02 true resid norm
3.414496482865e-02 ||Ae||/||Ax|| 7.885441854116e-02
6 KSP preconditioned resid norm 2.025183908090e-02 true resid norm
2.025183908092e-02 ||Ae||/||Ax|| 4.676961897981e-02
7 KSP preconditioned resid norm 1.023707179274e-02 true resid norm
1.023707179279e-02 ||Ae||/||Ax|| 2.364150462111e-02
8 KSP preconditioned resid norm 4.938281368004e-03 true resid norm
4.938281368031e-03 ||Ae||/||Ax|| 1.140447230867e-02
9 KSP preconditioned resid norm 2.373276511281e-03 true resid norm
2.373276511245e-03 ||Ae||/||Ax|| 5.480847330515e-03
10 KSP preconditioned resid norm 1.180493643594e-03 true resid norm
1.180493643635e-03 ||Ae||/||Ax|| 2.726233291716e-03
11 KSP preconditioned resid norm 7.142592937639e-04 true resid norm
7.142592936809e-04 ||Ae||/||Ax|| 1.649511181911e-03
12 KSP preconditioned resid norm 4.226036746778e-04 true resid norm
3) -ksp_type fgmres -ksp_gmres_restart 10 -ksp_monitor_true_residual
0 KSP preconditioned resid norm 8.770931682456e+05 true resid norm
8.770931682456e+05 ||Ae||/||Ax|| 2.025559907164e+06
1 KSP preconditioned resid norm 1.466659768959e+00 true resid norm
1.466659768959e+00 ||Ae||/||Ax|| 3.387105649673e+00
2 KSP preconditioned resid norm 2.135783619677e-01 true resid norm
2.135783619678e-01 ||Ae||/||Ax|| 4.932380991006e-01
3 KSP preconditioned resid norm 1.342354228152e-01 true resid norm
1.342354228151e-01 ||Ae||/||Ax|| 3.100034299884e-01
4 KSP preconditioned resid norm 7.872906297489e-02 true resid norm
7.872906297491e-02 ||Ae||/||Ax|| 1.818169828064e-01
5 KSP preconditioned resid norm 3.414496482864e-02 true resid norm
3.414496482865e-02 ||Ae||/||Ax|| 7.885441854116e-02
6 KSP preconditioned resid norm 2.025183908090e-02 true resid norm
2.025183908092e-02 ||Ae||/||Ax|| 4.676961897981e-02
7 KSP preconditioned resid norm 1.023707179274e-02 true resid norm
1.023707179279e-02 ||Ae||/||Ax|| 2.364150462111e-02
8 KSP preconditioned resid norm 4.938281368004e-03 true resid norm
4.938281368031e-03 ||Ae||/||Ax|| 1.140447230867e-02
9 KSP preconditioned resid norm 2.373276511281e-03 true resid norm
2.373276511245e-03 ||Ae||/||Ax|| 5.480847330515e-03
10 KSP preconditioned resid norm 1.180493643635e-03 true resid norm
1.180493643635e-03 ||Ae||/||Ax|| 2.726233291716e-03
11 KSP preconditioned resid norm 8.973593571160e-04 true resid norm
8.973593570913e-04 ||Ae||/||Ax|| 2.072362665506e-03
12 KSP preconditioned resid norm 8.954054180580e-04 true resid norm
More information about the petsc-users
mailing list