Move code to ITensors.ITensorMPS

Project.toml

@@ -6,16 +6,9 @@ version = "0.2.1"
 [deps]
 ITensors = "9136182c-28ba-11e9-034c-db9fb085ebd5"
 KrylovKit = "0b1a1467-8014-51b9-945f-bf0ae24f4b77"
-LinearAlgebra = "37e2e46d-f89d-539d-b4ee-838fcccc9c8e"
-Observers = "338f10d5-c7f1-4033-a7d1-f9dec39bcaa0"
-Printf = "de0858da-6303-5e67-8744-51eddeeeb8d7"
-TimerOutputs = "a759f4b9-e2f1-59dc-863e-4aeb61b1ea8f"
 
 [compat]
-ITensors = "0.3.58"
-KrylovKit = "0.6"
-Observers = "0.2"
-TimerOutputs = "0.5"
+ITensors = "0.3.60"
 julia = "1.6"
 
 [extras]

src/ITensorTDVP.jl

@@ -1,19 +1,25 @@
 module ITensorTDVP
-export TimeDependentSum, dmrg_x, linsolve, tdvp, to_vec
 
-include("defaults.jl")
-include("update_observer.jl")
+using ITensors: ITensors, MPO
+using ITensors.NDTensors: @Algorithm_str
+using ITensors.ITensorMPS: ITensorMPS
+using KrylovKit: KrylovKit
+
+tdvp(args...; kwargs...) = ITensorMPS.itensortdvp_tdvp(args...; kwargs...)
+
+function KrylovKit.linsolve(A::MPO, args...; kwargs...)
+  return ITensorMPS.itensortdvp_linsolve(A, args...; kwargs...)
+end
+
+dmrg(args...; kwargs...) = ITensorMPS.itensortdvp_dmrg(args...; kwargs...)
+
+dmrg_x(args...; kwargs...) = ITensorMPS.itensortdvp_dmrg_x(args...; kwargs...)
+
+function ITensors.contract(alg::Algorithm"fit", args...; kwargs...)
+  return ITensorMPS.itensortdvp_contract(alg, args...; kwargs...)
+end
+
 include("solver_utils.jl")
-include("tdvporder.jl")
-include("tdvpinfo.jl")
-include("sweep_update.jl")
-include("alternating_update.jl")
-include("tdvp.jl")
-include("dmrg.jl")
-include("dmrg_x.jl")
-include("projmpo_apply.jl")
-include("contract_mpo_mps.jl")
-include("projmps2.jl")
-include("projmpo_mps2.jl")
-include("linsolve.jl")
+export TimeDependentSum, to_vec
+
 end