Nothing
hat.transf <-
function(C22, transf, vc, w, k, N, tol.err = 1e-6, GPU = FALSE) {
qu <- numeric(k)
middle <- diag(N) - C22/vc
#eigen.test <- eigen(middle, only.values = TRUE)
#if (min(eigen.test$values) < tol.err) {
## cat("Matrix middle bended\n")
## cat("Minimum eigenvalue:\n", min(eigen.test$values))
# if (min(eigen.test$values) > (-tol.err)) middle <- middle + diag(N)*2*tol.err
# if (min(eigen.test$values) < (-tol.err)) middle <- middle - (min(eigen.test$values) - tol.err)*diag(N)
#}
#M <- chol(middle)
svdmid <- svd(middle)
M <- t(svdmid$u %*% diag(sqrt(svdmid$d)))
if (!GPU) {
qu <- 1 - colSums(tcrossprod(M, transf)**2)/w
}
#else {
# qu <- 1 - colSums(gpuMatMult(M, t(transf))**2)/w
#}
return(qu)
}
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.