it is still not very clear to me why I am not getting the same mass
matrix when doing this two different approaches when, as far as I
understand should be the same. In order to get a 2x2 block matrix
containing the mass matrix for the top left and bottom right I do the
following:
/Mass_matrix_tmp(i,j)+=((fe_values[real_component].value(i, q_point) *
fe_values[real_component].value(j, q_point)) * fe_values.JxW(q_point));/
/Mass_matrix_tmp(i,j)+=((fe_values[imag_component].value(i, q_point) *
fe_values[imag_component].value(j, q_point)) * fe_values.JxW(q_point));/
which gives different results compared to:
/if (fe.system_to_component_index(i).first ==
fe.system_to_component_index(j).first)/
/Mass_matrix_tmp(i,j)+= ((fe_values.shape_value(i, q_point) *
fe_values.shape_value(j, q_point)) * fe_values.JxW(q_point));/
/
/
Could you please clarify why this two approaches do not give the same
results?
It looks to me like these two methods should give the same matrix. In
which way are they different?
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/07b37563-5a0f-9fa7-f7c1-00e0a4a4243a%40colostate.edu.
