sysdeps/x86_64/fpu/test-libmvec-alias-mod.c

   1 /* Part of test to build shared library to ensure link against
   2    *_finite aliases from libmvec.
   3    Copyright (C) 2016 Free Software Foundation, Inc.
   4    This file is part of the GNU C Library.
   5
   6    The GNU C Library is free software; you can redistribute it and/or
   7    modify it under the terms of the GNU Lesser General Public
   8    License as published by the Free Software Foundation; either
   9    version 2.1 of the License, or (at your option) any later version.
  10
  11    The GNU C Library is distributed in the hope that it will be useful,
  12    but WITHOUT ANY WARRANTY; without even the implied warranty of
  13    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  14    Lesser General Public License for more details.
  15
  16    You should have received a copy of the GNU Lesser General Public
  17    License along with the GNU C Library; if not, see
  18    <http://www.gnu.org/licenses/>.  */
  19
  20 #define N 4000
  21 FLOAT log_arg[N];
  22 FLOAT exp_arg[N];
  23 FLOAT log_res[N];
  24 FLOAT exp_res[N];
  25 FLOAT pow_res[N];
  26 int arch_check = 1;
  27
  28 static void
  29 init_arg (void)
  30 {
  31   int i;
  32
  33   CHECK_ARCH_EXT;
  34
  35   arch_check = 0;
  36
  37   for (i = 0; i < N; i += 1)
  38     {
  39       log_arg[i] = 1.0;
  40       exp_arg[i] = 0.0;
  41     }
  42 }
  43
  44 int
  45 test_finite_alias (void)
  46 {
  47   int i;
  48
  49   init_arg ();
  50
  51   if (arch_check) return 77;
  52
  53 #pragma omp simd
  54   for (i = 0; i < N; i += 1)
  55     {
  56       log_res[i] = FUNC (log) (log_arg[i]);
  57       exp_res[i] = FUNC (exp) (exp_arg[i]);
  58       pow_res[i] = FUNC (pow) (log_arg[i], log_arg[i]);
  59     }
  60
  61   if (log_res[0] != 0.0) return 1;
  62   if (exp_res[0] != 1.0) return 1;
  63   if (pow_res[0] != 1.0) return 1;
  64
  65   return 0;
  66 }