src/OpenFOAM/primitives/Scalar/doubleFloat.H

   1 /*---------------------------------------------------------------------------*\
   2   =========                 |
   3   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   4    \\    /   O peration     |
   5     \\  /    A nd           | Copyright (C) 1991-2009 OpenCFD Ltd.
   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 the
  13     Free Software Foundation; either version 2 of the License, or (at your
  14     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, write to the Free Software Foundation,
  23     Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA
  24
  25 \*---------------------------------------------------------------------------*/
  26
  27 #ifndef doubleFloat_H
  28 #define doubleFloat_H
  29
  30 #include "label.H"
  31 #include "products.H"
  32
  33 #include <cmath>
  34
  35 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  36
  37 namespace Foam
  38 {
  39
  40 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  41
  42 template<class Cmpt>
  43 class typeOfRank<Cmpt, 0>
  44 {
  45 public:
  46
  47     typedef Cmpt type;
  48 };
  49
  50
  51 template<class Cmpt>
  52 class symmTypeOfRank<Cmpt, 0>
  53 {
  54 public:
  55
  56     typedef Cmpt type;
  57 };
  58
  59
  60 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  61
  62 template<class T>
  63 inline bool equal(const T& s1, const T& s2)
  64 {
  65     return s1 == s2;
  66 }
  67
  68
  69 #define MAXMINPOW(retType, type1, type2)          \
  70                                                   \
  71 MAXMIN(retType, type1, type2)                     \
  72                                                   \
  73 inline double pow(const type1 s, const type2 e)   \
  74 {                                                 \
  75     return ::pow(double(s), double(e));           \
  76 }
  77
  78
  79 MAXMINPOW(double, double, double)
  80 MAXMINPOW(double, double, float)
  81 MAXMINPOW(double, float, double)
  82 MAXMINPOW(float, float, float)
  83 MAXMINPOW(double, double, int)
  84 MAXMINPOW(double, int, double)
  85 MAXMINPOW(double, double, long)
  86 MAXMINPOW(double, long, double)
  87 MAXMINPOW(float, float, int)
  88 MAXMINPOW(float, int, float)
  89 MAXMINPOW(float, float, long)
  90 MAXMINPOW(float, long, float)
  91
  92 #undef MAXMINPOW
  93
  94
  95 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
  96
  97 } // End namespace Foam
  98
  99 // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
 100
 101 #endif
 102
 103 // ************************************************************************* //