-
Notifications
You must be signed in to change notification settings - Fork 1
/
main.R
37 lines (30 loc) · 930 Bytes
/
main.R
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
rm(list = ls())
data(Tr1)
data(V1)
ngrid.x=40
ngrid.y=20
ngrid.t=10
xx=seq(-0.89,3.39,length.out=ngrid.x)
yy=seq(-0.89,0.89,length.out=ngrid.y)
ss=expand.grid(xx,yy)
tt=(0:(ngrid.t-1))/(ngrid.t-1)
Data=data.frame(x=rep(ss[,1],ngrid.t),y=rep(ss[,2],ngrid.t),
t=rep(tt,each=dim(ss)[1]))
knots=c(0.2,0.4,0.6,0.8)
Boundary.knots=c(0,1)
# library(BPST)
# library(splines2)
# library(MGLM)
# library(Matrix)
library(TPST)
d <- 2
r <- 1
rho <- 3
Basis1 <- basis.tensor(ss = ss, tt = tt, V = V1, Tri = Tr1,
d = d, r = r, time.knots = knots, rho = rho,
time.bound = Boundary.knots, line.up = TRUE)
Basis2 <- basis.tensor(ss = Data[,1:2], tt = Data[,3],
V = V1, Tri = Tr1, d = d, r = r,
time.knots = knots, rho = rho,
time.bound = Boundary.knots, line.up = FALSE)
which(Basis1$Psi.Q2 != Basis2$Psi.Q2)