The Beylkin-Cramer Summation Rule and A New Fast Algorithm of Cosmic Statistics for Large Data Sets

Abstract

Based on the Beylkin-Cramer summation rule, we introduce a new fast algorithm that enable us to explore the high order statistics efficiently in large data sets. Central to this technique is to make decomposition both of fields and operators within the framework of multi-resolution analysis (MRA), and realize theirs discrete representations. Accordingly, a homogenous point process could be equivalently described by a operation of a Toeplitz matrix on a vector, which is accomplished by making use of fast Fourier transformation. The algorithm could be applied widely in the cosmic statistics to tackle large data sets. Especially, we demonstrate this novel technique using the spherical, cubic and cylinder counts in cells respectively. The numerical test shows that the algorithm produces an excellent agreement with the expected results. Moreover, the algorithm introduces naturally a sharp-filter, which is capable of suppressing shot noise in weak signals. In the numerical procedures, the algorithm is somewhat similar to particle-mesh (PM) methods in N-body simulations. As scaled with $O(N\log N)$, it is significantly faster than the current particle-based methods, and its computational cost does not relies on shape or size of sampling cells. In addition, based on this technique, we propose further a simple fast scheme to compute the second statistics for cosmic density fields and justify it using simulation samples. Hopefully, the technique developed here allows us to make a comprehensive study of non-Guassianity of the cosmic fields in high precision cosmology. A specific implementation of the algorithm is publicly available upon request to the author.

arXiv comment27 pages, 9 figures included. revised version, changes include (a) adding a new fast algorithm for 2nd statistics (b) more numerical tests including counts in asymmetric cells, the two-point correlation functions and 2nd variances (c) more discussions on technica

References in corpus (17)

• Simulating the joint evolution of quasars, galaxies and their large-scale distribution
• Stable clustering, the halo model and nonlinear cosmological power spectra
• Large-Scale Structure of the Universe and Cosmological Perturbation Theory
• Clustering of Galaxies in a Hierarchical Universe: I. Methods and Results at z=0
• The 2dF Galaxy Redshift Survey: The bias of galaxies and the density of the Universe
• The IRAS 1.2 Jy Survey: Redshift Data
• The 2dF Galaxy Redshift Survey: higher order galaxy correlation functions
• The 2dF Galaxy Redshift Survey: voids and hierarchical scaling models
• The Three-point Correlation Function of Galaxies Determined from the 2dF Galaxy Redshift Survey
• The Three-point Correlation Function of Galaxies: Comparing Halo Occupation Models with Observations
• A New Method for Calculating Counts in Cells
• Non-Gaussianity and the recovery of the mass power spectrum from the Ly$\alpha$ forest
• Fast Edge Corrected Measurement of the Two-Point Correlation Function and the Power Spectrum
• Breaking Degeneracy of Dark Matter Models by the Scale-Scale Correlations of Galaxies
• The intermittent behavior and hierarchical clustering of the cosmic mass field
• Conditional Cumulants in Weakly Non-linear Regime
• Quasi-Local Evolution of the Cosmic Gravitational Clustering in Halo Model

Cited by in corpus (4)

• Discreteness Effects in Lambda Cold Dark Matter Simulations: A Wavelet-Statistical View
• X-ray emission characteristics of two WR binaries : V444 Cyg and CD Cru
• Pair Counting without Binning -- A New Approach to Correlation Functions in Clustering Statistics
• An Optimal In-Situ Multipole Algorithm for the Isotropic Three-Point Correlation Functions