test up to L=10^4, normalize_eachmatrix! not necessary

stecrotti · May 14, 2024 · a0fffc4 · a0fffc4
1 parent fbd435c
commit a0fffc4
Show file tree

Hide file tree

Showing 3 changed files with 2 additions and 3 deletions.
diff --git a/src/abstract_tensor_train.jl b/src/abstract_tensor_train.jl
@@ -144,7 +144,6 @@ Compute the normalization of ``Z=\\sum_{x^1,\\ldots,x^L} A^1(x^1)\\cdots A^L(x^L
 Return the natural logarithm of the absolute normalization ``\\log|Z|``
 """
 function LinearAlgebra.normalize!(A::AbstractTensorTrain)
-    # normalize_eachmatrix!(A)
     absZ = abs(accumulate_L(A)[2])
     L = length(A)
     x = exp(1/L * log(absZ))

diff --git a/test/periodic_tensor_train.jl b/test/periodic_tensor_train.jl
@@ -34,7 +34,7 @@
     @testset "Long tensor trains" begin
         rng = MersenneTwister(0)
         qs = (2, 2)
-        L = 2000
+        L = 10000
 
         # overflow
         A = rand_periodic_tt( [1; rand(rng, 10:15, L-1); 1], qs... )

diff --git a/test/tensor_train.jl b/test/tensor_train.jl
@@ -122,7 +122,7 @@ end
     @testset "Long tensor trains" begin
         rng = MersenneTwister(0)
         qs = (2, 2)
-        L = 2000
+        L = 10000
 
         # overflow
         A = rand_tt( [1; rand(rng, 10:15, L-1); 1], qs... )