# [petsc-users] Supplying the RHS Vector to the DMMG solver from outside

khalid ashraf khalid_eee at yahoo.com
Sat Nov 27 22:02:56 CST 2010

```Here is the output with opt 1)  VecSet(b,h):

0 KSP preconditioned resid norm 5.656854249492e+00 true resid norm
5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00
1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm
8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16
Linear solve converged due to CONVERGED_RTOL iterations
Process [0]
0.666667
0.666667
Process [1]
0.666667
0.666667
Process [2]
0.666667
0.666667
Process [3]
0.666667
0.666667
0 KSP preconditioned resid norm 5.656854249492e+00 true resid norm
5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00
1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm
8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16
Linear solve converged due to CONVERGED_RTOL iterations 1
0 KSP preconditioned resid norm 5.656854249492e+00 true resid norm
5.656854249492e+00 ||Ae||/||Ax|| 1.000000000000e+00
1 KSP preconditioned resid norm 1.256073966947e-15 true resid norm
8.881784197001e-16 ||Ae||/||Ax|| 1.570092458684e-16
Linear solve converged due to CONVERGED_RTOL iterations 1

Here is the output with opt 2) Copied from another function:

Linear solve converged due to CONVERGED_ATOL iterations 0
Process [0]
0
0
Process [1]
0
0
Process [2]
0
0
Process [3]
0
0
0 KSP preconditioned resid norm 7.515606252360e+01 true resid norm
7.515606252360e+01 ||Ae||/||Ax|| 1.000000000000e+00
1 KSP preconditioned resid norm 2.014114995652e+01 true resid norm
2.014114995652e+01 ||Ae||/||Ax|| 2.679910213523e-01
2 KSP preconditioned resid norm 6.102231836870e+00 true resid norm
6.102231836870e+00 ||Ae||/||Ax|| 8.119413966044e-02
3 KSP preconditioned resid norm 3.582243925876e-01 true resid norm
3.582243925876e-01 ||Ae||/||Ax|| 4.766407134156e-03
4 KSP preconditioned resid norm 9.891345584245e-02 true resid norm
9.891345584245e-02 ||Ae||/||Ax|| 1.316107477176e-03
5 KSP preconditioned resid norm 3.103414625576e-14 true resid norm
2.787743265007e-14 ||Ae||/||Ax|| 3.709272640689e-16
Linear solve converged due to CONVERGED_RTOL iterations 5
0 KSP preconditioned resid norm 7.567218624719e+01 true resid norm
7.567218624719e+01 ||Ae||/||Ax|| 1.000000000000e+00
1 KSP preconditioned resid norm 2.027946648035e+01 true resid norm
2.027946648035e+01 ||Ae||/||Ax|| 2.679910213523e-01
2 KSP preconditioned resid norm 6.144138058565e+00 true resid norm
6.144138058565e+00 ||Ae||/||Ax|| 8.119413966044e-02
3 KSP preconditioned resid norm 3.606844483858e-01 true resid norm
3.606844483858e-01 ||Ae||/||Ax|| 4.766407134156e-03
4 KSP preconditioned resid norm 9.959273013421e-02 true resid norm
9.959273013421e-02 ||Ae||/||Ax|| 1.316107477176e-03
5 KSP preconditioned resid norm 2.724152110884e-14 true resid norm
3.587687035382e-14 ||Ae||/||Ax|| 4.741090767039e-16
Linear solve converged due to CONVERGED_RTOL iterations 5

Here the printed zeros are the print of DMMGGetx()

Thanks.

-------------- next part --------------
An HTML attachment was scrubbed...