Wigner symbols for large order

[Luis.Rivera at eost.u-strasbg.fr]

Hi,

I would like to make the computation described in the attached
figure (Wigner symbols or rotaion matrix elements) for very
large values of l ( l ~ 10000, abs(m) <= l, abs(k)<=l).

I coded it using mpz_t, mpq_t and mpf_t and it works properly
until about l ~ 70. But there, all of a sudden, it starts to
produce quite bad results. For sure I am making something wrong,
but I have hard time tracking the problem.

Some idea?

Thanks,

Luis




-------------- next part --------------
A non-text attachment was scrubbed...
Name: img92.gif
Type: image/gif
Size: 2198 bytes
Desc: not available
Url : http://gmplib.org/list-archives/gmp-discuss/attachments/20081227/a806b0bf/attachment.gif