diff --git a/Project.toml b/Project.toml index 8c2cee7..df4f77f 100644 --- a/Project.toml +++ b/Project.toml @@ -10,7 +10,7 @@ LinearOperators = "5c8ed15e-5a4c-59e4-a42b-c7e8811fb125" PartitionedStructures = "a6683cb1-4e99-4d03-b5b5-b79ed0acabcf" [compat] -Krylov = "=0.9.0" +Krylov = "0.9.8" LinearOperators = "2.4.0" PartitionedStructures = "0.1.4, 0.1.5" julia = "^1.6.0" diff --git a/docs/Project.toml b/docs/Project.toml index 089ea9e..49f0874 100644 --- a/docs/Project.toml +++ b/docs/Project.toml @@ -7,7 +7,7 @@ PartitionedVectors = "9d5e22db-ddc8-46c9-8309-91b035e7f733" [compat] Documenter = "~0.27" -Krylov = "=0.9.0" +Krylov = "0.9.8" LinearOperators = "2.4.0" PartitionedStructures = "0.1.4, 0.1.5" -julia = "^1.6.0" \ No newline at end of file +julia = "^1.6.0" diff --git a/src/krylov.jl b/src/krylov.jl index bbd2f3c..c07c7ca 100644 --- a/src/krylov.jl +++ b/src/krylov.jl @@ -138,7 +138,7 @@ function CgSolver(pv::PartitionedVector{T}) where {T} Ap .= (T)(0) z = similar(pv; simulate_vector = true) z .= (T)(0) - stats = Krylov.SimpleStats(0, false, false, T[], T[], T[], "unknown") + stats = Krylov.SimpleStats(0, false, false, T[], T[], T[], 0.0, "unknown") solver = Krylov.CgSolver{T, T, PartitionedVector{T}}(n, n, Δx, x, r, p, Ap, z, true, stats) return solver end diff --git a/test/krylov.jl b/test/krylov.jl index ce47bed..ca15904 100644 --- a/test/krylov.jl +++ b/test/krylov.jl @@ -99,18 +99,18 @@ end b2 = rand(T) c = rand(T) - Krylov.@kaxpy!(n, a, x, y) - Krylov.@kaxpy!(n, a2, x, y) - Krylov.@kaxpby!(n, a, x, b, y) - Krylov.@kaxpby!(n, a2, x, b, y) - Krylov.@kaxpby!(n, a, x, b2, y) - Krylov.@kaxpby!(n, a2, x, b2, y) - - @test (@allocated Krylov.@kaxpy!(n, a, x, y)) == 0 - @test (@allocated Krylov.@kaxpy!(n, a2, x, y)) == 0 - @test (@allocated Krylov.@kaxpby!(n, a, x, b, y)) == 0 - @test (@allocated Krylov.@kaxpby!(n, a2, x, b, y)) == 0 - @test (@allocated Krylov.@kaxpby!(n, a, x, b2, y)) == 0 - @test (@allocated Krylov.@kaxpby!(n, a2, x, b2, y)) == 0 + Krylov.kaxpy!(n, a, x, y) + Krylov.kaxpy!(n, a2, x, y) + Krylov.kaxpby!(n, a, x, b, y) + Krylov.kaxpby!(n, a2, x, b, y) + Krylov.kaxpby!(n, a, x, b2, y) + Krylov.kaxpby!(n, a2, x, b2, y) + + @test (@allocated Krylov.kaxpy!(n, a, x, y)) == 0 + @test (@allocated Krylov.kaxpy!(n, a2, x, y)) == 0 + @test (@allocated Krylov.kaxpby!(n, a, x, b, y)) == 0 + @test (@allocated Krylov.kaxpby!(n, a2, x, b, y)) == 0 + @test (@allocated Krylov.kaxpby!(n, a, x, b2, y)) == 0 + @test (@allocated Krylov.kaxpby!(n, a2, x, b2, y)) == 0 end end