On 10/24/22 13:16, Мария Бронзова wrote:
*** Caution: EXTERNAL Sender ***
Dear all,
I have been trying to implement orthotropic elasticity (similar as in
step-18) and in order to verify my results I decided to set all E and G
modulus values as in an isotropic case and to compare the outcome with
the one for the same isotropic elastic problem but with another
formulation (as in step-46 for the structure part).
I created the simplest 1 cell models and compared the matrices in both
cases (both LHS and RHS), they are absolutely identical. But weirdly
enough I receive different solution vectors.
Maria:
That is of course a contradiction. If you have the same linear system,
and you solve it the same way, you should get the same solution.
We are generally happy to help with advice on this forum, but we can't
debug other people's code. But I think you have everything you need to
debug your codes yourself. Just output the elements of the matrix and
rhs of the linear system for both codes. Make sure they are the exact
same. If they are not, then you know that the problem lies with the
assembly. If they are, then solve the linear system with the exact same
solver and output the solution element by element. It should be the
exact same. If it is not, you know that the problem lies with what must
be different ways you are solving the linear system. If it is, then the
problem lies in the part of the code that does something with the
solution later on.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups "deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/dad3cb12-9e95-33a3-1b1c-782ad6d68c29%40colostate.edu.
