Generated on Fri Oct 19 11:24:57 2018 for Gecode by doxygen 1.6.3

rounding.cpp

Go to the documentation of this file.
00001 /* -*- mode: C++; c-basic-offset: 2; indent-tabs-mode: nil -*- */
00002 /*
00003  *  Main authors:
00004  *     Vincent Barichard <Vincent.Barichard@univ-angers.fr>
00005  *
00006  *  Contributing authors:
00007  *     Christian Schulte <schulte@gecode.org>
00008  *
00009  *  Copyright:
00010  *     Christian Schulte, 2012
00011  *     Vincent Barichard, 2012
00012  *
00013  *  This file is part of Gecode, the generic constraint
00014  *  development environment:
00015  *     http://www.gecode.org
00016  *
00017  *  Permission is hereby granted, free of charge, to any person obtaining
00018  *  a copy of this software and associated documentation files (the
00019  *  "Software"), to deal in the Software without restriction, including
00020  *  without limitation the rights to use, copy, modify, merge, publish,
00021  *  distribute, sublicense, and/or sell copies of the Software, and to
00022  *  permit persons to whom the Software is furnished to do so, subject to
00023  *  the following conditions:
00024  *
00025  *  The above copyright notice and this permission notice shall be
00026  *  included in all copies or substantial portions of the Software.
00027  *
00028  *  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
00029  *  EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
00030  *  MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
00031  *  NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
00032  *  LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
00033  *  OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
00034  *  WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
00035  *
00036  */
00037 
00038 #include <gecode/float.hh>
00039 
00040 #ifdef GECODE_HAS_MPFR
00041 
00042 #include <gmp.h>
00043 #include <mpfr.h>
00044 
00045 namespace Gecode { namespace Float {
00046 
00048   typedef int mpfr_func(mpfr_t, const __mpfr_struct*, mp_rnd_t);
00049 
00051   forceinline double
00052   invoke_mpfr(FloatNum x, mpfr_func f, mp_rnd_t r) {
00053     mpfr_t xx;
00054     mpfr_init_set_d(xx, x, GMP_RNDN);
00055     f(xx, xx, r);
00056     FloatNum res = mpfr_get_d(xx, r);
00057     mpfr_clear(xx);
00058     return res;
00059   }
00060 
00062 #define GECODE_GENR_FUNC(name) \
00063   FloatNum Rounding::name##_down(FloatNum x) { \
00064     return invoke_mpfr(x, mpfr_##name, GMP_RNDD); \
00065   } \
00066   FloatNum Rounding::name##_up(FloatNum x) { \
00067     return invoke_mpfr(x, mpfr_##name, GMP_RNDU); \
00068   }
00069   GECODE_GENR_FUNC(exp)
00070   GECODE_GENR_FUNC(log)
00071   GECODE_GENR_FUNC(sin)
00072   GECODE_GENR_FUNC(cos)
00073   GECODE_GENR_FUNC(tan)
00074   GECODE_GENR_FUNC(asin)
00075   GECODE_GENR_FUNC(acos)
00076   GECODE_GENR_FUNC(atan)
00077   GECODE_GENR_FUNC(sinh)
00078   GECODE_GENR_FUNC(cosh)
00079   GECODE_GENR_FUNC(tanh)
00080   GECODE_GENR_FUNC(asinh)
00081   GECODE_GENR_FUNC(acosh)
00082   GECODE_GENR_FUNC(atanh)
00083 
00084 #undef GECODE_GENR_FUNC
00085 
00086 }}
00087 
00088 #endif
00089 
00090 // STATISTICS: float-var
00091