ENH: add field operations for complex (#1365)

applications/test/complex/Test-complex.C

@@ -108,6 +108,7 @@ int main(int argc, char *argv[])
     Info<< "sin: " << sin(fld1) << nl;
 
     Info<< "operator + : " << (fld1 + fld2) << nl;
+    Info<< "operator + : " << (fld1 + fld2 + complex(1,0)) << nl;
 
     // Some operators are still incomplete
 

src/OpenFOAM/fields/Fields/complex/complexField.C

@@ -144,6 +144,15 @@ complex sumProd(const UList<complex>& f1, const UList<complex>& f2)
 }
 
 
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+BINARY_TYPE_OPERATOR(complex, complex, complex, +, add)
+BINARY_TYPE_OPERATOR(complex, complex, complex, -, subtract)
+
+BINARY_OPERATOR(complex, complex, complex, *, multiply)
+BINARY_OPERATOR(complex, complex, complex, /, divide)
+
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 UNARY_FUNCTION(complex, complex, pow3)

src/OpenFOAM/fields/Fields/complex/complexField.H

@@ -81,6 +81,15 @@ template<>
 complex sumProd(const UList<complex>& f1, const UList<complex>& f2);
 
 
+// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
+
+BINARY_TYPE_OPERATOR(complex, complex, complex, +, add)
+BINARY_TYPE_OPERATOR(complex, complex, complex, -, subtract)
+
+BINARY_OPERATOR(complex, complex, complex, *, multiply)
+BINARY_OPERATOR(complex, complex, complex, /, divide)
+
+
 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
 UNARY_FUNCTION(complex, complex, pow3)