Hi
I am having trouble understanding the documentation for
PDL::GSLSF::GAMMA. For example, consider gsl_sf_gamma_inc_Q
Normalized Incomplete Gamma Function
Q(a,x) = 1/Gamma(a) Integral[ t^(a-1) e^(-t), {t,x,Infinity} ]
I see a signature
Signature: (double x(); double [o]y(); double [o]e(); double a)
If I refer to the GSL documentation, I see functions
double gsl_sf_gamma_inc_Q( double a, double x );
int gsl_sf_gamma_inc_Q_e( double a, double x, gsl_sf_result* result);
I don't see how to map from the PDL documentation to the GSL function.
For the GSL, in each case we have 3 doubles: 2 input parameters and
one result. In the case of the "_e" version of the function, we also
have an integer status return value.
The PDL signature has *4* doubles. 2 would map to the function arguments,
and one would hold the result. So what is the additional
double [o]e() output argument for? This holds for all the routines
in this module - an extra output double argument, with no explanation
of what it means.
Thanks.
Terry Gaetz
_______________________________________________
Perldl mailing list
[email protected]
http://mailman.jach.hawaii.edu/mailman/listinfo/perldl
