[Matrix-SIG] Re: cephesmodule
Les Schaffer
godzilla@netmeg.net (Les Schaffer)
Wed, 10 Mar 1999 16:01:45 -0500 (EST)
> Do you know of some code that computes them? Isn't there someway to
> write them as a hypergeometric functions?
i wouldnt be surprised, but i doubt its worth it :-)
Numerical Recipes has some suggestions for using recurrence relations,
which i think i read, somewhere, in this particular case are very
fast. i could see coding it up in such a way that it would work great
with Python. given all the normalization issues, etc.
one thing, a spherical harmonic generally takes two scaler indicies,
and two angular positions, which ideally we could use ufuncs on. would
ufuncs handle such a thing? say we have a 2D matrix each of theta and
phi values; can we call
Y( l, m, theta, phi) l,m = integers, theta,phi = 2D matrix of floats
and get back a 2D matrix of values? that would be very very cool.
les
--
____ Les Schaffer ___| --->> Engineering R&D <<---
Theoretical & Applied Mechanics | Designspring, Inc.
Center for Radiophysics & Space Research | http://www.designspring.com/
Cornell Univ. schaffer@tam.cornell.edu | les@designspring.com