Remove the kaxis function

Author: amontoison

docs/src/reference.md

@@ -11,7 +11,6 @@ Krylov.FloatOrComplex
 Krylov.niterations
 Krylov.Aprod
 Krylov.Atprod
-Krylov.kaxis
 Krylov.kaxes
 Base.show
 ```

src/fgmres.jl

@@ -111,7 +111,6 @@ function fgmres!(solver :: FgmresSolver{T,FC,S}, A, b :: AbstractVector{FC};
                  callback = solver -> false, iostream :: IO=kstdout) where {T <: AbstractFloat, FC <: FloatOrComplex{T}, S <: AbstractVector{FC}}
 
   m, n = size(A)
-  ixn = kaxis(solver.x)
   m == n || error("System must be square")
   length(b) == m || error("Inconsistent problem size")
   (verbose > 0) && @printf(iostream, "FGMRES: system of size %d\n", n)
@@ -229,7 +228,7 @@ function fgmres!(solver :: FgmresSolver{T,FC,S}, A, b :: AbstractVector{FC};
         end
         push!(s, zero(FC))
         push!(c, zero(T))
-        push!(Z, similar(S, ixn))
+        push!(Z, similar(x))
       end
 
       # Continue the process.
@@ -296,7 +295,7 @@ function fgmres!(solver :: FgmresSolver{T,FC,S}, A, b :: AbstractVector{FC};
       # Compute vₖ₊₁
       if !(solved || inner_tired || breakdown)
         if !restart && (inner_iter ≥ mem)
-          push!(V, S(undef, n))
+          push!(V, similar(x))
           push!(z, zero(FC))
         end
         @. V[inner_iter+1] = q / Hbis  # hₖ₊₁.ₖvₖ₊₁ = q

src/fom.jl

@@ -104,7 +104,6 @@ function fom!(solver :: FomSolver{T,FC,S}, A, b :: AbstractVector{FC};
               callback = solver -> false, iostream :: IO=kstdout) where {T <: AbstractFloat, FC <: FloatOrComplex{T}, S <: AbstractVector{FC}}
 
   m, n = size(A)
-  ixn = kaxis(solver.x)
   m == n || error("System must be square")
   length(b) == m || error("Inconsistent problem size")
   (verbose > 0) && @printf(iostream, "FOM: system of size %d\n", n)
@@ -279,7 +278,7 @@ function fom!(solver :: FomSolver{T,FC,S}, A, b :: AbstractVector{FC};
       # Compute vₖ₊₁.
       if !(solved || inner_tired || breakdown)
         if !restart && (inner_iter ≥ mem)
-          push!(V, similar(S, ixn))
+          push!(V, similar(x))
         end
         @. V[inner_iter+1] = q / Hbis  # hₖ₊₁.ₖvₖ₊₁ = q
       end

src/gmres.jl

@@ -104,7 +104,6 @@ function gmres!(solver :: GmresSolver{T,FC,S}, A, b :: AbstractVector{FC};
                 callback = solver -> false, iostream :: IO=kstdout) where {T <: AbstractFloat, FC <: FloatOrComplex{T}, S <: AbstractVector{FC}}
 
   m, n = size(A)
-  ixn = kaxis(solver.x)
   m == n || error("System must be square")
   length(b) == m || error("Inconsistent problem size")
   (verbose > 0) && @printf(iostream, "GMRES: system of size %d\n", n)
@@ -289,7 +288,7 @@ function gmres!(solver :: GmresSolver{T,FC,S}, A, b :: AbstractVector{FC};
       # Compute vₖ₊₁
       if !(solved || inner_tired || breakdown)
         if !restart && (inner_iter ≥ mem)
-          push!(V, similar(S, ixn))
+          push!(V, similar(x))
           push!(z, zero(FC))
         end
         @. V[inner_iter+1] = q / Hbis  # hₖ₊₁.ₖvₖ₊₁ = q

src/gpmr.jl

@@ -144,8 +144,6 @@ function gpmr!(solver :: GpmrSolver{T,FC,S}, A, B, b :: AbstractVector{FC}, c ::
 
   m, n = size(A)
   s, t = size(B)
-  ixm = kaxis(solver.x)
-  ixn = kaxis(solver.y)
   m == t         || error("Inconsistent problem size")
   s == n         || error("Inconsistent problem size")
   length(b) == m || error("Inconsistent problem size")
@@ -437,8 +435,8 @@ function gpmr!(solver :: GpmrSolver{T,FC,S}, A, B, b :: AbstractVector{FC}, c ::
     # Compute vₖ₊₁ and uₖ₊₁
     if !(solved || tired || breakdown || user_requested_exit)
       if iter ≥ mem
-        push!(V, similar(S, ixm))
-        push!(U, similar(S, ixn))
+        push!(V, similar(x))
+        push!(U, similar(y))
         push!(zt, zero(FC), zero(FC))
       end
 

src/krylov_utils.jl

@@ -198,17 +198,6 @@ function vec2str(x :: AbstractVector{T}; ndisp :: Int=7) where T <: Union{Abstra
   return s
 end
 
-"""
-  ixm = kaxis(b)
-
-Return the length of `b` if `b` is a `DenseVector`.
-Otherwise, it returns the axis of `b` along the first dimension.
-"""
-function kaxis end
-
-kaxis(v::S) where S <: DenseVector = length(v)
-kaxis(v::S) where S <: AbstractVector = axes(v, 1)
-
 """
   ixm, ixn = kaxes(S, A)