Skip to content

Commit

Permalink
switching to IterativeSolvers 0.3
Browse files Browse the repository at this point in the history
  • Loading branch information
montyvesselinov committed May 15, 2017
1 parent 33011b2 commit ed7f787
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 3 deletions.
2 changes: 1 addition & 1 deletion REQUIRE
@@ -1,4 +1,4 @@
julia 0.5
IterativeSolvers 0.0- 0.3.0
IterativeSolvers 0.3.0
RobustPmap
Grid
3 changes: 1 addition & 2 deletions src/lsqr.jl
Expand Up @@ -51,10 +51,9 @@ function pcgalsqriteration(forwardmodel::Function, s::Vector, X::Vector, xis::Ar
Hs = (results[length(xis)+2] - hs) / delta
b = [y - hs + Hs; zeros(p)];
bigA = PCGALowRankMatrix(etas, HX, R)
x = IterativeSolvers.lsqr(bigA, b)[1]
x = IterativeSolvers.lsqr(bigA, b)
beta_bar = x[end]
xi_bar = x[1:end-1]
@show x[1:end-1]
s = X * beta_bar
for i = 1:length(xis)#add HQ' * xi_bar to s
etai = (results[i] - hs) / delta
Expand Down

0 comments on commit ed7f787

Please sign in to comment.