Dear Wolfgang,

Thanks for your reply.
My expectation stems from the fact that I am using the exact same mesh, but 
instead of using "GridGenerator" function of dealii, this time I import a 
mesh from Cubit. Then, I mark the boundaries in the same order: left: 1, 
right: 2, bottom: 3, up: 4. So, I think I did the exact same thing.

What happens is the order of DoF numbering changes (as shown in the 
pictures). But since my mesh is the same, and so are the boundary IDs, I 
expect the results to be as I described since 
"DoFTools::make_periodicity_constraints" 
function gets the boundary IDs as inputs.

Please let me know if I am still not clear.

Masoud

On Friday, 2 December 2022 at 23:43:08 UTC Wolfgang Bangerth wrote:

> On 12/2/22 10:20, Masoud Ahmadi wrote:
> > 
> > The problem is, when I import the mesh from Cubit (Lines 68-71 in the 
> uploaded 
> > code), I will have different numbering for the nodes (Case <II> in the 
> below 
> > image). While I expect the constraints to be as follows:
> > U4 = U6
> > U5 = U7
> > U0 = U2
> > U1 = U3
> > for the first function, and
> > U2 = U6
> > U3 = U7
> > they are the same as case <I>!
>
> Masoud -- you don't say in your question, but are the constraints wrong? 
> You 
> say that you *expect* them to be different, but not where your expectation 
> stems from.
>
> 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].
To view this discussion on the web visit 
https://groups.google.com/d/msgid/dealii/37a17f88-f95e-4fcc-8425-c8e2fa08ded9n%40googlegroups.com.

Reply via email to