zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

ENH: cref for grad(Ua) in the adjoint library

Browse Source 
Necessary when grad(Ua) is cached
This commit is contained in:
Vaggelis Papoutsis
2022-06-06 12:04:26 +03:00
committed by Andrew Heather
parent a5f59cd903
commit f6897a595f
3 changed files with 3 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

2
src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurface/sensitivitySurfaceIncompressible.C
View File

@ -700,7 +700,7 @@ void sensitivitySurface::accumulateIntegrand(const scalar dt)
{
tmp<volTensorField> tgradUa = fvc::grad(Ua);
const volTensorField::Boundary& gradUabf =
tgradUa.ref().boundaryField();
tgradUa.cref().boundaryField();
for (const label patchI : sensitivityPatchIDs_)
{

2
src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/sensitivitySurfacePoints/sensitivitySurfacePointsIncompressible.C
View File

@ -513,7 +513,7 @@ void sensitivitySurfacePoints::accumulateIntegrand(const scalar dt)
{
tmp<volTensorField> tgradUa = fvc::grad(Ua);
const volTensorField::Boundary& gradUabf =
tgradUa.ref().boundaryField();
tgradUa.cref().boundaryField();
for (const label patchI : sensitivityPatchIDs_)
{
const fvPatch& patch = mesh_.boundary()[patchI];

2
src/optimisation/adjointOptimisation/adjoint/optimisation/adjointSensitivity/incompressible/shapeSensitivities/shapeSensitivitiesIncompressible.C
View File

@ -105,7 +105,7 @@ tmp<boundaryVectorField> shapeSensitivities::dvdbMult() const
const singlePhaseTransportModel& lamTransp = primalVars_.laminarTransport();
volScalarField nuEff(lamTransp.nu() + turbVars->nutRef());
tmp<volTensorField> tgradUa = fvc::grad(Ua);
const volTensorField::Boundary& gradUabf = tgradUa.ref().boundaryField();
const volTensorField::Boundary& gradUabf = tgradUa.cref().boundaryField();
for (const label patchI : sensitivityPatchIDs_)
{