[MOAB-dev] Edges at the ends of a periodic mesh?

Robert Jacob jacob at mcs.anl.gov
Wed May 23 21:10:06 CDT 2012



On 5/23/12 12:34 PM, Iulian Grindeanu wrote:
>
>
> ------------------------------------------------------------------------
>
>
>     Hi,
>
>     If I have a 2D rectangular mesh which is periodic in one direction,
>     lets
>     say the j direction, then the leftmost and rightmost edge of the mesh
>     are really the same edge. Do they both need to be declared in MOAB? Or
>     does the periodic property mean you only need to describe one of the
>     edges and the other end will know its adjacent to it?
>
>     Thanks,
>     Rob
>
> Hi,
> I assume that you use ScdInterface;
> If the box is periodic in j direction, the nodes in that direction will
> be "periodic".
> Nodes are not duplicated, the edges that you get on the leftmost and
> rightmost sides are really pointing to the same nodes, they will not be
> duplicated, either.
>
> How do you get edges after you create quads?
> You can use "get_adjacency" methods, with "create_if_missing" flag true.
That's what we're going to do.

> Only one edge would be created.
> The problem is that the order in which these are created is not very
> "structured"
> Or do you create them with new scd boxes, of dimension 1, in a loop over
> i (and then j?)? In that case, you have complete control over "order"
> they are created in.

We're just making one big call to create the box.  What do you mean 
about the order not being very structured?

Rob

> I hope this helps.
> Iulian
>
>
> Iulian


More information about the moab-dev mailing list