I would like this to work on 2d arrays and other dimensions if possible
using Distributions, LinearAlgebra
dist = MvNormal(I(3))
u = ones(3,4)
gradlogpdf(dist, u) # intenteded behvavior is reshape(hcat([gradlogpdf(dist, @view u[:,i]) for i in axes(u, 2)]...), size(u))```
I would like this to work on 2d arrays and other dimensions if possible