-
Notifications
You must be signed in to change notification settings - Fork 21
Open
Description
ToeplitzMatrices v0.8.2 gives the following error message on tall matrices:
julia> A = Toeplitz(1:4, [1.0,0,0])
4×3 Toeplitz{Float64, Vector{Float64}, Vector{Float64}}:
1.0 0.0 0.0
2.0 1.0 0.0
3.0 2.0 1.0
4.0 3.0 2.0
julia> A\randn(4)
ERROR: DimensionMismatch: second dimension of A, 3, does not match length of x, 4
Stacktrace:
[1] mul!(y::Vector{Float64}, A::Toeplitz{Float64, Vector{Float64}, Vector{Float64}}, x::Vector{Float64}, α::Float64, β::Float64)
@ ToeplitzMatrices C:\Users\michele.zaffalon\.julia\packages\ToeplitzMatrices\uvasQ\src\linearalgebra.jl:14
[2] cgs(A::Toeplitz{Float64, Vector{Float64}, Vector{Float64}}, x::Vector{Float64}, b::Vector{Float64}, M::ToeplitzMatrices.ToeplitzFactorization{Float64, ToeplitzMatrices.Circulant{Float64, Vector{Float64}}, ComplexF64, FFTW.cFFTWPlan{ComplexF64, -1, true, 1, Tuple{Int64}}}, max_it::Int64, tol::Float64)
@ ToeplitzMatrices.IterativeLinearSolvers C:\Users\michele.zaffalon\.julia\packages\ToeplitzMatrices\uvasQ\src\iterativeLinearSolvers.jl:150
<snip>
julia> Float64.(A)\randn(4)
3-element Vector{Float64}:
0.47506079964203446
-0.6631620446434979
-0.12761857223813952
Maybe a not implemented error method or an implicit conversion to Float64 would be better.
Metadata
Metadata
Assignees
Labels
No labels