1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
|
#include <errno.h>
#include "headers/truncd2.h"
#include "headers/tgammad2.h"
static __inline double _tgamma(double x)
{
double res;
vector double vx;
vector double truncx;
vector double vc = { 0.0, 0.0 };
vector unsigned long long cmpres;
vector signed int verrno, ferrno;
vector signed int fail = { EDOM, EDOM, EDOM, EDOM };
vx = spu_promote(x, 0);
res = spu_extract(_tgammad2(vx), 0);
#ifndef _IEEE_LIBM
/*
* use vector truncd2 rather than splat x, and splat truncx.
*/
truncx = _truncd2(vx);
cmpres = spu_cmpeq(truncx, vx);
verrno = spu_splats(errno);
ferrno = spu_sel(verrno, fail, (vector unsigned int) cmpres);
cmpres = spu_cmpgt(vc, vx);
errno = spu_extract(spu_sel(verrno, ferrno, (vector unsigned int) cmpres), 0);
#endif
return res;
}
|
