From 26b99d66005fcaf235968278e38bb1d336fb4cba Mon Sep 17 00:00:00 2001 From: Anthony Royer <anthony.royer@uliege.be> Date: Wed, 23 Nov 2022 11:11:20 +0100 Subject: [PATCH] up --- examples/helmholtz/crossPoints/Subproblem3D.cpp | 2 +- examples/helmholtz/crossPoints/ddm3D.cpp | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/examples/helmholtz/crossPoints/Subproblem3D.cpp b/examples/helmholtz/crossPoints/Subproblem3D.cpp index 9d9a162e..bee5fc19 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 2e44cebb..93818b90 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); } } -- GitLab