Add an error message if trying to directsum two tensors with `EmptyNu…

…mber`
ITensor · May 15, 2024 · d58d354 · d58d354
1 parent fa43d64
commit d58d354
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/src/tensor_operations/tensor_algebra.jl b/src/tensor_operations/tensor_algebra.jl
@@ -261,6 +261,12 @@ function directsum_projectors(
   return D1, D2
 end
 
+function directsum_projectors(::Type{<:EmptyNumber}, ::Type{<:EmptyNumber}, ::Index, ::Index, ::Index)
+  error("It is not possible to call directsum on two tensors with element type EmptyNumber.
+  If you are inputting ITensors constructor like ITensor(i, j), try specifying the element type, 
+  e.g. ITensor(Float64, i, j), or filling them with zero value, e.g. ITensor(0.0, i, j).")
+end
+
 function check_directsum_inds(A::ITensor, I, B::ITensor, J)
   a = uniqueinds(A, I)
   b = uniqueinds(B, J)