dmaldona/RLinearAlgebra.jl

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