You are not providing full information to determine the problem. Why are there three linear solves? You never said anything about three linear solves. Why are you printing the solution after the first linear solve? Likely that "first" linear solve where you got CONVERGED_ATOL the right hand side is never set or it is set to zero.
Barry On Nov 27, 2010, at 10:02 PM, khalid ashraf wrote: > 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. > >
