up

examples/helmholtz/crossPoints/ddm2D.cpp

@@ -318,16 +318,15 @@ namespace D2 {
       gamma = gmshfem::domain::Domain("gamma");
     }
     else if (benchmark == "marmousi") {
-//      gamma = corner[0](3) |
-//              corner[0](7) |
-//              corner[0](11) |
-//              corner[0](15) |
-//              corner[0](19) |
-//              corner[0](23) |
-//              corner[0](27) |
-//              corner[0](31) |
-//              corner[0](35);
-      gamma = corner[0](3);
+      gamma = corner[0](3) |
+              corner[0](7) |
+              corner[0](11) |
+              corner[0](15) |
+              corner[0](19) |
+              corner[0](23) |
+              corner[0](27) |
+              corner[0](31) |
+              corner[0](35);
     }
 
     // Kappa definition
@@ -425,14 +424,14 @@ namespace D2 {
     }
     else if (benchmark == "marmousi") {
       u(3).addConstraint(corner[0](3), 1.);
-//      u(7).addConstraint(corner[0](5), 1.);
-//      u(11).addConstraint(corner[0](11), 1.);
-//      u(15).addConstraint(corner[0](15), 1.);
-//      u(19).addConstraint(corner[0](19), 1.);
-//      u(23).addConstraint(corner[0](23), 1.);
-//      u(27).addConstraint(corner[0](27), 1.);
-//      u(31).addConstraint(corner[0](31), 1.);
-//      u(35).addConstraint(corner[0](35), 1.);
+      u(7).addConstraint(corner[0](5), 1.);
+      u(11).addConstraint(corner[0](11), 1.);
+      u(15).addConstraint(corner[0](15), 1.);
+      u(19).addConstraint(corner[0](19), 1.);
+      u(23).addConstraint(corner[0](23), 1.);
+      u(27).addConstraint(corner[0](27), 1.);
+      u(31).addConstraint(corner[0](31), 1.);
+      u(35).addConstraint(corner[0](35), 1.);
     }
     
     std::vector< std::vector< Subproblem * > > subproblem(nDomX);