Function2.H

Go to the documentation of this file.

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Copyright (C) 2020-2024 OpenFOAM Foundation
     \\/     M anipulation  |
-------------------------------------------------------------------------------
License
    This file is part of OpenFOAM.
 
    OpenFOAM is free software: you can redistribute it and/or modify it
    under the terms of the GNU General Public License as published by
    the Free Software Foundation, either version 3 of the License, or
    (at your option) any later version.
 
    OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
    ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
    FITNESS FOR A PARTICULAR PURPOSE.  See the GNU General Public License
    for more details.
 
    You should have received a copy of the GNU General Public License
    along with OpenFOAM.  If not, see <http://www.gnu.org/licenses/>.
 
Class
    Foam::Function2
 
Description
    Run-time selectable function of two variables
 
    with many options provided from simple constant values to complex
    functions, interpolated tabulated data etc. etc.
 
SourceFiles
    Function2.C
    Function2New.C
 
\*---------------------------------------------------------------------------*/
 
#ifndef Function2_H
#define Function2_H
 
#include "dictionary.H"
#include "unitConversion.H"
#include "tmp.H"
#include "typeInfo.H"
#include "Field.H"
#include "addToRunTimeSelectionTable.H"
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
namespace Foam
{
 
// Forward declaration of friend functions and operators
template<class Type> class Function2;
template<class Type> Ostream& operator<<(Ostream&, const Function2<Type>&);
 
/*---------------------------------------------------------------------------*\
                 Class Function2s::unitConversions Declaration
\*---------------------------------------------------------------------------*/
 
namespace Function2s
{
struct unitConversions
{
    // Public Data
 
        //- Unit conversion for x-axis values
        unitConversion x;
 
        //- Unit conversion for y-axis values
        unitConversion y;
 
        //- Unit conversion for result values
        unitConversion value;
};
}
 
 
// Global Functions
 
//- Generate an error in an context where unit conversions are not supported
void assertNoConvertUnits
(
    const word& typeName,
    const Function2s::unitConversions& units,
    const dictionary& dict
);
 
 
/*---------------------------------------------------------------------------*\
                          Class Function2 Declaration
\*---------------------------------------------------------------------------*/
 
template<class Type>
class Function2
:
    public tmp<Function2<Type>>::refCount
{
 
protected:
 
    // Protected data
 
        //- Name of entry
        const word name_;
 
 
public:
 
    typedef Type returnType;
 
 
    //- Runtime type information
    TypeName("Function2");
 
 
    // Declare runtime constructor selection tables
 
        declareRunTimeSelectionTable
        (
            autoPtr,
            Function2,
            dictionary,
            (
                const word& name,
                const Function2s::unitConversions& units,
                const dictionary& dict
            ),
            (name, units, dict)
        );
 
        declareRunTimeSelectionTable
        (
            autoPtr,
            Function2,
            Istream,
            (
                const word& name,
                const Function2s::unitConversions& units,
                Istream& is
            ),
            (name, units, is)
        );
 
 
    // Constructors
 
        //- Construct from name
        Function2(const word& name);
 
        //- Copy constructor
        Function2(const Function2<Type>& f2);
 
        //- Construct and return a clone
        virtual tmp<Function2<Type>> clone() const = 0;
 
 
    // Selectors
 
        //- Select from dictionary
        static autoPtr<Function2<Type>> New
        (
            const word& name,
            const Function2s::unitConversions& units,
            const dictionary& dict
        );
 
        //- Select from dictionary
        static autoPtr<Function2<Type>> New
        (
            const word& name,
            const unitConversion& xUnits,
            const unitConversion& yUnits,
            const unitConversion& valueUnits,
            const dictionary& dict
        );
 
        //- Select from Istream
        static autoPtr<Function2<Type>> New
        (
            const word& name,
            const Function2s::unitConversions& units,
            const word& Function2Type,
            Istream& is
        );
 
        //- Select from Istream
        static autoPtr<Function2<Type>> New
        (
            const word& name,
            const unitConversion& xUnits,
            const unitConversion& yUnits,
            const unitConversion& valueUnits,
            const word& Function2Type,
            Istream& is
        );
 
 
    //- Destructor
    virtual ~Function2();
 
 
    // Member Functions
 
        //- Return the name of the entry
        const word& name() const;
 
        //- Return value as a function of two scalars
        virtual Type value(const scalar x, const scalar y) const = 0;
 
        //- Return value as a function of two scalar fields
        virtual tmp<Field<Type>> value
        (
            const scalarField& x,
            const scalarField& y
        ) const = 0;
 
        //- Write in dictionary format
        virtual void write
        (
            Ostream& os,
            const Function2s::unitConversions&
        ) const = 0;
 
 
    // Member Operators
 
        //- Assignment
        void operator=(const Function2<Type>&);
 
 
    // IOstream Operators
 
        //- Ostream Operator
        friend Ostream& operator<< <Type>
        (
            Ostream& os,
            const Function2<Type>& func
        );
};
 
 
template<class Type>
void writeEntry(Ostream& os, const Function2<Type>& f2);
 
template<class Type>
void writeEntry
(
    Ostream& os,
    const Function2s::unitConversions& units,
    const Function2<Type>& f2
);
 
template<class Type>
void writeEntry
(
    Ostream& os,
    const unitConversion& xUnits,
    const unitConversion& yUnits,
    const unitConversion& valueUnits,
    const Function2<Type>& f2
);
 
 
/*---------------------------------------------------------------------------*\
                       Class FieldFunction2 Declaration
\*---------------------------------------------------------------------------*/
 
template<class Type, class Function2Type>
class FieldFunction2
:
    public Function2<Type>
{
 
public:
 
    // Constructors
 
        //- Construct from name
        FieldFunction2(const word& name);
 
        //- Construct and return a clone
        virtual tmp<Function2<Type>> clone() const;
 
 
    //- Destructor
    virtual ~FieldFunction2();
 
 
    // Member Functions
 
        //- Return value as a function of two scalars
        virtual Type value(const scalar x, const scalar y) const = 0;
 
        //- Return value as a function of two scalar fields
        virtual tmp<Field<Type>> value
        (
            const scalarField& x,
            const scalarField& y
        ) const;
};
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
} // End namespace Foam
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#define defineFunction2(Type)                                                  \
                                                                               \
    defineNamedTemplateTypeNameAndDebug(Function2<Type>, 0);                   \
    defineTemplateRunTimeSelectionTable(Function2<Type>, dictionary);          \
    defineTemplateRunTimeSelectionTable(Function2<Type>, Istream);
 
 
#define addFunction2(SS, Type)                                                 \
                                                                               \
    defineNamedTemplateTypeNameAndDebug(SS<Type>, 0);                          \
    typedef Function2<Type> Type##Function2;                                   \
    typedef SS<Type> Type##SS##Function2;                                      \
    addToRunTimeSelectionTable                                                 \
    (                                                                          \
        Type##Function2,                                                       \
        Type##SS##Function2,                                                   \
        dictionary                                                             \
    )
 
 
#define addStreamConstructableFunction2(SS, Type)                              \
                                                                               \
    addFunction2(SS, Type);                                                    \
    addToRunTimeSelectionTable                                                 \
    (                                                                          \
        Type##Function2,                                                       \
        Type##SS##Function2,                                                   \
        Istream                                                                \
    )
 
 
#define addScalarFunction2(SS)                                                 \
                                                                               \
    defineTypeNameAndDebug(SS, 0);                                             \
    typedef Function2<scalar> scalarFunction2;                                 \
    addToRunTimeSelectionTable(scalarFunction2, SS, dictionary)
 
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#ifdef NoRepository
    #include "Function2.C"
    #include "Constant2.H"
#endif
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 
#endif
 
// ************************************************************************* //