FieldFieldFunctions.C
Go to the documentation of this file.
1 /*---------------------------------------------------------------------------*\
2  ========= |
3  \\ / F ield | OpenFOAM: The Open Source CFD Toolbox
4  \\ / O peration | Website: https://openfoam.org
5  \\ / A nd | Copyright (C) 2011-2020 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 "PstreamReduceOps.H"
27 #include "FieldFieldReuseFunctions.H"
28 
29 #define TEMPLATE template<template<class> class Field, class Type>
30 #include "FieldFieldFunctionsM.C"
31 
32 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
33 
34 namespace Foam
35 {
36 
37 /* * * * * * * * * * * * * * * * Global functions * * * * * * * * * * * * * */
38 
39 template<template<class> class Field, class Type>
40 void component
41 (
42  FieldField<Field, typename FieldField<Field, Type>::cmptType>& sf,
43  const FieldField<Field, Type>& f,
44  const direction d
45 )
46 {
47  forAll(sf, i)
48  {
49  component(sf[i], f[i], d);
50  }
51 }
52 
53 
54 template<template<class> class Field, class Type>
55 void T(FieldField<Field, Type>& f1, const FieldField<Field, Type>& f2)
56 {
57  forAll(f1, i)
58  {
59  T(f1[i], f2[i]);
60  }
61 }
62 
63 
64 template<template<class> class Field, class Type, direction r>
65 void pow
66 (
67  FieldField<Field, typename powProduct<Type, r>::type>& f,
68  const FieldField<Field, Type>& vf
69 )
70 {
71  forAll(f, i)
72  {
73  pow(f[i], vf[i]);
74  }
75 }
76 
77 template<template<class> class Field, class Type, direction r>
78 tmp<FieldField<Field, typename powProduct<Type, r>::type>>
79 pow
80 (
81  const FieldField<Field, Type>& f, typename powProduct<Type, r>::type
82 )
83 {
84  typedef typename powProduct<Type, r>::type powProductType;
85  tmp<FieldField<Field, powProductType>> tRes
86  (
87  FieldField<Field, powProductType>::NewCalculatedType(f)
88  );
89  pow<Type, r>(tRes.ref(), f);
90  return tRes;
91 }
92 
93 template<template<class> class Field, class Type, direction r>
94 tmp<FieldField<Field, typename powProduct<Type, r>::type>>
95 pow
96 (
97  const tmp<FieldField<Field, Type>>& tf, typename powProduct<Type, r>::type
98 )
99 {
100  typedef typename powProduct<Type, r>::type powProductType;
101  tmp<FieldField<Field, powProductType>> tRes
102  (
103  reuseTmpFieldField<Field, powProductType, Type>::New(tf)
104  );
105  pow<Type, r>(tRes.ref(), tf());
106  tf.clear();
107  return tRes;
108 }
109 
110 
111 template<template<class> class Field, class Type>
112 void sqr
113 (
114  FieldField<Field, typename outerProduct<Type, Type>::type>& f,
115  const FieldField<Field, Type>& vf
116 )
117 {
118  forAll(f, i)
119  {
120  sqr(f[i], vf[i]);
121  }
122 }
123 
124 template<template<class> class Field, class Type>
125 tmp<FieldField<Field, typename outerProduct<Type, Type>::type>>
126 sqr(const FieldField<Field, Type>& f)
127 {
128  typedef typename outerProduct<Type, Type>::type outerProductType;
129  tmp<FieldField<Field, outerProductType>> tRes
130  (
131  FieldField<Field, outerProductType>::NewCalculatedType(f)
132  );
133  sqr(tRes.ref(), f);
134  return tRes;
135 }
136 
137 template<template<class> class Field, class Type>
138 tmp<FieldField<Field, typename outerProduct<Type, Type>::type>>
139 sqr(const tmp<FieldField<Field, Type>>& tf)
140 {
141  typedef typename outerProduct<Type, Type>::type outerProductType;
142  tmp<FieldField<Field, outerProductType>> tRes
143  (
144  reuseTmpFieldField<Field, outerProductType, Type>::New(tf)
145  );
146  sqr(tRes.ref(), tf());
147  tf.clear();
148  return tRes;
149 }
150 
151 
152 template<template<class> class Field, class Type>
153 void magSqr(FieldField<Field, scalar>& sf, const FieldField<Field, Type>& f)
154 {
155  forAll(sf, i)
156  {
157  magSqr(sf[i], f[i]);
158  }
159 }
160 
161 template<template<class> class Field, class Type>
162 tmp<FieldField<Field, scalar>> magSqr(const FieldField<Field, Type>& f)
163 {
164  tmp<FieldField<Field, scalar>> tRes
165  (
166  FieldField<Field, scalar>::NewCalculatedType(f)
167  );
168 
169  magSqr(tRes.ref(), f);
170  return tRes;
171 }
172 
173 template<template<class> class Field, class Type>
174 tmp<FieldField<Field, scalar>> magSqr(const tmp<FieldField<Field, Type>>& tf)
175 {
176  tmp<FieldField<Field, scalar>> tRes
177  (
178  reuseTmpFieldField<Field, scalar, Type>::New(tf)
179  );
180 
181  magSqr(tRes.ref(), tf());
182  tf.clear();
183  return tRes;
184 }
185 
186 
187 template<template<class> class Field, class Type>
188 void mag(FieldField<Field, scalar>& sf, const FieldField<Field, Type>& f)
189 {
190  forAll(sf, i)
191  {
192  mag(sf[i], f[i]);
193  }
194 }
195 
196 template<template<class> class Field, class Type>
197 tmp<FieldField<Field, scalar>> mag(const FieldField<Field, Type>& f)
198 {
199  tmp<FieldField<Field, scalar>> tRes
200  (
201  FieldField<Field, scalar>::NewCalculatedType(f)
202  );
203 
204  mag(tRes.ref(), f);
205  return tRes;
206 }
207 
208 template<template<class> class Field, class Type>
209 tmp<FieldField<Field, scalar>> mag(const tmp<FieldField<Field, Type>>& tf)
210 {
211  tmp<FieldField<Field, scalar>> tRes
212  (
213  reuseTmpFieldField<Field, scalar, Type>::New(tf)
214  );
215 
216  mag(tRes.ref(), tf());
217  tf.clear();
218  return tRes;
219 }
220 
221 
222 template<template<class> class Field, class Type>
223 void cmptMax
224 (
225  FieldField<Field, typename FieldField<Field, Type>::cmptType>& cf,
226  const FieldField<Field, Type>& f
227 )
228 {
229  forAll(cf, i)
230  {
231  cmptMax(cf[i], f[i]);
232  }
233 }
234 
235 template<template<class> class Field, class Type>
236 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptMax
237 (
238  const FieldField<Field, Type>& f
239 )
240 {
241  typedef typename FieldField<Field, Type>::cmptType cmptType;
242  tmp<FieldField<Field, cmptType>> tRes
243  (
244  FieldField<Field, cmptType>::NewCalculatedType(f)
245  );
246  cmptMax(tRes.ref(), f);
247  return tRes;
248 }
249 
250 template<template<class> class Field, class Type>
251 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptMax
252 (
253  const tmp<FieldField<Field, Type>>& tf
254 )
255 {
256  typedef typename FieldField<Field, Type>::cmptType cmptType;
257  tmp<FieldField<Field, cmptType>> tRes
258  (
259  reuseTmpFieldField<Field, cmptType, Type>::New(tf)
260  );
261  cmptMax(tRes.ref(), tf());
262  tf.clear();
263  return tRes;
264 }
265 
266 
267 template<template<class> class Field, class Type>
268 void cmptMin
269 (
270  FieldField<Field, typename FieldField<Field, Type>::cmptType>& cf,
271  const FieldField<Field, Type>& f
272 )
273 {
274  forAll(cf, i)
275  {
276  cmptMin(cf[i], f[i]);
277  }
278 }
279 
280 template<template<class> class Field, class Type>
281 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptMin
282 (
283  const FieldField<Field, Type>& f
284 )
285 {
286  typedef typename FieldField<Field, Type>::cmptType cmptType;
287  tmp<FieldField<Field, cmptType>> tRes
288  (
289  FieldField<Field, cmptType>::NewCalculatedType(f)
290  );
291  cmptMin(tRes.ref(), f);
292  return tRes;
293 }
294 
295 template<template<class> class Field, class Type>
296 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptMin
297 (
298  const tmp<FieldField<Field, Type>>& tf
299 )
300 {
301  typedef typename FieldField<Field, Type>::cmptType cmptType;
302  tmp<FieldField<Field, cmptType>> tRes
303  (
304  reuseTmpFieldField<Field, cmptType, Type>::New(tf)
305  );
306  cmptMin(tRes.ref(), tf());
307  tf.clear();
308  return tRes;
309 }
310 
311 
312 template<template<class> class Field, class Type>
313 void cmptAv
314 (
315  FieldField<Field, typename FieldField<Field, Type>::cmptType>& cf,
316  const FieldField<Field, Type>& f
317 )
318 {
319  forAll(cf, i)
320  {
321  cmptAv(cf[i], f[i]);
322  }
323 }
324 
325 template<template<class> class Field, class Type>
326 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptAv
327 (
328  const FieldField<Field, Type>& f
329 )
330 {
331  typedef typename FieldField<Field, Type>::cmptType cmptType;
332  tmp<FieldField<Field, cmptType>> tRes
333  (
334  FieldField<Field, cmptType>::NewCalculatedType(f)
335  );
336  cmptAv(tRes.ref(), f);
337  return tRes;
338 }
339 
340 template<template<class> class Field, class Type>
341 tmp<FieldField<Field, typename FieldField<Field, Type>::cmptType>> cmptAv
342 (
343  const tmp<FieldField<Field, Type>>& tf
344 )
345 {
346  typedef typename FieldField<Field, Type>::cmptType cmptType;
347  tmp<FieldField<Field, cmptType>> tRes
348  (
349  reuseTmpFieldField<Field, cmptType, Type>::New(tf)
350  );
351  cmptAv(tRes.ref(), tf());
352  tf.clear();
353  return tRes;
354 }
355 
356 
357 template<template<class> class Field, class Type>
358 void cmptMag
359 (
360  FieldField<Field, Type>& cf,
361  const FieldField<Field, Type>& f
362 )
363 {
364  forAll(cf, i)
365  {
366  cmptMag(cf[i], f[i]);
367  }
368 }
369 
370 template<template<class> class Field, class Type>
371 tmp<FieldField<Field, Type>> cmptMag
372 (
373  const FieldField<Field, Type>& f
374 )
375 {
376  tmp<FieldField<Field, Type>> tRes
377  (
378  FieldField<Field, Type>::NewCalculatedType(f)
379  );
380  cmptMag(tRes.ref(), f);
381  return tRes;
382 }
383 
384 template<template<class> class Field, class Type>
385 tmp<FieldField<Field, Type>> cmptMag
386 (
387  const tmp<FieldField<Field, Type>>& tf
388 )
389 {
390  tmp<FieldField<Field, Type>> tRes(New(tf));
391  cmptMag(tRes.ref(), tf());
392  tf.clear();
393  return tRes;
394 }
395 
396 
397 #define TMP_UNARY_FUNCTION(returnType, func) \
398  \
399 template<template<class> class Field, class Type> \
400 returnType func(const tmp<FieldField<Field, Type>>& tf1) \
401 { \
402  returnType res = func(tf1()); \
403  tf1.clear(); \
404  return res; \
405 }
406 
407 template<template<class> class Field, class Type>
408 Type max(const FieldField<Field, Type>& f)
409 {
410  label i = 0;
411  while (i < f.size() && !f[i].size()) i++;
412 
413  if (i < f.size())
414  {
415  Type Max(max(f[i]));
416 
417  for (label j=i+1; j<f.size(); j++)
418  {
419  if (f[j].size())
420  {
421  Max = max(max(f[j]), Max);
422  }
423  }
424 
425  return Max;
426  }
427  else
428  {
429  return pTraits<Type>::min;
430  }
431 }
432 
433 TMP_UNARY_FUNCTION(Type, max)
434 
435 template<template<class> class Field, class Type>
436 Type min(const FieldField<Field, Type>& f)
437 {
438  label i = 0;
439  while (i < f.size() && !f[i].size()) i++;
440 
441  if (i < f.size())
442  {
443  label i = 0;
444  while (!f[i].size()) i++;
445 
446  Type Min(min(f[i]));
447 
448  for (label j=i+1; j<f.size(); j++)
449  {
450  if (f[j].size())
451  {
452  Min = min(min(f[j]), Min);
453  }
454  }
455 
456  return Min;
457  }
458  else
459  {
460  return pTraits<Type>::max;
461  }
462 }
463 
464 TMP_UNARY_FUNCTION(Type, min)
465 
466 template<template<class> class Field, class Type>
467 Type sum(const FieldField<Field, Type>& f)
468 {
469  if (f.size())
470  {
471  Type Sum = Zero;
472 
473  forAll(f, i)
474  {
475  Sum += sum(f[i]);
476  }
477 
478  return Sum;
479  }
480  else
481  {
482  return Zero;
483  }
484 }
485 
486 TMP_UNARY_FUNCTION(Type, sum)
487 
488 template<template<class> class Field, class Type>
489 scalar sumMag(const FieldField<Field, Type>& f)
490 {
491  if (f.size())
492  {
493  scalar SumMag = 0.0;
494 
495  forAll(f, i)
496  {
497  SumMag += sumMag(f[i]);
498  }
499 
500  return SumMag;
501  }
502  else
503  {
504  return 0.0;
505  }
506 }
507 
508 TMP_UNARY_FUNCTION(scalar, sumMag)
509 
510 template<template<class> class Field, class Type>
511 Type average(const FieldField<Field, Type>& f)
512 {
513  if (f.size())
514  {
515  label n = 0;
516 
517  forAll(f, i)
518  {
519  n += f[i].size();
520  }
521 
522  if (n == 0)
523  {
524  WarningInFunction
525  << "empty fieldField, returning zero" << endl;
526 
527  return Zero;
528  }
529 
530  Type avrg = sum(f)/n;
531 
532  return avrg;
533  }
534  else
535  {
536  WarningInFunction
537  << "empty fieldField, returning zero" << endl;
538 
539  return Zero;
540  }
541 }
542 
543 TMP_UNARY_FUNCTION(Type, average)
544 
545 
546 #define G_UNARY_FUNCTION(returnType, gFunc, func, rFunc) \
547  \
548 template<template<class> class Field, class Type> \
549 returnType gFunc(const FieldField<Field, Type>& f) \
550 { \
551  returnType res = func(f); \
552  reduce(res, rFunc##Op<Type>()); \
553  return res; \
554 } \
555 TMP_UNARY_FUNCTION(returnType, gFunc)
556 
557 G_UNARY_FUNCTION(Type, gMax, max, max)
558 G_UNARY_FUNCTION(Type, gMin, min, min)
559 G_UNARY_FUNCTION(Type, gSum, sum, sum)
560 G_UNARY_FUNCTION(scalar, gSumMag, sumMag, sum)
561 
562 #undef G_UNARY_FUNCTION
563 
564 
565 template<template<class> class Field, class Type>
566 Type gAverage(const FieldField<Field, Type>& f)
567 {
568  label n = 0;
569 
570  forAll(f, i)
571  {
572  n += f[i].size();
573  }
574 
575  reduce(n, sumOp<label>());
576 
577  if (n > 0)
578  {
579  Type avrg = gSum(f)/n;
580 
581  return avrg;
582  }
583  else
584  {
585  WarningInFunction
586  << "empty fieldField, returning zero" << endl;
587 
588  return Zero;
589  }
590 }
591 
592 TMP_UNARY_FUNCTION(Type, gAverage)
593 
594 #undef TMP_UNARY_FUNCTION
595 
596 
597 BINARY_FUNCTION(Type, Type, Type, max)
598 BINARY_FUNCTION(Type, Type, Type, min)
599 BINARY_FUNCTION(Type, Type, Type, cmptMultiply)
600 BINARY_FUNCTION(Type, Type, Type, cmptDivide)
601 
602 BINARY_TYPE_FUNCTION(Type, Type, Type, max)
603 BINARY_TYPE_FUNCTION(Type, Type, Type, min)
604 BINARY_TYPE_FUNCTION(Type, Type, Type, cmptMultiply)
605 BINARY_TYPE_FUNCTION(Type, Type, Type, cmptDivide)
606 
607 
608 /* * * * * * * * * * * * * * * * Global operators * * * * * * * * * * * * * */
609 
610 UNARY_OPERATOR(Type, Type, -, negate)
611 
612 BINARY_OPERATOR(Type, Type, scalar, *, multiply)
613 BINARY_OPERATOR(Type, scalar, Type, *, multiply)
614 BINARY_OPERATOR(Type, Type, scalar, /, divide)
615 
616 BINARY_TYPE_OPERATOR_SF(Type, scalar, Type, *, multiply)
617 BINARY_TYPE_OPERATOR_FS(Type, Type, scalar, *, multiply)
618 
619 BINARY_TYPE_OPERATOR_FS(Type, Type, scalar, /, divide)
620 
621 
622 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
623 
624 #define PRODUCT_OPERATOR(product, op, opFunc) \
625  \
626 template \
627 < \
628  template<class> class Field1, \
629  template<class> class Field2, \
630  class Type1, \
631  class Type2 \
632 > \
633 void opFunc \
634 ( \
635  FieldField<Field1, typename product<Type1, Type2>::type>& f, \
636  const FieldField<Field1, Type1>& f1, \
637  const FieldField<Field2, Type2>& f2 \
638 ) \
639 { \
640  forAll(f, i) \
641  { \
642  opFunc(f[i], f1[i], f2[i]); \
643  } \
644 } \
645  \
646 template \
647 < \
648  template<class> class Field1, \
649  template<class> class Field2, \
650  class Type1, \
651  class Type2 \
652 > \
653 tmp<FieldField<Field1, typename product<Type1, Type2>::type>> \
654 operator op \
655 ( \
656  const FieldField<Field1, Type1>& f1, \
657  const FieldField<Field2, Type2>& f2 \
658 ) \
659 { \
660  typedef typename product<Type1, Type2>::type productType; \
661  tmp<FieldField<Field1, productType>> tRes \
662  ( \
663  FieldField<Field1, productType>::NewCalculatedType(f1) \
664  ); \
665  opFunc(tRes.ref(), f1, f2); \
666  return tRes; \
667 } \
668  \
669 template<template<class> class Field, class Type1, class Type2> \
670 tmp<FieldField<Field, typename product<Type1, Type2>::type>> \
671 operator op \
672 ( \
673  const FieldField<Field, Type1>& f1, \
674  const tmp<FieldField<Field, Type2>>& tf2 \
675 ) \
676 { \
677  typedef typename product<Type1, Type2>::type productType; \
678  tmp<FieldField<Field, productType>> tRes \
679  ( \
680  reuseTmpFieldField<Field, productType, Type2>::New(tf2) \
681  ); \
682  opFunc(tRes.ref(), f1, tf2()); \
683  tf2.clear(); \
684  return tRes; \
685 } \
686  \
687 template \
688 < \
689  template<class> class Field1, \
690  template<class> class Field2, \
691  class Type1, \
692  class Type2 \
693 > \
694 tmp<FieldField<Field, typename product<Type1, Type2>::type>> \
695 operator op \
696 ( \
697  const FieldField<Field1, Type1>& f1, \
698  const tmp<FieldField<Field2, Type2>>& tf2 \
699 ) \
700 { \
701  typedef typename product<Type1, Type2>::type productType; \
702  tmp<FieldField<Field1, productType>> tRes \
703  ( \
704  FieldField<Field1, productType>::NewCalculatedType(f1) \
705  ); \
706  opFunc(tRes.ref(), f1, tf2()); \
707  tf2.clear(); \
708  return tRes; \
709 } \
710  \
711 template \
712 < \
713  template<class> class Field1, \
714  template<class> class Field2, \
715  class Type1, \
716  class Type2 \
717 > \
718 tmp<FieldField<Field1, typename product<Type1, Type2>::type>> \
719 operator op \
720 ( \
721  const tmp<FieldField<Field1, Type1>>& tf1, \
722  const FieldField<Field2, Type2>& f2 \
723 ) \
724 { \
725  typedef typename product<Type1, Type2>::type productType; \
726  tmp<FieldField<Field1, productType>> tRes \
727  ( \
728  reuseTmpFieldField<Field1, productType, Type1>::New(tf1) \
729  ); \
730  opFunc(tRes.ref(), tf1(), f2); \
731  tf1.clear(); \
732  return tRes; \
733 } \
734  \
735 template \
736 < \
737  template<class> class Field1, \
738  template<class> class Field2, \
739  class Type1, \
740  class Type2 \
741 > \
742 tmp<FieldField<Field1, typename product<Type1, Type2>::type>> \
743 operator op \
744 ( \
745  const tmp<FieldField<Field1, Type1>>& tf1, \
746  const tmp<FieldField<Field2, Type2>>& tf2 \
747 ) \
748 { \
749  typedef typename product<Type1, Type2>::type productType; \
750  tmp<FieldField<Field1, productType>> tRes \
751  ( \
752  reuseTmpTmpFieldField<Field1, productType, Type1, Type2>::New(tf1, tf2)\
753  ); \
754  opFunc(tRes.ref(), tf1(), tf2()); \
755  tf1.clear(); \
756  tf2.clear(); \
757  return tRes; \
758 } \
759  \
760 template \
761 < \
762  template<class> class Field, \
763  class Type, \
764  class Form, \
765  class Cmpt, \
766  direction nCmpt \
767 > \
768 void opFunc \
769 ( \
770  FieldField<Field, typename product<Type, Form>::type>& f, \
771  const FieldField<Field, Type>& f1, \
772  const VectorSpace<Form,Cmpt,nCmpt>& vs \
773 ) \
774 { \
775  forAll(f, i) \
776  { \
777  opFunc(f[i], f1[i], vs); \
778  } \
779 } \
780  \
781 template \
782 < \
783  template<class> class Field, \
784  class Type, \
785  class Form, \
786  class Cmpt, \
787  direction nCmpt \
788 > \
789 tmp<FieldField<Field, typename product<Type, Form>::type>> \
790 operator op \
791 ( \
792  const FieldField<Field, Type>& f1, \
793  const VectorSpace<Form,Cmpt,nCmpt>& vs \
794 ) \
795 { \
796  typedef typename product<Type, Form>::type productType; \
797  tmp<FieldField<Field, productType>> tRes \
798  ( \
799  FieldField<Field, productType>::NewCalculatedType(f1) \
800  ); \
801  opFunc(tRes.ref(), f1, static_cast<const Form&>(vs)); \
802  return tRes; \
803 } \
804  \
805 template \
806 < \
807  template<class> class Field, \
808  class Type, \
809  class Form, \
810  class Cmpt, \
811  direction nCmpt \
812 > \
813 tmp<FieldField<Field, typename product<Type, Form>::type>> \
814 operator op \
815 ( \
816  const tmp<FieldField<Field, Type>>& tf1, \
817  const VectorSpace<Form,Cmpt,nCmpt>& vs \
818 ) \
819 { \
820  typedef typename product<Type, Form>::type productType; \
821  tmp<FieldField<Field, productType>> tRes \
822  ( \
823  reuseTmpFieldField<Field, productType, Type>::New(tf1) \
824  ); \
825  opFunc(tRes.ref(), tf1(), static_cast<const Form&>(vs)); \
826  tf1.clear(); \
827  return tRes; \
828 } \
829  \
830 template \
831 < \
832  template<class> class Field, \
833  class Type, \
834  class Form, \
835  class Cmpt, \
836  direction nCmpt \
837 > \
838 void opFunc \
839 ( \
840  FieldField<Field, typename product<Form, Type>::type>& f, \
841  const VectorSpace<Form,Cmpt,nCmpt>& vs, \
842  const FieldField<Field, Type>& f1 \
843 ) \
844 { \
845  forAll(f, i) \
846  { \
847  opFunc(f[i], vs, f1[i]); \
848  } \
849 } \
850  \
851 template \
852 < \
853  template<class> class Field, \
854  class Type, \
855  class Form, \
856  class Cmpt, \
857  direction nCmpt \
858 > \
859 tmp<FieldField<Field, typename product<Form, Type>::type>> \
860 operator op \
861 ( \
862  const VectorSpace<Form,Cmpt,nCmpt>& vs, \
863  const FieldField<Field, Type>& f1 \
864 ) \
865 { \
866  typedef typename product<Form, Type>::type productType; \
867  tmp<FieldField<Field, productType>> tRes \
868  ( \
869  FieldField<Field, productType>::NewCalculatedType(f1) \
870  ); \
871  opFunc(tRes.ref(), static_cast<const Form&>(vs), f1); \
872  return tRes; \
873 } \
874  \
875 template \
876 < \
877  template<class> class Field, \
878  class Type, \
879  class Form, \
880  class Cmpt, \
881  direction nCmpt \
882 > \
883 tmp<FieldField<Field, typename product<Form, Type>::type>> \
884 operator op \
885 ( \
886  const VectorSpace<Form,Cmpt,nCmpt>& vs, \
887  const tmp<FieldField<Field, Type>>& tf1 \
888 ) \
889 { \
890  typedef typename product<Form, Type>::type productType; \
891  tmp<FieldField<Field, productType>> tRes \
892  ( \
893  reuseTmpFieldField<Field, productType, Type>::New(tf1) \
894  ); \
895  opFunc(tRes.ref(), static_cast<const Form&>(vs), tf1()); \
896  tf1.clear(); \
897  return tRes; \
898 }
899 
900 PRODUCT_OPERATOR(typeOfSum, +, add)
901 PRODUCT_OPERATOR(typeOfSum, -, subtract)
902 
903 PRODUCT_OPERATOR(outerProduct, *, outer)
904 PRODUCT_OPERATOR(crossProduct, ^, cross)
905 PRODUCT_OPERATOR(innerProduct, &, dot)
906 PRODUCT_OPERATOR(scalarProduct, &&, dotdot)
907 
908 #undef PRODUCT_OPERATOR
909 
910 
911 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
912 
913 } // End namespace Foam
914 
915 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
916 
917 #include "undefFieldFunctionsM.H"
918 
919 // ************************************************************************* //
BINARY_FUNCTION
#define BINARY_FUNCTION(ReturnType, Type1, Type2, Func)
Definition: DimensionedFieldFunctionsM.C:132
UNARY_OPERATOR
#define UNARY_OPERATOR(ReturnType, Type1, Op, OpFunc, Dfunc)
Definition: DimensionedFieldFunctionsM.C:81
BINARY_OPERATOR
#define BINARY_OPERATOR(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:411
BINARY_TYPE_OPERATOR_SF
#define BINARY_TYPE_OPERATOR_SF(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:539
BINARY_TYPE_OPERATOR_FS
#define BINARY_TYPE_OPERATOR_FS(ReturnType, Type1, Type2, Op, OpName, OpFunc)
Definition: DimensionedFieldFunctionsM.C:621
BINARY_TYPE_FUNCTION
#define BINARY_TYPE_FUNCTION(ReturnType, Type1, Type2, Func)
Definition: DimensionedFieldFunctionsM.C:404
G_UNARY_FUNCTION
#define G_UNARY_FUNCTION(returnType, gFunc, func, rFunc)
Definition: FieldFieldFunctions.C:546
TMP_UNARY_FUNCTION
#define TMP_UNARY_FUNCTION(returnType, func)
Definition: FieldFieldFunctions.C:397
PRODUCT_OPERATOR
#define PRODUCT_OPERATOR(product, op, opFunc)
Definition: FieldFieldFunctions.C:624
PstreamReduceOps.H
Inter-processor communication reduction functions.
n
label n
Definition: TABSMDCalcMethod2.H:31
forAll
#define forAll(list, i)
Loop across all elements in list.
Definition: UList.H:434
Foam::FieldField
Generic field type.
Definition: FieldField.H:77
Foam::FieldField::cmptType
pTraits< Type >::cmptType cmptType
Component type.
Definition: FieldField.H:82
Foam::Field
Pre-declare SubField and related Field type.
Definition: Field.H:83
Foam::List::size
void size(const label)
Override size to be inconsistent with allocated storage.
Definition: ListI.H:164
Foam::List::clear
void clear()
Clear the list, i.e. set size to zero.
Definition: ListI.H:125
Foam::outerProduct::type
typeOfRank< typename pTraits< arg1 >::cmptType, direction(pTraits< arg1 >::rank)+direction(pTraits< arg2 >::rank) >::type type
Definition: products.H:90
Foam::pTraits
Traits class for primitives.
Definition: pTraits.H:53
Foam::powProduct::type
symmTypeOfRank< typename pTraits< arg1 >::cmptType, arg2 *direction(pTraits< arg1 >::rank) >::type type
Definition: products.H:136
Foam::tmp
A class for managing temporary objects.
Definition: tmp.H:55
Foam::tmp::ref
T & ref() const
Return non-const reference or generate a fatal error.
Definition: tmpI.H:181
sf
volScalarField sf(fieldObject, mesh)
tf
const tensorField & tf
Definition: getPatchFieldTensor.H:36
WarningInFunction
#define WarningInFunction
Report a warning using Foam::Warning.
Definition: messageStream.H:251
Foam
Namespace for OpenFOAM.
Definition: atmBoundaryLayer.H:214
Foam::negate
void negate(FieldField< Field, Type > &res, const FieldField< Field, Type > &f)
Foam::Zero
static const zero Zero
Definition: zero.H:97
Foam::gSum
Type gSum(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:559
Foam::label
intWM_LABEL_SIZE_t label
A label is an int32_t or int64_t as specified by the pre-processor macro WM_LABEL_SIZE.
Definition: label.H:59
Foam::gSumMag
scalar gSumMag(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:560
Foam::sqr
dimensionedSymmTensor sqr(const dimensionedVector &dv)
Definition: dimensionedSymmTensor.C:49
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:257
Foam::component
void component(FieldField< Field, typename FieldField< Field, Type >::cmptType > &sf, const FieldField< Field, Type > &f, const direction d)
Definition: FieldFieldFunctions.C:41
Foam::sum
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
Definition: DimensionedFieldFunctions.C:311
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:901
Foam::cmptMin
void cmptMin(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:269
Foam::average
dimensioned< Type > average(const DimensionedField< Type, GeoMesh > &df)
Definition: DimensionedFieldFunctions.C:313
Foam::dot
void dot(FieldField< Field1, typename innerProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:905
Foam::min
layerAndWeight min(const layerAndWeight &a, const layerAndWeight &b)
Definition: fvMeshStitchersMoving.C:110
Foam::divide
void divide(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
Foam::pow
dimensionedScalar pow(const dimensionedScalar &ds, const dimensionedScalar &expt)
Definition: dimensionedScalar.C:73
Foam::cmptAv
tmp< DimensionedField< typename DimensionedField< Type, GeoMesh >::cmptType, GeoMesh >> cmptAv(const DimensionedField< Type, GeoMesh > &df)
Definition: DimensionedFieldFunctions.C:232
Foam::reduce
void reduce(const List< UPstream::commsStruct > &comms, T &Value, const BinaryOp &bop, const int tag, const label comm)
Definition: PstreamReduceOps.H:49
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:900
Foam::mag
dimensioned< scalar > mag(const dimensioned< Type > &)
Foam::New
tmp< DimensionedField< TypeR, GeoMesh > > New(const tmp< DimensionedField< TypeR, GeoMesh >> &tdf1, const word &name, const dimensionSet &dimensions)
Definition: DimensionedFieldReuseFunctions.H:38
Foam::max
layerAndWeight max(const layerAndWeight &a, const layerAndWeight &b)
Definition: fvMeshStitchersMoving.C:104
Foam::cmptMultiply
dimensioned< Type > cmptMultiply(const dimensioned< Type > &, const dimensioned< Type > &)
Foam::cmptMag
dimensionSet cmptMag(const dimensionSet &)
Definition: dimensionSet.C:296
Foam::cmptMax
void cmptMax(FieldField< Field, typename FieldField< Field, Type >::cmptType > &cf, const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:224
Foam::cmptDivide
dimensioned< Type > cmptDivide(const dimensioned< Type > &, const dimensioned< Type > &)
Foam::gAverage
Type gAverage(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:566
Foam::dotdot
void dotdot(FieldField< Field1, typename scalarProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:906
Foam::gMin
Type gMin(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:558
Foam::cross
void cross(FieldField< Field1, typename crossProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:904
Foam::T
void T(FieldField< Field, Type > &f1, const FieldField< Field, Type > &f2)
Definition: FieldFieldFunctions.C:55
Foam::magSqr
dimensioned< scalar > magSqr(const dimensioned< Type > &)
Foam::multiply
void multiply(FieldField< Field, Type > &f, const FieldField< Field, Type > &f1, const FieldField< Field, scalar > &f2)
Foam::gMax
Type gMax(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:557
Foam::outer
void outer(FieldField< Field1, typename outerProduct< Type1, Type2 >::type > &f, const FieldField< Field1, Type1 > &f1, const FieldField< Field2, Type2 > &f2)
Definition: FieldFieldFunctions.C:903
Foam::direction
uint8_t direction
Definition: direction.H:45
Foam::sumMag
dimensioned< scalar > sumMag(const DimensionedField< Type, GeoMesh > &df)
Definition: DimensionedFieldFunctions.C:312
f
labelList f(nPoints)