Commit 6402edd5 authored by Corson N. Areshenkoff's avatar Corson N. Areshenkoff

Remove kl distance

parent 27cd4ca6
......@@ -16,7 +16,3 @@ spd.dist.riemannian <- function(x,y,...){
z <- x.inv.sqrt %*% y %*% x.inv.sqrt
return(norm(logm2(x.inv.sqrt %*% y %*% x.inv.sqrt), type = 'F'))
}
spd.dist.kl <- function(x,y,...){
sqrt(.5 * sum(diag(x %*% solve(y) + solve(x) %*% y - 2 * diag(rep(1, dim(x)[1])))))
}
......@@ -12,7 +12,6 @@
#' \item{"cholesky": }{The Frobenius norm of the difference between the cholesky factors
#' of \code{x} and \code{y}. Not affinely invariant.}
#' \item{"riemannian": }{The Riemmanian distance proposed by Barachant, et al. (2013)}
#' \item{"kl": }{The symmetric KL divergence between \code{x} and \code{y}}
#' }
spd.dist <- function(x, y, method = 'euclidean', ...){
......@@ -31,7 +30,6 @@ spd.dist <- function(x, y, method = 'euclidean', ...){
euclidean = spd.dist.euclidean(x, y, ...),
logeuclidean = spd.dist.logeuclidean(x, y, ...),
cholesky = spd.dist.cholesky(x, y, ...),
riemannian = spd.dist.riemannian(x, y, ...),
kl = spd.dist.kl(x, y, ...))
riemannian = spd.dist.riemannian(x, y, ...))
return(d)
}
......@@ -10,7 +10,6 @@ test_that("All distance methods run without error", {
expect_error(spd.dist(x, y, method = 'logeuclidean'), NA)
expect_error(spd.dist(x, y, method = 'cholesky'), NA)
expect_error(spd.dist(x, y, method = 'riemannian'), NA)
expect_error(spd.dist(x, y, method = 'kl'), NA)
})
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment