[stBern]

I would like to create a new unary_function in OpenFOAM

inline scalar newSign(const scalar s)
{
    return (s > 0)? 1: -1;
}

How can I add this function to OpenFOAM to use with geometricFields and Fields?

[olesen]

Originally Posted by stBern

I would like to create a new unary_function in OpenFOAM

inline scalar newSign(const scalar s)
{
    return (s > 0)? 1: -1;
}

How can I add this function to OpenFOAM to use with geometricFields and Fields?

You may want something something similar for vectors ? I don't know.
To add in everywhere you will need to touch the following files (header + source):

• src/OpenFOAM/fields/Fields/Field/FieldFunctions.H
• src/OpenFOAM/fields/DimensionedFields/DimensionedField/DimensionedFieldFunctions.H
• src/OpenFOAM/fields/FieldFields/FieldField/FieldFieldFunctions.H
• src/OpenFOAM/fields/GeometricFields/GeometricField/GeometricFieldFunctions.H
• src/OpenFOAM/fields/pointPatchFields/pointPatchField/pointPatchFieldFunctions.H (need to double-check this one)

Essentially you have a lot of boilerplate to add - FieldFunctions for regular fields, DimensionedFieldFunctions to wrap for internal fields, FieldFieldFunctions to wrap all of the boundary field handling and finally GeometricFieldFunctions for the volume fields.


In addition, you will also need to define how the dimensions change (or don't) - some code in dimensionSet.{C,H}




However, how useful is it to have -1 as the sign for a zero value?