Erreur en lançant unions_ellipsoids_kplm.jl pour certains paramètres

[clairebrecheteau]

Mais c'est pas tout le temps....

ERROR: TaskFailedException

nested task error: ArgumentError: matrix contains Infs or NaNs
Stacktrace:
 [1] chkfinite
   @ /opt/julias/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/lapack.jl:97 [inlined]
 [2] geevx!(balanc::Char, jobvl::Char, jobvr::Char, sense::Char, A::Matrix{Float64})
   @ LinearAlgebra.LAPACK /opt/julias/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/lapack.jl:2056
 [3] eigen!(A::Matrix{Float64}; permute::Bool, scale::Bool, sortby::typeof(LinearAlgebra.eigsortby))
   @ LinearAlgebra /opt/julias/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/eigen.jl:152
 [4] eigen(A::Matrix{Float64}; permute::Bool, scale::Bool, sortby::typeof(LinearAlgebra.eigsortby))
   @ LinearAlgebra /opt/julias/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/eigen.jl:237
 [5] eigen
   @ /opt/julias/julia-1.7.2/share/julia/stdlib/v1.7/LinearAlgebra/src/eigen.jl:235 [inlined]
 [6] aux_dim_d(Σ::Matrix{Float64}, s2min::Float64, s2max::Float64, λmin::Float64, d_prim::Int64)
   @ Main ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/examples/unions_ellipsoids_kplm.jl:28
 [7] f_Σ_dim_d(Σ::Matrix{Float64})
   @ Main ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/examples/unions_ellipsoids_kplm.jl:52
 [8] macro expansion
   @ ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/src/kplm.jl:197 [inlined]
 [9] (::GeometricClusterAnalysis.var"#25#37"{Matrix{Float64}, Int64, typeof(f_Σ_dim_d), Vector{Int64}, Vector{Vector{Int64}}, Vector{Vector{Float64}}, BitVector, UnitRange{Int64}})()
   @ GeometricClusterAnalysis ./threadingconstructs.jl:178

Stacktrace: [1] sync_end(c::Channel{Any}) @ Base ./task.jl:381 [2] macro expansion @ ./task.jl:400 [inlined] [3] kplm(rng::MersenneTwister, points::Matrix{Float64}, k::Int64, n_centers::Int64, signal::Int64, iter_max::Int64, nstart::Int64, f_Σ!::typeof(f_Σ_dim_d), first_centers::Vector{Int64}) @ GeometricClusterAnalysis ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/src/kplm.jl:174 [4] kplm(rng::MersenneTwister, points::Matrix{Float64}, k::Int64, n_centers::Int64, signal::Int64, iter_max::Int64, nstart::Int64, f_Σ!::Function) @ GeometricClusterAnalysis ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/src/kplm.jl:57 [5] top-level scope @ ~/Bureau/Codes/Codes Julia/GeometricClusterAnalysis/examples/unions_ellipsoids_kplm.jl:56

[pnavaro]

Ce fichier n'existe pas sur la branche master, j'essaye de voir ce qui ne va pas dans la branche dev_Claire.

[pnavaro]

La branche dev_Claire n'est pas à jour, je n'arrive pas à reproduire le bug :-(