=
(imag * scalar_product(J_a, phi_i)) *
fe_values.JxW(q_point);
cell_rhs(i) += rhs_value.real();
Thanks & Regards,
Pratyusa Kar
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?h
and proper refined mesh.
Thank you.
Regards,
Pratyusa
On Sun, Dec 15, 2024, 17:41 Wolfgang Bangerth
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
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.
Thank you.
On Saturday, December 14, 2024 at 6:41:36 PM UTC-6 Pratyusa Kar wrote:
> Hi,
>
> I have a simple 2D linear elasticit