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/
