Test block_minres on GPUs

JuliaSmoothOptimizers · Nov 4, 2024 · 1ab7105 · 1ab7105
1 parent 60ffac6
commit 1ab7105
Show file tree

Hide file tree

Showing 3 changed files with 28 additions and 1 deletion.
diff --git a/test/gpu/amd.jl b/test/gpu/amd.jl
@@ -103,6 +103,15 @@ include("gpu.jl")
       @test norm(b - A * x) ≤ atol + rtol * norm(b)
     end
 
+    @testset "block-MINRES -- $FC" begin
+      A, b = symmetric_indefinite(FC=FC)
+      B = hcat(b, -b)
+      A = M(A)
+      B = M(B)
+      X, stats = block_minres(A, B)
+      @test norm(B - A * X) ≤ atol + rtol * norm(B)
+    end
+
     @testset "block-GMRES -- $FC" begin
       A, b = nonsymmetric_indefinite(FC=FC)
       B = hcat(b, -b)

diff --git a/test/gpu/intel.jl b/test/gpu/intel.jl
@@ -86,8 +86,17 @@ include("gpu.jl")
       @test norm(b - A * x) ≤ atol + rtol * norm(b)
     end
 
+    @testset "block-MINRES -- $FC" begin
+      A, b = symmetric_indefinite(FC=FC)
+      B = hcat(b, -b)
+      A = M(A)
+      B = M(B)
+      X, stats = block_minres(A, B)
+      @test norm(B - A * X) ≤ atol + rtol * norm(B)
+    end
+
     @testset "block-GMRES -- $FC" begin
-      A, b = symmetric_definite(FC=FC)
+      A, b = nonsymmetric_definite(FC=FC)
       B = hcat(b, -b)
       A = M(A)
       B = M(B)

diff --git a/test/gpu/nvidia.jl b/test/gpu/nvidia.jl
@@ -186,6 +186,15 @@ include("gpu.jl")
       @test norm(b - A * x) ≤ atol + rtol * norm(b)
     end
 
+    @testset "block-MINRES -- $FC" begin
+      A, b = symmetric_indefinite(FC=FC)
+      B = hcat(b, -b)
+      A = M(A)
+      B = M(B)
+      X, stats = block_minres(A, B)
+      @test norm(B - A * X) ≤ atol + rtol * norm(B)
+    end
+
     @testset "block-GMRES -- $FC" begin
       A, b = nonsymmetric_indefinite(FC=FC)
       B = hcat(b, -b)