Modular parameter for transfer matrices

Project.toml

@@ -11,6 +11,7 @@ KrylovKit = "0b1a1467-8014-51b9-945f-bf0ae24f4b77"
 LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
 LoggingExtras = "e6f89c97-d47a-5376-807f-9c37f3926c36"
 MatrixAlgebraKit = "6c742aac-3347-4629-af66-fc926824e5e4"
+NonlinearSolve = "8913a72c-1f9b-4ce2-8d82-65094dcecaec"
 OptimKit = "77e91f04-9b3b-57a6-a776-40b61faaebe0"
 Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
 QuadGK = "1fd47b50-473d-5c70-9696-f719f8f3bcdc"
@@ -28,6 +29,7 @@ FastGaussQuadrature = "1.1.0"
 KrylovKit = "0.10"
 LoggingExtras = "~1.0"
 MatrixAlgebraKit = "0.6.1"
+NonlinearSolve = "4.12.0"
 OptimKit = "0.4.0"
 QuadGK = "2.11.2"
 Roots = "3.0.0"

src/TNRKit.jl

@@ -11,6 +11,7 @@ using SpecialFunctions
 using FastGaussQuadrature
 using QuadGK
 using Roots
+using NonlinearSolve
 using Base.Threads
 using Combinatorics: permutations
 

src/schemes/looptnr.jl

@@ -433,7 +433,9 @@ function loop_opt(
 end
 
 """
-Coarse-grain `ΨB` to renormalized `TA`, `TB` tensors
+Coarse-grain `ΨB` to renormalized `TA`, `TB` tensors.
+The lattice is rotated by 135 degrees in counter clockwise direction.
+The elementary modular parameter `τ₀ ↦ (1 + τ₀) / (1 - τ₀)`.
 """
 function ΨB_to_TATB(psiB::Vector{T}) where {T <: AbstractTensorMap{<:Any, <:Any, 1, 2}}
     #= 

src/schemes/trg.jl

@@ -9,7 +9,9 @@ Tensor Renormalization Group
 # Running the algorithm
     run!(::TRG, trunc::TruncationStrategy, stop::Stopcrit[, finalizer=default_Finalizer, finalize_beginning=true, verbosity=1])
 
-Each step rescales the lattice by a (linear) factor of √2
+Each step rescales the lattice by a (linear) factor of √2,
+and rotate the lattice by 45 degrees in counter clockwise direction.
+The elementary modular parameter `τ₀ ↦ (τ₀ - 1) / (τ₀ + 1)`.
 
 !!! info "verbosity levels"
     - 0: No output