[857] | 1 | /* Boost interval/detail/sparc_rounding_control.hpp file
|
---|
| 2 | *
|
---|
| 3 | * Copyright 2000 Jens Maurer
|
---|
| 4 | * Copyright 2002 Hervé Brönnimann, Guillaume Melquiond, Sylvain Pion
|
---|
| 5 | *
|
---|
| 6 | * Distributed under the Boost Software License, Version 1.0.
|
---|
| 7 | * (See accompanying file LICENSE_1_0.txt or
|
---|
| 8 | * copy at http://www.boost.org/LICENSE_1_0.txt)
|
---|
| 9 | *
|
---|
| 10 | * The basic code in this file was kindly provided by Jeremy Siek.
|
---|
| 11 | */
|
---|
| 12 |
|
---|
| 13 | #ifndef BOOST_NUMERIC_INTERVAL_DETAIL_SPARC_ROUNDING_CONTROL_HPP
|
---|
| 14 | #define BOOST_NUMERIC_INTERVAL_DETAIL_SPARC_ROUNDING_CONTROL_HPP
|
---|
| 15 |
|
---|
| 16 | #if !defined(sparc) && !defined(__sparc__)
|
---|
| 17 | # error This header is only intended for SPARC CPUs.
|
---|
| 18 | #endif
|
---|
| 19 |
|
---|
| 20 | #ifdef __SUNPRO_CC
|
---|
| 21 | # include <ieeefp.h>
|
---|
| 22 | #endif
|
---|
| 23 |
|
---|
| 24 |
|
---|
| 25 | namespace boost {
|
---|
| 26 | namespace numeric {
|
---|
| 27 | namespace interval_lib {
|
---|
| 28 | namespace detail {
|
---|
| 29 |
|
---|
| 30 | struct sparc_rounding_control
|
---|
| 31 | {
|
---|
| 32 | typedef unsigned int rounding_mode;
|
---|
| 33 |
|
---|
| 34 | static void set_rounding_mode(const rounding_mode& mode)
|
---|
| 35 | {
|
---|
| 36 | # if defined(__GNUC__)
|
---|
| 37 | __asm__ __volatile__("ld %0, %%fsr" : : "m"(mode));
|
---|
| 38 | # elif defined (__SUNPRO_CC)
|
---|
| 39 | fpsetround(fp_rnd(mode));
|
---|
| 40 | # elif defined(__KCC)
|
---|
| 41 | asm("sethi %hi(mode), %o1");
|
---|
| 42 | asm("ld [%o1+%lo(mode)], %fsr");
|
---|
| 43 | # else
|
---|
| 44 | # error Unsupported compiler for Sparc rounding control.
|
---|
| 45 | # endif
|
---|
| 46 | }
|
---|
| 47 |
|
---|
| 48 | static void get_rounding_mode(rounding_mode& mode)
|
---|
| 49 | {
|
---|
| 50 | # if defined(__GNUC__)
|
---|
| 51 | __asm__ __volatile__("st %%fsr, %0" : "=m"(mode));
|
---|
| 52 | # elif defined (__SUNPRO_CC)
|
---|
| 53 | mode = fpgetround();
|
---|
| 54 | # elif defined(__KCC)
|
---|
| 55 | # error KCC on Sun SPARC get_round_mode: please fix me
|
---|
| 56 | asm("st %fsr, [mode]");
|
---|
| 57 | # else
|
---|
| 58 | # error Unsupported compiler for Sparc rounding control.
|
---|
| 59 | # endif
|
---|
| 60 | }
|
---|
| 61 |
|
---|
| 62 | #if defined(__SUNPRO_CC)
|
---|
| 63 | static void downward() { set_rounding_mode(FP_RM); }
|
---|
| 64 | static void upward() { set_rounding_mode(FP_RP); }
|
---|
| 65 | static void to_nearest() { set_rounding_mode(FP_RN); }
|
---|
| 66 | static void toward_zero() { set_rounding_mode(FP_RZ); }
|
---|
| 67 | #else
|
---|
| 68 | static void downward() { set_rounding_mode(0xc0000000); }
|
---|
| 69 | static void upward() { set_rounding_mode(0x80000000); }
|
---|
| 70 | static void to_nearest() { set_rounding_mode(0x00000000); }
|
---|
| 71 | static void toward_zero() { set_rounding_mode(0x40000000); }
|
---|
| 72 | #endif
|
---|
| 73 | };
|
---|
| 74 |
|
---|
| 75 | } // namespace detail
|
---|
| 76 |
|
---|
| 77 | extern "C" {
|
---|
| 78 | float rintf(float);
|
---|
| 79 | double rint(double);
|
---|
| 80 | }
|
---|
| 81 |
|
---|
| 82 | template<>
|
---|
| 83 | struct rounding_control<float>:
|
---|
| 84 | detail::sparc_rounding_control
|
---|
| 85 | {
|
---|
| 86 | static const float& force_rounding(const float& x) { return x; }
|
---|
| 87 | static float to_int(const float& x) { return rintf(x); }
|
---|
| 88 | };
|
---|
| 89 |
|
---|
| 90 | template<>
|
---|
| 91 | struct rounding_control<double>:
|
---|
| 92 | detail::sparc_rounding_control
|
---|
| 93 | {
|
---|
| 94 | static const double& force_rounding(const double& x) { return x; }
|
---|
| 95 | static double to_int(const double& x) { return rint(x); }
|
---|
| 96 | };
|
---|
| 97 |
|
---|
| 98 | template<>
|
---|
| 99 | struct rounding_control<long double>:
|
---|
| 100 | detail::sparc_rounding_control
|
---|
| 101 | {
|
---|
| 102 | static const long double& force_rounding(const long double& x) { return x; }
|
---|
| 103 | static long double to_int(const long double& x) { return rint(x); }
|
---|
| 104 | };
|
---|
| 105 |
|
---|
| 106 | } // namespace interval_lib
|
---|
| 107 | } // namespace numeric
|
---|
| 108 | } // namespace boost
|
---|
| 109 |
|
---|
| 110 | #undef BOOST_NUMERIC_INTERVAL_NO_HARDWARE
|
---|
| 111 |
|
---|
| 112 | #endif /* BOOST_NUMERIC_INTERVAL_DETAIL_SPARC_ROUNDING_CONTROL_HPP */
|
---|