Partition Meshes 2D

[zhang-alvin]

Mandatory Checklist

Please ensure that the following criteria are met:

• [x] Title of pull request describes the changes/features
• [x] Request at least 2 reviewers
• [x] If new files are being added, the files are no larger than 100kB. Post the file sizes.
• [x] Code coverage did not decrease. If this is a bug fix, a test should cover that bug fix. If a new feature is added, a test should be made to cover that feature.
• [x] New features have appropriate documentation strings (readable by sphinx)
• [x] Contributor has read and agreed with CONTRIBUTING.md and has added themselves to CONTRIBUTORS.md

As a general rule of thumb, try to follow PEP8 guidelines.

Description

Try to get the partition logic for 2D triangle meshes working

[codecov[bot]]

Codecov Report

Merging #1239 (4bcdf11) into master (78b6d07) will decrease coverage by 0.01%. The diff coverage is 23.33%.

@@            Coverage Diff             @@
##           master    #1239      +/-   ##
==========================================
- Coverage   52.82%   52.81%   -0.02%     
==========================================
  Files         531      531              
  Lines      109745   109772      +27     
==========================================
+ Hits        57971    57974       +3     
- Misses      51774    51798      +24     

Impacted Files	Coverage Δ
proteus/MeshTools.py	55.13% <23.33%> (-0.24%)	:arrow_down:
proteus/tests/test_boundaryconditions.py	97.46% <0.00%> (-0.13%)	:arrow_down:

---

Continue to review full report at Codecov.

Legend - Click here to learn more Δ = absolute <relative> (impact), ø = not affected, ? = missing data Powered by Codecov. Last update 78b6d07...8d3e1a8. Read the comment docs.

[zhang-alvin]

When running a serial case with the -F option on a mesh that has already been generated, the analysis proceeds with identical results (top is serial case w/o -F and bottom is with -F):

[zhang-alvin]

Running the same case with two cores yields the following assertion error:

python: proteus/partitioning.cpp:3961: int proteus::partitionNodesFromTriangleFiles(const MPI_Comm&, const char*, int, Mesh&, int): Assertion `edges_overlap.size() + nEdges_subdomain_new[rank] == edgeNodesMap.size()' failed.

A print statement for edges_overlap.size(), nEdges_subdomain_new, and edgeNodesMap.size shows the following:

rank 0 
edges_overlap=6 
nEdges_subdomain_new[0]=235 
edgeNodesMap.size()=241

rank 1
edges_overlap=18
nEdges_subdomain_new[1]=227
edgeNodesMap.size()=244

The assertion fails for rank 1 because the sum of the number of overlapping edges and number of owned edges is off by one.