Constrained Realizations of Gaussian Fields: A Simple Algorithm
Abstract
A straightforward method is presented for the construction of constrained realizations of Gaussian fields, considering a Gaussian random field and its ensemble mean field, given a set of constraints. The residual of the field from its mean is statistically independent of the actual numerical value of the constraints. The algorithm is based on a simple construction of this residual field, which is then added to the analytically calculated mean field. This algorithm is exact and involves no iterations. The algorithm facilitates the generation of the initial conditions for N-body simulations which obey a few hundreds of constraints imposed by the observable universe.
- Publication:
-
The Astrophysical Journal
- Pub Date:
- October 1991
- DOI:
- 10.1086/186160
- Bibcode:
- 1991ApJ...380L...5H
- Keywords:
-
- Astronomical Models;
- Computational Astrophysics;
- Computerized Simulation;
- Cosmology;
- Algorithms;
- Dark Matter;
- Fast Fourier Transformations;
- Galactic Clusters;
- Many Body Problem;
- Power Spectra;
- Astrophysics;
- COSMOLOGY;
- GALAXIES: CLUSTERING