21#include "rheolef/compiler.h"
26 if (x ==
T(0))
return 1;
29 for (
size_t n =
static_cast<int>(x)-1;
n > 0;
n--)
34 static T sqrt_pi = sqrt(acos(
T(-1.)));
35 if (x == -
T(0.5))
return -2*sqrt_pi;
36 if (x-floor(x) ==
T(0.5)) {
38 for (
size_t n =
static_cast<int>(x);
n != 0;
n--)
42 std::cerr <<
"gamma: " << x <<
" is not of integer of half order" << std::endl;
This file is part of Rheolef.