From e9009cdd21415eee6102875c45f6aa1d5ee50752 Mon Sep 17 00:00:00 2001 From: Aaron Watry Date: Thu, 15 Sep 2016 00:17:31 +0000 Subject: [PATCH] math: Implement lgamma Just use lgamma_r and ignore the value returned in the second argument Signed-off-by: Aaron Watry Reviewed-by: Tom Stellard llvm-svn: 281565 --- libclc/generic/include/clc/clc.h | 1 + libclc/generic/include/clc/math/lgamma.h | 2 ++ libclc/generic/include/clc/math/lgamma.inc | 1 + libclc/generic/lib/SOURCES | 1 + libclc/generic/lib/math/lgamma.cl | 44 ++++++++++++++++++++++++++++++ 5 files changed, 49 insertions(+) create mode 100644 libclc/generic/include/clc/math/lgamma.h create mode 100644 libclc/generic/include/clc/math/lgamma.inc create mode 100644 libclc/generic/lib/math/lgamma.cl diff --git a/libclc/generic/include/clc/clc.h b/libclc/generic/include/clc/clc.h index 934606c..fecbd59 100644 --- a/libclc/generic/include/clc/clc.h +++ b/libclc/generic/include/clc/clc.h @@ -69,6 +69,7 @@ #include #include #include +#include #include #include #include diff --git a/libclc/generic/include/clc/math/lgamma.h b/libclc/generic/include/clc/math/lgamma.h new file mode 100644 index 0000000..9181856 --- /dev/null +++ b/libclc/generic/include/clc/math/lgamma.h @@ -0,0 +1,2 @@ +#define __CLC_BODY +#include diff --git a/libclc/generic/include/clc/math/lgamma.inc b/libclc/generic/include/clc/math/lgamma.inc new file mode 100644 index 0000000..95e4079 --- /dev/null +++ b/libclc/generic/include/clc/math/lgamma.inc @@ -0,0 +1 @@ +_CLC_OVERLOAD _CLC_DECL __CLC_GENTYPE lgamma(__CLC_GENTYPE a); diff --git a/libclc/generic/lib/SOURCES b/libclc/generic/lib/SOURCES index 6ddae4d..bcd0991 100644 --- a/libclc/generic/lib/SOURCES +++ b/libclc/generic/lib/SOURCES @@ -97,6 +97,7 @@ math/hypot.cl math/ilogb.cl math/clc_ldexp.cl math/ldexp.cl +math/lgamma.cl math/lgamma_r.cl math/log.cl math/log10.cl diff --git a/libclc/generic/lib/math/lgamma.cl b/libclc/generic/lib/math/lgamma.cl new file mode 100644 index 0000000..26cd20e --- /dev/null +++ b/libclc/generic/lib/math/lgamma.cl @@ -0,0 +1,44 @@ +/* + * Copyright (c) 2016 Aaron Watry + * Copyright (c) 2014 Advanced Micro Devices, Inc. + * + * Permission is hereby granted, free of charge, to any person obtaining a copy + * of this software and associated documentation files (the "Software"), to deal + * in the Software without restriction, including without limitation the rights + * to use, copy, modify, merge, publish, distribute, sublicense, and/or sell + * copies of the Software, and to permit persons to whom the Software is + * furnished to do so, subject to the following conditions: + * + * The above copyright notice and this permission notice shall be included in + * all copies or substantial portions of the Software. + * + * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR + * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, + * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE + * AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER + * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, + * OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN + * THE SOFTWARE. + */ + +#include +#include "../clcmacro.h" + +_CLC_OVERLOAD _CLC_DEF float lgamma(float x) { + int s; + return lgamma_r(x, &s); +} + +_CLC_UNARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, float, lgamma, float) + +#ifdef cl_khr_fp64 +#pragma OPENCL EXTENSION cl_khr_fp64 : enable + +_CLC_OVERLOAD _CLC_DEF double lgamma(double x) { + int s; + return lgamma_r(x, &s); +} + +_CLC_UNARY_VECTORIZE(_CLC_OVERLOAD _CLC_DEF, double, lgamma, double) + +#endif \ No newline at end of file -- 2.7.4