NumSBT: A subroutine for calculating spherical Bessel transforms numerically

Published: 1 February 2009| Version 1 | DOI: 10.17632/y294ttxyw4.1
J.D. Talman


Abstract A previous subroutine, LSFBTR, for computing numerical spherical Bessel (Hankel) transforms is updated with several improvements and modifications. The procedure is applicable if the input radial function and the output transform are defined on logarithmic meshes and if the input function satisfies reasonable smoothness conditions. Important aspects of the procedure are that it is simply implemented with two successive applications of the fast Fourier transform, and it yields accurate results... Title of program: NumSBT Catalogue Id: AANZ_v2_0 Nature of problem This program is a subroutine which, for a function defined numerically on a logarithmic mesh in the radial coordinate, generates the spherical Bessel, or Hankel, transform on a logarithmic mesh in the transform variable. Accurate results for large values of the transform variable are obtained, that would not be otherwise obtainable. Versions of this program held in the CPC repository in Mendeley Data AANZ_v1_0; LSFBTR; 10.1016/0010-4655(83)90057-7 AANZ_v2_0; NumSBT; 10.1016/j.cpc.2008.10.003 AANZ_v3_0; SBT, version number 3; 10.1016/j.cpc.2010.08.024 This program has been imported from the CPC Program Library held at Queen's University Belfast (1969-2019)



Computational Physics, Computational Method