Another approach I tried recently: FullMatrix<double> B (dofs_per_cell, dofs_per_cell);
fe.get_interpolation_matrix (fe, B); But it seems that I receive a matrix full of 1.0 on the diagonal only like the identity matrix. What exactly is the purpose of this interpolation matrix? I thought maybe it is the matrix I seek. They usually call it differential operator B(xi) or interpolation matrix. Then the stiffness matrix is constructed by means of: integral K = transpose(B)*C*B*det(J)*h dxi1 dxi2 I just need the B-Operator somehow. BR, S. A. Mohseni -- 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]. For more options, visit https://groups.google.com/d/optout.
