Merge pull request #367 from emstoudenmire/localop_bugfix

localop bugfix
ITensor · Oct 12, 2020 · 5d36dbe · 5d36dbe
2 parents 7af4030 + 51f4a5b
commit 5d36dbe
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/itensor/mps/localop.h b/itensor/mps/localop.h
@@ -231,7 +231,7 @@ LocalOp(const ITensor& Op1,
     R_(nullptr),
     size_(-1)
     {
-    nc_ = args.getInt("NumCenter",2);
+    nc_ = args.getInt("NumCenter",1);
     if(nc_ == 1)
       update(Op1,L,R);
     else

diff --git a/itensor/mps/mps.cc b/itensor/mps/mps.cc
@@ -1260,7 +1260,7 @@ template <typename MPSType>
 IndexSet
 linkInds(MPSType const& x, int i)
     {
-    if( i == 1 ) return IndexSet(rightLinkIndex(x,i));
+    if(i == 1) return IndexSet(rightLinkIndex(x,i));
     else if( i == length(x) ) return IndexSet(leftLinkIndex(x,i));
     return IndexSet(leftLinkIndex(x,i),rightLinkIndex(x,i));
     }
@@ -1276,6 +1276,7 @@ linkInds(MPSType const& x)
     for( auto n : range1(N-1) )
       {
       auto s = linkIndex(x,n);
+      if(!s) Error("MPS or MPO has missing/null link index");
       inds.nextIndex(std::move(s));
       }
     return inds.build();