[MOAB-dev] ghost exchange with triangles?
Ryan O'Kuinghttons
ryan.okuinghttons at noaa.gov
Mon Nov 12 10:45:21 CST 2018
Hi Iulian,
After looking into this further I realize that the issue has nothing to
do with meshes composed of triangles, I do actually get the same issue
with quadrilaterals (my previous mail said that quads were working). In
fact, I even see the same issue with the last test case we were
debugging. The difference between now and then is that the previous test
case was only calling exchange_tags with the elem_coords tag (as a
minimal reproducer), which now works in both versions. Now I am seeing
that every other tag on this test mesh is not being updated by the
exchange_tags call, even though the debug code says that they are:
1 ParallelComm(0.04 s) Unpacking tag GLOBAL_ID
1 ParallelComm(0.04 s) Unpacking tag orig_pos
1 ParallelComm(0.04 s) Unpacking tag owner
1 ParallelComm(0.04 s) Done unpacking tags.
The only difference I can see between the elem_coords and the GLOBAL_ID
or owner tags is that elem_coords is initialized as a double, and the
others are integers. Is there some reason that a call like the following
would not work with exchange_tags?
int int_def_val = 0;
merr=moab_mesh->tag_get_handle(GLOBAL_ID_TAG_NAME, 1, MB_TYPE_INTEGER,
mbmp->gid_tag, MB_TAG_DENSE, &int_def_val);
Ryan
On 11/9/18 10:56, Ryan O'Kuinghttons wrote:
> Hi Iulian, I ran into an issue with the ghost element exchange this
> week. I'm using a small test mesh with one quadrilateral, two
> triangles and two pentagon/hexagons that are broken down into
> triangles. The ghost exchange seems to work, but the tag exchange does
> not. I modified the test we were using before, if you want to have a
> look. When I plot the mesh after the exchange_tags call I see 0s at
> the the node and element ids locations, which is strange because this
> works just fine with a mesh of quadrilaterals. The test is on the
> moab4 branch, at
> esmf/src/Infrastructure/Mesh/tests/ESMC_MeshMOABGhostUTest. You should
> call 'make distclean' before rebuilding esmf, then just "make" in test
> directory and you can run it with:
>
> mpirun -np 4 $ESMF_DIR/tests/*/*/ESMC_MeshMOABGhostUTest
>
> I am traveling today, and I will be in and out most of the day, but
> will check email when I can. I have a feeling this is another issue
> with how we are setting the tags, but hopefully it's a quick fix.
> Thanks again,
>
More information about the moab-dev
mailing list