# define the harmonic acceleration operator used in the
# analysis of the daily temperature data
Lcoef <- c(0,(2*pi/365)^2,0)
harmaccelLfd <- vec2Lfd(Lcoef, c(0,365))
hmat <- vec2Lfd(matrix(Lcoef, 1), c(0, 365))
stopifnot(
all.equal(harmaccelLfd, hmat)
)
Run the code above in your browser using DataLab