https://github.com/cran/fields
Tip revision: a81d415efd7dfef5f862ed54900428b54512e5e9 authored by Doug Nychka on 24 May 2001, 16:33:19 UTC
version 1.0
version 1.0
Tip revision: a81d415
rad.cov.r
"rad.cov" <-
function (x1, x2, p = 1, with.log = T, with.constant = T, C = NULL)
{
if (!is.loaded(symbol.For("radbas"))) {
temp <- dyn.load(paste(FIELDS.BIN, "fields.o", sep = ""),
2)
}
if (is.null(C)) {
if (!is.matrix(x1))
x1 <- as.matrix(x1)
if (!is.matrix(x2))
x2 <- as.matrix(x2)
d <- ncol(x1)
n1 <- nrow(x1)
n2 <- nrow(x2)
m <- (d + p)/2
par <- c(p/2, ifelse((d%%2 == 0) & (with.log), 1, 0))
temp <- .Fortran("radbas", nd = as.integer(d), x1 = as.double(x1),
n1 = as.integer(n1), x2 = as.double(x2), n2 = as.integer(n2),
par = as.double(par), k = as.double(rep(0, n1 * n2)))
if (with.constant) {
Amd <- radbas.constant(m, d)
}
else {
Amd <- 1
}
Amd * matrix(temp$k, ncol = n2, nrow = n1)
}
else {
rad.covC(x1, x2, p = p, with.log = with.log, with.constant = with.constant,
C = C)
}
}