Hi Dr. Bangerth, Thank you for the explanation.
Yes I just did the problem with 1 element and figured out that because I was using vector to output the strain the errors come up. I should have checked this before. Now with scalar definition my problem is converging even with multiple elements and proper refined mesh. Thank you. Regards, Pratyusa On Sun, Dec 15, 2024, 17:41 Wolfgang Bangerth <bange...@colostate.edu> wrote: > 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 > . > -- 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/CAG0QpcbMaMA6aye13YDeL%3D5CTcAajFkN8pYesD5qxwBySi28yA%40mail.gmail.com.
