28 #define TEMPLATE template<template<class> class PatchField, class GeoMesh>
38 template<
template<
class>
class PatchField,
class GeoMesh>
51 template<
template<
class>
class PatchField,
class GeoMesh>
52 tmp<GeometricField<scalar, PatchField, GeoMesh> >
stabilise
64 "stabilise(" + gsf.
name() +
',' + ds.
name() +
')',
81 template<
template<
class>
class PatchField,
class GeoMesh>
82 tmp<GeometricField<scalar, PatchField, GeoMesh> >
stabilise
95 "stabilise(" + gsf.
name() +
',' + ds.
name() +
')',
120 template<template<class> class PatchField, class GeoMesh>
128 pow(Pow.internalField(), gsf1.internalField(), gsf2.internalField());
129 pow(Pow.boundaryField(), gsf1.boundaryField(), gsf2.boundaryField());
133 template<
template<
class>
class PatchField,
class GeoMesh>
134 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
146 "pow(" + gsf1.
name() +
',' + gsf2.
name() +
')',
161 pow(tPow(), gsf1, gsf2);
167 template<
template<
class>
class PatchField,
class GeoMesh>
168 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
181 "pow(" + gsf1.
name() +
',' + gsf2.
name() +
')',
190 pow(tPow(), gsf1, gsf2);
198 template<
template<
class>
class PatchField,
class GeoMesh>
199 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
212 "pow(" + gsf1.
name() +
',' + gsf2.
name() +
')',
221 pow(tPow(), gsf1, gsf2);
228 template<
template<
class>
class PatchField,
class GeoMesh>
229 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
241 <scalar, scalar, scalar, scalar, PatchField,
GeoMesh>::New
245 "pow(" + gsf1.
name() +
',' + gsf2.
name() +
')',
254 pow(tPow(), gsf1, gsf2);
257 <scalar, scalar, scalar, scalar, PatchField,
GeoMesh>
264 template<
template<
class>
class PatchField,
class GeoMesh>
277 template<
template<
class>
class PatchField,
class GeoMesh>
278 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
290 "pow(" + gsf.
name() +
',' + ds.
name() +
')',
301 pow(tPow(), gsf, ds);
306 template<
template<
class>
class PatchField,
class GeoMesh>
307 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
320 "pow(" + gsf.
name() +
',' + ds.
name() +
')',
325 pow(tPow(), gsf, ds);
332 template<
template<
class>
class PatchField,
class GeoMesh>
333 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
342 template<
template<
class>
class PatchField,
class GeoMesh>
343 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
353 template<
template<
class>
class PatchField,
class GeoMesh>
366 template<
template<
class>
class PatchField,
class GeoMesh>
367 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
379 "pow(" + ds.
name() +
',' + gsf.
name() +
')',
390 pow(tPow(), ds, gsf);
396 template<
template<
class>
class PatchField,
class GeoMesh>
397 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
410 "pow(" + ds.
name() +
',' + gsf.
name() +
')',
415 pow(tPow(), ds, gsf);
422 template<
template<
class>
class PatchField,
class GeoMesh>
423 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
432 template<
template<
class>
class PatchField,
class GeoMesh>
433 tmp<GeometricField<scalar, PatchField, GeoMesh> >
pow
480 #define BesselFunc(func) \
482 template<template<class> class PatchField, class GeoMesh> \
485 GeometricField<scalar, PatchField, GeoMesh>& gsf, \
487 const GeometricField<scalar, PatchField, GeoMesh>& gsf1 \
490 func(gsf.internalField(), n, gsf1.internalField()); \
491 func(gsf.boundaryField(), n, gsf1.boundaryField()); \
494 template<template<class> class PatchField, class GeoMesh> \
495 tmp<GeometricField<scalar, PatchField, GeoMesh> > func \
498 const GeometricField<scalar, PatchField, GeoMesh>& gsf \
501 if (!gsf.dimensions().dimensionless()) \
505 #func"(const int n, " \
506 "const GeometricField<scalar, PatchField, GeoMesh>& gsf)" \
507 ) << "gsf not dimensionless" \
508 << abort(FatalError); \
511 tmp<GeometricField<scalar, PatchField, GeoMesh> > tFunc \
513 new GeometricField<scalar, PatchField, GeoMesh> \
517 #func "(" + gsf.name() + ')', \
528 func(tFunc(), n, gsf); \
533 template<template<class> class PatchField, class GeoMesh> \
534 tmp<GeometricField<scalar, PatchField, GeoMesh> > func \
537 const tmp<GeometricField<scalar, PatchField, GeoMesh> >& tgsf \
540 const GeometricField<scalar, PatchField, GeoMesh>& gsf = tgsf(); \
542 if (!gsf.dimensions().dimensionless()) \
546 #func"(const int n, " \
547 "const tmp<GeometricField<scalar, PatchField, GeoMesh> >& gsf)" \
548 ) << " : gsf not dimensionless" \
549 << abort(FatalError); \
552 tmp<GeometricField<scalar, PatchField, GeoMesh> > tFunc \
554 reuseTmpGeometricField<scalar, scalar, PatchField, GeoMesh>::New \
557 #func "(" + gsf.name() + ')', \
562 func(tFunc(), n, gsf); \
564 reuseTmpGeometricField<scalar, scalar, PatchField, GeoMesh> \