# NOT RUN {
Ns = (1:100)*10
ns = 1:100
fulltime = outer(Ns,ns,FUN = function(N,n) multivariance.timing(N,n))
contour(Ns,ns,fulltime,xlab = "N",ylab = "n",
main = "computation time of multivariance in secs",
sub = "using default parameters -
use 'determine.parameters = TRUE' to compute machine specific values")
# Run to determine the parameters of your system:
# multivariance.timing(determine.parameters = TRUE)
# }
Run the code above in your browser using DataLab