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

[khalid ashraf]

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...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20101127/adbf9090/attachment.htm>