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

Weird behavior of transposes and TriangularToeplitz #34

Closed
riccardosven opened this issue Jan 11, 2018 · 3 comments
Closed

Weird behavior of transposes and TriangularToeplitz #34

riccardosven opened this issue Jan 11, 2018 · 3 comments

Comments

@riccardosven
Copy link

Multiplying a vector by a lower triangular Toeplitz matrix gives weird results when associated with transposition:

g = [1.; 2.; 3.]
u = TriangularToeplitz([1.; 1.; 1.], :L)

then

g'*u
1×3 RowVector{Float64,Array{Float64,1}}:
 6.0  5.0  3.0

but

g'*u
1×3 RowVector{Float64,Array{Float64,1}}:
 3.0  2.0  0.0

Similarly

u.'*g
3-element Array{Float64,1}:
 3.0
 2.0
 0.0

u*g
3-element Array{Float64,1}:
 1.0
 3.0
 6.0
@andreasnoack
Copy link
Member

Thanks for reporting. The methods here predates the RowVector so I'm not sure what this dispatches to. I'll take a look.

@riccardosven
Copy link
Author

I've tried following the chain of calls, but I'm just a Julia beginner :)
I believe it is something in A_mul_B! as it also happens in conjunction to normal vector-matrix multiplication (not only with RowVector.

@dkarrasch
Copy link
Member

I couldn't find an issue even with the currently released version, even more so since I don't understand where the issue should come from: two different results are reported for the exact same call g'*u. So I think this can be closed.

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