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