Also another thing which gives me a headache is that why is the jacobian value
twice the value it should be. For a rectangular bilinear lagrange element we
can compute simply J by a/2 and b/2, but we would receive 1 and 1 for the
above case where a=b=2.
This doesn't correspond to what deal.II outputs.
It's not clear to me what you mean by a and b, and "1 and 1". Can you say what
you expect the Jacobian matrix to be, and what deal.II outputs?
Additionally, using
std::vector<DerivativeForm<1, dim, dim> > J = fe_values.get_jacobians ();
gives a vector containing n x n matrices?
Correct.
Since I can just output J[0][i][j], J[1][i][j] or J[10][i][j]. All work, but
why? Shouldn't it be somehow limited? (Stupid question maybe...)
What do you mean by "limited"? And do the J[0][i][j] match what you expect
them to be?
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].
For more options, visit https://groups.google.com/d/optout.
