Move code to ITensors.ITensorMPS

Project.toml

@@ -5,17 +5,11 @@ 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"
+Reexport = "189a3867-3050-52da-a836-e630ba90ab69"
 
 [compat]
-ITensors = "0.3.58"
-KrylovKit = "0.6"
-Observers = "0.2"
-TimerOutputs = "0.5"
+ITensors = "0.3.59"
+Reexport = "1"
 julia = "1.6"
 
 [extras]

examples/05_utils.jl

@@ -1,5 +1,4 @@
-using ITensors: MPS, maxlinkdim
-using ITensorTDVP: ITensorTDVP
+using ITensors.ITensorMPS: ITensorMPS, MPS, maxlinkdim
 using Observers: observer, update!
 using Printf: @printf
 
@@ -12,22 +11,22 @@ function tdvp_nonuniform_timesteps(
   time_start=0.0,
   order=2,
   (step_observer!)=observer(),
-  maxdim=ITensorTDVP.default_maxdim(),
-  mindim=ITensorTDVP.default_mindim(),
-  cutoff=ITensorTDVP.default_cutoff(),
-  noise=ITensorTDVP.default_noise(),
-  outputlevel=ITensorTDVP.default_outputlevel(),
+  maxdim=ITensorMPS.default_maxdim(),
+  mindim=ITensorMPS.default_mindim(),
+  cutoff=ITensorMPS.default_cutoff(),
+  noise=ITensorMPS.default_noise(),
+  outputlevel=ITensorMPS.default_outputlevel(),
   kwargs...,
 )
   nsweeps = length(time_steps)
-  maxdim, mindim, cutoff, noise = ITensorTDVP.process_sweeps(;
+  maxdim, mindim, cutoff, noise = ITensorMPS.process_sweeps(;
     nsweeps, maxdim, mindim, cutoff, noise
   )
-  tdvp_order = ITensorTDVP.TDVPOrder(order, Base.Forward)
+  tdvp_order = ITensorMPS.TDVPOrder(order, Base.Forward)
   current_time = time_start
   for sw in 1:nsweeps
     sw_time = @elapsed begin
-      psi, PH, info = ITensorTDVP.sweep_update(
+      psi, PH, info = ITensorMPS.sweep_update(
         tdvp_order,
         solver,
         PH,
@@ -61,16 +60,16 @@ end
 function tdvp_nonuniform_timesteps(
   H,
   psi::MPS;
-  ishermitian=ITensorTDVP.default_ishermitian(),
-  issymmetric=ITensorTDVP.default_issymmetric(),
-  solver_tol=ITensorTDVP.default_solver_tol(exponentiate),
-  solver_krylovdim=ITensorTDVP.default_solver_krylovdim(exponentiate),
-  solver_maxiter=ITensorTDVP.default_solver_maxiter(exponentiate),
-  solver_outputlevel=ITensorTDVP.default_solver_outputlevel(exponentiate),
+  ishermitian=ITensorMPS.default_ishermitian(),
+  issymmetric=ITensorMPS.default_issymmetric(),
+  solver_tol=ITensorMPS.default_solver_tol(exponentiate),
+  solver_krylovdim=ITensorMPS.default_solver_krylovdim(exponentiate),
+  solver_maxiter=ITensorMPS.default_solver_maxiter(exponentiate),
+  solver_outputlevel=ITensorMPS.default_solver_outputlevel(exponentiate),
   kwargs...,
 )
   return tdvp_nonuniform_timesteps(
-    ITensorTDVP.tdvp_solver(
+    ITensorMPS.tdvp_solver(
       exponentiate;
       ishermitian,
       issymmetric,

src/ITensorTDVP.jl

@@ -1,19 +1,9 @@
 module ITensorTDVP
-export TimeDependentSum, dmrg_x, linsolve, tdvp, to_vec
 
-include("defaults.jl")
-include("update_observer.jl")
-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")
+using Reexport: @reexport
+@reexport using ITensors.ITensorMPS: tdvp, dmrg_x, to_vec, TimeDependentSum, linsolve
+
+using ITensors.ITensorMPS: ITensorMPS
+dmrg(args...; kwargs...) = ITensorMPS.alternating_update_dmrg(args...; kwargs...)
+
 end