On 4/11/19 9:17 PM, Phạm Ngọc Kiên wrote:
> 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 parallel
> to the edge.
> However, if I create a mesh with GMSH and check the shape function vectors,
> some of them are: (-0.0868383, -0.428, 1), (-0.129933, -0.402772, 1),
> (-0.052103 -0.2568 0.6),..,
> which are not parallel to the edge.
>
> In my limited opinion, I will get the wrong solution as the shape functions
> created by my code are not parallel to the edge in my model.
> Could you please tell me why this happen and how to fix it?
Good question. This does sound wrong. To figure out what is happening exactly,
we typically construct small test cases that try to illustrate one particular
issue. In your case, I would try to create a mesh with exactly one cell --
either the one you have above, or maybe even simpler a cube that is just
rotated. Then create a DoFHandler on it and output the values of the shape
functions. This way, you really have only one thing that could go wrong, and
it's easy to demonstrate what you see and how it differs from expectations.
Do you think you could come up with a small program that does this?
Best
WB
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
