Populate ITensor with zeros if empty after contracting with ProjTTN.

ITensor · Jan 25, 2024 · 1ef03c5 · 1ef03c5
1 parent 56fcb93
commit 1ef03c5
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/src/treetensornetworks/projttns/abstractprojttn.jl b/src/treetensornetworks/projttns/abstractprojttn.jl
@@ -80,6 +80,10 @@ function contract(P::AbstractProjTTN, v::ITensor)::ITensor
   for it in itensor_map
     Hv *= it
   end
+  if isempty(Hv)
+    Hv = similar(v)
+    Hv .= 0
+  end
   return Hv
 end