Get QNs working for the new FormConstructors

Author: JoeyT1994

src/ITensorNetworks.jl

@@ -31,7 +31,6 @@ include("formnetworks/abstractformnetwork.jl")
 include("formnetworks/bilinearformnetwork.jl")
 include("formnetworks/linearformnetwork.jl")
 include("formnetworks/quadraticformnetwork.jl")
-include("contraction_tree_to_graph.jl")
 include("gauging.jl")
 include("utils.jl")
 include("update_observer.jl")

src/contraction_tree_to_graph.jl

@@ -65,7 +65,7 @@ function BilinearFormNetwork(
   operator_inds = union_all_inds(s, s_mapped)
   constructor_f =
     v -> inds -> if !isempty(inds)
-      reduce(*, [delta(s, sm) for (s, sm) in zip(s[v], s_mapped[v])])
+      reduce(*, [delta(dag(s), sm) for (s, sm) in zip(s[v], s_mapped[v])])
     else
       ITensor(one(Bool))
     end

src/formnetworks/bilinearformnetwork.jl

@@ -40,10 +40,11 @@ using Test: @test, @testset
   sz_exact = expect(ψ, "Sz"; alg="exact")
   @test sz_bp ≈ sz_exact
 
-  #Test with QNS, product state so should be immediately exact
+  #Test with Quantum Numbers, product state so BP should be exact
   L, χ = 2, 2
   g = named_grid((L, L))
   s = siteinds("S=1/2", g; conserve_qns=true)
+
   ψ = ITensorNetwork(v -> isodd(sum(v)) ? "↑" : "↓", s)
 
   sz_bp = expect(ψ, "Sz"; alg="bp")