spheresimgrid {FieldSim}R Documentation

Random spherical field simulation by the mifieldsim method on user grid

Description

The function spheresimgrid yields discretization of sample path of a Gaussian spherical field following the procedure described in Brouste et al. (2009).

Usage

spheresimgrid(R,Ne=100,Nr=100,nbNeighbor=4,S)

Arguments

R a covariance function (defined on the sphere) of a Random spherical field to simulate.
Ne a positive integer. Ne is the number of simulation points associated with the uniform distributed discretization of the sphere for the first step of the algorithm (Accurate simulation step)
Nr a positive integer. Nr is the number of simulation points associated with the uniform distributed discretization of the sphere for the second step of the algorithm (Refined simulation step).
nbNeighbor a positive integer. nbNeighbor must be between 1 and 32. nbNeighbor is the number of neighbors to use in the second step of the algorithm.
S a matrix with 3 rows and N columns. (S[1,n],S[2,n],S[3,n]) is the point where the field must be simulated.

Details

The function spheresimgrid yields discretization of sample path of a Gaussian spherical field associated with the covariance function given by R and at each point of the grid given by S.

Value

A list with the following components:

S a matrix with 3 rows and N columns. (S[1,n],S[2,n],S[3,n]) is the point where the field has been simulated.
Wg the vector of length at more N containing the value of the simulated field at point (S[1,n],S[2,n],S[3,n])
time the CPU time

Author(s)

Alexandre Brouste (http://subaru.univ-lemans.fr/sciences/statist/pages_persos/Brouste/) and Sophie Lambert-Lacroix (http://ljk.imag.fr/membres/Sophie.Lambert).

References

A. Brouste, J. Istas and S. Lambert-Lacroix (2009). On simulation of manifold indexed fractional Gaussian fields.

See Also

spheresim

Examples

# load FieldSim library
library(FieldSim)

## Define the grid
n <- 30
phi <- pi/n*(1:(n-1))
theta <- 2*pi/n*(1:(n-1))
S <- numeric(0)
for (kk in 1:(n-1)){
S <- cbind(S,rbind(cos(theta[kk])*sin(phi),sin(theta[kk])*sin(phi),cos(phi)))
}
## Define the autocovariance function
d<-function(x){
u <- x[1]*x[4]+x[2]*x[5]+x[3]*x[6]
if (u<(-1))
u<- -1
if (u>1)
u<-1
acos(u)
}
R<-function(x){
H<-0.1 
1/2*(d(c(1,0,0,x[1:3]))^{2*H}+d(c(1,0,0,x[4:6]))^{2*H}-d(x)^{2*H})
}

##Simulate the path
Wg <- spheresimgrid(R,Ne=50,Nr=50,nbNeighbor=4,S)$Wg 


[Package FieldSim version 2.1 Index]