[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Bug-gsl] [bug #39466] errors in gsl_sf_coupling_6j_e and gsl_sf_couplin
From: |
Patrick Alken |
Subject: |
[Bug-gsl] [bug #39466] errors in gsl_sf_coupling_6j_e and gsl_sf_coupling_9j_e |
Date: |
Thu, 11 Jul 2013 17:52:27 +0000 |
User-agent: |
Mozilla/5.0 (X11; Linux i686 on x86_64; rv:22.0) Gecko/20100101 Firefox/22.0 |
Follow-up Comment #2, bug #39466 (project gsl):
(From Alexey Illarionov)
Hi, the function should be
static
int
triangle_selection_fails(int two_ja, int two_jb, int two_jc)
{
return ((two_jb < abs(two_ja - two_jc)) || (two_jb > two_ja +
two_jc)) || GSL_IS_ODD(two_jb + two_ja + two_jc);
}
This should also solve the problem of 6j and 9j.
Btw, the the algorithm of evaluation of 3j, 6j, 9j in GSL is the worst
possible. It just follows the Racah formula and does not use any
symmetry properties. See [Tuzun et all, Comp Phys Comm, 1998] for 3j,
6j algorithms.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/bugs/?39466>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/