diff --git a/math/Makefile b/math/Makefile index 8a89308..a2240d0 100644 --- a/math/Makefile +++ b/math/Makefile @@ -97,8 +97,9 @@ tests-static = test-fpucw-static test-fpucw-ieee-static test-longdouble-yes = test-ldouble test-ildoubl ifneq (no,$(PERL)) +libm-vec-tests = $(addprefix test-,$(libmvec-tests)) libm-tests = test-float test-double $(test-longdouble-$(long-double-fcts)) \ - test-ifloat test-idouble + test-ifloat test-idouble $(libm-vec-tests) libm-tests.o = $(addsuffix .o,$(libm-tests)) tests += $(libm-tests) @@ -125,8 +126,22 @@ $(objpfx)test-double.o: $(objpfx)libm-test.stmp $(objpfx)test-idouble.o: $(objpfx)libm-test.stmp $(objpfx)test-ldouble.o: $(objpfx)libm-test.stmp $(objpfx)test-ildoubl.o: $(objpfx)libm-test.stmp + +$(objpfx)test-double-vlen2.o: $(objpfx)libm-test.stmp +$(objpfx)test-double-vlen4.o: $(objpfx)libm-test.stmp + +$(objpfx)test-double-vlen2: $(common-objpfx)mathvec/libmvec.so \ + $(objpfx)init-arch.o +$(objpfx)test-double-vlen4: $(common-objpfx)mathvec/libmvec.so \ + $(objpfx)init-arch.o endif +CFLAGS-test-double-vlen2.c = -fno-inline -ffloat-store -fno-builtin -frounding-math \ + -D__FAST_MATH__ -DTEST_FAST_MATH -D_OPENMP=201307 \ + -Wno-unknown-pragmas +CFLAGS-test-double-vlen4.c = -fno-inline -ffloat-store -fno-builtin -frounding-math \ + -D__FAST_MATH__ -DTEST_FAST_MATH -D_OPENMP=201307 \ + -Wno-unknown-pragmas $(arch-ext-cflags) CFLAGS-test-float.c = -fno-inline -ffloat-store -fno-builtin -frounding-math CFLAGS-test-double.c = -fno-inline -ffloat-store -fno-builtin -frounding-math CFLAGS-test-ldouble.c = -fno-inline -ffloat-store -fno-builtin -frounding-math diff --git a/math/test-double-vlen2.h b/math/test-double-vlen2.h new file mode 100644 index 0000000..6ce47eb --- /dev/null +++ b/math/test-double-vlen2.h @@ -0,0 +1,42 @@ +/* Copyright (C) 2014 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#define FLOAT double +#define FUNC(function) function +#define TEST_MSG "testing double vector math (without inline functions)\n" +#define MATHCONST(x) x +#define CHOOSE(Clongdouble,Cdouble,Cfloat,Cinlinelongdouble,Cinlinedouble,Cinlinefloat) Cdouble +#define PRINTF_EXPR "e" +#define PRINTF_XEXPR "a" +#define PRINTF_NEXPR "f" +#define TEST_DOUBLE 1 +#define TEST_MATHVEC 1 + +#ifndef __NO_MATH_INLINES +# define __NO_MATH_INLINES +#endif + +#define EXCEPTION_TESTS_double 0 +#define ROUNDING_TESTS_double(MODE) ((MODE) == FE_TONEAREST) + +#define VEC_PREFIX vlen2_ + +#define CONCAT(prefix,func) __CONCAT (prefix,func) + +#define WRAPPER_NAME(function) CONCAT (VEC_PREFIX,function) + +#define FUNC_TEST(function) VEC_PREFIX_ ## function diff --git a/math/test-double-vlen4.h b/math/test-double-vlen4.h new file mode 100644 index 0000000..4e26f97 --- /dev/null +++ b/math/test-double-vlen4.h @@ -0,0 +1,40 @@ +/* Copyright (C) 2014 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#define FLOAT double +#define FUNC(function) function +#define TEST_MSG "testing double vector math (without inline functions)\n" +#define MATHCONST(x) x +#define CHOOSE(Clongdouble,Cdouble,Cfloat,Cinlinelongdouble,Cinlinedouble,Cinlinefloat) Cdouble +#define PRINTF_EXPR "e" +#define PRINTF_XEXPR "a" +#define PRINTF_NEXPR "f" +#define TEST_DOUBLE 1 +#define TEST_MATHVEC 1 + +#ifndef __NO_MATH_INLINES +# define __NO_MATH_INLINES +#endif + +#define EXCEPTION_TESTS_double 0 +#define ROUNDING_TESTS_double(MODE) ((MODE) == FE_TONEAREST) + +#define CONCAT(prefix,func) __CONCAT (prefix,func) + +#define WRAPPER_NAME(function) CONCAT (VEC_PREFIX,function) + +#define FUNC_TEST(function) VEC_PREFIX_ ## function diff --git a/sysdeps/x86_64/fpu/Makefile b/sysdeps/x86_64/fpu/Makefile index 25fe0d4..e6add6d 100644 --- a/sysdeps/x86_64/fpu/Makefile +++ b/sysdeps/x86_64/fpu/Makefile @@ -2,3 +2,20 @@ ifeq ($(subdir),mathvec) libmvec-support += svml_d_cos2_core svml_d_cos4_core_avx \ svml_d_cos4_core_avx2 svml_d_cos_data endif + +# Rules for libmvec tests +ifeq ($(subdir),math) +ifeq ($(build-mathvec),yes) +libmvec-tests += double-vlen2 double-vlen4 double-vlen4-avx2 +arch-ext-cflags = -mavx + +$(objpfx)test-double-vlen4-avx2.o: $(objpfx)libm-test.stmp + +$(objpfx)test-double-vlen4-avx2: $(common-objpfx)mathvec/libmvec.so \ + $(objpfx)init-arch.o + +CFLAGS-test-double-vlen4-avx2.c = -fno-inline -ffloat-store -fno-builtin -frounding-math \ + -D__FAST_MATH__ -DTEST_FAST_MATH -D_OPENMP=201307 \ + -Wno-unknown-pragmas -mavx2 +endif +endif diff --git a/sysdeps/x86_64/fpu/libm-test-ulps b/sysdeps/x86_64/fpu/libm-test-ulps index 36e1b76..b5c88d4 100644 --- a/sysdeps/x86_64/fpu/libm-test-ulps +++ b/sysdeps/x86_64/fpu/libm-test-ulps @@ -1961,6 +1961,15 @@ ifloat: 3 ildouble: 4 ldouble: 4 +Function: "vlen2_cos": +double: 1 + +Function: "vlen4_avx2_cos": +double: 1 + +Function: "vlen4_avx_cos": +double: 1 + Function: "y0": double: 2 float: 1 diff --git a/sysdeps/x86_64/fpu/test-double-vlen2.c b/sysdeps/x86_64/fpu/test-double-vlen2.c new file mode 100644 index 0000000..674c5de --- /dev/null +++ b/sysdeps/x86_64/fpu/test-double-vlen2.c @@ -0,0 +1,44 @@ +/* Tests for SSE4 ISA versions of vector math functions. + Copyright (C) 2014 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "test-double-vlen2.h" +#include + +// Wrapper from scalar to vector function implemented in SSE4. +#define VECTOR_WRAPPER(scalar_func,vector_func) \ +extern __m128d vector_func ( __m128d);\ +FLOAT scalar_func (FLOAT x)\ +{\ + int i;\ + __m128d mx = _mm_set1_pd (x);\ + __m128d mr = vector_func (mx);\ + for (i=1;i<2;i++)\ + {\ + if (((FLOAT*)&mr)[0]!=((FLOAT*)&mr)[i])\ + {\ + return ((FLOAT*)&mr)[0]+0.1;\ + }\ + }\ + return ((FLOAT*)&mr)[0];\ +} + +VECTOR_WRAPPER (WRAPPER_NAME (cos),_ZGVbN2v_cos) + +#define TEST_VECTOR_cos 1 + +#include "libm-test.c" diff --git a/sysdeps/x86_64/fpu/test-double-vlen4-avx2.c b/sysdeps/x86_64/fpu/test-double-vlen4-avx2.c new file mode 100644 index 0000000..15b7930 --- /dev/null +++ b/sysdeps/x86_64/fpu/test-double-vlen4-avx2.c @@ -0,0 +1,48 @@ +/* Tests for AVX2 ISA versions of vector math functions. + Copyright (C) 2014 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "test-double-vlen4.h" +#include + +// Wrapper from scalar to vector function implemented in AVX2. +#define VECTOR_WRAPPER(scalar_func,vector_func) \ +extern __m256d vector_func (__m256d);\ +FLOAT scalar_func (FLOAT x)\ +{\ + int i;\ + __m256d mx = _mm256_set1_pd (x);\ + __m256d mr = vector_func (mx);\ + for (i=1;i<4;i++)\ + {\ + if (((FLOAT*)&mr)[0]!=((FLOAT*)&mr)[i])\ + {\ + return ((FLOAT*)&mr)[0]+0.1;\ + }\ + }\ + return ((FLOAT*)&mr)[0];\ +} + +#define VEC_PREFIX vlen4_avx2_ + +VECTOR_WRAPPER (WRAPPER_NAME (cos),_ZGVdN4v_cos) + +#define TEST_VECTOR_cos 1 + +#define REQUIRE_AVX2 + +#include "libm-test.c" diff --git a/sysdeps/x86_64/fpu/test-double-vlen4.c b/sysdeps/x86_64/fpu/test-double-vlen4.c new file mode 100644 index 0000000..5f68af5 --- /dev/null +++ b/sysdeps/x86_64/fpu/test-double-vlen4.c @@ -0,0 +1,46 @@ +/* Tests for AVX ISA versions of vector math functions. + Copyright (C) 2014 Free Software Foundation, Inc. + This file is part of the GNU C Library. + + The GNU C Library is free software; you can redistribute it and/or + modify it under the terms of the GNU Lesser General Public + License as published by the Free Software Foundation; either + version 2.1 of the License, or (at your option) any later version. + + The GNU C Library is distributed in the hope that it will be useful, + but WITHOUT ANY WARRANTY; without even the implied warranty of + MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + Lesser General Public License for more details. + + You should have received a copy of the GNU Lesser General Public + License along with the GNU C Library; if not, see + . */ + +#include "test-double-vlen4.h" +#include + +// Wrapper from scalar to vector function implemented in AVX. +#define VECTOR_WRAPPER(scalar_func,vector_func) \ +extern __m256d vector_func (__m256d);\ +FLOAT scalar_func (FLOAT x)\ +{\ + int i;\ + __m256d mx = _mm256_set1_pd (x);\ + __m256d mr = vector_func (mx);\ + for (i=1;i<4;i++)\ + {\ + if (((FLOAT*)&mr)[0]!=((FLOAT*)&mr)[i])\ + {\ + return ((FLOAT*)&mr)[0]+0.1;\ + }\ + }\ + return ((FLOAT*)&mr)[0];\ +} + +#define VEC_PREFIX vlen4_avx_ + +VECTOR_WRAPPER (WRAPPER_NAME (cos),_ZGVcN4v_cos) + +#define TEST_VECTOR_cos 1 + +#include "libm-test.c"