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

Iulian Grindeanu iulian at mcs.anl.gov
Fri May 25 15:24:48 CDT 2012


So, first you decide if the global mesh is periodic. Let's say it is periodic in i direction. | When you SCDInterface::compute_partition, you pass gdims [3]-gdims[0]
| = N+1, where N is the number of cells, even if it is periodic in i.
That was wrong, it is gdims [3]-gdims[0] = N Iulian | If (ldims[3]-ldims[0] == gdims[3]-gdims[0]) it means that your
| partition was not broken in i direction, so when you create the box,
| you have to do periodic in i, so it will create only N nodes in i
| direction.
| It is a bad case (small number of processors, or number of data in j
| direction is very large compared to i). In "most" of the cases, you
| should have a break in both directions, then your periodicity will be
| broken.
| But you have to account for these "worse case scenarios" too.
| Iulian
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/moab-dev/attachments/20120525/c1af8ac1/attachment.html>


More information about the moab-dev mailing list