[Getfem-commits] (no subject)

[Tetsuo Koyama]

branch: devel-tetsuo-outer_faces_in_ball
commit 2ebd0a91a55caf28e0b7867abe936e035f991886
Author: Tetsuo Koyama <address@hidden>
AuthorDate: Sat Feb 22 14:56:05 2020 +0900

    :heavy_plus_sign: select_faces_in_ball
---
 tests/thermo_elasticity_electrical_coupling.cc | 21 ++++++++++++++++++++-
 1 file changed, 20 insertions(+), 1 deletion(-)

diff --git a/tests/thermo_elasticity_electrical_coupling.cc 
b/tests/thermo_elasticity_electrical_coupling.cc
index eb38cda..a9968ce 100644
--- a/tests/thermo_elasticity_electrical_coupling.cc
+++ b/tests/thermo_elasticity_electrical_coupling.cc
@@ -150,13 +150,32 @@ int main(int argc, char *argv[]) {
   getfem::mesh_region fb5
     = getfem::select_faces_of_normal(mesh, border_faces,
                                      base_small_vector(0., -1.), 0.01);
+  getfem::mesh_region fb6
+    = getfem::select_faces_in_ball(mesh, border_faces, base_node(25., 12.5),
+                                   8.+0.01*h);
+  getfem::mesh_region fb7
+    = getfem::select_faces_in_ball(mesh, border_faces, base_node(50., 12.5),
+                                   8.+0.01*h);
+  getfem::mesh_region fb8
+    = getfem::select_faces_in_ball(mesh, border_faces, base_node(75., 12.5),
+                                   8.+0.01*h);
 
   size_type RIGHT_BOUND = 1, LEFT_BOUND = 2, TOP_BOUND = 3, BOTTOM_BOUND = 4;
+  size_type HOLE_BOUND = 5, HOLE1_BOUND = 6, HOLE2_BOUND = 7, HOLE3_BOUND = 8;
   mesh.region( RIGHT_BOUND) = getfem::mesh_region::subtract(fb2, fb1);
   mesh.region(  LEFT_BOUND) = getfem::mesh_region::subtract(fb3, fb1);
   mesh.region(   TOP_BOUND) = getfem::mesh_region::subtract(fb4, fb1);
   mesh.region(BOTTOM_BOUND) = getfem::mesh_region::subtract(fb5, fb1);
- 
+  mesh.region(  HOLE_BOUND) = fb1;
+  mesh.region( HOLE1_BOUND) = fb6;
+  mesh.region( HOLE2_BOUND) = fb7;
+  mesh.region( HOLE3_BOUND) = fb8;
+
+  mesh.region_merge(HOLE1_BOUND, HOLE2_BOUND);
+  mesh.region_merge(HOLE1_BOUND, HOLE3_BOUND);
+
+  GMM_ASSERT1(mesh.region(HOLE_BOUND) == mesh.region(HOLE1_BOUND), "Region is 
not equal.");
+
   if (export_mesh) {
     getfem::vtk_export exp("mesh.vtk", false);
     exp.exporting(mesh);