scalarFieldField.H
Go to the documentation of this file.
1/*---------------------------------------------------------------------------*\
2 ========= |
3 \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4 \\ / O peration |
5 \\ / A nd | www.openfoam.com
6 \\/ M anipulation |
7-------------------------------------------------------------------------------
8 Copyright (C) 2011-2017 OpenFOAM Foundation
9 Copyright (C) 2019 OpenCFD Ltd.
10-------------------------------------------------------------------------------
11License
12 This file is part of OpenFOAM.
13
14 OpenFOAM is free software: you can redistribute it and/or modify it
15 under the terms of the GNU General Public License as published by
16 the Free Software Foundation, either version 3 of the License, or
17 (at your option) any later version.
18
19 OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
20 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
21 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
22 for more details.
23
24 You should have received a copy of the GNU General Public License
25 along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
26
27InClass
28 Foam::scalarFieldField
29
30Description
31 Specialisation of FieldField<T> for scalar.
32
33SourceFiles
34 scalarFieldField.C
35
36\*---------------------------------------------------------------------------*/
37
38#ifndef scalarFieldField_H
39#define scalarFieldField_H
40
41#include "scalarField.H"
42#include "FieldField.H"
43
44#define TEMPLATE template<template<class> class Field>
46
47// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
48
49namespace Foam
50{
51
52// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
53
54template<template<class> class Field>
55void stabilise
56(
57 FieldField<Field, scalar>& Res,
58 const FieldField<Field, scalar>& sf,
59 const scalar s
60);
61
62
63template<template<class> class Field>
64tmp<FieldField<Field, scalar>> stabilise
65(
66 const FieldField<Field, scalar>&,
67 const scalar s
68);
69
70
71template<template<class> class Field>
72tmp<FieldField<Field, scalar>> stabilise
73(
74 const tmp<FieldField<Field, scalar>>&,
75 const scalar s
76);
77
78
79// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
80
81BINARY_TYPE_OPERATOR(scalar, scalar, scalar, +, add)
82BINARY_TYPE_OPERATOR(scalar, scalar, scalar, -, subtract)
83
84BINARY_OPERATOR(scalar, scalar, scalar, *, multiply)
85BINARY_OPERATOR(scalar, scalar, scalar, /, divide)
86
87BINARY_TYPE_OPERATOR_SF(scalar, scalar, scalar, /, divide)
88
89BINARY_FUNCTION(scalar, scalar, scalar, pow)
90BINARY_TYPE_FUNCTION(scalar, scalar, scalar, pow)
91
92BINARY_FUNCTION(scalar, scalar, scalar, atan2)
93BINARY_TYPE_FUNCTION(scalar, scalar, scalar, atan2)
94
95BINARY_FUNCTION(scalar, scalar, scalar, hypot)
96BINARY_TYPE_FUNCTION(scalar, scalar, scalar, hypot)
97
98
99// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
100
101UNARY_FUNCTION(scalar, scalar, pow3)
102UNARY_FUNCTION(scalar, scalar, pow4)
103UNARY_FUNCTION(scalar, scalar, pow5)
104UNARY_FUNCTION(scalar, scalar, pow6)
105UNARY_FUNCTION(scalar, scalar, pow025)
106UNARY_FUNCTION(scalar, scalar, sqrt)
107UNARY_FUNCTION(scalar, scalar, cbrt)
108UNARY_FUNCTION(scalar, scalar, sign)
109UNARY_FUNCTION(scalar, scalar, pos)
110UNARY_FUNCTION(scalar, scalar, pos0)
111UNARY_FUNCTION(scalar, scalar, neg)
112UNARY_FUNCTION(scalar, scalar, neg0)
113UNARY_FUNCTION(scalar, scalar, posPart)
114UNARY_FUNCTION(scalar, scalar, negPart)
115UNARY_FUNCTION(scalar, scalar, exp)
116UNARY_FUNCTION(scalar, scalar, log)
117UNARY_FUNCTION(scalar, scalar, log10)
118UNARY_FUNCTION(scalar, scalar, sin)
119UNARY_FUNCTION(scalar, scalar, cos)
120UNARY_FUNCTION(scalar, scalar, tan)
121UNARY_FUNCTION(scalar, scalar, asin)
122UNARY_FUNCTION(scalar, scalar, acos)
123UNARY_FUNCTION(scalar, scalar, atan)
124UNARY_FUNCTION(scalar, scalar, sinh)
125UNARY_FUNCTION(scalar, scalar, cosh)
126UNARY_FUNCTION(scalar, scalar, tanh)
127UNARY_FUNCTION(scalar, scalar, asinh)
128UNARY_FUNCTION(scalar, scalar, acosh)
129UNARY_FUNCTION(scalar, scalar, atanh)
130UNARY_FUNCTION(scalar, scalar, erf)
131UNARY_FUNCTION(scalar, scalar, erfc)
132UNARY_FUNCTION(scalar, scalar, lgamma)
133UNARY_FUNCTION(scalar, scalar, j0)
134UNARY_FUNCTION(scalar, scalar, j1)
135UNARY_FUNCTION(scalar, scalar, y0)
136UNARY_FUNCTION(scalar, scalar, y1)
137
138
139#define BesselFunc(func) \
140void func \
141( \
142 FieldField<Field, scalar>& Res, \
143 const int n, \
144 const FieldField<Field, scalar>& sf \
145); \
146tmp<scalarField> func(const int n, const FieldField<Field, scalar>&); \
147tmp<scalarField> func(const int n, const tmp<FieldField<Field, scalar>>&);
148
151
152#undef BesselFunc
153
154
155// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
156
157} // End namespace Foam
158
159// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
160
161#include "undefFieldFunctionsM.H"
162
163#ifdef NoRepository
164 #include "scalarFieldField.C"
165#endif
166
167// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
168
169#endif
170
171// ************************************************************************* //
#define BINARY_FUNCTION(ReturnType, Type1, Type2, Func)
#define BINARY_TYPE_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
#define BINARY_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
#define BINARY_TYPE_OPERATOR_SF(ReturnType, Type1, Type2, Op, OpName, OpFunc)
#define UNARY_FUNCTION(ReturnType, Type1, Func, Dfunc)
#define BINARY_TYPE_FUNCTION(ReturnType, Type1, Type2, Func)
#define BesselFunc(func)
High performance macro functions for Field<Type> algebra. These expand using either array element acc...
gmvFile<< "tracers "<< particles.size()<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().x()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().y()<< " ";}gmvFile<< nl;for(const passiveParticle &p :particles){ gmvFile<< p.position().z()<< " ";}gmvFile<< nl;forAll(lagrangianScalarNames, i){ word name=lagrangianScalarNames[i];IOField< scalar > s(IOobject(name, runTime.timeName(), cloud::prefix, mesh, IOobject::MUST_READ, IOobject::NO_WRITE))
Namespace for OpenFOAM.
dimensionedScalar pow6(const dimensionedScalar &ds)
dimensionedScalar pos(const dimensionedScalar &ds)
dimensionedScalar pow5(const dimensionedScalar &ds)
dimensionedScalar erfc(const dimensionedScalar &ds)
dimensionedScalar asin(const dimensionedScalar &ds)
dimensionedScalar exp(const dimensionedScalar &ds)
dimensionedScalar tan(const dimensionedScalar &ds)
dimensionedScalar pos0(const dimensionedScalar &ds)
dimensionedScalar sign(const dimensionedScalar &ds)
dimensionedScalar lgamma(const dimensionedScalar &ds)
dimensionedScalar j1(const dimensionedScalar &ds)
dimensionedScalar pow3(const dimensionedScalar &ds)
dimensionedScalar y0(const dimensionedScalar &ds)
dimensionedScalar cosh(const dimensionedScalar &ds)
void subtract(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
dimensionedScalar sin(const dimensionedScalar &ds)
dimensionedScalar tanh(const dimensionedScalar &ds)
dimensionedScalar erf(const dimensionedScalar &ds)
dimensionedScalar sinh(const dimensionedScalar &ds)
dimensionedScalar log10(const dimensionedScalar &ds)
dimensionedScalar yn(const int n, const dimensionedScalar &ds)
dimensionedScalar log(const dimensionedScalar &ds)
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
dimensionedScalar atan2(const dimensionedScalar &x, const dimensionedScalar &y)
dimensionedScalar y1(const dimensionedScalar &ds)
dimensionedScalar negPart(const dimensionedScalar &ds)
dimensionedScalar acosh(const dimensionedScalar &ds)
dimensionedScalar sqrt(const dimensionedScalar &ds)
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
dimensionedScalar hypot(const dimensionedScalar &x, const dimensionedScalar &y)
dimensionedScalar pow4(const dimensionedScalar &ds)
dimensionedScalar jn(const int n, const dimensionedScalar &ds)
void multiply(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
dimensionedScalar neg(const dimensionedScalar &ds)
dimensionedScalar atanh(const dimensionedScalar &ds)
dimensionedScalar stabilise(const dimensionedScalar &x, const dimensionedScalar &y)
dimensionedScalar neg0(const dimensionedScalar &ds)
dimensionedScalar cbrt(const dimensionedScalar &ds)
dimensionedScalar atan(const dimensionedScalar &ds)
void divide(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
dimensionedScalar cos(const dimensionedScalar &ds)
dimensionedScalar posPart(const dimensionedScalar &ds)
dimensionedScalar acos(const dimensionedScalar &ds)
dimensionedScalar j0(const dimensionedScalar &ds)
dimensionedScalar pow025(const dimensionedScalar &ds)
dimensionedScalar asinh(const dimensionedScalar &ds)
Specialisation of FieldField<T> for scalar.