Deepak,
By the way, I am constructing the cell matrices using B'DB, where B is
the strain-displacement matrix and D is the constitutive relation. Then
I add it to the system matrix using distribute_local_to_gobal function.
Since B and D are the same for global and adaptive case, I do not expect
any difference in the cell matrices. The system_rhs is correct since I
have only a point load and I can check it by printing the vector. Thus,
the possible issue can be only in the assembly (for the local refinement
one). My problem is vector-values. In case this information can lead to
any other suggestion from someone, it would be great.
I don't have any experience with this approach, but at the end of the
day, you get a local matrix A_K = B_K' D_K B_K for each cell K, and
that should be all that matters. I assume you mean
ConstraintMatrix::distribute_local_to_gobal, not
cell->distribute_local_to_gobal, right?
Meanwhile I am going to extend the simple example for a basic version of
the elastic problem I am trying to solve. Hope then I can figure out the
error.
Good plan.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
