fix completely wrong umfpack test

nep-pack · Jun 28, 2023 · 333ab38 · 333ab38
1 parent 2ed9f55
commit 333ab38
Show file tree

Hide file tree

Showing 2 changed files with 7 additions and 6 deletions.
diff --git a/test/cd_player_native.jl b/test/cd_player_native.jl
@@ -9,11 +9,12 @@ using LinearAlgebra
 
     (Λ,V)=polyeig(nep)
 
-    (λ,v)=quasinewton(nep,λ=Λ[1],v=V[:,1],logger=displaylevel, armijo_factor=0.5,armijo_max=10, tol=1e-10)
-    verify_lambdas(1, nep, λ, v, 1e-10)
+    (λ,v)=quasinewton(nep,λ=Λ[1],v=V[:,1],logger=displaylevel, armijo_factor=0.2,armijo_max=10, tol=1e-17)
 
-    (λ,v)=resinv(nep,λ=Λ[2],v=V[:,2],logger=displaylevel, armijo_factor=0.5,armijo_max=10, tol=1e-10)
-    verify_lambdas(1, nep, λ, v, 1e-10)
+    verify_lambdas(1, nep, λ, v, 1e-9)
+
+    (λ,v)=resinv(nep,λ=Λ[1],v=V[:,1],logger=displaylevel, armijo_factor=0.2,armijo_max=10, tol=1e-17)
+    verify_lambdas(1, nep, λ, v, 1e-9)
 
     @testset "Errors thrown" begin
         @test_throws MethodError nep_gallery("nlevp_native_cd_player", 15)

diff --git a/test/linsolver.jl b/test/linsolver.jl
@@ -87,8 +87,8 @@ using LinearAlgebra
         r2 = @timed lin_solve(c2, x)
 
         # test that the norm is at least 10 times bigger with no iterative refinements
-        norm0 = norm(M*r0[1] - x)
-        norm2 = norm(M*r2[1] - x)
+        norm0 = norm(M*r0 - x)
+        norm2 = norm(M*r2 - x)
         @test norm0 > 10 * norm2
 
     end