## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
    collapse = TRUE,
    comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(mcgf)

h1 <- array(
    c(0, -1, 1, 0),
    dim = c(2, 2, 3)
)
h2 <- array(
    c(0, -2, 2, 0),
    dim = c(2, 2, 3)
)
h <- sqrt(h1^2 + h2^2)
u <- array(
    rep(0:2, each = 4),
    dim = c(2, 2, 3)
)

## ----exp----------------------------------------------------------------------
c_exp <- cor_exp(
    h,
    c = 0.2,
    gamma = 0.5,
    nugget = 0.05,
    is.dist = TRUE
)

## ----cauchy-------------------------------------------------------------------
c_time <- cor_cauchy(
    u,
    a = 0.4,
    alpha = 0.5
)

## ----sep----------------------------------------------------------------------
par_s <- list(
    nugget = 0.05,
    c = 0.2,
    gamma = 0.5
)

par_t <- list(
    a = 0.4,
    alpha = 0.5
)

c_sep <- cor_sep(
    spatial = "exp",
    temporal = "cauchy",
    par_s = par_s,
    par_t = par_t,
    h = h,
    u = u
)

## ----fs-----------------------------------------------------------------------
c_fs <- cor_fs(
    nugget = 0.05,
    c = 0.2,
    gamma = 0.5,
    a = 0.4,
    alpha = 0.5,
    beta = 0.5,
    h = h,
    u = u
)

## ----tri----------------------------------------------------------------------
c_tri <- cor_lagr_tri(
    v1 = 2,
    v2 = 1,
    k = 3,
    h1 = h1,
    h2 = h2,
    u = u
)

## ----askey--------------------------------------------------------------------
c_askey <- cor_lagr_askey(
    v1 = 2,
    v2 = 1,
    k = 3,
    h1 = h1,
    h2 = h2,
    u = u
)

## ----lagr-exp-----------------------------------------------------------------
c_lagr_exp <- cor_lagr_exp(
    v1 = 2,
    v2 = 1,
    k = 3,
    h1 = h1,
    h2 = h2,
    u = u
)

## ----stat---------------------------------------------------------------------
par_base_fs <- list(
    nugget = 0.05,
    c = 0.2,
    gamma = 0.5,
    a = 0.4,
    alpha = 0.5,
    beta = 0.3
)

par_lagr <- list(
    v1 = 2,
    v2 = 1,
    k = 3
)

c_stat <- cor_stat(
    base = "fs",
    lagrangian = "lagr_exp",
    par_base = par_base_fs,
    par_lagr = par_lagr,
    lambda = 0.2,
    h = h,
    h1 = h1,
    h2 = h2,
    u = u
)