src/foam/primitives/Scalar/doubleFloat.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | foam-extend: Open Source CFD
   4    \\    /   O peration     | Version:     3.2
   5     \\  /    A nd           | Web:         http://www.foam-extend.org
   6      \\/     M anipulation  | For copyright notice see file Copyright
   7 -------------------------------------------------------------------------------
   8 License
   9     This file is part of foam-extend.
  10
  11     foam-extend is free software: you can redistribute it and/or modify it
  12     under the terms of the GNU General Public License as published by the
  13     Free Software Foundation, either version 3 of the License, or (at your
  14     option) any later version.
  15
  16     foam-extend is distributed in the hope that it will be useful, but
  17     WITHOUT ANY WARRANTY; without even the implied warranty of
  18     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  19     General Public License for more details.
  20
  21     You should have received a copy of the GNU General Public License
  22     along with foam-extend.  If not, see <http://www.gnu.org/licenses/>.
  23
  24 \*---------------------------------------------------------------------------*/
  25
  26 #ifndef doubleFloat_H
  27 #define doubleFloat_H
  28
  29 #include "label.H"
  30 #include "products.H"
  31
  32 #include <cmath>
  33
  34 #ifndef DUMMY_SCALAR_FUNCTIONS
  35 #define DUMMY_SCALAR_FUNCTIONS
  36 inline float j0f(float x) { return float(j0(double(x)));}
  37 inline float j1f(float x) { return float(j1(double(x)));}
  38 inline float y0f(float x) { return float(y0(double(x)));}
  39 inline float y1f(float x) { return float(y1(double(x)));}
  40 inline float jnf(const int n, const float s) { return float(jn(n, double(s))); }
  41 inline float ynf(const int n, const float s) { return float(yn(n, double(s))); }
  42 #endif
  43
  44 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  45
  46 namespace Foam
  47 {
  48
  49 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  50
  51 template<class Cmpt>
  52 class typeOfRank<Cmpt, 0>
  53 {
  54 public:
  55
  56     typedef Cmpt type;
  57 };
  58
  59
  60 template<class Cmpt>
  61 class symmTypeOfRank<Cmpt, 0>
  62 {
  63 public:
  64
  65     typedef Cmpt type;
  66 };
  67
  68
  69 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  70
  71 template<class T>
  72 inline bool equal(const T& s1, const T& s2)
  73 {
  74     return s1 == s2;
  75 }
  76
  77
  78 #define MAXMINPOW(retType, type1, type2)          \
  79                                                   \
  80 MAXMIN(retType, type1, type2)                     \
  81                                                   \
  82 inline double pow(const type1 s, const type2 e)   \
  83 {                                                 \
  84     return ::pow(double(s), double(e));           \
  85 }
  86
  87 MAXMINPOW(long double, long double, long double)
  88 MAXMINPOW(long double, long double, double)
  89 MAXMINPOW(long double, long double, float)
  90 MAXMINPOW(long double, long double, int)
  91 MAXMINPOW(long double, double, long double)
  92 MAXMINPOW(long double, float, long double)
  93 MAXMINPOW(long double, int, long double)
  94
  95 MAXMINPOW(double, double, double)
  96 MAXMINPOW(double, double, float)
  97 MAXMINPOW(double, float, double)
  98 MAXMINPOW(float, float, float)
  99 MAXMINPOW(double, double, int)
 100 MAXMINPOW(double, int, double)
 101 MAXMINPOW(double, double, long)
 102 MAXMINPOW(double, long, double)
 103 MAXMINPOW(float, float, int)
 104 MAXMINPOW(float, int, float)
 105 MAXMINPOW(float, float, long)
 106 MAXMINPOW(float, long, float)
 107
 108 #undef MAXMINPOW
 109
 110
 111 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 112
 113 } // End namespace Foam
 114
 115 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 116
 117 #endif
 118
 119 // ************************************************************************* //