Fix some issues with op

Author: mtfishman

src/SiteTypes/sitetype.jl

@@ -337,14 +337,14 @@ function op(name::AbstractString, s::Index...; adjoint::Bool=false, kwargs...)
   # otherwise try calling a function of the form:
   #    op!(::ITensor, ::OpName, ::SiteType, ::Index...; kwargs...)
   #
-  Op = ITensor(prime.(s)..., dag.(s)...)
-  for st in common_stypes
-    op!(Op, opn, st, s...; kwargs...)
-    if !isempty(Op)
-      adjoint && return swapprime(dag(Op), 0 => 1)
-      return Op
-    end
-  end
+  ## Op = ITensor(prime.(s)..., dag.(s)...)
+  ## for st in common_stypes
+  ##   op!(Op, opn, st, s...; kwargs...)
+  ##   if !isempty(Op)
+  ##     adjoint && return swapprime(dag(Op), 0 => 1)
+  ##     return Op
+  ##   end
+  ## end
 
   if length(s) > 1
     # No overloads for common tags found. It might be a

src/SiteTypes/sitetypes/qubit.jl

@@ -90,6 +90,11 @@ state(::StateName"Tetra4", t::SiteType"Qubit") = [
   √2 / √3 * exp(im * 4π / 3)
 ]
 
+op(::OpName"Id", ::SiteType"Qubit") = [
+  1 0
+  0 1
+]
+
 #
 # 1-Qubit gates
 #

src/SiteTypes/sitetypes/spinone.jl

@@ -54,6 +54,12 @@ state(::StateName"Y+", ::SiteType"S=1") = [-1 / 2, -im / sqrt(2), 1 / 2]
 state(::StateName"Y0", ::SiteType"S=1") = [1 / sqrt(2), 0, 1 / sqrt(2)]
 state(::StateName"Y-", ::SiteType"S=1") = [-1 / 2, im / sqrt(2), 1 / 2]
 
+op(::OpName"Id", ::SiteType"S=1") = [
+  1.0 0.0 0.0
+  0.0 1.0 0.0
+  0.0 0.0 1.0
+]
+
 op(::OpName"Sz", ::SiteType"S=1") = [
   1.0 0.0 0.0
   0.0 0.0 0.0