From 07268961f44104ed797e8ebd5f1a01ee2c1c10c2 Mon Sep 17 00:00:00 2001 From: mtfishman Date: Wed, 13 Aug 2025 16:20:35 -0400 Subject: [PATCH] Loosen some truncated SVD test bounds --- Project.toml | 2 +- test/test_matrixalgebra.jl | 10 +++++----- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/Project.toml b/Project.toml index 570e05e..1c77028 100644 --- a/Project.toml +++ b/Project.toml @@ -1,7 +1,7 @@ name = "TensorAlgebra" uuid = "68bd88dc-f39d-4e12-b2ca-f046b68fcc6a" authors = ["ITensor developers and contributors"] -version = "0.3.10" +version = "0.3.11" [deps] ArrayLayouts = "4c555306-a7a7-4459-81d9-ec55ddd5c99a" diff --git a/test/test_matrixalgebra.jl b/test/test_matrixalgebra.jl index d95e51b..36ee71c 100644 --- a/test/test_matrixalgebra.jl +++ b/test/test_matrixalgebra.jl @@ -162,7 +162,7 @@ elts = (Float32, Float64, ComplexF32, ComplexF64) # p = 2, relative = true ũ, s̃, ṽ = svd_trunc( - a; trunc=truncerr(; rtol=norm([0.3, 0.2, 0.01]) / norm(diag(s)) + eps(real(elt))) + a; trunc=truncerr(; rtol=norm([0.3, 0.2, 0.01]) / norm(diag(s)) + 10eps(real(elt))) ) @test size(ũ) == (n, 2) @test size(s̃) == (2, 2) @@ -188,7 +188,7 @@ elts = (Float32, Float64, ComplexF32, ComplexF64) # p = 2, relative = false ũ, s̃, ṽ = svd_trunc( - a; trunc=truncerr(; atol=norm([0.3, 0.2, 0.01]) + eps(real(elt))) + a; trunc=truncerr(; atol=norm([0.3, 0.2, 0.01]) + 10eps(real(elt))) ) @test size(ũ) == (n, 2) @test size(s̃) == (2, 2) @@ -207,7 +207,7 @@ elts = (Float32, Float64, ComplexF32, ComplexF64) @test size(ṽ) == (n, n) @test ũ * s̃ * ṽ ≈ a ũ, s̃, ṽ = svd_trunc( - a; trunc=truncerr(; atol=(norm(diag(s)) * (one(real(elt)) + eps(real(elt))))) + a; trunc=truncerr(; atol=(norm(diag(s)) * (one(real(elt)) + 10eps(real(elt))))) ) @test size(ũ) == (n, 0) @test size(s̃) == (0, 0) @@ -218,7 +218,7 @@ elts = (Float32, Float64, ComplexF32, ComplexF64) ũ, s̃, ṽ = svd_trunc( a; trunc=truncerr(; - rtol=(norm([0.3, 0.2, 0.01], 1) / norm(diag(s), 1) + eps(real(elt))), p=1 + rtol=(norm([0.3, 0.2, 0.01], 1) / norm(diag(s), 1) + 10eps(real(elt))), p=1 ), ) @test size(ũ) == (n, 2) @@ -228,7 +228,7 @@ elts = (Float32, Float64, ComplexF32, ComplexF64) ũ, s̃, ṽ = svd_trunc( a; trunc=truncerr(; - rtol=(norm([0.3, 0.2, 0.01], 1) / norm(diag(s), 1) - eps(real(elt))), p=1 + rtol=(norm([0.3, 0.2, 0.01], 1) / norm(diag(s), 1) - 10eps(real(elt))), p=1 ), ) @test size(ũ) == (n, 3)