[SciPy-dev] Hankel transform

Tom Grydeland tom.grydeland@gmail....
Wed Dec 3 01:57:59 CST 2008


Hi developers,

When I recently needed a Hankel transform I was unable to find one in
Scipy.  I found one in MATLAB though[1], written by Prof. Brian
Borchers at New Mexico Tech. The code is based on a previous FORTRAN
implementation by W. L. Anderson [2], and the MATLAB code is not
marked with any copyright statements.  Hankel transforms of the first
and second order can be computed through digital filtering.

I have rewritten the code in Python/numpy, complete with tests and
acknowledgements of the origins, and my employer has agreed that I can
donate the code to Scipy.  I believe this should be of interest.
Hankel transforms arise often in acoustics and other systems with
cylinder symmetry.  If you want it, I would like a suggestion where it
belongs (with other integral transforms, probably) and how I should
shape the tests to make them suitable for inclusion.

The tests I currently have use Hankel transforms of five different
functions with known analytical transforms and compares the
transformed values to the numerically evaluated analytical
expressions.  Currently plots are generated, but for automated testing
I suppose something else would be better.  Pointing me at an example
is sufficient.

[1] http://infohost.nmt.edu/~borchers/hankel.html
[2] Anderson, W. L., 1979, Computer Program Numerical Integration of
Related Hankel Transforms of Orders 0 and 1 by Adaptive Digital
Filtering. Geophysics, 44(7):1287-1305.

Best regards,

-- 
Tom Grydeland
  <Tom.Grydeland@(discoverpetroleum.com|gmail.com)>


More information about the Scipy-dev mailing list