[deal.II] The shape function on physical cell in hexahedral is not parallel with the edge (FeNedelecSZ and FENedelec)

Hi colleagues, When using FE_NedelecSZ in my code. Testing for an edge whose global vertices located from (0,0,0) to (0,0,1) in real coordinates. With a cube I get the shape function vectors at the dof related to the edge, for examples, (0,0,0), (0,0,-0.25), (0,0,-0.5), (0,0,-1), which are

Re: [deal.II] Mass matrix for a distributed vector problem

Dear Luca, Thanks for your suggestion. Unfortunately, it did not solve the problem. I am sending the modified version, as well as the output of the program. Out of curiosity, what is the reason to use (2*fe_degree + 1)? Checking step-8, I notice that there a quadrature degree one larger than the

Re: [deal.II] Mass matrix for a distributed vector problem

You are integrating using two quadrature points per direction. Can you raise that to (2*fe.degree+1)? L. > On 11 Apr 2019, at 10:58, bobspar...@gmail.com wrote: > > Dear all, > > First of all, thanks for the awesome library. I am starting to learn deal.II, > and it is great! The

[deal.II] Mass matrix for a distributed vector problem

Dear all, First of all, thanks for the awesome library. I am starting to learn deal.II, and it is great! The documentation has been extremely helpful so far. I am trying to solve the time dependent elastic equation, and there we naturally have a mass matrix-like for the velocities. To build