Faster symmetry planes

2011-08-12 16:02:50 +01:00 · 2011-08-12 16:02:50 +01:00 · 3a98907f8f
commit 3a98907f8f
parent 8cb18ad453
1 changed files with 3 additions and 4 deletions
--- a/src/finiteVolume/fields/fvPatchFields/basic/basicSymmetry/basicSymmetryFvPatchField.C
+++ b/src/finiteVolume/fields/fvPatchFields/basic/basicSymmetry/basicSymmetryFvPatchField.C
@ -118,12 +118,11 @@ void basicSymmetryFvPatchField<Type>::evaluate(const Pstream::commsTypes)
    }

    vectorField nHat = this->patch().nf();
+
+    Field<Type> pif = this->patchInternalField();
    Field<Type>::operator=
    (
-        (
-            this->patchInternalField()
-          + transform(I - 2.0*sqr(nHat), this->patchInternalField())
-        )/2.0
+        0.5*(pif + transform(I - 2.0*sqr(nHat), pif))
    );

    transformFvPatchField<Type>::evaluate();