Skip to content

zerocorr + inflationfactor => crash / not defined for strided matrices #37

Description

@behinger

There will be more bugreports coming, but so I dont forget

zerocorr + inflationfactor => crash / not defined for strided matrices

cov(λmle):3×3 LinearAlgebra.Diagonal{Float64, Vector{Float64}}:
 1.40284    ⋅       ⋅ 
  ⋅       14.1031   ⋅ 
  ⋅         ⋅      0.00451851
cov(λemp):3×3 Matrix{Float64}:
 14.0876      1.16592      0.0133924
  1.16592     1.38775     -0.00259049
  0.0133924  -0.00259049   0.00103344

in inflationfactor we do: mle / emp => results in following error

MethodError: no method matching rdiv!(::SparseArrays.SparseMatrixCSC{Float64, Int64}, ::LinearAlgebra.LowerTriangular{Float64, Matrix{Float64}})

Closest candidates are:

rdiv!(!Matched::StridedMatrix{T}, ::LinearAlgebra.LowerTriangular{T, <:StridedMatrix{T} where T}) where T<:Union{Float32, Float64, ComplexF32, ComplexF64} at /opt/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/triangular.jl:731

rdiv!(!Matched::StridedMatrix{T} where T, ::LinearAlgebra.LowerTriangular) at /opt/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/triangular.jl:1361

rdiv!(!Matched::ChainRulesCore.AbstractThunk, ::Any) at ~/.julia/packages/ChainRulesCore/sHMAp/src/tangent_types/thunks.jl:86

...

    /(::LinearAlgebra.Diagonal{Float64, Vector{Float64}}, ::LinearAlgebra.LowerTriangular{Float64, Matrix{Float64}})@triangular.jl:1677
    (::Main.workspace#31.var"#1#2"{Float64})(::Tuple{MixedModels.ReMat{Float64, 3}, Matrix{Float64}})@[Local: 38](http://localhost:1234/edit?id=13a303e0-b69d-11ec-289a-a56f888854a6#)
    iterate@generator.jl:47[inlined]
    collect(::Base.Generator{Base.Iterators.Zip{Tuple{Vector{MixedModels.AbstractReMat{Float64}}, Vector{Matrix{Float64}}}}, Main.workspace#31.var"#1#2"{Float64}})@array.jl:724
    map(::Function, ::Base.Iterators.Zip{Tuple{Vector{MixedModels.AbstractReMat{Float64}}, Vector{Matrix{Float64}}}})@abstractarray.jl:2878
    top-level scope@[Local: 5](http://localhost:1234/edit?id=13a303e0-b69d-11ec-289a-a56f888854a6#)

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions