diff options
| author | eschnett <eschnett@105869f7-3296-0410-a4ea-f4349344b45a> | 2012-04-02 22:00:40 +0000 |
|---|---|---|
| committer | eschnett <eschnett@105869f7-3296-0410-a4ea-f4349344b45a> | 2012-04-02 22:00:40 +0000 |
| commit | e07a351d694a99f412ecf505311e9d71dbb1332f (patch) | |
| tree | be3a2f0ccd69da5b94cb180c328db14180c211e7 /src/vectors-8-DoubleHummer.h | |
| parent | b21da6f8ea2a34554be2f87b27eef895c888ccf4 (diff) | |
Implement asin, sinh, asinh, and friends
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/Vectors/trunk@55 105869f7-3296-0410-a4ea-f4349344b45a
Diffstat (limited to 'src/vectors-8-DoubleHummer.h')
| -rw-r--r-- | src/vectors-8-DoubleHummer.h | 35 |
1 files changed, 27 insertions, 8 deletions
diff --git a/src/vectors-8-DoubleHummer.h b/src/vectors-8-DoubleHummer.h index fdc4be7..3366a0a 100644 --- a/src/vectors-8-DoubleHummer.h +++ b/src/vectors-8-DoubleHummer.h @@ -224,21 +224,40 @@ vec8_set(f(vec8_elt0(x)), \ f(vec8_elt1(x))); \ }) -#define K8REPL2(f,x_,a_) \ +#define K8REPL2S(f,x_,a_) \ ({ \ CCTK_REAL8_VEC const x__=(x_); \ - CCTK_REAL8_VEC const x=x__; \ CCTK_REAL8 const a__=(a_); \ + CCTK_REAL8_VEC const x=x__; \ CCTK_REAL8 const a=a__; \ vec8_set(f(vec8_elt0(x),a), \ f(vec8_elt1(x),a)); \ }) +#define K8REPL2(f,x_,y_) \ + ({ \ + CCTK_REAL8_VEC const x__=(x_); \ + CCTK_REAL8_VEC const y__=(y_); \ + CCTK_REAL8_VEC const x=x__; \ + CCTK_REAL8_VEC const y=y__; \ + vec8_set(f(vec8_elt0(x),vec8_elt0(y)), \ + f(vec8_elt1(x),vec8_elt1(y))); \ + }) -#define k8cos(x) K8REPL(cos,x) -#define k8exp(x) K8REPL(exp,x) -#define k8log(x) K8REPL(log,x) -#define k8pow(x,a) K8REPL2(pow,x,a) -#define k8sin(x) K8REPL(sin,x) -#define k8tan(x) K8REPL(tan,x) +#define k8acos(x) K8REPL(acos,x) +#define k8acosh(x) K8REPL(acosh,x) +#define k8asin(x) K8REPL(asin,x) +#define k8asinh(x) K8REPL(asinh,x) +#define k8atan(x) K8REPL(atan,x) +#define k8atan2(x,y) K8REPL2(atan2,x,y) +#define k8atanh(x) K8REPL(atanh,x) +#define k8cos(x) K8REPL(cos,x) +#define k8cosh(x) K8REPL(cosh,x) +#define k8exp(x) K8REPL(exp,x) +#define k8log(x) K8REPL(log,x) +#define k8pow(x,a) K8REPL2S(pow,x,a) +#define k8sin(x) K8REPL(sin,x) +#define k8sinh(x) K8REPL(sinh,x) +#define k8tan(x) K8REPL(tan,x) +#define k8tanh(x) K8REPL(tanh,x) #define k8ifmsb(x,y,z) (__fpsel(x,y,z)) |