Thank you so much Sir, appreciate your help.
Regards
Shahin
On Mon, May 24, 2021 at 5:46 PM Wolfgang Bangerth <bange...@colostate.edu>
wrote:
> On 5/21/21 1:52 PM, Shahin Heydari wrote:
> > m_i = sigma (m_ij) = int(phi_i dx)
>
> Right. You compute the integral as a sum over the integrals on each cell.
> So
> when you do
>
> for (cell = ...)
> for (i=...)
> {
> lump = 0;
> for (q=...)
> lump += ...
>
> lump_cell_matrix(i,i) = lump;
> }
>
> then this is correct if you then add lump_cell_matrix(i,i) to your global
> matrix.
>
> You mention in your original question that you are unsure about what to do
> with Q2 (or higher) elements. This is a difficult question, but there is a
> good amount of literature on the topic. Just search for "mass lumping
> higher
> order elements".
>
> 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/1b0908b3-e43c-2d6a-b2a0-09a7a51b4636%40colostate.edu
> .
>
--
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/CALUbCFwyACZFLqF2ZEXztJ3H4b2iEhcVGYJ90beAMNPnUVz6aQ%40mail.gmail.com.
