From c9ca6b9f19602c4770598b9892f595952937aa19 Mon Sep 17 00:00:00 2001
From: Vaggelis Papoutsis <vaggelisp@gmail.com>
Date: Wed, 2 Jun 2021 14:27:45 +0300
Subject: [PATCH] ENH: the adjoint eikonal equation grabs the epsilon value

from the equivalent entry of the primal eikonal equation solver, unless
specified in the corresponding adjoint dictionary.
---
 .../adjointEikonalSolverIncompressible.C                 | 9 +++++++--
 1 file changed, 7 insertions(+), 2 deletions(-)

diff --git a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointEikonalSolver/adjointEikonalSolverIncompressible.C b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointEikonalSolver/adjointEikonalSolverIncompressible.C
index 90ca968de8..41348f1bf6 100644
--- a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointEikonalSolver/adjointEikonalSolverIncompressible.C
+++ b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointEikonalSolver/adjointEikonalSolverIncompressible.C
@@ -5,8 +5,8 @@
     \\  /    A nd           | www.openfoam.com
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
-    Copyright (C) 2007-2020 PCOpt/NTUA
-    Copyright (C) 2013-2020 FOSS GP
+    Copyright (C) 2007-2021 PCOpt/NTUA
+    Copyright (C) 2013-2021 FOSS GP
     Copyright (C) 2019-2020 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
@@ -67,6 +67,11 @@ void adjointEikonalSolver::read()
     nEikonalIters_ = dict_.getOrDefault<label>("iters", 1000);
     tolerance_ = dict_.getOrDefault<scalar>("tolerance", 1e-6);
     epsilon_ = dict_.getOrDefault<scalar>("epsilon", 0.1);
+    const scalar defaultEps =
+        mesh_.schemesDict().subDict("wallDist").
+            subOrEmptyDict("advectionDiffusionCoeffs").
+                getOrDefault<scalar>("epsilon", 0.1);
+    epsilon_ = dict_.getOrDefault<scalar>("epsilon", defaultEps);
 }