replace allocator with TensorOperations implementation

Author: lkdvos

Project.toml

@@ -53,7 +53,7 @@ RecipesBase = "1.1"
 TensorKit = "0.16.3"
 TensorKitManifolds = "0.7"
 TensorKitTensors = "0.2"
-TensorOperations = "5"
+TensorOperations = "5.5"
 Test = "1"
 TestExtras = "0.3"
 VectorInterface = "0.2, 0.3, 0.4, 0.5"

src/MPSKit.jl

@@ -66,7 +66,7 @@ using MatrixAlgebraKit: TruncationStrategy, PolarViaSVD, LAPACK_SVDAlgorithm
 using BlockTensorKit
 using BlockTensorKit: TensorMapSumSpace
 using TensorOperations
-using TensorOperations: AbstractBackend, DefaultBackend, DefaultAllocator
+using TensorOperations: AbstractBackend, DefaultBackend, DefaultAllocator, BufferAllocator
 using KrylovKit
 using KrylovKit: KrylovAlgorithm
 using OptimKit
@@ -99,7 +99,6 @@ include("utility/logging.jl")
 using .IterativeLoggers
 include("utility/iterativesolvers.jl")
 
-include("utility/allocator.jl")
 include("utility/styles.jl")
 include("utility/periodicarray.jl")
 include("utility/windowarray.jl")

src/algorithms/derivatives/derivatives.jl

@@ -222,6 +222,8 @@ Base.:*(h::LazySum{<:Union{DerivativeOrMultiplied}}, v) = h(v)
 """
     prepare_operator!!(O, [backend], [allocator]) -> O′
 
-Given an operator, try to construct a more efficient representation of that operator for repeated application.
+Given an operator, try to construct a more efficient representation of that operator.
+This typically consists of precomputing some parts of the application,
+and is expected to only pay off for repeated applications.
 """
-prepare_operator!!(O, backend::AbstractBackend = DefaultBackend(), allocator = GrowingBuffer()) = O
+prepare_operator!!(O, backend::AbstractBackend = DefaultBackend(), allocator = BufferAllocator()) = O