Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error displaying L factor of Tridiagonal matrix #29644

Closed
ranjanan opened this issue Oct 15, 2018 · 0 comments
Closed

Error displaying L factor of Tridiagonal matrix #29644

ranjanan opened this issue Oct 15, 2018 · 0 comments
Labels
linear algebra Linear algebra

Comments

@ranjanan
Copy link
Contributor

julia> a = sprand(10,10,0.1);

julia> a = a  * a';

julia> a = a  * a' + I;

julia> lu(a)
UMFPACK LU Factorization of a (10, 10) sparse matrix
Ptr{Nothing} @0x00007ffc56963190

julia> lu(Tridiagonal(a))
LU{Float64,Tridiagonal{Float64,SparseVector{Float64,Int64}}}
L factor:
Error showing value of type LU{Float64,Tridiagonal{Float64,SparseVector{Float64,Int64}}}:
ERROR: MethodError: no method matching Bidiagonal(::Array{Float64,1}, ::SparseVector{Float64,Int64}, ::Symbol)
Closest candidates are:
  Bidiagonal(::V<:AbstractArray{T,1}, ::V<:AbstractArray{T,1}, ::Symbol) where {T, V<:AbstractArray{T,1}} at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/LinearAlgebra/src/bidiag.jl:65
  Bidiagonal(::V<:AbstractArray{T,1}, ::V<:AbstractArray{T,1}, ::AbstractChar) where {T, V<:AbstractArray{T,1}} at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/LinearAlgebra/src/bidiag.jl:68
Stacktrace:
 [1] getproperty(::LU{Float64,Tridiagonal{Float64,SparseVector{Float64,Int64}}}, ::Symbol) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/LinearAlgebra/src/lu.jl:499
 [2] show(::IOContext{REPL.Terminals.TTYTerminal}, ::MIME{Symbol("text/plain")}, ::LU{Float64,Tridiagonal{Float64,SparseVector{Float64,Int64}}}) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/LinearAlgebra/src/lu.jl:310
 [3] display(::REPL.REPLDisplay, ::MIME{Symbol("text/plain")}, ::Any) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:131
 [4] display(::REPL.REPLDisplay, ::Any) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:135
 [5] display(::Any) at ./multimedia.jl:287
 [6] #invokelatest#1 at ./essentials.jl:697 [inlined]
 [7] invokelatest at ./essentials.jl:696 [inlined]
 [8] print_response(::IO, ::Any, ::Any, ::Bool, ::Bool, ::Any) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:154
 [9] print_response(::REPL.AbstractREPL, ::Any, ::Any, ::Bool, ::Bool) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:139
 [10] (::getfield(REPL, Symbol("#do_respond#40")){Bool,getfield(REPL, Symbol("##50#59")){REPL.LineEditREPL,REPL.REPLHistoryProvider},REPL.LineEditREPL,REPL.LineEdit.Prompt})(::Any, ::Any, ::Any) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:708
 [11] #invokelatest#1 at ./essentials.jl:697 [inlined]
 [12] invokelatest at ./essentials.jl:696 [inlined]
 [13] run_interface(::REPL.Terminals.TextTerminal, ::REPL.LineEdit.ModalInterface, ::REPL.LineEdit.MIState) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/LineEdit.jl:2261
 [14] run_frontend(::REPL.LineEditREPL, ::REPL.REPLBackendRef) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:1029
 [15] run_repl(::REPL.AbstractREPL, ::Any) at /Users/osx/buildbot/slave/package_osx64/build/usr/share/julia/stdlib/v1.0/REPL/src/REPL.jl:191
 [16] (::getfield(Base, Symbol("##719#721")){Bool,Bool,Bool,Bool})(::Module) at ./logging.jl:311
 [17] #invokelatest#1 at ./essentials.jl:697 [inlined]
 [18] invokelatest at ./essentials.jl:696 [inlined]
 [19] macro expansion at ./logging.jl:308 [inlined]
 [20] run_main_repl(::Bool, ::Bool, ::Bool, ::Bool, ::Bool) at ./client.jl:330
 [21] exec_options(::Base.JLOptions) at ./client.jl:242
 [22] _start() at ./client.jl:421
julia> versioninfo()
Julia Version 1.0.1
Commit 0d713926f8 (2018-09-29 19:05 UTC)
Platform Info:
  OS: macOS (x86_64-apple-darwin14.5.0)
  CPU: Intel(R) Core(TM) i5-5257U CPU @ 2.70GHz
  WORD_SIZE: 64
  LIBM: libopenlibm
  LLVM: libLLVM-6.0.0 (ORCJIT, broadwell)

On a separate note, /Users/osx/buildbot/slave/package_osx64/ seems wrong too.

@fredrikekre fredrikekre added the linear algebra Linear algebra label Oct 16, 2018
KristofferC pushed a commit that referenced this issue Oct 19, 2018
* fix #29644: getfield of Tridiagonal{SparseVector}

* Fix tests.

(cherry picked from commit 0902d58)
KristofferC pushed a commit that referenced this issue Feb 11, 2019
* fix #29644: getfield of Tridiagonal{SparseVector}

* Fix tests.

(cherry picked from commit 0902d58)
KristofferC pushed a commit that referenced this issue Feb 20, 2020
* fix #29644: getfield of Tridiagonal{SparseVector}

* Fix tests.

(cherry picked from commit 0902d58)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
linear algebra Linear algebra
Projects
None yet
Development

No branches or pull requests

2 participants