From 7192cc924c499c258ab2aeb9ecd4468aafe1cf68 Mon Sep 17 00:00:00 2001
From: Vaggelis Papoutsis <vaggelisp@gmail.com>
Date: Thu, 30 Jan 2020 19:19:21 +0200
Subject: [PATCH] ENH: added function just returning the sensitivities

without recomputing them. Use with caution!
---
 .../adjointSensitivity/adjointSensitivityIncompressible.C   | 6 ++++++
 .../adjointSensitivity/adjointSensitivityIncompressible.H   | 4 ++++
 2 files changed, 10 insertions(+)

diff --git a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.C b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.C
index b904c4032e3..5ed9ec2fa1e 100644
--- a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.C
+++ b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.C
@@ -121,6 +121,12 @@ const scalarField& adjointSensitivity::calculateSensitivities()
 }
 
 
+const scalarField& adjointSensitivity::getSensitivities() const
+{
+    return derivatives_;
+}
+
+
 void adjointSensitivity::clearSensitivities()
 {
     derivatives_ = scalar(0);
diff --git a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.H b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.H
index 71b0f3b1bd0..f0727a17eac 100644
--- a/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.H
+++ b/src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/adjointSensitivity/adjointSensitivityIncompressible.H
@@ -177,6 +177,10 @@ public:
        //  Used with optimisation libraries
        virtual const scalarField& calculateSensitivities();
 
+       //- Returns the sensitivity fields
+       //  Assumes it has already been updated/computed
+       const scalarField& getSensitivities() const;
+
        //- Zero sensitivity fields and their constituents
        virtual void clearSensitivities();
 
-- 
GitLab