Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression

Details Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression

2011-10-21

arxiv.org/abs/1110.4874v2

Astronomy and Astrophysics

Astrophysics

Instrumentation and Methods for Astrophysics

13 pages, 6 figures, accepted by ApJS

Scientific paper

10.1088/0067-0049/199/1/5

We present Wavemoth, an experimental open source code for computing scalar spherical harmonic transforms (SHTs). Such transforms are ubiquitous in astronomical data analysis. Our code performs substantially better than existing publicly available codes due to improvements on two fronts. First, the computational core is made more efficient by using small amounts of precomputed data, as well as paying attention to CPU instruction pipelining and cache usage. Second, Wavemoth makes use of a fast and numerically stable algorithm based on compressing a set of linear operators in a precomputation step. The resulting SHT scales as O(L^2 (log L)^2) for the resolution range of practical interest, where L denotes the spherical harmonic truncation degree. For low and medium-range resolutions, Wavemoth tends to be twice as fast as libpsht, which is the current state of the art implementation for the HEALPix grid. At the resolution of the Planck experiment, L ~ 4000, Wavemoth is between three and six times faster than libpsht, depending on the computer architecture and the required precision. Due to the experimental nature of the project, only spherical harmonic synthesis is currently supported, although adding support or spherical harmonic analysis should be trivial.

Affiliated with

Also associated with

No associations

Rating

Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.

If you have personal experience with Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Wavemoth -- Fast spherical harmonic transforms by butterfly matrix compression will most certainly appreciate the feedback.

Rate now

Comments { 0 }

Profile ID: LFWR-SCP-O-85933