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
paral
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
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 documentation
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 th