Fix bug in QuadraticFormNetwork constructor (#169)

ITensor · May 7, 2024 · d18a61f · d18a61f
1 parent df09ba5
commit d18a61f
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 2 deletions.
diff --git a/src/formnetworks/quadraticformnetwork.jl b/src/formnetworks/quadraticformnetwork.jl
@@ -59,7 +59,7 @@ function QuadraticFormNetwork(
   kwargs...,
 )
   blf = BilinearFormNetwork(
-    bra,
+    ket,
     ket;
     dual_site_index_map=dual_index_map,
     dual_link_index_map=dual_index_map,

diff --git a/test/test_forms.jl b/test/test_forms.jl
@@ -43,8 +43,12 @@ using Random: Random
   @test underlying_graph(operator_network(blf)) == underlying_graph(A)
   @test underlying_graph(bra_network(blf)) == underlying_graph(ψbra)
 
+  qf = QuadraticFormNetwork(ψket)
+  @test nv(qf) == 3 * nv(ψket)
+  @test isempty(flatten_siteinds(qf))
+
   qf = QuadraticFormNetwork(A, ψket)
-  @test nv(qf) == 2 * nv(ψbra) + nv(A)
+  @test nv(qf) == 2 * nv(ψket) + nv(A)
   @test isempty(flatten_siteinds(qf))
 
   v = (1, 1)