Hello,
Wolfgang Bangerth wrote:
AMESOS ERROR -22, ../../../../packages/amesos/src/Amesos_Klu.cpp, line
656
AMESOS ERROR -22, ../../../../packages/amesos/src/Amesos_Klu.cpp, line
696
2.175e-01 0000 0.00e+00
.
I think I am creating a matrix and a right hand side that are not the
actual ones I am looking for since I want the Jacobian matrix of the
residual with respect to velocity vector.
Yes, it looks like you're doing something wrong. I think you'll have to
find out for yourself what is going wrong, but if you take a look at the
Trilinos sources, you'll recognize that error -22 means that the matrix is
singular. I would also be suspicious to see that the residual is 0.00+e00
after 0 iterations: it may be that either the matrix is empty or the right
hand side vector is empty or both. You may want to check this.
W.
I checked everything and I found some mistakes. I changed what was badly
coded. And now, using the direct solver Amesos and following some points
of tutorial step-33, I obtain:
NonLin Res Lin Iter Lin Res
_____________________________________
8.103e-02 0000 0.00e+00
1.294e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
1.280e-01 0000 0.00e+00
.
.
.
So, it doesn't converge. Perhaps, you would have more ideas where the
problem could be.
Thanks in advance.
Best.
Isa
_______________________________________________
