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

Support TriangularToeplitz(Float64[],:U) #22

Closed
dlfivefifty opened this issue Oct 26, 2016 · 4 comments
Closed

Support TriangularToeplitz(Float64[],:U) #22

dlfivefifty opened this issue Oct 26, 2016 · 4 comments

Comments

@dlfivefifty
Copy link
Member

TriangularToeplitz(Float64[],:U) currently gives an unreadable error. I think it should be supported (to represent a 0 x 0 matrix), or otherwise given a better error.

Supporting it will take some thought about how to construct tmp.

@andreasnoack
Copy link
Member

I think the main challenge is the dft since we could just condition https://github.com/JuliaMatrices/ToeplitzMatrices.jl/blob/master/src/ToeplitzMatrices.jl#L366-L373 on n>0. We might want to avoid the cost of creating the dft for small n>0 so maybe a way forward is to make that field Nullable{Base.DFT.FFTW.cFFTWPlan{Complex{T},-1,true,1}} or something like that.

@dlfivefifty
Copy link
Member Author

Yes I think that's a good idea. (I forgot about Nullable...)

@dlfivefifty
Copy link
Member Author

For my usage, I think I can work around this by using Nullable{TriangularToeplitz}.

@dkarrasch
Copy link
Member

The following works, so I believe the issue has been solved some time ago.

julia> U = TriangularToeplitz(Float64[],:U)
0×0 TriangularToeplitz{Float64}

julia> x = Float64[]
Float64[]

julia> U*x
Float64[]

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants