31 #define TEMPLATE template<class Type, class GeoMesh> 41 template<
class Type,
class GeoMesh, direction r>
56 "pow(" + df.
name() +
',' +
name(r) +
')',
64 pow<Type, r, GeoMesh>(tres.ref().field(), df.
field());
70 template<
class Type,
class GeoMesh, direction r>
74 const tmp<DimensionedField<Type, GeoMesh>>& tdf,
80 const DimensionedField<Type, GeoMesh>& df = tdf();
86 "pow(" + df.name() +
',' +
name(r) +
')',
87 pow(df.dimensions(), r)
90 pow<Type, r, GeoMesh>(tres.ref().field(), df.field());
97 template<
class Type,
class GeoMesh>
99 sqr(
const DimensionedField<Type, GeoMesh>& df)
104 tmp<DimensionedField<outerProductType, GeoMesh>>
::New 108 "sqr(" + df.name() +
')',
116 sqr(tres.ref().field(), df.field());
121 template<
class Type,
class GeoMesh>
123 sqr(
const tmp<DimensionedField<Type, GeoMesh>>& tdf)
127 const DimensionedField<Type, GeoMesh>& df = tdf();
133 "sqr(" + df.name() +
')',
137 sqr(tres.ref().field(), df.field());
144 template<
class Type,
class GeoMesh>
146 magSqr(
const DimensionedField<Type, GeoMesh>& df)
151 tmp<DimensionedField<magType, GeoMesh>>
::New 155 "magSqr(" + df.name() +
')',
168 template<
class Type,
class GeoMesh>
170 magSqr(
const tmp<DimensionedField<Type, GeoMesh>>& tdf)
174 const DimensionedField<Type, GeoMesh>& df = tdf();
180 "magSqr(" + df.name() +
')',
184 magSqr(tres.ref().field(), df.field());
191 template<
class Type,
class GeoMesh>
193 mag(
const DimensionedField<Type, GeoMesh>& df)
198 tmp<DimensionedField<magType, GeoMesh>>
::New 202 "mag(" + df.name() +
')',
210 mag(tres.ref().field(), df.field());
215 template<
class Type,
class GeoMesh>
217 mag(
const tmp<DimensionedField<Type, GeoMesh>>& tdf)
221 const DimensionedField<Type, GeoMesh>& df = tdf();
227 "mag(" + df.name() +
')',
231 mag(tres.ref().field(), df.field());
238 template<
class Type,
class GeoMesh>
255 "cmptAv(" + df.
name() +
')',
269 template<
class Type,
class GeoMesh>
287 "cmptAv(" + df.
name() +
')',
298 #define UNARY_REDUCTION_FUNCTION(returnType, func, dfunc) \ 300 template<class Type, class GeoMesh> \ 301 dimensioned<returnType> func \ 303 const DimensionedField<Type, GeoMesh>& df \ 306 return dimensioned<returnType> \ 308 #func "(" + df.name() + ')', \ 314 template<class Type, class GeoMesh> \ 315 dimensioned<returnType> func \ 317 const tmp<DimensionedField<Type, GeoMesh>>& tdf1 \ 320 dimensioned<returnType> res = func(tdf1()); \ 334 #undef UNARY_REDUCTION_FUNCTION 366 #define PRODUCT_OPERATOR(product, op, opFunc) \ 368 template<class Type1, class Type2, class GeoMesh> \ 369 tmp<DimensionedField<typename product<Type1, Type2>::type, GeoMesh>> \ 372 const DimensionedField<Type1, GeoMesh>& df1, \ 373 const DimensionedField<Type2, GeoMesh>& df2 \ 376 typedef typename product<Type1, Type2>::type productType; \ 379 tmp<DimensionedField<productType, GeoMesh>>::New \ 383 '(' + df1.name() + #op + df2.name() + ')', \ 388 df1.dimensions() op df2.dimensions() \ 391 Foam::opFunc(tres.ref().field(), df1.field(), df2.field()); \ 397 template<class Type1, class Type2, class GeoMesh> \ 398 tmp<DimensionedField<typename product<Type1, Type2>::type, GeoMesh>> \ 401 const DimensionedField<Type1, GeoMesh>& df1, \ 402 const tmp<DimensionedField<Type2, GeoMesh>>& tdf2 \ 405 typedef typename product<Type1, Type2>::type productType; \ 407 const DimensionedField<Type2, GeoMesh>& df2 = tdf2(); \ 410 reuseTmpDimensionedField<productType, Type2, GeoMesh>::New \ 413 '(' + df1.name() + #op + df2.name() + ')', \ 414 df1.dimensions() op df2.dimensions() \ 417 Foam::opFunc(tres.ref().field(), df1.field(), df2.field()); \ 424 template<class Type1, class Type2, class GeoMesh> \ 425 tmp<DimensionedField<typename product<Type1, Type2>::type, GeoMesh>> \ 428 const tmp<DimensionedField<Type1, GeoMesh>>& tdf1, \ 429 const DimensionedField<Type2, GeoMesh>& df2 \ 432 typedef typename product<Type1, Type2>::type productType; \ 434 const DimensionedField<Type1, GeoMesh>& df1 = tdf1(); \ 437 reuseTmpDimensionedField<productType, Type1, GeoMesh>::New \ 440 '(' + df1.name() + #op + df2.name() + ')', \ 441 df1.dimensions() op df2.dimensions() \ 444 Foam::opFunc(tres.ref().field(), df1.field(), df2.field()); \ 451 template<class Type1, class Type2, class GeoMesh> \ 452 tmp<DimensionedField<typename product<Type1, Type2>::type, GeoMesh>> \ 455 const tmp<DimensionedField<Type1, GeoMesh>>& tdf1, \ 456 const tmp<DimensionedField<Type2, GeoMesh>>& tdf2 \ 459 typedef typename product<Type1, Type2>::type productType; \ 461 const DimensionedField<Type1, GeoMesh>& df1 = tdf1(); \ 462 const DimensionedField<Type2, GeoMesh>& df2 = tdf2(); \ 465 reuseTmpTmpDimensionedField \ 466 <productType, Type1, Type1, Type2, GeoMesh>::New \ 470 '(' + df1.name() + #op + df2.name() + ')', \ 471 df1.dimensions() op df2.dimensions() \ 474 Foam::opFunc(tres.ref().field(), df1.field(), df2.field()); \ 482 template<class Form, class Type, class GeoMesh> \ 483 tmp<DimensionedField<typename product<Type, Form>::type, GeoMesh>> \ 486 const DimensionedField<Type, GeoMesh>& df1, \ 487 const dimensioned<Form>& dvs \ 490 typedef typename product<Type, Form>::type productType; \ 493 tmp<DimensionedField<productType, GeoMesh>>::New \ 497 '(' + df1.name() + #op + dvs.name() + ')', \ 502 df1.dimensions() op dvs.dimensions() \ 505 Foam::opFunc(tres.ref().field(), df1.field(), dvs.value()); \ 511 template<class Form, class Cmpt, direction nCmpt, class Type, class GeoMesh> \ 512 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 515 const DimensionedField<Type, GeoMesh>& df1, \ 516 const VectorSpace<Form,Cmpt,nCmpt>& vs \ 519 return df1 op dimensioned<Form>(static_cast<const Form&>(vs)); \ 523 template<class Form, class Type, class GeoMesh> \ 524 tmp<DimensionedField<typename product<Type, Form>::type, GeoMesh>> \ 527 const tmp<DimensionedField<Type, GeoMesh>>& tdf1, \ 528 const dimensioned<Form>& dvs \ 531 typedef typename product<Type, Form>::type productType; \ 533 const DimensionedField<Type, GeoMesh>& df1 = tdf1(); \ 536 reuseTmpDimensionedField<productType, Type, GeoMesh>::New \ 539 '(' + df1.name() + #op + dvs.name() + ')', \ 540 df1.dimensions() op dvs.dimensions() \ 543 Foam::opFunc(tres.ref().field(), df1.field(), dvs.value()); \ 550 template<class Form, class Cmpt, direction nCmpt, class Type, class GeoMesh> \ 551 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 554 const tmp<DimensionedField<Type, GeoMesh>>& tdf1, \ 555 const VectorSpace<Form,Cmpt,nCmpt>& vs \ 558 return tdf1 op dimensioned<Form>(static_cast<const Form&>(vs)); \ 562 template<class Form, class Type, class GeoMesh> \ 563 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 566 const dimensioned<Form>& dvs, \ 567 const DimensionedField<Type, GeoMesh>& df1 \ 570 typedef typename product<Form, Type>::type productType; \ 573 tmp<DimensionedField<productType, GeoMesh>>::New \ 577 '(' + dvs.name() + #op + df1.name() + ')', \ 582 dvs.dimensions() op df1.dimensions() \ 585 Foam::opFunc(tres.ref().field(), dvs.value(), df1.field()); \ 591 template<class Form, class Cmpt, direction nCmpt, class Type, class GeoMesh> \ 592 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 595 const VectorSpace<Form,Cmpt,nCmpt>& vs, \ 596 const DimensionedField<Type, GeoMesh>& df1 \ 599 return dimensioned<Form>(static_cast<const Form&>(vs)) op df1; \ 603 template<class Form, class Type, class GeoMesh> \ 604 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 607 const dimensioned<Form>& dvs, \ 608 const tmp<DimensionedField<Type, GeoMesh>>& tdf1 \ 611 typedef typename product<Form, Type>::type productType; \ 613 const DimensionedField<Type, GeoMesh>& df1 = tdf1(); \ 616 reuseTmpDimensionedField<productType, Type, GeoMesh>::New \ 619 '(' + dvs.name() + #op + df1.name() + ')', \ 620 dvs.dimensions() op df1.dimensions() \ 623 Foam::opFunc(tres.ref().field(), dvs.value(), df1.field()); \ 630 template<class Form, class Cmpt, direction nCmpt, class Type, class GeoMesh> \ 631 tmp<DimensionedField<typename product<Form, Type>::type, GeoMesh>> \ 634 const VectorSpace<Form,Cmpt,nCmpt>& vs, \ 635 const tmp<DimensionedField<Type, GeoMesh>>& tdf1 \ 638 return dimensioned<Form>(static_cast<const Form&>(vs)) op tdf1; \ 650 #undef PRODUCT_OPERATOR void divide(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
#define PRODUCT_OPERATOR(product, op, opFunc)
#define BINARY_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
void multiply(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
static tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< Type1, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
Type gMin(const FieldField< Field, Type > &f)
const word & name() const noexcept
Return the object name.
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
A min/max value pair with additional methods. In addition to conveniently storing values...
#define BINARY_TYPE_OPERATOR_SF(ReturnType, Type1, Type2, Op, OpName, OpFunc)
#define BINARY_TYPE_OPERATOR_FS(ReturnType, Type1, Type2, Op, OpName, OpFunc)
void cross(FieldField< Field1, typename crossProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) >::type type
#define BINARY_TYPE_FUNCTION(ReturnType, Type1, Type2, Func)
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions, const bool initCopy=false)
Global function forwards to reuseTmpDimensionedField::New.
Field< Type >::cmptType cmptType
Component type of the field elements.
void dotdot(FieldField< Field1, typename scalarProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
void subtract(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
void outer(FieldField< Field1, typename outerProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
void dot(FieldField< Field1, typename innerProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
dimensioned< Type > cmptDivide(const dimensioned< Type > &, const dimensioned< Type > &)
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
MinMax< label > minMax(const labelHashSet &set)
Find the min/max values of labelHashSet.
#define BINARY_FUNCTION(ReturnType, Type1, Type2, Func)
Type gSum(const FieldField< Field, Type > &f)
word name(const expressions::valueTypeCode typeCode)
A word representation of a valueTypeCode. Empty for INVALID.
dimensionSet clip(const dimensionSet &ds1, const dimensionSet &ds2)
#define BINARY_TYPE_FUNCTION_FS(ReturnType, Type1, Type2, Func)
const objectRegistry & db() const noexcept
Return the local objectRegistry.
MinMax< Type > gMinMax(const FieldField< Field, Type > &f)
#define UNARY_REDUCTION_FUNCTION(returnType, func, dfunc)
dimensioned< Type > average(const DimensionedField< Type, GeoMesh > &df)
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
dimensioned< Type > cmptMultiply(const dimensioned< Type > &, const dimensioned< Type > &)
const Mesh & mesh() const noexcept
Return mesh.
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Type gMax(const FieldField< Field, Type > &f)
scalarMinMax gMinMaxMag(const FieldField< Field, Type > &f)
const fileName & instance() const noexcept
Read access to instance path component.
typeOfMag< Type >::type gSumMag(const FieldField< Field, Type > &f)
dimensioned< typename typeOfMag< Type >::type > sumMag(const DimensionedField< Type, GeoMesh > &df)
const Field< Type > & field() const noexcept
Return const-reference to the field values.
symmTypeOfRank< typename pTraits< arg1 >::cmptType, arg2 *direction(pTraits< arg1 >::rank) >::type type
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
void negate(FieldField< Field, Type > &res, const FieldField< Field, Type > &f)
Type gAverage(const FieldField< Field, Type > &f)
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
#define UNARY_OPERATOR(ReturnType, Type1, Op, OpFunc, Dfunc)
Generic mesh wrapper used by volMesh, surfaceMesh, pointMesh etc.
A class for managing temporary objects.
dimensionSet transform(const dimensionSet &ds)
Return the argument; transformations do not change the dimensions.
Defines the attributes of an object for which implicit objectRegistry management is supported...
tmp< DimensionedField< typename DimensionedField< Type, GeoMesh >::cmptType, GeoMesh >> cmptAv(const DimensionedField< Type, GeoMesh > &df)
dimensioned< typename typeOfMag< Type >::type > magSqr(const dimensioned< Type > &dt)
dimensioned< scalarMinMax > minMaxMag(const DimensionedField< Type, GeoMesh > &df)
const dimensionSet & dimensions() const noexcept
Return dimensions.
pTraits< typename pTraits< arg1 >::cmptType >::magType type