/usr/share/doc/r-cran-fields/tests/evlpoly.test.R is in r-cran-fields 9.0-3.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 | # fields is a package for analysis of spatial data written for
# the R software environment .
# Copyright (C) 2017
# University Corporation for Atmospheric Research (UCAR)
# Contact: Douglas Nychka, nychka@ucar.edu,
# National Center for Atmospheric Research,
# PO Box 3000, Boulder, CO 80307-3000
#
# 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.
library( fields)
options( echo=FALSE)
test.for.zero.flag<-1
set.seed( 245)
x<- runif(3)
coef<- runif( 5)
temp<- fields.evlpoly( x, coef)
temp2<- coef[1]
for( k in (2:5) ){
temp2<- temp2 + coef[k]*x**(k-1)
}
test.for.zero( temp, temp2)
set.seed( 124)
x<- matrix( runif(12), ncol=3)
fields.mkpoly(x, m=3)-> out
attr( out, "ptab")-> ptab
J<- nrow( ptab)
coef<- runif( J)
temp<- fields.evlpoly2( x, coef, ptab)
temp2<-out%*% coef
test.for.zero( temp,temp2)
fields.derivative.poly( x, m=3, coef)-> temp
fields.mkpoly( cbind( x[,1:2], x[,3]+1e-6), m=3)%*% coef-> temp2
fields.mkpoly( cbind( x[,1:2], x[,3]-1e-6), m=3)%*% coef-> temp3
temp2<- (temp2- temp3)/ 2e-6
test.for.zero( temp[,3], temp2)
cat("Done testing polynomial evaluation",fill=TRUE)
options( echo=FALSE)
|