Nothing
GetGHmatrix <- function(quadpoint = quadpoint, p1a = p1a) {
gq_vals <- statmod::gauss.quad(n = quadpoint, kind = "hermite")
xs <- gq_vals$nodes
ws <- gq_vals$weights
if (p1a == 3) {
xsmatrix <- matrix(0, nrow = 4, ncol = quadpoint^4)
wsmatrix <- xsmatrix
xsmatrix[4, ] <- rep(xs, quadpoint^3)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint)
Total <- c(Total, sub)
}
xsmatrix[3, ] <- rep(Total, quadpoint^2)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint^2)
Total <- c(Total, sub)
}
xsmatrix[2, ] <- rep(Total, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint^3)
Total <- c(Total, sub)
}
xsmatrix[1, ] <- Total
xsmatrix <- t(xsmatrix)
wsmatrix[4, ] <- rep(ws, quadpoint^3)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint)
Total <- c(Total, sub)
}
wsmatrix[3, ] <- rep(Total, quadpoint^2)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint^2)
Total <- c(Total, sub)
}
wsmatrix[2, ] <- rep(Total, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint^3)
Total <- c(Total, sub)
}
wsmatrix[1, ] <- Total
wsmatrix <- t(wsmatrix)
} else if (p1a == 2) {
xsmatrix <- matrix(0, nrow = 3, ncol = quadpoint^3)
wsmatrix <- xsmatrix
xsmatrix[3, ] <- rep(xs, quadpoint^2)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint)
Total <- c(Total, sub)
}
xsmatrix[2, ] <- rep(Total, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint^2)
Total <- c(Total, sub)
}
xsmatrix[1, ] <- Total
xsmatrix <- t(xsmatrix)
wsmatrix[3, ] <- rep(ws, quadpoint^2)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint)
Total <- c(Total, sub)
}
wsmatrix[2, ] <- rep(Total, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint^2)
Total <- c(Total, sub)
}
wsmatrix[1, ] <- Total
wsmatrix <- t(wsmatrix)
} else {
xsmatrix <- matrix(0, nrow = 2, ncol = quadpoint^2)
wsmatrix <- xsmatrix
xsmatrix[2, ] <- rep(xs, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(xs[i], quadpoint)
Total <- c(Total, sub)
}
xsmatrix[1, ] <- Total
xsmatrix <- t(xsmatrix)
wsmatrix[2, ] <- rep(ws, quadpoint)
Total <- NULL
for (i in 1:quadpoint) {
sub <- rep(ws[i], quadpoint)
Total <- c(Total, sub)
}
wsmatrix[1, ] <- Total
wsmatrix <- t(wsmatrix)
}
return(list(xsmatrix = xsmatrix, wsmatrix = wsmatrix))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.