libstdc++-v3/testsuite/26_numerics/headers/cmath/hypot.cc

   1 // Copyright (C) 2016-2018 Free Software Foundation, Inc.
   2 //
   3 // This file is part of the GNU ISO C++ Library.  This library is free
   4 // software; you can redistribute it and/or modify it under the
   5 // terms of the GNU General Public License as published by the
   6 // Free Software Foundation; either version 3, or (at your option)
   7 // any later version.
   8
   9 // This library is distributed in the hope that it will be useful,
  10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
  11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  12 // GNU General Public License for more details.
  13
  14 // You should have received a copy of the GNU General Public License along
  15 // with this library; see the file COPYING3.  If not see
  16 // <http://www.gnu.org/licenses/>.
  17
  18 // { dg-options "-std=gnu++17" }
  19 // { dg-do run { target c++17 } }
  20
  21 #include <cmath>
  22 #include <type_traits>
  23 #if defined(__TEST_DEBUG)
  24 #include <iostream>
  25 #define VERIFY(A) \
  26 if (!(A)) \
  27   { \
  28     std::cout << "line " << __LINE__ \
  29       << "  max_abs_frac = " << max_abs_frac \
  30       << "  tolerance = " << toler \
  31       << std::endl; \
  32   }
  33 #else
  34 #include <testsuite_hooks.h>
  35 #endif
  36
  37 using std::is_same_v;
  38 static_assert(is_same_v<double, decltype(std::hypot(0.0, 0.0, 0.0))>);
  39 static_assert(is_same_v<double, decltype(std::hypot(0.0f, 0.0, 0.0))>);
  40 static_assert(is_same_v<double, decltype(std::hypot(0.0, 0.0f, 0.0))>);
  41 static_assert(is_same_v<double, decltype(std::hypot(0.0, 0.0, 0.0f))>);
  42 static_assert(is_same_v<double, decltype(std::hypot(0.0f, 0.0f, 0.0))>);
  43 static_assert(is_same_v<double, decltype(std::hypot(0.0f, 0.0, 0))>);
  44 static_assert(is_same_v<long double, decltype(std::hypot(0.0f, 0.0, 0.0l))>);
  45 static_assert(is_same_v<long double, decltype(std::hypot(0, 0.0, 0.0l))>);
  46
  47 template<typename T> struct testcase_hypot { T x, y, z, f0; };
  48
  49 template<typename Tp, unsigned int Num>
  50   void
  51   test(const testcase_hypot<Tp> (&data)[Num], Tp toler)
  52   {
  53     const Tp eps = std::numeric_limits<Tp>::epsilon();
  54     Tp max_abs_diff = -Tp(1);
  55     Tp max_abs_frac = -Tp(1);
  56     unsigned int num_datum = Num;
  57     for (unsigned int i = 0; i < num_datum; ++i)
  58       {
  59         const Tp f = std::hypot(data[i].x, data[i].y, data[i].z);
  60         const Tp f0 = data[i].f0;
  61         const Tp diff = f - f0;
  62         if (std::abs(diff) > max_abs_diff)
  63           max_abs_diff = std::abs(diff);
  64         if (std::abs(f0) > Tp(10) * eps && std::abs(f) > Tp(10) * eps)
  65           {
  66             const Tp frac = diff / f0;
  67             if (std::abs(frac) > max_abs_frac)
  68               max_abs_frac = std::abs(frac);
  69           }
  70       }
  71     VERIFY(max_abs_frac < toler);
  72   }
  73
  74 const testcase_hypot<double> data1[] = {
  75   { 0.0, 0.0, 0.0, 0.0 },
  76   { 0.0, 1.0, 1.0, std::sqrt(2.0) },
  77   { 1.0, 1.0, 1.0, std::sqrt(3.0) },
  78   { 1.0, 2.0, 2.0, 3.0 },
  79   { 2.0, 3.0, 6.0, 7.0 },
  80   { 1.0, 4.0, 8.0, 9.0 },
  81   { 4.0, 4.0, 7.0, 9.0 },
  82   { 12.0, 16.0, 21.0, 29.0 },
  83   { 1e8, 1., 1e-8, 1e8 },
  84   { 1., 1e8, 1e-8, 1e8 },
  85   { 1e-8, 1., 1e8, 1e8 },
  86   { 1e-2, 1e-4, 1e-4, 0.01000099995 },
  87   { 214748364., 214748364., 214748364., 371955077.2902952 }
  88 };
  89 const double toler1 = 1e-12;
  90
  91 const testcase_hypot<float> data2[] = {
  92   { 0.0f, 0.0f, 0.0f, 0.0f },
  93   { 0.0f, 1.0f, 1.0f, std::sqrt(2.0f) },
  94   { 1.0f, 1.0f, 1.0f, std::sqrt(3.0f) },
  95   { 1.0f, 2.0f, 2.0f, 3.0f },
  96   { 2.0f, 3.0f, 6.0f, 7.0f },
  97   { 1.0f, 4.0f, 8.0f, 9.0f },
  98   { 4.0f, 4.0f, 7.0f, 9.0f },
  99   { 12.0f, 16.0f, 21.0f, 29.0f },
 100   { 1e8f, 1.f, 1e-8f, 1e8f },
 101   { 1.f, 1e8f, 1e-8f, 1e8f },
 102   { 1e-8f, 1.f, 1e8f, 1e8f },
 103   { 1e-2f, 1e-4f, 1e-4f, 0.010001f },
 104   { 214748364.f, 214748364.f, 214748364.f, 371955072.f }
 105 };
 106 const float toler2 = 1e-7f;
 107
 108 const testcase_hypot<long double> data3[] = {
 109   { 0.0l, 0.0l, 0.0l, 0.0l },
 110   { 0.0l, 1.0l, 1.0l, std::sqrt(2.0l) },
 111   { 1.0l, 1.0l, 1.0l, std::sqrt(3.0l) },
 112   { 1.0l, 2.0l, 2.0l, 3.0l },
 113   { 2.0l, 3.0l, 6.0l, 7.0l },
 114   { 1.0l, 4.0l, 8.0l, 9.0l },
 115   { 4.0l, 4.0l, 7.0l, 9.0l },
 116   { 12.0l, 16.0l, 21.0l, 29.0l },
 117   { 1e8l, 1.l, 1e-8l, 1e8l },
 118   { 1.l, 1e8l, 1e-8l, 1e8l },
 119   { 1e-8l, 1.l, 1e8l, 1e8l },
 120   { 1e-2l, 1e-4l, 1e-4l, 0.010000999950004999375l },
 121   { 2147483647.l, 2147483647.l, 2147483647.l, 3719550785.027307813987l }
 122 };
 123 const long double toler3 = 1e-16l;
 124
 125 void
 126 test01()
 127 {
 128   // See hypot-long-double.cc for this macro
 129 #ifndef TEST_HYPOT_LONG_DOUBLE
 130   test(data1, toler1);
 131   test(data2, toler2);
 132 #else
 133   test(data3, toler3);
 134 #endif
 135 }
 136
 137 int
 138 main()
 139 {
 140   test01();
 141 }