Commit 65bea38b authored by Corson N. Areshenkoff's avatar Corson N. Areshenkoff

Fix heatkern "stat"

parent a8ae5c33
......@@ -44,7 +44,7 @@ heatkern.heatcontent <- function(x, m = 1, ...) {
# Compute heat content invariance
qm <- sapply(1:m, function(i) {
sum(colSums(V^2) * (-lambda)^i / factorial(i))
sum(apply(V, 2, function(j) sum(j %o% j)) * (-lambda)^i / factorial(i))
})
return(qm)
......
......@@ -38,7 +38,7 @@ heatkern <- function(x, stat = 'kernel', ...){
}
# Wrapper
h <- switch(method,
h <- switch(stat,
kernel = heatkern.kernel(x, ...),
trace = heatkern.trace(x, ...),
heatcontent = heatkern.heatcontent(x, ...),
......
No preview for this file type
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