DimensionedScalarField.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration |
5  \\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
6  \\/ M anipulation |
7 -------------------------------------------------------------------------------
8 License
9  This file is part of OpenFOAM.
10 
11  OpenFOAM is free software: you can redistribute it and/or modify it
12  under the terms of the GNU General Public License as published by
13  the Free Software Foundation, either version 3 of the License, or
14  (at your option) any later version.
15 
16  OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
17  ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
18  FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
19  for more details.
20 
21  You should have received a copy of the GNU General Public License
22  along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
23 
24 \*---------------------------------------------------------------------------*/
25 
26 #include "DimensionedScalarField.H"
27 
28 #define TEMPLATE template<class GeoMesh>
29 #include "DimensionedFieldFunctionsM.C"
30 
31 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
32 
33 namespace Foam
34 {
35 
36 // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
37 
38 template<class GeoMesh>
39 tmp<DimensionedField<scalar, GeoMesh> > stabilise
40 (
41  const DimensionedField<scalar, GeoMesh>& dsf,
42  const dimensioned<scalar>& ds
43 )
44 {
45  tmp<DimensionedField<scalar, GeoMesh> > tRes
46  (
47  new DimensionedField<scalar, GeoMesh>
48  (
49  IOobject
50  (
51  "stabilise(" + dsf.name() + ',' + ds.name() + ')',
52  dsf.instance(),
53  dsf.db()
54  ),
55  dsf.mesh(),
56  dsf.dimensions() + ds.dimensions()
57  )
58  );
59 
60  stabilise(tRes().field(), dsf.field(), ds.value());
61 
62  return tRes;
63 }
64 
65 
66 template<class GeoMesh>
67 tmp<DimensionedField<scalar, GeoMesh> > stabilise
68 (
69  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf,
70  const dimensioned<scalar>& ds
71 )
72 {
73  const DimensionedField<scalar, GeoMesh>& dsf = tdsf();
74 
75  tmp<DimensionedField<scalar, GeoMesh> > tRes =
76  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
77  (
78  tdsf,
79  "stabilise(" + dsf.name() + ',' + ds.name() + ')',
80  dsf.dimensions() + ds.dimensions()
81  );
82 
83  stabilise(tRes().field(), dsf.field(), ds.value());
84 
85  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf);
86 
87  return tRes;
88 }
89 
90 
91 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
92 
93 BINARY_TYPE_OPERATOR(scalar, scalar, scalar, +, '+', add)
94 BINARY_TYPE_OPERATOR(scalar, scalar, scalar, -, '-', subtract)
95 
96 BINARY_OPERATOR(scalar, scalar, scalar, *, '*', multiply)
97 BINARY_OPERATOR(scalar, scalar, scalar, /, '|', divide)
98 
99 BINARY_TYPE_OPERATOR_SF(scalar, scalar, scalar, /, '|', divide)
100 
101 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
102 
103 template<class GeoMesh>
104 tmp<DimensionedField<scalar, GeoMesh> > pow
105 (
106  const DimensionedField<scalar, GeoMesh>& dsf1,
107  const DimensionedField<scalar, GeoMesh>& dsf2
108 )
109 {
110  tmp<DimensionedField<scalar, GeoMesh> > tPow
111  (
112  new DimensionedField<scalar, GeoMesh>
113  (
114  IOobject
115  (
116  "pow(" + dsf1.name() + ',' + dsf2.name() + ')',
117  dsf1.instance(),
118  dsf1.db()
119  ),
120  dsf1.mesh(),
121  pow
122  (
123  dsf1.dimensions(),
124  dimensionedScalar("1", 1.0, dsf2.dimensions())
125  )
126  )
127  );
128 
129  pow(tPow().field(), dsf1.field(), dsf2.field());
130 
131  return tPow;
132 }
133 
134 
135 template<class GeoMesh>
136 tmp<DimensionedField<scalar, GeoMesh> > pow
137 (
138  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf1,
139  const DimensionedField<scalar, GeoMesh>& dsf2
140 )
141 {
142  const DimensionedField<scalar, GeoMesh>& dsf1 = tdsf1();
143 
144  tmp<DimensionedField<scalar, GeoMesh> > tPow =
145  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
146  (
147  tdsf1,
148  "pow(" + dsf1.name() + ',' + dsf2.name() + ')',
149  pow
150  (
151  dsf1.dimensions(),
152  dimensionedScalar("1", 1.0, dsf2.dimensions())
153  )
154  );
155 
156  pow(tPow().field(), dsf1.field(), dsf2.field());
157 
158  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf1);
159 
160  return tPow;
161 }
162 
163 
164 template<class GeoMesh>
165 tmp<DimensionedField<scalar, GeoMesh> > pow
166 (
167  const DimensionedField<scalar, GeoMesh>& dsf1,
168  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf2
169 )
170 {
171  const DimensionedField<scalar, GeoMesh>& dsf2 = tdsf2();
172 
173  tmp<DimensionedField<scalar, GeoMesh> > tPow =
174  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
175  (
176  tdsf2,
177  "pow(" + dsf1.name() + ',' + dsf2.name() + ')',
178  pow
179  (
180  dsf1.dimensions(),
181  dimensionedScalar("1", 1.0, dsf2.dimensions())
182  )
183  );
184 
185  pow(tPow().field(), dsf1.field(), dsf2.field());
186 
187  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf2);
188 
189  return tPow;
190 }
191 
192 template<class GeoMesh>
193 tmp<DimensionedField<scalar, GeoMesh> > pow
194 (
195  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf1,
196  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf2
197 )
198 {
199  const DimensionedField<scalar, GeoMesh>& dsf1 = tdsf1();
200  const DimensionedField<scalar, GeoMesh>& dsf2 = tdsf2();
201 
202  tmp<DimensionedField<scalar, GeoMesh> > tPow =
203  reuseTmpTmpDimensionedField<scalar, scalar, scalar, scalar, GeoMesh>::
204  New
205  (
206  tdsf1,
207  tdsf2,
208  "pow(" + dsf1.name() + ',' + dsf2.name() + ')',
209  pow
210  (
211  dsf1.dimensions(),
212  dimensionedScalar("1", 1.0, dsf2.dimensions())
213  )
214  );
215 
216  pow(tPow().field(), dsf1.field(), dsf2.field());
217 
218  reuseTmpTmpDimensionedField<scalar, scalar, scalar, scalar, GeoMesh>::clear
219  (
220  tdsf1,
221  tdsf2
222  );
223 
224  return tPow;
225 }
226 
227 
228 template<class GeoMesh>
229 tmp<DimensionedField<scalar, GeoMesh> > pow
230 (
231  const DimensionedField<scalar, GeoMesh>& dsf,
232  const dimensionedScalar& ds
233 )
234 {
235  tmp<DimensionedField<scalar, GeoMesh> > tPow
236  (
237  new DimensionedField<scalar, GeoMesh>
238  (
239  IOobject
240  (
241  "pow(" + dsf.name() + ',' + ds.name() + ')',
242  dsf.instance(),
243  dsf.db()
244  ),
245  dsf.mesh(),
246  pow(dsf.dimensions(), ds)
247  )
248  );
249 
250  pow(tPow().field(), dsf.field(), ds.value());
251 
252  return tPow;
253 }
254 
255 template<class GeoMesh>
256 tmp<DimensionedField<scalar, GeoMesh> > pow
257 (
258  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf,
259  const dimensionedScalar& ds
260 )
261 {
262  const DimensionedField<scalar, GeoMesh>& dsf = tdsf();
263 
264  tmp<DimensionedField<scalar, GeoMesh> > tPow =
265  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
266  (
267  tdsf,
268  "pow(" + dsf.name() + ',' + ds.name() + ')',
269  pow(dsf.dimensions(), ds)
270  );
271 
272  pow(tPow().field(), dsf.field(), ds.value());
273 
274  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf);
275 
276  return tPow;
277 }
278 
279 template<class GeoMesh>
280 tmp<DimensionedField<scalar, GeoMesh> > pow
281 (
282  const DimensionedField<scalar, GeoMesh>& dsf,
283  const scalar& s
284 )
285 {
286  return pow(dsf, dimensionedScalar(s));
287 }
288 
289 template<class GeoMesh>
290 tmp<DimensionedField<scalar, GeoMesh> > pow
291 (
292  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf,
293  const scalar& s
294 )
295 {
296  return pow(tdsf, dimensionedScalar(s));
297 }
298 
299 
300 template<class GeoMesh>
301 tmp<DimensionedField<scalar, GeoMesh> > pow
302 (
303  const dimensionedScalar& ds,
304  const DimensionedField<scalar, GeoMesh>& dsf
305 )
306 {
307  tmp<DimensionedField<scalar, GeoMesh> > tPow
308  (
309  new DimensionedField<scalar, GeoMesh>
310  (
311  IOobject
312  (
313  "pow(" + ds.name() + ',' + dsf.name() + ')',
314  dsf.instance(),
315  dsf.db()
316  ),
317  dsf.mesh(),
318  pow(ds, dsf.dimensions())
319  )
320  );
321 
322  pow(tPow().field(), ds.value(), dsf.field());
323 
324  return tPow;
325 }
326 
327 
328 template<class GeoMesh>
329 tmp<DimensionedField<scalar, GeoMesh> > pow
330 (
331  const dimensionedScalar& ds,
332  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf
333 )
334 {
335  const DimensionedField<scalar, GeoMesh>& dsf = tdsf();
336 
337  tmp<DimensionedField<scalar, GeoMesh> > tPow =
338  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
339  (
340  tdsf,
341  "pow(" + ds.name() + ',' + dsf.name() + ')',
342  pow(ds, dsf.dimensions())
343  );
344 
345  pow(tPow().field(), ds.value(), dsf.field());
346 
347  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf);
348 
349  return tPow;
350 }
351 
352 template<class GeoMesh>
353 tmp<DimensionedField<scalar, GeoMesh> > pow
354 (
355  const scalar& s,
356  const DimensionedField<scalar, GeoMesh>& dsf
357 )
358 {
359  return pow(dimensionedScalar(s), dsf);
360 }
361 
362 template<class GeoMesh>
363 tmp<DimensionedField<scalar, GeoMesh> > pow
364 (
365  const scalar& s,
366  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf
367 )
368 {
369  return pow(dimensionedScalar(s), tdsf);
370 }
371 
372 
373 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
374 
375 template<class GeoMesh>
376 tmp<DimensionedField<scalar, GeoMesh> > atan2
377 (
378  const DimensionedField<scalar, GeoMesh>& dsf1,
379  const DimensionedField<scalar, GeoMesh>& dsf2
380 )
381 {
382  tmp<DimensionedField<scalar, GeoMesh> > tAtan2
383  (
384  new DimensionedField<scalar, GeoMesh>
385  (
386  IOobject
387  (
388  "atan2(" + dsf1.name() + ',' + dsf2.name() + ')',
389  dsf1.instance(),
390  dsf1.db()
391  ),
392  dsf1.mesh(),
393  atan2(dsf1.dimensions(), dsf2.dimensions())
394  )
395  );
396 
397  atan2(tAtan2().field(), dsf1.field(), dsf2.field());
398 
399  return tAtan2;
400 }
401 
402 
403 template<class GeoMesh>
404 tmp<DimensionedField<scalar, GeoMesh> > atan2
405 (
406  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf1,
407  const DimensionedField<scalar, GeoMesh>& dsf2
408 )
409 {
410  const DimensionedField<scalar, GeoMesh>& dsf1 = tdsf1();
411 
412  tmp<DimensionedField<scalar, GeoMesh> > tAtan2 =
413  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
414  (
415  tdsf1,
416  "atan2(" + dsf1.name() + ',' + dsf2.name() + ')',
417  atan2(dsf1.dimensions(), dsf2.dimensions())
418  );
419 
420  atan2(tAtan2().field(), dsf1.field(), dsf2.field());
421 
422  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf1);
423 
424  return tAtan2;
425 }
426 
427 
428 template<class GeoMesh>
429 tmp<DimensionedField<scalar, GeoMesh> > atan2
430 (
431  const DimensionedField<scalar, GeoMesh>& dsf1,
432  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf2
433 )
434 {
435  const DimensionedField<scalar, GeoMesh>& dsf2 = tdsf2();
436 
437  tmp<DimensionedField<scalar, GeoMesh> > tAtan2 =
438  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
439  (
440  tdsf2,
441  "atan2(" + dsf1.name() + ',' + dsf2.name() + ')',
442  atan2(dsf1.dimensions(), dsf2.dimensions())
443  );
444 
445  atan2(tAtan2().field(), dsf1.field(), dsf2.field());
446 
447  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf2);
448 
449  return tAtan2;
450 }
451 
452 template<class GeoMesh>
453 tmp<DimensionedField<scalar, GeoMesh> > atan2
454 (
455  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf1,
456  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf2
457 )
458 {
459  const DimensionedField<scalar, GeoMesh>& dsf1 = tdsf1();
460  const DimensionedField<scalar, GeoMesh>& dsf2 = tdsf2();
461 
462  tmp<DimensionedField<scalar, GeoMesh> > tAtan2 =
463  reuseTmpTmpDimensionedField<scalar, scalar, scalar, scalar, GeoMesh>::
464  New
465  (
466  tdsf1,
467  tdsf2,
468  "atan2(" + dsf1.name() + ',' + dsf2.name() + ')',
469  atan2(dsf1.dimensions(), dsf2.dimensions())
470  );
471 
472  atan2(tAtan2().field(), dsf1.field(), dsf2.field());
473 
474  reuseTmpTmpDimensionedField<scalar, scalar, scalar, scalar, GeoMesh>::clear
475  (
476  tdsf1,
477  tdsf2
478  );
479 
480  return tAtan2;
481 }
482 
483 
484 template<class GeoMesh>
485 tmp<DimensionedField<scalar, GeoMesh> > atan2
486 (
487  const DimensionedField<scalar, GeoMesh>& dsf,
488  const dimensionedScalar& ds
489 )
490 {
491  tmp<DimensionedField<scalar, GeoMesh> > tAtan2
492  (
493  new DimensionedField<scalar, GeoMesh>
494  (
495  IOobject
496  (
497  "atan2(" + dsf.name() + ',' + ds.name() + ')',
498  dsf.instance(),
499  dsf.db()
500  ),
501  dsf.mesh(),
502  atan2(dsf.dimensions(), ds)
503  )
504  );
505 
506  atan2(tAtan2().field(), dsf.field(), ds.value());
507 
508  return tAtan2;
509 }
510 
511 template<class GeoMesh>
512 tmp<DimensionedField<scalar, GeoMesh> > atan2
513 (
514  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf,
515  const dimensionedScalar& ds
516 )
517 {
518  const DimensionedField<scalar, GeoMesh>& dsf = tdsf();
519 
520  tmp<DimensionedField<scalar, GeoMesh> > tAtan2 =
521  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
522  (
523  tdsf,
524  "atan2(" + dsf.name() + ',' + ds.name() + ')',
525  atan2(dsf.dimensions(), ds)
526  );
527 
528  atan2(tAtan2().field(), dsf.field(), ds.value());
529 
530  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf);
531 
532  return tAtan2;
533 }
534 
535 template<class GeoMesh>
536 tmp<DimensionedField<scalar, GeoMesh> > atan2
537 (
538  const DimensionedField<scalar, GeoMesh>& dsf,
539  const scalar& s
540 )
541 {
542  return atan2(dsf, dimensionedScalar(s));
543 }
544 
545 template<class GeoMesh>
546 tmp<DimensionedField<scalar, GeoMesh> > atan2
547 (
548  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf,
549  const scalar& s
550 )
551 {
552  return atan2(tdsf, dimensionedScalar(s));
553 }
554 
555 
556 template<class GeoMesh>
557 tmp<DimensionedField<scalar, GeoMesh> > atan2
558 (
559  const dimensionedScalar& ds,
560  const DimensionedField<scalar, GeoMesh>& dsf
561 )
562 {
563  tmp<DimensionedField<scalar, GeoMesh> > tAtan2
564  (
565  new DimensionedField<scalar, GeoMesh>
566  (
567  IOobject
568  (
569  "atan2(" + ds.name() + ',' + dsf.name() + ')',
570  dsf.instance(),
571  dsf.db()
572  ),
573  dsf.mesh(),
574  atan2(ds, dsf.dimensions())
575  )
576  );
577 
578  atan2(tAtan2().field(), ds.value(), dsf.field());
579 
580  return tAtan2;
581 }
582 
583 
584 template<class GeoMesh>
585 tmp<DimensionedField<scalar, GeoMesh> > atan2
586 (
587  const dimensionedScalar& ds,
588  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf
589 )
590 {
591  const DimensionedField<scalar, GeoMesh>& dsf = tdsf();
592 
593  tmp<DimensionedField<scalar, GeoMesh> > tAtan2 =
594  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New
595  (
596  tdsf,
597  "atan2(" + ds.name() + ',' + dsf.name() + ')',
598  atan2(ds, dsf.dimensions())
599  );
600 
601  atan2(tAtan2().field(), ds.value(), dsf.field());
602 
603  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf);
604 
605  return tAtan2;
606 }
607 
608 template<class GeoMesh>
609 tmp<DimensionedField<scalar, GeoMesh> > atan2
610 (
611  const scalar& s,
612  const DimensionedField<scalar, GeoMesh>& dsf
613 )
614 {
615  return atan2(dimensionedScalar(s), dsf);
616 }
617 
618 template<class GeoMesh>
619 tmp<DimensionedField<scalar, GeoMesh> > atan2
620 (
621  const scalar& s,
622  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf
623 )
624 {
625  return atan2(dimensionedScalar(s), tdsf);
626 }
627 
628 
629 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
630 
631 UNARY_FUNCTION(scalar, scalar, pow3, pow3)
632 UNARY_FUNCTION(scalar, scalar, pow4, pow4)
633 UNARY_FUNCTION(scalar, scalar, pow5, pow5)
634 UNARY_FUNCTION(scalar, scalar, pow6, pow6)
635 UNARY_FUNCTION(scalar, scalar, pow025, pow025)
636 UNARY_FUNCTION(scalar, scalar, sqrt, sqrt)
637 UNARY_FUNCTION(scalar, scalar, cbrt, cbrt)
638 UNARY_FUNCTION(scalar, scalar, sign, sign)
639 UNARY_FUNCTION(scalar, scalar, pos, pos)
640 UNARY_FUNCTION(scalar, scalar, neg, neg)
641 UNARY_FUNCTION(scalar, scalar, posPart, posPart)
642 UNARY_FUNCTION(scalar, scalar, negPart, negPart)
643 
644 UNARY_FUNCTION(scalar, scalar, exp, trans)
645 UNARY_FUNCTION(scalar, scalar, log, trans)
646 UNARY_FUNCTION(scalar, scalar, log10, trans)
647 UNARY_FUNCTION(scalar, scalar, sin, trans)
648 UNARY_FUNCTION(scalar, scalar, cos, trans)
649 UNARY_FUNCTION(scalar, scalar, tan, trans)
650 UNARY_FUNCTION(scalar, scalar, asin, trans)
651 UNARY_FUNCTION(scalar, scalar, acos, trans)
652 UNARY_FUNCTION(scalar, scalar, atan, trans)
653 UNARY_FUNCTION(scalar, scalar, sinh, trans)
654 UNARY_FUNCTION(scalar, scalar, cosh, trans)
655 UNARY_FUNCTION(scalar, scalar, tanh, trans)
656 UNARY_FUNCTION(scalar, scalar, asinh, trans)
657 UNARY_FUNCTION(scalar, scalar, acosh, trans)
658 UNARY_FUNCTION(scalar, scalar, atanh, trans)
659 UNARY_FUNCTION(scalar, scalar, erf, trans)
660 UNARY_FUNCTION(scalar, scalar, erfc, trans)
661 UNARY_FUNCTION(scalar, scalar, lgamma, trans)
662 UNARY_FUNCTION(scalar, scalar, j0, trans)
663 UNARY_FUNCTION(scalar, scalar, j1, trans)
664 UNARY_FUNCTION(scalar, scalar, y0, trans)
665 UNARY_FUNCTION(scalar, scalar, y1, trans)
666 
667 
668 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
669 
670 #define BesselFunc(func) \
671  \
672 template<class GeoMesh> \
673 tmp<DimensionedField<scalar, GeoMesh> > func \
674 ( \
675  const int n, \
676  const DimensionedField<scalar, GeoMesh>& dsf \
677 ) \
678 { \
679  if (!dsf.dimensions().dimensionless()) \
680  { \
681  FatalErrorInFunction \
682  << "dsf not dimensionless" \
683  << abort(FatalError); \
684  } \
685  \
686  tmp<DimensionedField<scalar, GeoMesh> > tFunc \
687  ( \
688  new DimensionedField<scalar, GeoMesh> \
689  ( \
690  IOobject \
691  ( \
692  #func "(" + name(n) + ',' + dsf.name() + ')', \
693  dsf.instance(), \
694  dsf.db() \
695  ), \
696  dsf.mesh(), \
697  dimless \
698  ) \
699  ); \
700  \
701  func(tFunc().field(), n, dsf.field()); \
702  \
703  return tFunc; \
704 } \
705  \
706 template<class GeoMesh> \
707 tmp<DimensionedField<scalar, GeoMesh> > func \
708 ( \
709  const int n, \
710  const tmp<DimensionedField<scalar, GeoMesh> >& tdsf \
711 ) \
712 { \
713  const DimensionedField<scalar, GeoMesh>& dsf = tdsf(); \
714  \
715  if (!dsf.dimensions().dimensionless()) \
716  { \
717  FatalErrorInFunction \
718  << " : dsf not dimensionless" \
719  << abort(FatalError); \
720  } \
721  \
722  tmp<DimensionedField<scalar, GeoMesh> > tFunc \
723  ( \
724  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::New \
725  ( \
726  tdsf, \
727  #func "(" + name(n) + ',' + dsf.name() + ')', \
728  dimless \
729  ) \
730  ); \
731  \
732  func(tFunc().field(), n, dsf.field()); \
733  \
734  reuseTmpDimensionedField<scalar, scalar, GeoMesh>::clear(tdsf); \
735  \
736  return tFunc; \
737 }
738 
739 BesselFunc(jn)
740 BesselFunc(yn)
741 
742 #undef BesselFunc
743 
744 
745 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
746 
747 } // End namespace Foam
748 
749 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
750 
751 #include "undefFieldFunctionsM.H"
752 
753 // ************************************************************************* //
BINARY_TYPE_OPERATOR_SF
#define BINARY_TYPE_OPERATOR_SF(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:531
BINARY_OPERATOR
#define BINARY_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:417
Foam::IOobject
IOobject defines the attributes of an object for which implicit objectRegistry management is supporte...
Definition: IOobject.H:91
Foam::subtract
void subtract(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:871
Foam::tan
dimensionedScalar tan(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:257
UNARY_FUNCTION
#define UNARY_FUNCTION(ReturnType, Type1, Func, Dfunc)
Definition: DimensionedFieldFunctionsM.C:30
Foam::cosh
dimensionedScalar cosh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:262
Foam::y1
dimensionedScalar y1(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:273
Foam::tmp
A class for managing temporary objects.
Definition: PtrList.H:118
Foam::DimensionedField::field
const Field< Type > & field() const
Definition: DimensionedFieldI.H:61
Foam::sin
dimensionedScalar sin(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:255
Foam::jn
dimensionedScalar jn(const int n, const dimensionedScalar &ds)
Definition: dimensionedScalar.C:296
Foam::posPart
dimensionedScalar posPart(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:212
Foam::atan2
dimensionedScalar atan2(const dimensionedScalar &x, const dimensionedScalar &y)
Definition: dimensionedScalar.C:303
BINARY_TYPE_OPERATOR
#define BINARY_TYPE_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:684
Foam::dimensioned::name
const word & name() const
Return const reference to name.
Definition: dimensionedType.C:235
undefFieldFunctionsM.H
Foam::dimensioned::value
const Type & value() const
Return const reference to value.
Definition: dimensionedType.C:261
Foam::exp
dimensionedScalar exp(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:252
Foam::sign
dimensionedScalar sign(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:179
Foam::erf
dimensionedScalar erf(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:267
DimensionedFieldFunctionsM.C
Foam::divide
void divide(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
Foam::pow025
dimensionedScalar pow025(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:131
Foam::atanh
dimensionedScalar atanh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:266
DimensionedScalarField.H
Scalar specific part of the implementation of DimensionedField.
BesselFunc
#define BesselFunc(func)
Definition: DimensionedScalarField.C:670
Foam::lgamma
dimensionedScalar lgamma(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:269
Foam::pow4
dimensionedScalar pow4(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:98
Foam::pow6
dimensionedScalar pow6(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:120
Foam::DimensionedField::mesh
const Mesh & mesh() const
Return mesh.
Definition: DimensionedFieldI.H:38
Foam::tanh
dimensionedScalar tanh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:263
Foam::pow3
dimensionedScalar pow3(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:87
Foam::stabilise
tmp< DimensionedField< scalar, GeoMesh > > stabilise(const DimensionedField< scalar, GeoMesh > &dsf, const dimensioned< scalar > &ds)
Definition: DimensionedScalarField.C:40
Foam::log10
dimensionedScalar log10(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:254
Foam::y0
dimensionedScalar y0(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:272
Foam::erfc
dimensionedScalar erfc(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:268
Foam::dimensionedScalar
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Definition: dimensionedScalarFwd.H:41
Foam::asinh
dimensionedScalar asinh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:264
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:73
Foam::add
void add(FieldField< Field1, typename typeOfSum< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:870
Foam::pow5
dimensionedScalar pow5(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:109
Foam::log
dimensionedScalar log(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:253
Foam::dimensioned< scalar >
Foam
Namespace for OpenFOAM.
Definition: combustionModel.C:30
s
gmvFile<< "tracers "<< particles.size()<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().x()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().position().y()<< " ";}gmvFile<< nl;forAllConstIter(Cloud< passiveParticle >, particles, iter){ gmvFile<< iter().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))
Foam::DimensionedField::dimensions
const dimensionSet & dimensions() const
Return dimensions.
Definition: DimensionedFieldI.H:46
Foam::yn
dimensionedScalar yn(const int n, const dimensionedScalar &ds)
Definition: dimensionedScalar.C:297
Foam::acosh
dimensionedScalar acosh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:265
Foam::negPart
dimensionedScalar negPart(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:223
Foam::sqrt
dimensionedScalar sqrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:142
Foam::j0
dimensionedScalar j0(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:270
Foam::acos
dimensionedScalar acos(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:259
Foam::atan
dimensionedScalar atan(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:260
Foam::reuseTmpTmpDimensionedField::clear
static void clear(const tmp< DimensionedField< Type1, GeoMesh > > &tdf1, const tmp< DimensionedField< Type2, GeoMesh > > &tdf2)
Definition: DimensionedFieldReuseFunctions.H:155
Foam::reuseTmpDimensionedField::clear
static void clear(const tmp< DimensionedField< Type1, GeoMesh > > &tdf1)
Definition: DimensionedFieldReuseFunctions.H:66
Foam::multiply
void multiply(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
Foam::dimensioned::dimensions
const dimensionSet & dimensions() const
Return const reference to dimensions.
Definition: dimensionedType.C:248
Foam::cbrt
dimensionedScalar cbrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:153
Foam::trans
dimensionSet trans(const dimensionSet &)
Definition: dimensionSet.C:438
Foam::j1
dimensionedScalar j1(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:271
Foam::neg
dimensionedScalar neg(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:201
Foam::asin
dimensionedScalar asin(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:258
Foam::reuseTmpTmpDimensionedField::New
static tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< Type1, GeoMesh > > &tdf1, const tmp< DimensionedField< Type2, GeoMesh > > &tdf2, const word &name, const dimensionSet &dimensions)
Definition: DimensionedFieldReuseFunctions.H:129
Foam::reuseTmpDimensionedField::New
static tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< Type1, GeoMesh > > &tdf1, const word &name, const dimensionSet &dimensions)
Definition: DimensionedFieldReuseFunctions.H:42
Foam::DimensionedField
Field with dimensions and associated with geometry type GeoMesh which is used to size the field and a...
Definition: DimensionedField.H:51
Foam::cos
dimensionedScalar cos(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:256
Foam::pos
dimensionedScalar pos(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:190
Foam::sinh
dimensionedScalar sinh(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:261

Copyright © 2011-2018  OpenFOAM | OPENFOAM® is a registered trademark of OpenCFD Ltd.

Generated by   doxygen 1.8.17