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.