From 86e5b360dba97646f2c7bb785ac44ca1196ed227 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Thu, 17 Mar 2016 11:33:54 +0000
Subject: [PATCH] Added cmptSign

---
 src/foam/fields/Fields/Field/FieldFunctions.C | 24 +++++++++++++++++++
 src/foam/fields/Fields/Field/FieldFunctions.H | 10 ++++++++
 2 files changed, 34 insertions(+)

diff --git a/src/foam/fields/Fields/Field/FieldFunctions.C b/src/foam/fields/Fields/Field/FieldFunctions.C
index 8719f3cf6..1ad527657 100644
--- a/src/foam/fields/Fields/Field/FieldFunctions.C
+++ b/src/foam/fields/Fields/Field/FieldFunctions.C
@@ -297,6 +297,30 @@ tmp<Field<Type> > cmptMag(const tmp<Field<Type> >& tf)
 }
 
 
+template<class Type>
+void cmptSign(Field<Type>& res, const UList<Type>& f)
+{
+    TFOR_ALL_F_OP_FUNC_F(Type, res, =, cmptSign, Type, f)
+}
+
+template<class Type>
+tmp<Field<Type> > cmptSign(const UList<Type>& f)
+{
+    tmp<Field<Type> > tRes(new Field<Type>(f.size()));
+    cmptSign(tRes(), f);
+    return tRes;
+}
+
+template<class Type>
+tmp<Field<Type> > cmptSign(const tmp<Field<Type> >& tf)
+{
+    tmp<Field<Type> > tRes = reuseTmp<Type, Type>::New(tf);
+    cmptSign(tRes(), tf());
+    reuseTmp<Type, Type>::clear(tf);
+    return tRes;
+}
+
+
 #define TMP_UNARY_FUNCTION(ReturnType, Func)                                  \
                                                                               \
 template<class Type>                                                          \
diff --git a/src/foam/fields/Fields/Field/FieldFunctions.H b/src/foam/fields/Fields/Field/FieldFunctions.H
index b4a6eb515..7400083d1 100644
--- a/src/foam/fields/Fields/Field/FieldFunctions.H
+++ b/src/foam/fields/Fields/Field/FieldFunctions.H
@@ -152,6 +152,16 @@ tmp<Field<Type> > cmptMag(const UList<Type>& f);
 template<class Type>
 tmp<Field<Type> > cmptMag(const tmp<Field<Type> >& tf);
 
+
+template<class Type>
+void cmptSign(Field<Type>& res, const UList<Type>& f);
+
+template<class Type>
+tmp<Field<Type> > cmptSign(const UList<Type>& f);
+
+template<class Type>
+tmp<Field<Type> > cmptSign(const tmp<Field<Type> >& tf);
+
 #define TMP_UNARY_FUNCTION(ReturnType, Func)                                  \
                                                                               \
 template<class Type>                                                          \