Hello all!
I write pdl program to calculate hyperfine splitting of atomic levels
and I would use module PDL::GSL::COUPLING for calculations of intensity.
However, I faced with a problem. Function gsl_sf_coupling_6j returns
zero values for all my data.
___________________________________
use PDL;
use PDL::GSLSF::COUPLING;
my ($y,$e)=gsl_sf_coupling_6j(7, 6, 5, 4, 5, 2);
print $y;
___________________________________
OUTPUT: 0
When I use pure C-library:
___________________________
#include<stdio.h>
#include<gsl_sf.h>
int main(void){
double ret=gsl_sf_coupling_6j(7, 6, 5, 4, 5, 2);
printf("%f",ret);
return 1;
}
____________________________
I see value -0.123718 which is right. And similar using fortran code.
I use Slackware Linux x86_64 14.0, PDL-2.016, gsl-1.16, perl 5.16.1
Such error occurs at all computers with other versions of libraries.
How to fix this problem?
With respect, Yury
------------------------------------------------------------------------------
_______________________________________________
pdl-general mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/pdl-general
