Remove unused non-mutating iterable constructors

Author: haampie

src/bicgstabl.jl

@@ -24,8 +24,6 @@ mutable struct BiCGStabIterable{precT, matT, solT, vecT <: AbstractVector, small
     M::smallMatT
 end
 
-bicgstabl_iterator(A, b, l; kwargs...) = bicgstabl_iterator!(zerox(A, b), A, b, l; initial_zero = true, kwargs...)
-
 function bicgstabl_iterator!(x, A, b, l::Int = 2;
     Pl = Identity(),
     max_mv_products = size(A, 2),

src/cg.jl

@@ -1,6 +1,6 @@
 import Base: start, next, done
 
-export cg, cg!, CGIterable, PCGIterable, cg_iterator, cg_iterator!
+export cg, cg!, CGIterable, PCGIterable, cg_iterator!
 
 mutable struct CGIterable{matT, solT, vecT, numT <: Real}
     A::matT
@@ -92,8 +92,6 @@ end
 
 # Utility functions
 
-@inline cg_iterator(A, b, Pl = Identity(); kwargs...) = cg_iterator!(zerox(A, b), A, b, Pl; initially_zero = true, kwargs...)
-
 function cg_iterator!(x, A, b, Pl = Identity();
     tol = sqrt(eps(real(eltype(b)))),
     maxiter::Int = size(A, 2),

src/chebyshev.jl

@@ -54,9 +54,6 @@ function next(cheb::ChebyshevIterable, iteration::Int)
     cheb.resnorm, iteration + 1
 end
 
-chebyshev_iterable(A, b, λmin::Real, λmax::Real; kwargs...) =
-    chebyshev_iterable!(zerox(A, b), A, b, λmin, λmax; kwargs...)
-
 function chebyshev_iterable!(x, A, b, λmin::Real, λmax::Real;
     tol = sqrt(eps(real(eltype(b)))), maxiter = size(A, 2), Pl = Identity(), initially_zero = false)
 

src/gmres.jl

@@ -98,8 +98,6 @@ function next(g::GMRESIterable, iteration::Int)
     g.residual.current, iteration + 1
 end
 
-gmres_iterable(A, b; kwargs...) = gmres_iterable!(zerox(A, b), A, b; initially_zero = true, kwargs...)
-
 function gmres_iterable!(x, A, b;
     Pl = Identity(),
     Pr = Identity(),

src/minres.jl

@@ -36,8 +36,6 @@ mutable struct MINRESIterable{matT, solT, vecT <: DenseVector, smallVecT <: Dens
     resnorm::realT
 end
 
-minres_iterable(A, b; kwargs...) = minres_iterable!(zerox(A, b), A, b; initially_zero = true, kwargs...)
-
 function minres_iterable!(x, A, b; 
     initially_zero::Bool = false, 
     skew_hermitian::Bool = false, 

src/simple.jl

@@ -54,12 +54,6 @@ function powm_iterable!(A, x; tol = eps(real(eltype(A))) * size(A, 2) ^ 3, maxit
     PowerMethodIterable(A, x, tol, maxiter, zero(T), similar(x), similar(x), realmax(real(T)))
 end
 
-function powm_iterable(A; kwargs...)
-    x0 = rand(Complex{real(eltype(A))}, size(A, 1))
-    scale!(x0, one(eltype(A)) / norm(x0))
-    powm_iterable!(A, x0; kwargs...)
-end
-
 """
     powm(B; kwargs...) -> λ, x, [history]