examples/kronSolve.R
In jmhewitt/telefit: Estimation and Prediction for Remote Effects Spatial Process Models
set.seed(2018)
coord.s = matrix(runif(100), ncol=2)
coord.r = matrix(runif(50), ncol=2)
d.s = as.matrix(dist(coord.s))
d.r = as.matrix(dist(coord.r))
S1 = exp(-d.s)
S2 = exp(-d.r)
A = t(chol(S1))
B = t(chol(S2))
s = 15
x = matrix(runif(nrow(S1)*nrow(S2)*s), ncol=s)
y = kronecker(A,B) %*% x
x.solved = forwardsolve.kron(A, B, y)
max(abs(x - x.solved))
jmhewitt/telefit documentation built on Feb. 9, 2020, 7:15 p.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
set.seed(2018)
coord.s = matrix(runif(100), ncol=2)
coord.r = matrix(runif(50), ncol=2)
d.s = as.matrix(dist(coord.s))
d.r = as.matrix(dist(coord.r))
S1 = exp(-d.s)
S2 = exp(-d.r)
A = t(chol(S1))
B = t(chol(S2))
s = 15
x = matrix(runif(nrow(S1)*nrow(S2)*s), ncol=s)
y = kronecker(A,B) %*% x
x.solved = forwardsolve.kron(A, B, y)
max(abs(x - x.solved))
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.