R Under development (unstable) (2024-06-26 r86840 ucrt) -- "Unsuffered Consequences"
Copyright (C) 2024 The R Foundation for Statistical Computing
Platform: x86_64-w64-mingw32/x64

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> #
> # fields  is a package for analysis of spatial data written for
> # the R software environment.
> # Copyright (C) 2022 Colorado School of Mines
> # 1500 Illinois St., Golden, CO 80401
> # Contact: Douglas Nychka,  douglasnychka@gmail.edu,
> #
> # This program is free software; you can redistribute it and/or modify
> # it under the terms of the GNU General Public License as published by
> # the Free Software Foundation; either version 2 of the License, or
> # (at your option) any later version.
> # This program is distributed in the hope that it will be useful,
> # but WITHOUT ANY WARRANTY; without even the implied warranty of
> # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
> # GNU General Public License for more details.
> #
> # You should have received a copy of the GNU General Public License
> # along with the R software environment if not, write to the Free Software
> # Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301  USA
> # or see http://www.r-project.org/Licenses/GPL-2
> ##END HEADER
> ##END HEADER
> 
> 
> suppressMessages(library(fields))
> options(echo=FALSE)
Testing:  Tps radial basis function 2d
PASSED test at tolerance  1e-08
Testing:  Tps radial basis function 2d
PASSED test at tolerance  1e-08
Testing:  Tps radial basis function 2d
PASSED test at tolerance  1e-08
Testing:  Tps radial basis function 3d
PASSED test at tolerance  1e-08
Testing:  Tps radial basis function 3d
PASSED test at tolerance  1e-08
Testing:  Tps radial basis function 3d
PASSED test at tolerance  1e-08
PASSED test at tolerance  1e-10
PASSED test at tolerance  1e-10
Testing:   der of Rad.cov
PASSED test at tolerance  1e-06
Testing:   der of Rad.cov
PASSED test at tolerance  1e-06
PASSED test at tolerance  1e-08
Testing:  Tps 2-d m=2 sanity check
PASSED test at tolerance  1e-08
Testing:  Tps 2-d m=2 sanity predict
PASSED test at tolerance  1e-08
Testing:  Tps vs.  Krig w/ GCV
PASSED test at tolerance  1e-08
Testing:  Tps vs. radial basis w Krig
PASSED test at tolerance  1e-08
Warning: 
Grid searches over lambda (nugget and sill variances) with  minima at the endpoints: 
  (GCV) Generalized Cross-Validation 
   minimum at  right endpoint  lambda  =  6.735369e-05 (eff. df= 475 )
PASSED test at tolerance  0.001
PASSED test at tolerance  0.001
Testing:  Tps.cov marginal
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess tau
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess default
PASSED test at tolerance  1e-08
Testing:  Eff.df Tps vs spatialProcess
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess SE
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess predictions on a grid
PASSED test at tolerance  1e-08
Testing:  Tps vs spatialProcess SE predictions on a grid
PASSED test at tolerance  1e-08
> cat("all done testing Tps and spatialProcess with Tps.cov", fill=TRUE)
all done testing Tps and spatialProcess with Tps.cov
> # 
> # lambda0<- out0$lambda.est["REML", "lambda"]
> # 
> # Krig.flplike(out0, lambda0)
> 
> 
> proc.time()
   user  system elapsed 
   2.07    0.34    2.42