gamma
double gamma(double);
y = gamma(y);
y = -x*gamma(-x);
y = gamma(x);
DFUNC_DBL(gamma)
extern void stbi_hdr_to_ldr_gamma(float gamma);
extern void stbi_ldr_to_hdr_gamma(float gamma);
void stbi_hdr_to_ldr_gamma(float gamma) { h2l_gamma_i = 1/gamma; }
void stbi_ldr_to_hdr_gamma(float gamma) { l2h_gamma = gamma; }
TEST(gamma(5.0));