[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-06 Thread jyr
Just to say that I have opened a ticket now with Dan's changes: http://trac.sagemath.org/sage_trac/ticket/5996 Jens --~--~-~--~~~---~--~~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-un

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-05 Thread jyr
Hi Dan, Thanks for reviewing the code. I have included your suggestions and I have done a spell check on the file now and will submit the updated version to trac once I get an account. (Michael?) As a side remark concerning the prec argument, it is passed straight into the sqrt() functions used

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-04 Thread Ondrej Certik
On Mon, May 4, 2009 at 4:52 PM, Dan Christensen wrote: > > Jens, > > I now see that you've written an article on this topic and tested > exactly the case I've also tested: floating point 6j symbols in a > compiled language.  Since your conclusions are different from mine, I'm > curious whether yo

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-04 Thread Dan Christensen
Jens, I now see that you've written an article on this topic and tested exactly the case I've also tested: floating point 6j symbols in a compiled language. Since your conclusions are different from mine, I'm curious whether your storage system is faster or your 6j routine is slower. Is your co

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-04 Thread Dan Christensen
jyr writes: > def test_calc_factlist(nn): > r""" > Function calculates a list of precomputed factorials in order to > massively accelerate consequetive calculations of the various Typo. And maybe say "future" instead of "consecutive"? >

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-04 Thread Dan Christensen
jyr writes: > The thought has occured to me that one could use the index functions > for the 3j, 6j, and Gaunt coefficients for a much simpler storage > scheme in python by using the index as a key for a dictionary of > stored symbols. I could then extend the above published routines with > an o

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-03 Thread jyr
On May 3, 11:29 am, mabshoff wrote: > On May 3, 3:07 am, jyr wrote: > > > If there is interest I could post it here or get a trac account and > > open a ticket for it for review. > > Yes, do both. > Hi Michael, Could you open a trac account for me please, so that I can add the code to a tick

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-03 Thread jyr
Hi Ondrej, No, I have not come round to publishing it yet in Comput. Phys. Commun. or some similar journal. The thought has occured to me that one could use the index functions for the 3j, 6j, and Gaunt coefficients for a much simpler storage scheme in python by using the index as a key for a d

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-03 Thread Ondrej Certik
On Sun, May 3, 2009 at 3:59 AM, jyr wrote: > > Ok, here it is: > > --- > r""" > Calculate Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt > coefficients > > Collection of functions for calculating er 3j

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-03 Thread jyr
Ok, here it is: --- r""" Calculate Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients Collection of functions for calculating er 3j, 6j, 9j, Clebsch-Gordan, Racah as well as Gaunt coefficient

[sage-devel] Re: Wigner 3j, 6j, 9j, Clebsch-Gordan, Racah and Gaunt coefficients

2009-05-03 Thread mabshoff
On May 3, 3:07 am, jyr wrote: > Hi, Hi Jens, > I have coded some routines that calculate Wigner 3j, 6j, 9j, Clebsch- > Gordan, Racah and Gaunt coefficients (integrals over 3 spherical > harmonics) exactly. It is all in a single python file with doc tests. > If I am not mistaken then Sage curr