Mostly a todo list

Author: kmp5VT

ext/ITensorsVectorInterfaceExt/ITensorsVectorInterfaceExt.jl

@@ -38,7 +38,12 @@ function VectorInterface.add(a::ITensor, b::ITensor, α::Number, β::Number)
   return a * β + b * α
 end
 function VectorInterface.add!(a::ITensor, b::ITensor, α::Number, β::Number)
-  a .= a .* β .+ b .* α
+  ## TODO Coding it this way skips code which fails on GPU in broadcast.jl
+  ## Not a permanant fix, just here until the broadcast function can be
+  ## updated
+  a .*= β
+  a .+= b .* α
+  # a .= a .* β .+ b .* α
   return a
 end
 function VectorInterface.add!!(a::ITensor, b::ITensor, α::Number, β::Number)

src/broadcast.jl

@@ -414,6 +414,10 @@ function Base.copyto!(
     A, C = C, A
   end
   if !isnothing(A) && !isnothing(C) && !isnothing(α) && !isnothing(β)
+    ## TODO this line fails, because the function (r, t) -> β * r + α * t 
+    ## is an anonymous ITensor function which the GPU compiler cannot process
+    ## This code is being called in ITensorsVectorInterfaceExt but checked was missed because
+    
     map!((r, t) -> β * r + α * t, T, T, A)
   else
     bc_bc_α = find_type(Broadcasted, bc_α.args)