https://github.com/cran/fields
Tip revision: e41cdb90fbcc54695ce1724284e5b3f3abcf643a authored by Doug Nychka on 21 October 2008, 16:58:57 UTC
version 5.01
version 5.01
Tip revision: e41cdb9
fields.mkpoly.R
# fields, Tools for spatial data
# Copyright 2004-2007, Institute for Mathematics Applied Geosciences
# University Corporation for Atmospheric Research
# Licensed under the GPL -- www.gpl.org/licenses/gpl.html
"fields.mkpoly" <-
function (x, m = 2)
{
if (m < 1) stop("'m' has to be larger than zero.")
if (!is.matrix(x)) x <- as.matrix(x)
d <- ncol(x)
n <- nrow(x)
nterms <- .Fortran("mkpoly", as.integer(m), as.integer(d),
nterms = as.integer(0), PACKAGE = "fields")$nterms
temp <- .Fortran("dmaket", m = as.integer(m), n = as.integer(n),
dim = as.integer(d), des = as.double(x), lddes = as.integer(n),
npoly = as.integer(nterms), tmatrix = as.double(rep(0,
n * (nterms))), ldt = as.integer(n), wptr = as.integer(rep(0,
d * m)), info = as.integer(0), ptab = as.integer(rep(0,
nterms * d)), ldptab = as.integer(nterms), PACKAGE = "fields")
temp2 <- matrix(temp$tmatrix, nrow = n)
attr(temp2, "ptab") <- matrix(temp$ptab, nrow = nterms, ncol = d)
temp2
}