Summary
This paper proposes a new scheme for the generation of Gaussian random fields over large domains (domain size much larger than the correlation length). The scheme decomposes the simulation domain into overlapping subdomains and essentially generates independent random fields over each of them before merging them on the overlaps. It is naturally suited for simulation over clusters of computers. With this approach, the number of operations for each processor depends only on the number of local degrees of freedom and not on the total number over all processors. Hence, weak scalability is perfectly met. This paper describes the general scheme and introduces two error estimates for comparison with classical sampling schemes. Improvements in terms of scalability are demonstrated both theoretically and through numerical examples. The behavior in the overlap is studied in detail. Simulations using the localized approach were performed in up to 512 processors and allowed to generate in 41 seconds a realization of a random field over a cube of side 300 correlation lengths (close to 2‐billion sampling points), much more efficiently and rapidly than with classical methods.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.