This is the mail archive of the libstdc++@gcc.gnu.org mailing list for the libstdc++ project.


Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]
Other format: [Raw text]

18/numeric_limits fail on solaris 2.8



from include/std/std_limits.h:
#define __glibcpp_f128_min 3.362103143112093506262677817321752603E-4932L
#define __glibcpp_f128_max 1.189731495357231765085759326628007016E+4932L
define __glibcpp_f128_epsilon 1.925929944387235853055977942584927319E-34L

from usr/include/float.h:
/* Follows IEEE standards for 128-bit floating point */
#define LDBL_MANT_DIG   113
#define LDBL_EPSILON    1.925929944387235853055977942584927319E-34L
#define LDBL_DIG        33
#define LDBL_MIN_EXP    (-16381)
#define LDBL_MIN        3.362103143112093506262677817321752603E-4932L
#define LDBL_MIN_10_EXP (-4931)
#define LDBL_MAX_EXP    (+16384)
#define LDBL_MAX        1.189731495357231765085759326628007016E+4932L
#define LDBL_MAX_10_EXP (+4932)

#define DEFINE_EXTREMA(T, m, M) \
  template<> T extrema<T>::min = m; \
  template<> T extrema<T>::max = M

DEFINE_EXTREMA(long double, LDBL_MIN, LDBL_MAX);

the applicable parts of numeric_limits.cc:
  T limits_min = std::numeric_limits<T>::min();
  T limits_max = std::numeric_limits<T>::max();
  T extrema_min = extrema<T>::min;
  T extrema_max = extrema<T>::max;
  T epsilon = std::numeric_limits<T>::epsilon();

  VERIFY( (extrema_min - limits_min) < epsilon );
  VERIFY( (limits_min - extrema_min) < epsilon );
  VERIFY( (extrema_max - limits_max) < (1 + epsilon) );
  VERIFY( (limits_max - extrema_max) < (1 + epsilon) );

% a.out
Assertion failed: (limits_max - extrema_max) < (1 + epsilon), file 18_support/numeric_limits.cc, line 91
Abort (core dumped)

debugging on solaris:

(gdb) s
void test_extrema<long double>() () at 18_support/numeric_limits.cc:81
81        bool test = true;
(gdb) n
82        T limits_min = std::numeric_limits<T>::min();
(gdb) n
83        T limits_max = std::numeric_limits<T>::max();
(gdb) n
84        T extrema_min = extrema<T>::min;
(gdb) n
85        T extrema_max = extrema<T>::max;
(gdb) n
86        T epsilon = std::numeric_limits<T>::epsilon();
(gdb) n
88        VERIFY( (extrema_min - limits_min) < epsilon );
(gdb) p limits_min
$1 = 1.3906711615670008644313954483327525e-309
(gdb) p extrema_min
$2 = 1.1519648082658484611329186009243131e-19
(gdb) p extrema_max
$3 = 35740566642812252160
(gdb) p limits_max
$4 = nan(0xeffffffffffff)
(gdb) ptype limits_min
type = long double
(gdb) p sizeof(limits_min)
$5 = 16
(gdb) p epsilon
$6 = 0.01513671875

?

Thoughts?

-benjamin



Index Nav: [Date Index] [Subject Index] [Author Index] [Thread Index]
Message Nav: [Date Prev] [Date Next] [Thread Prev] [Thread Next]