Sahel Iqbal

I have the same issue.

It's interesting to note that this difference does not exist if `samples` is a `Vector{Vector}` instead of a matrix. ```julia using LinearAlgebra, Random, Distributions, BenchmarkTools, Test samples = [randn(5) for...

## For `logpdf` ```julia-repl julia> Profile.clear() julia> @profile (for _ in 1:1000; logpdf(dist, samples); end) julia> Profile.print() Overhead ╎ [+additional indent] Count File:Line; Function ========================================================= ╎86 @Base/task.jl:514; (::VSCodeServer.var"#62#63")() ╎ 86...

I'm assuming the generated images are what you asked for (I'm currently working on my first project in Julia, so there's a lot to learn). ## logpdf  ## mod_logpdf...

Sorry, my bad XD. This is the profile view plot for `logpdf`:  Most of the time seems to be taken up by `sqmahal`. Is it possible that computing the...

Yes I can. But just a question, is it problematic if `logpdf` calls the mutating version underneath? I don't know your design principles behind this package, but if we want...

> No, ideally we would not mix both paths, also eg for better compatibility with static arrays. Even though probably currently many methods don't work (in an optimized way) with...

Ok, thanks for the clarification. Then I'm afraid I don't know a solution to this.