From 312c7a1c32d8ddb17bdf9d82e077096ad9507d0c Mon Sep 17 00:00:00 2001
From: Vaggelis Papoutsis <vaggelisp@gmail.com>
Date: Wed, 24 Jan 2024 14:39:20 +0200
Subject: [PATCH] FIX: replaced temp internalField with DimensionedField::null
(fixes #3082)
---
.../adjointOptimisation/adjoint/global/boundaryFieldsFwd.H | 6 +++---
.../adjointOptimisation/adjoint/global/createZeroField.H | 4 ++--
2 files changed, 5 insertions(+), 5 deletions(-)
diff --git a/src/optimisation/adjointOptimisation/adjoint/global/boundaryFieldsFwd.H b/src/optimisation/adjointOptimisation/adjoint/global/boundaryFieldsFwd.H
index 4a4292e4e18..0b1feeb72fc 100644
--- a/src/optimisation/adjointOptimisation/adjoint/global/boundaryFieldsFwd.H
+++ b/src/optimisation/adjointOptimisation/adjoint/global/boundaryFieldsFwd.H
@@ -50,13 +50,13 @@ namespace Foam
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
-//volFields
+//- volFields
typedef typename volScalarField::Boundary boundaryScalarField;
typedef typename volVectorField::Boundary boundaryVectorField;
typedef typename volTensorField::Boundary boundaryTensorField;
-//pointFields - actually a plain list of fields with dimension equal to the
-// number of points per patch
+//- pointFields - actually a plain list of fields with dimension equal to the
+//- number of points per patch
typedef List<Field<scalar>> pointBoundaryScalarField;
typedef List<Field<vector>> pointBoundaryVectorField;
typedef List<Field<tensor>> pointBoundaryTensorField;
diff --git a/src/optimisation/adjointOptimisation/adjoint/global/createZeroField.H b/src/optimisation/adjointOptimisation/adjoint/global/createZeroField.H
index 93035d286f4..8f9d7734563 100644
--- a/src/optimisation/adjointOptimisation/adjoint/global/createZeroField.H
+++ b/src/optimisation/adjointOptimisation/adjoint/global/createZeroField.H
@@ -89,7 +89,7 @@ createZeroBoundaryPtr
Boundary;
// Make sure that the patchFields to be generated will be of type
- // calculated, even if they are of constraint type
+ // calculated, even if they are of constraint type.
// Necessary to avoid unexpected behaviour when computing sensitivities
// on symmetry patches (not a good practice either way)
const fvBoundaryMesh& bm = mesh.boundary();
@@ -113,7 +113,7 @@ createZeroBoundaryPtr
new Boundary
(
mesh.boundary(),
- mesh.V()*pTraits<Type>::zero, // Dummy internal field,
+ DimensionedField<Type, volMesh>::null(),
wordList(bm.size(), fvPatchFieldBase::calculatedType()),
actualPatchTypes
)
--
GitLab