Unable to produce solutions from rsolve
Closed this issue · 2 comments
vp314 commented
Minimum Working Example
Run this interactively or in a script
using LinearAlgebra, RLinearAlgebra
A = RLinearAlgebra.generate_matrix(100, 20, 20, 8.0)
x_oracle = randn(20) * 100
b = A * x_oracle
sol = LinearSolver(TypeRPM())
x = rsolve(sol, A, b)
println("Solution when calling rsolve")
println(x)
println()
z = randn(20);
rsolve!(z, sol, A, b)
println("Solution when calling rsolve!")
println(z)
Output
The print statement of x
generates all zeros for me, while z
generates a solution. This is the case whether I run it interactively or as a script.
vp314 commented
When I run examples/rpm.jl
it works fine, so I am not sure what is going on.
vp314 commented
New merge resolves the issue