From 65c394743b54405102c2da5cdb124e9d3b44bf98 Mon Sep 17 00:00:00 2001 From: Gaspard Jankowiak Date: Wed, 12 Feb 2025 11:04:31 +0100 Subject: [PATCH] clean up --- test/erf.jl | 2 -- test/gauss.jl | 17 ----------------- test/tensors.jl | 17 ++--------------- 3 files changed, 2 insertions(+), 34 deletions(-) diff --git a/test/erf.jl b/test/erf.jl index 98e50f3..e8be465 100644 --- a/test/erf.jl +++ b/test/erf.jl @@ -132,8 +132,6 @@ end x = σ * (2 * rand() - 0.5) + μ p = [L, μ, σ] - @show diff_of_erf(p; x=x) - @test TaylorTest.check(diff_of_erf, J_diff_of_erf, p; x=x) @test TaylorTest.check(J_diff_of_erf, H_diff_of_erf, p, [2, 3]; x=x) end diff --git a/test/gauss.jl b/test/gauss.jl index edef5a5..f174504 100644 --- a/test/gauss.jl +++ b/test/gauss.jl @@ -38,23 +38,6 @@ end @test TaylorTest.check(J_gauss, H_gauss, p; x=x) end -function check_symmetry(a) - R = true - if ndims(a) == 2 - R = norm(a - a') < 1e-10 - else - for i in axes(a, 1) - r = a[i, :, :]' == a[i, :, :] - R &= r - if !r - n = maximum(abs, a[i, :, :]' - a[i, :, :]) - @warn "Tensor not symmetric at i=$i: $n" - end - end - end - return R -end - @testset "Composed Gauss function" begin p = 1 .+ 5rand(6) diff --git a/test/tensors.jl b/test/tensors.jl index 47e4406..7badd00 100644 --- a/test/tensors.jl +++ b/test/tensors.jl @@ -29,19 +29,6 @@ end TP = TO.tensorproduct -function check_symmetry(a) - R = true - for i in axes(a, 1) - r = a[i,:,:]' == a[i,:,:] - R &= r - if !r - n = maximum(abs, a[i,:,:]' - a[i,:,:]) - @warn "Tensor not symmetric at i=$i: $n" - end - end - return R -end - function prod(A, B) return TP([1, 2, 3], A, [1, 2], B, [3]) end @@ -76,14 +63,14 @@ end end φ = x -> x[1] * x[2]^2 - ∇φ = x -> [x[2]^2; 2x[1]*x[2]] + ∇φ = x -> [x[2]^2; 2x[1] * x[2]] Hφ = x -> [ 0 2x[2]; 2x[2] 2*x[1] ] f = x -> v(x) * φ(x) - Jf = x -> φ(x) * Jv(x) + v(x)*∇φ(x)' + Jf = x -> φ(x) * Jv(x) + v(x) * ∇φ(x)' Hf = x -> φ(x) * Hv(x) + symprod(Jv(x), ∇φ(x)) + TP(v(x), [1], Hφ(x), [2, 3]) x = 2rand(2) .- 0.5