diff --git a/examples/helmholtz/crossPoints/Subproblem3D.cpp b/examples/helmholtz/crossPoints/Subproblem3D.cpp
index 9d9a162e57931fcee057c306ca09f7b9b9c39fba..bee5fc196f4d886f1a7512b57558e9b10dd0fe99 100755
--- a/examples/helmholtz/crossPoints/Subproblem3D.cpp
+++ b/examples/helmholtz/crossPoints/Subproblem3D.cpp
@@ -493,7 +493,7 @@ namespace D3 {
       sigma = sigmaStar * distSigma * distSigma / (_size * _size);
     }
     
-    gmshfem::function::ScalarFunction< std::complex< double > > kMod = 1. + im * sigma/* / kappaMod*/;
+    gmshfem::function::ScalarFunction< std::complex< double > > kMod = 1. + im * sigma / kappaMod;
     if(orientation() == "E" || orientation() == "W") {
       D = gmshfem::function::tensorDiag< std::complex< double > >(1./kMod, kMod, kMod);
       E = kMod;
diff --git a/examples/helmholtz/crossPoints/ddm3D.cpp b/examples/helmholtz/crossPoints/ddm3D.cpp
index 2e44cebb50a578e91899bd361e127ea3cea8d670..93818b90bc0df0582a033705d63f91c3d6eebddd 100755
--- a/examples/helmholtz/crossPoints/ddm3D.cpp
+++ b/examples/helmholtz/crossPoints/ddm3D.cpp
@@ -764,7 +764,7 @@ namespace D3 {
             if(index == 0) {
               formulation(index).integral(dof(lambda), tf(u(index)), gamma, gauss);
               formulation(index).integral(dof(u(index)), tf(lambda), gamma, gauss);
-              formulation(index).integral(formulation.physicalSource(- 1.), tf(lambda), gamma, gauss);
+              formulation(index).integral(formulation.physicalSource(- fAnalytic), tf(lambda), gamma, gauss);
             }
           }