On 12/15/24 08:40, Pratyusa Kar wrote:
Also I checked with a 4 element mesh and 1 quadrature point. However, I see a
noninvertible mass matrix. I do not understand what is missing in the code.
Pratyusa:
OK, so you already determined that the matrix is not invertible. The question
is why. What happens if you have only one cell? Or you look at the cell matrix
of a single cell?
I don't understand what you mean by "and 1 quadrature point". I checked with
your code, and if I understand correctly, you are using a Q1 element to
project the strain onto, and that you are using a 2x2 Gauss formula. That
seems appropriate.
I looked at your code, and I do not understand why fe_L2 is vector valued. You
are trying to project individual components of the strain onto a finite
element. The individual components are scalars. Why is fe_L2 not also a scalar?
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 visit
https://groups.google.com/d/msgid/dealii/b4075a0d-7dd5-4253-a695-0d3b104a7d1f%40colostate.edu.
