tests/exact_mle_loaded_data_test2.R

In ecrc/exageostatR: R Package Demonstrates the R / C Language Interface for Exageostat

#
#
# Copyright (c) 2017-2023, King Abdullah University of Science and Technology
# All rights reserved.
#
# ExaGeoStat-R is a software package provided by KAUST
#
#
#
# @file exact_mle_loaded_data_test2.R
# ExaGeoStat R wrapper test example
#
# @version 1.2.0

# @author Kesen Wang
# @date 2021-12-08

# Load ExaGeoStat-R lib
library("exageostatr")                                           # Load ExaGeoStat-R lib.

# 0 --> Euclidean distance, 1--> great circle distance.
dmetric         = "euclidean"                                   #0 --> Euclidean distance, 1--> great circle distance.

WORKDIR <- Sys.getenv("PWD")
FILE_PATH <- paste(WORKDIR, "/../exageostatr/data/Data.rda", sep = "")

load(file=FILE_PATH)                                  #load data


# Convert data for input
data=list("x" = as.numeric(Data1Large$x), "y" = as.numeric(Data1Large$y),
          "z" = as.numeric(Data1Large$z))

# Initiate exageostat instance
exageostat_init(hardware = list (ncores = 2, ngpus = 0, ts = 320, lts = 600,  pgrid = 1, qgrid = 1))

# Estimate MLE parameters with tile low rank method
result       = exact_mle(data, "ugsm-s", dmetric, optimization = list(clb = c(0.001, 0.001, 0.001), cub = c(5, 5, 5), tol = 1e-4, max_iters = 13))
print(result)

# Finalize exageostat instance
exageostat_finalize()