@Denis Davydov & Jean-Paul Pelteret:
Thanks. I will check it out.
@Prof. Bangerth:
For a rectangular bilinear lagrange element we define two parameters "a"
and "b". They define the width b and height a of the rectangle within
physical coordinates. The Jacobian matrix, B-operator, stiffness matrix
etc. can all be computed then using this special case formulas. 1 and 1
means a = b = 1, so we obtain a square with same width and height in an
amount of 1.
I expect the Jacobian matrix to be
b/2 0
0 a/2
For the previous cube where I used a=b=2, hence a square with width and
length 2, we should receive the following:
1 0
0 1
But if you compare this to my results above, I receive 2 instead of 1? Why?
What do you mean by "limited"? And do the J[0][i][j] match what you expect
> them to be?
Limited means, if you try to access data from a vector/matrix, either you
receive an error or just 0 since the matrix has only values where you store
them. In this case I can just put 10 as first index, it still gives me
values, e.g. -3.9853 something.
This makes no sense to me. For shape functions I checked. Only values from
n_q_points x dofs_per_cell are stored there, rest is 0. But here it seems
different. Or am I overseeing something?
Best regards,
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.
