StackOverflowError: infinite call loop between LazyArrays and ArrayLayouts
Los mantenedores suelen responder en 1 día
Nadie ha tomado este issue todavía.
Evaluación
- Dificultad
- 4/5
- Tiempo estimado
- 3-5 días
- Aptitud para principiantes
- 35/100
Línea de trabajo
Reproduce el ejemplo e inspecciona LazyArrays/src/linalg/mul.jl junto con ArrayLayouts/src/mul.jl y ArrayLayouts/src/ArrayLayouts.jl, donde el stack trace entra en la ruta de materialización recursiva. Determina por qué esta multiplicación anidada sigue haciendo dispatch entre los dos paquetes; el trabajo estará terminado cuando Matrix(A) se materialice sin un StackOverflowError para la expresión indicada.
Escrito por el modelo de indexación a partir del texto del issue.
Descripción
julia> A = zeros(1,1) * BandedMatrix(0=>Zeros(∞))[1:1,1:1] * zeros(1,1)
(1×1 Array{Float64,2}) * (1×1 BandedMatrix{Float64} with bandwidths (0, 0) with data (1×0 view(::Array{Float64,2}, 1:1, 1:0) with eltype Float64) * (0×1 view(::LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}}, 1:0, 1:1) with eltype Float64)) * (1×1 Array{Float64,2}):
0.0
julia> Matrix(A)
ERROR: StackOverflowError:
Stacktrace:
[1] simplify at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\linalg\mul.jl:290 [inlined]
[2] simplify(::ArrayLayouts.Mul{LazyArrays.ApplyLayout{typeof(*)},ArrayLayouts.DenseColumnMajor,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}}}},Array{Float64,2}}) at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\linalg\mul.jl:300
[3] copy at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\linalg\mul.jl:305 [inlined]
[4] materialize at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\mul.jl:105 [inlined]
[5] mul at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\mul.jl:106 [inlined]
[6] * at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\mul.jl:170 [inlined]
[7] * at .\operators.jl:538 [inlined]
[8] _default_materialize at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:86 [inlined]
[9] copy at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:89 [inlined]
[10] materialize at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:81 [inlined]
[11] copyto! at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:92 [inlined]
[12] _copyto! at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:334 [inlined]
[13] _copyto!(::Array{Float64,2}, ::LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}},Array{Float64,2}}}) at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\ArrayLayouts.jl:169
[14] copyto!(::Array{Float64,2}, ::LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}},Array{Float64,2}}}) at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\ArrayLayouts.jl:171
... (the last 4 lines are repeated 21758 more times)
[87047] copyto! at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:92 [inlined]
[87048] _copyto! at C:\Users\pty\.julia\packages\LazyArrays\6zNYO\src\lazyapplying.jl:334 [inlined]
[87049] _copyto!(::Array{Float64,2}, ::LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}},Array{Float64,2}}}) at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\ArrayLayouts.jl:169
[87050] copyto! at C:\Users\pty\.julia\packages\ArrayLayouts\ATRic\src\ArrayLayouts.jl:171 [inlined]
[87051] copyto_axcheck! at .\abstractarray.jl:946 [inlined]
[87052] Array{Float64,2}(::LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}},Array{Float64,2}}}) at .\array.jl:562
[87053] Array{T,2} where T(::LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{Array{Float64,2},BandedMatrix{Float64,LazyArrays.ApplyArray{Float64,2,typeof(*),Tuple{SubArray{Float64,2,Array{Float64,2},Tuple{UnitRange{Int64},UnitRange{Int64}},false},SubArray{Float64,2,LazyArrays.ApplyArray{Float64,2,typeof(vcat),Tuple{Zeros{Float64,2,Tuple{Base.OneTo{Int64},InfiniteArrays.OneToInf{Int64}}}}},Tuple{UnitRange{Int64},UnitRange{Int64}},false}}},Base.OneTo{Int64}},Array{Float64,2}}}) at .\boot.jl:428
In the above example, the error won't happen with one of the following modifications:
- Delete
zeros(1,1), the first or the last or both. - Change
BandedMatrix(0=>Zeros(∞))[1:1,1:1]to `BandedMatrix(0=>Zeros(∞)[1:1]).
- Lenguaje dominante
- Julia
- Estrellas
- 57
- Forks
- 10
- Merge medio
- 1 d 9 h
- PR fusionados (30 d)
- 11
Preparar el entorno
Este proyecto no incluye contenedor de desarrollo, Dockerfile ni guía de contribución, así que la configuración corre por tu cuenta: empieza por su README y consulta nuestra guía para la primera contribución para los pasos generales.
Primeros pasos
- Lee el issue completo y luego la guía de contribución del proyecto.
- Comenta en el issue que vas a ocuparte — evita que dos personas hagan lo mismo.
- Haz un fork del repositorio y trabaja en una rama.
- Abre un pull request que haga referencia al número del issue.
Más de JuliaLinearAlgebra/ArrayLayouts.jl
-
Dificultad 4/5 3-5 días Aptitud para principiantes 45/100
JuliaLinearAlgebra/ArrayLayouts.jl#277 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 4/5 3-5 días Aptitud para principiantes 38/100
JuliaLinearAlgebra/ArrayLayouts.jl#253 · 5 comentarios ·
Los mantenedores suelen responder en 1 día
-
Dificultad 3/5 1-2 días Aptitud para principiantes 48/100
JuliaLinearAlgebra/ArrayLayouts.jl#243 ·
Los mantenedores suelen responder en 1 día
-
Add RectDiagonalLayoutAbierto
Dificultad 5/5 Más de una semana Aptitud para principiantes 15/100
JuliaLinearAlgebra/ArrayLayouts.jl#185 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 35/100
JuliaLinearAlgebra/ArrayLayouts.jl#157 · 2 comentarios ·
Los mantenedores suelen responder en 1 día
Todos los issues de JuliaLinearAlgebra/ArrayLayouts.jl
Issues similares
-
documentation invalid
Dificultad 1/5 1-3 horas Aptitud para principiantes 88/100
JuliaPhysics/BeamletOptics.jl#190 ·
Los mantenedores suelen responder en 2 días
-
mul! with a Symmetric or Hermitian factor throws for a matrix view whose columns are not contiguousAbierto
Dificultad 2/5 1-3 horas Aptitud para principiantes 72/100
JuliaLang/LinearAlgebra.jl#1829 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
MagneticResonanceImaging/MRIReco.jl#320 ·
Los mantenedores suelen responder en 1 día
-
Dificultad 2/5 1-3 horas Aptitud para principiantes 78/100
rfordatascience/tidytuesday#1098 ·
-
documentation
Dificultad 2/5 1-3 horas Aptitud para principiantes 70/100
trixi-framework/TrixiParticles.jl#1364 ·
Los mantenedores suelen responder en 2 días