It looks like there is a typo in the gsl_sf_coupling.pd file for the
gsl_sf_coupling_6j where the call to the underlying GSL routine is to
gsl_sf_coupling_3j_e which should be changed to gsl_sf_coupling_6j_e on
line 41 of the current git version:
36 pp_def('gsl_sf_coupling_6j',
37 GenericTypes => [L],
38 Pars=>'ja(); jb(); jc(); jd(); je(); jf(); double [o]y();
double [o]e()',
39 Code =>'
40 gsl_sf_result r;
41 GSLERR(gsl_sf_coupling_3j_e,($ja(),$jb(),$jc(),$jd(),$je(),$jf(),&r))
42 $y() = r.val;
43 $e() = r.err;
44 ',
45 Doc =>'6j Symbols: (ja jb jc) over (jd je jf).'
46 );
Apologies for not having time to make the fix at the moment. I should
be able to get to this by the weekend unless another developer gets to
it before then.
Pakhomov, if you are comfortable building PDL yourself, you could make
the fix yourself. Again, apologies for not being able to resolve this
at this time.
--Chris
On 9/13/2016 16:22, Pakhomov Yury wrote:
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
------------------------------------------------------------------------------
_______________________________________________
pdl-general mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/pdl-general
