BUG: fan: bc value not updated. See #3211

src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricField.C

@@ -1092,6 +1092,10 @@ void Foam::GeometricField<Type, PatchField, GeoMesh>::relax(const scalar alpha)
         << "Relaxing" << nl << this->info() << " by " << alpha << endl;
 
     operator==(prevIter() + alpha*(*this - prevIter()));
+
+    // Make sure any e.g. jump-cyclic are updated. Note: unfortunately also
+    // triggers e.g. cyclic interpolation.
+    this->correctLocalBoundaryConditions();
 }
 
 

src/finiteArea/fields/faPatchFields/constraint/processor/processorFaPatchField.H

@@ -175,7 +175,7 @@ public:
         virtual bool ready() const;
 
         //- Return neighbour field given internal field
-        tmp<Field<Type>> patchNeighbourField() const;
+        virtual tmp<Field<Type>> patchNeighbourField() const;
 
 
     // Evaluation

src/finiteVolume/fields/fvPatchFields/constraint/cyclic/cyclicFvPatchField.H

@@ -170,7 +170,7 @@ public:
         // Evaluation functions
 
             //- Return neighbour coupled internal cell data
-            tmp<Field<Type>> patchNeighbourField() const;
+            virtual tmp<Field<Type>> patchNeighbourField() const;
 
             //- Return reference to neighbour patchField
             const cyclicFvPatchField<Type>& neighbourPatchField() const;

src/finiteVolume/fields/fvPatchFields/constraint/jumpCyclic/jumpCyclicFvPatchField.H

@@ -128,7 +128,7 @@ public:
         // Evaluation functions
 
             //- Return neighbour coupled given internal cell data
-            tmp<Field<Type>> patchNeighbourField() const;
+            virtual tmp<Field<Type>> patchNeighbourField() const;
 
             //- Update result field based on interface functionality
             virtual void updateInterfaceMatrix

src/finiteVolume/fields/fvPatchFields/constraint/jumpCyclicAMI/jumpCyclicAMIFvPatchField.H

@@ -137,7 +137,7 @@ public:
             virtual void evaluate(const Pstream::commsTypes commsType);
 
             //- Initialise interface functionality
-            void initInterfaceMatrixUpdate
+            virtual void initInterfaceMatrixUpdate
             (
                 solveScalarField& result,
                 const bool add,