Skip to content

Network solvers redesign #242

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 25 commits into
base: main
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
25 commits
Select commit Hold shift + click to select a range
eef7d0f
Delete and move some files to make space for new code
emstoudenmire Apr 1, 2025
b15d045
Add redesigned solver codes
emstoudenmire Apr 2, 2025
3edcc3d
Update solvers codes to latest versions
emstoudenmire Jul 1, 2025
8930d80
Remove previous solver includes
emstoudenmire Jul 1, 2025
9235349
Add ConstructionBase dependency
emstoudenmire Jul 1, 2025
d70bf03
Adapt new solver codes into ITensorNetworks module
emstoudenmire Jul 1, 2025
c081dfb
Continue adapting code and improve DMRG test
emstoudenmire Jul 2, 2025
e5047fc
Rename tdvp to time_evolve. Add tests.
emstoudenmire Jul 2, 2025
3f7442f
Change default outputlevel in tests
emstoudenmire Jul 2, 2025
8ef353f
Rename and reorder some arguments
emstoudenmire Jul 2, 2025
7391c0c
Rename test files
emstoudenmire Jul 2, 2025
735fa7e
Bump version number to 0.14.0
emstoudenmire Jul 3, 2025
1d1dbc0
Rename permute_indices to align_indices
emstoudenmire Jul 4, 2025
b918731
Rename file
emstoudenmire Jul 4, 2025
8928162
Implement many style changes
emstoudenmire Jul 4, 2025
0a4db77
Improve time_steps calculation code
emstoudenmire Jul 5, 2025
9abc732
Change name of callbacks and use overloading
emstoudenmire Jul 6, 2025
7719cc2
Improve import/using code style
emstoudenmire Jul 6, 2025
751902c
Rename region_iterator_action to region_step
emstoudenmire Jul 6, 2025
ddafd9f
Change callbacks and printer to take problem
emstoudenmire Jul 6, 2025
31a86d3
Remove older region_plans.jl file
emstoudenmire Jul 6, 2025
3882692
Change extracter etc to extract, update, insert. Change related kwarg…
emstoudenmire Jul 6, 2025
ad874a5
Fix line of code in time_evolve
emstoudenmire Jul 6, 2025
4a5f686
Merge branch 'network_solvers' of https://github.com/ITensor/ITensorN…
emstoudenmire Jul 6, 2025
ba9c1ef
Commit fitting code and test
emstoudenmire Jul 7, 2025
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 3 additions & 1 deletion Project.toml
Original file line number Diff line number Diff line change
@@ -1,12 +1,13 @@
name = "ITensorNetworks"
uuid = "2919e153-833c-4bdc-8836-1ea460a35fc7"
authors = ["Matthew Fishman <[email protected]>, Joseph Tindall <[email protected]> and contributors"]
version = "0.13.12"
version = "0.14.0"

[deps]
AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"
Combinatorics = "861a8166-3701-5b0c-9a16-15d98fcdc6aa"
Compat = "34da2185-b29b-5c13-b0c7-acf172513d20"
ConstructionBase = "187b0558-2788-49d3-abe0-74a17ed4e7c9"
DataGraphs = "b5a273c3-7e6c-41f6-98bd-8d7f1525a36a"
DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"
Dictionaries = "85a47980-9c8c-11e8-2b9f-f7ca1fa99fb4"
Expand Down Expand Up @@ -53,6 +54,7 @@ AbstractTrees = "0.4.4"
Adapt = "4"
Combinatorics = "1"
Compat = "3, 4"
ConstructionBase = "1.6.0"
DataGraphs = "0.2.3"
DataStructures = "0.18"
Dictionaries = "0.4"
Expand Down
40 changes: 23 additions & 17 deletions src/ITensorNetworks.jl
Original file line number Diff line number Diff line change
Expand Up @@ -34,11 +34,7 @@ include("formnetworks/quadraticformnetwork.jl")
include("gauging.jl")
include("utils.jl")
include("update_observer.jl")
include("solvers/local_solvers/eigsolve.jl")
include("solvers/local_solvers/exponentiate.jl")
include("solvers/local_solvers/dmrg_x.jl")
include("solvers/local_solvers/contract.jl")
include("solvers/local_solvers/linsolve.jl")

include("treetensornetworks/abstracttreetensornetwork.jl")
include("treetensornetworks/treetensornetwork.jl")
include("treetensornetworks/opsum_to_ttn/matelem.jl")
Expand All @@ -48,18 +44,28 @@ include("treetensornetworks/projttns/abstractprojttn.jl")
include("treetensornetworks/projttns/projttn.jl")
include("treetensornetworks/projttns/projttnsum.jl")
include("treetensornetworks/projttns/projouterprodttn.jl")
include("solvers/solver_utils.jl")
include("solvers/defaults.jl")
include("solvers/insert/insert.jl")
include("solvers/extract/extract.jl")
include("solvers/alternating_update/alternating_update.jl")
include("solvers/alternating_update/region_update.jl")
include("solvers/tdvp.jl")
include("solvers/dmrg.jl")
include("solvers/dmrg_x.jl")
include("solvers/contract.jl")
include("solvers/linsolve.jl")
include("solvers/sweep_plans/sweep_plans.jl")

include("solvers/local_solvers/eigsolve.jl")
include("solvers/local_solvers/exponentiate.jl")
include("solvers/local_solvers/runge_kutta.jl")
include("solvers/truncation_parameters.jl")
include("solvers/iterators.jl")
include("solvers/adapters.jl")
include("solvers/sweep_solve.jl")
include("solvers/region_plans/dfs_plans.jl")
include("solvers/region_plans/euler_tour.jl")
include("solvers/region_plans/euler_plans.jl")
include("solvers/region_plans/tdvp_region_plans.jl")
include("solvers/extracter.jl")
include("solvers/inserter.jl")
include("solvers/subspace/subspace.jl")
include("solvers/subspace/densitymatrix.jl")
include("solvers/align_indices.jl")
include("solvers/operator_map.jl")
include("solvers/eigsolve.jl")
include("solvers/applyexp.jl")
include("solvers/fitting.jl")

include("apply.jl")
include("inner.jl")
include("normalize.jl")
Expand Down
2 changes: 1 addition & 1 deletion src/apply.jl
Original file line number Diff line number Diff line change
Expand Up @@ -186,7 +186,7 @@ function simple_update_bp(o, ψ, v⃗; envs, callback=Returns(nothing), apply_kw
end

function ITensors.apply(
o,
o::Union{ITensor,NamedEdge},
ψ::AbstractITensorNetwork;
envs=ITensor[],
normalize=false,
Expand Down
32 changes: 32 additions & 0 deletions src/solvers/adapters.jl
Original file line number Diff line number Diff line change
@@ -0,0 +1,32 @@

#
# TupleRegionIterator
#
# Adapts outputs to be (region, region_kwargs) tuples
#
# More generic design? maybe just assuming RegionIterator
# or its outputs implement some interface function that
# generates each tuple?
#

mutable struct TupleRegionIterator{RegionIter}
Copy link
Member

@mtfishman mtfishman Jul 2, 2025

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Maybe we can come up with a name that is more descriptive, like RegionIteratorWithKwargs. Tuple is a bit vague.

region_iterator::RegionIter
end

region_iterator(T::TupleRegionIterator) = T.region_iterator

function Base.iterate(T::TupleRegionIterator, which=1)
state = iterate(region_iterator(T), which)
isnothing(state) && return nothing
(current_region, region_kwargs) = current_region_plan(region_iterator(T))
return (current_region, region_kwargs), last(state)
end

"""
region_tuples(R::RegionIterator)
The `region_tuples` adapter converts a RegionIterator into an
iterator which outputs a tuple of the form (current_region, current_region_kwargs)
at each step.
"""
region_tuples(R::RegionIterator) = TupleRegionIterator(R)
16 changes: 16 additions & 0 deletions src/solvers/align_indices.jl
Original file line number Diff line number Diff line change
@@ -0,0 +1,16 @@

function align_indices(tn)
si = siteinds(tn)
ptn = copy(tn)
for v in vertices(tn)
is = inds(tn[v])
ls = setdiff(is, si[v])
isempty(ls) && continue
new_is = [first(ls), si[v]...]
if length(ls) >= 2
new_is = vcat(new_is, ls[2:end])
end
ptn[v] = permute(tn[v], new_is)
end
return ptn
end
143 changes: 0 additions & 143 deletions src/solvers/alternating_update/alternating_update.jl

This file was deleted.

77 changes: 0 additions & 77 deletions src/solvers/alternating_update/region_update.jl

This file was deleted.

Loading
Loading