因为涉及的维度非常小n <5,并且出于效率考虑,我想将Distributions.jl与StaticArrays.jl向量和矩阵一起使用。
但是,我有种感觉,我们不能一起使用这些软件包:
可以使用Julia“ native”向量/矩阵类型:
julia> MvNormal(zeros(3),reshape(Float64[2,0,1,0,2,0,1,0,2],3,3))
FullNormal(
dim: 3
μ: [0.0, 0.0, 0.0]
Σ: [2.0 0.0 1.0; 0.0 2.0 0.0; 1.0 0.0 2.0]
)
使用StaticArrays.jl调用时出错:
julia> MvNormal(zeros(SVector{3}),SMatrix{3,3}(Float64[2,0,1,0,2,0,1,0,2]))
ERROR: MethodError: no method matching MvNormal(::SArray{Tuple{3},Float64,1,3}, ::SArray{Tuple{3,3},Float64,2,9})
Stacktrace:
[1] top-level scope at none:0
问题:这是真的吗,还是有一种方法可以使用StaticArrays.jl定义分布?