Set concrete truncation strategy directly into decomposition alg

Author: leburgel

src/algorithms/ctmrg/c4v.jl

@@ -191,7 +191,7 @@ Compute the C₄ᵥ projector from `eigh` decomposing the Hermitian `enlarged_co
 Also return the normalized eigenvalues as the new corner tensor.
 """
 function c4v_projector!(enlarged_corner, alg::C4vEighProjector)
-    alg = @set alg.trunc = truncation_strategy(alg, enlarged_corner)
+    alg = _set_truncation(alg, truncation_strategy(alg, enlarged_corner))
     eigh_alg = decomposition_algorithm(alg)
 
     D, V, truncation_error = eigh_trunc!(enlarged_corner, eigh_alg)

src/algorithms/ctmrg/projectors.jl

@@ -55,11 +55,7 @@ end
 
 Return the tensor decomposition algorithm of the `alg` projector algorithm.
 """
-function decomposition_algorithm(alg::ProjectorAlgorithm)
-    decomposition_alg = alg.decomposition_alg
-    decomposition_alg = @set decomposition_alg.trunc = alg.trunc
-    return decomposition_alg
-end
+decomposition_algorithm(alg::ProjectorAlgorithm) = alg.decomposition_alg
 
 function truncation_strategy(alg::ProjectorAlgorithm, edge)
     if alg.trunc isa FixedSpaceTruncation
@@ -77,7 +73,7 @@ Update the truncation strategy of a given projector algorithm, keeping all other
 the same.
 """
 function _set_truncation(alg::ProjectorAlgorithm, trunc::TruncationStrategy)
-    alg´ = @set alg.trunc = trunc
+    alg´ = @set alg.decomposition_alg.trunc = trunc
     return alg´
 end
 

src/algorithms/ctmrg/sequential.jl

@@ -95,10 +95,10 @@ function sequential_projectors(
     _, r, c = coordinate
     r′ = _prev(r, size(env, 2))
     trunc = truncation_strategy(alg, env.edges[WEST, r′, c])
-    alg′ = @set alg.trunc = trunc
+    alg´ = _set_truncation(alg, trunc)
     Q1 = TensorMap(EnlargedCorner(network, env, (SOUTHWEST, r, c)))
     Q2 = TensorMap(EnlargedCorner(network, env, (NORTHWEST, r′, c)))
-    return compute_projector((Q1, Q2), alg′)
+    return compute_projector((Q1, Q2), alg´)
 end
 function sequential_projectors(
         coordinate::NTuple{3, Int}, network, env::CTMRGEnv, alg::FullInfiniteProjector
@@ -108,14 +108,14 @@ function sequential_projectors(
     coordinate_ne = _next_coordinate(coordinate_nw, rowsize, colsize)
     coordinate_se = _next_coordinate(coordinate_ne, rowsize, colsize)
     trunc = truncation_strategy(alg, env.edges[WEST, coordinate_nw[2:3]...])
-    alg′ = @set alg.trunc = trunc
+    alg´ = _set_truncation(alg, trunc)
     ec = (
         TensorMap(EnlargedCorner(network, env, coordinate_se)),
         TensorMap(EnlargedCorner(network, env, coordinate)),
         TensorMap(EnlargedCorner(network, env, coordinate_nw)),
         TensorMap(EnlargedCorner(network, env, coordinate_ne)),
     )
-    return compute_projector(ec, alg′)
+    return compute_projector(ec, alg´)
 end
 
 """

src/algorithms/ctmrg/simultaneous.jl

@@ -95,7 +95,7 @@ function simultaneous_projectors(
     ) where {E}
     coordinate′ = _next_coordinate(coordinate, size(env)[2:3]...)
     trunc = truncation_strategy(alg, env.edges[coordinate[1], coordinate′[2:3]...])
-    alg′ = @set alg.trunc = trunc
+    alg′ = _set_truncation(alg, trunc)
     ec = (enlarged_corners[coordinate...], enlarged_corners[coordinate′...])
     return compute_projector(ec, alg′)
 end
@@ -107,7 +107,7 @@ function simultaneous_projectors(
     coordinate3 = _next_coordinate(coordinate2, rowsize, colsize)
     coordinate4 = _next_coordinate(coordinate3, rowsize, colsize)
     trunc = truncation_strategy(alg, env.edges[coordinate[1], coordinate2[2:3]...])
-    alg′ = @set alg.trunc = trunc
+    alg′ = _set_truncation(alg, trunc)
     ec = (
         enlarged_corners[coordinate4...],
         enlarged_corners[coordinate...],