Merge pull request #78 from issp-center-dev/fix_tenes_simple_Sy

issp-center-dev · Sep 17, 2023 · deb9d83 · deb9d83
2 parents ec83a73 + bfa44f9
commit deb9d83
Show file tree

Hide file tree

Showing 4 changed files with 130 additions and 88 deletions.
diff --git a/NEWS.md b/NEWS.md
@@ -8,3 +8,7 @@
 - Multisite observables (https://github.com/issp-center-dev/TeNeS/pull/69)
 - Finite temperature calculation (https://github.com/issp-center-dev/TeNeS/pull/75)
 - The bond dimension of CTM is automatically reduced in updating CTM (https://github.com/issp-center-dev/TeNeS/pull/77)
+
+### Bug fixes
+
+- The sign of Sy in `tenes_simple` is fixed (https://github.com/issp-center-dev/TeNeS/pull/78)
diff --git a/test/data/output_simple_mode.toml b/test/data/output_simple_mode.toml
@@ -27,18 +27,18 @@ bonds = """
 3 1 0
 """
 elements = """
-0 0 0 0 -0.25 0.0
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
-0 0 1 0 -0.0125 0.00125
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.5 0.0
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.5 0.0
-1 1 0 1 -0.0125 -0.00125
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.25 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[hamiltonian]]
@@ -50,38 +50,22 @@ bonds = """
 3 0 1
 """
 elements = """
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
+0 0 0 0 0.16666666666666669 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.25 0.0
-0 0 1 0 -0.0125 0.00125
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.25 0.0
 0 1 1 0 0.25 0.0
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.25 0.0
 0 1 0 1 -0.25 0.0
-1 1 0 1 -0.0125 -0.00125
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.25 0.0
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.5 0.0
-"""
-
-[[hamiltonian]]
-dim = [2, 2]
-bonds = """
-"""
-elements = """
-0 0 0 0 -0.25 0.0
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
-0 0 1 0 -0.0125 0.00125
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
-1 1 0 1 -0.0125 -0.00125
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.25 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.3333333333333333 0.0
 """
 
 [[hamiltonian]]
@@ -93,8 +77,18 @@ bonds = """
 3 1 1
 """
 elements = """
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.25 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.25 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[hamiltonian]]
@@ -106,14 +100,22 @@ bonds = """
 3 -1 1
 """
 elements = """
-0 0 0 0 0.125 0.0
+0 0 0 0 0.04166666666666667 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.125 0.0
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.125 0.0
 0 1 1 0 0.125 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.125 0.0
 0 1 0 1 -0.125 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.125 0.0
-1 1 1 1 0.125 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.20833333333333331 0.0
 """
 
 [[hamiltonian]]
@@ -125,8 +127,18 @@ bonds = """
 3 2 0
 """
 elements = """
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.125 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.125 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[hamiltonian]]
@@ -138,14 +150,22 @@ bonds = """
 3 0 2
 """
 elements = """
-0 0 0 0 0.0625 0.0
+0 0 0 0 -0.02083333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.0625 0.0
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.0625 0.0
 0 1 1 0 0.0625 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.0625 0.0
 0 1 0 1 -0.0625 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.0625 0.0
-1 1 1 1 0.0625 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.14583333333333331 0.0
 """
 
 [observable]
@@ -175,8 +195,8 @@ group = 2
 sites = []
 dim = 2
 elements = """
-1 0 0.0 0.5
-0 1 0.0 -0.5
+1 0 0.0 -0.5
+0 1 0.0 0.5
 """
 
 [[observable.twosite]]
@@ -190,18 +210,18 @@ bonds = """
 3 1 0
 """
 elements = """
-0 0 0 0 -0.25 0.0
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
-0 0 1 0 -0.0125 0.00125
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.5 0.0
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.5 0.0
-1 1 0 1 -0.0125 -0.00125
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.25 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[observable.twosite]]
@@ -215,40 +235,22 @@ bonds = """
 3 0 1
 """
 elements = """
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
+0 0 0 0 0.16666666666666669 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.25 0.0
-0 0 1 0 -0.0125 0.00125
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.25 0.0
 0 1 1 0 0.25 0.0
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.25 0.0
 0 1 0 1 -0.25 0.0
-1 1 0 1 -0.0125 -0.00125
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.25 0.0
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.5 0.0
-"""
-
-[[observable.twosite]]
-name = "bond_hamiltonian"
-group = 0
-dim = [2, 2]
-bonds = """
-"""
-elements = """
-0 0 0 0 -0.25 0.0
-1 0 0 0 -0.0125 -0.00125
-0 1 0 0 -0.0125 -0.00125
-0 0 1 0 -0.0125 0.00125
-1 1 1 0 -0.0125 -0.00125
-0 0 0 1 -0.0125 0.00125
-1 1 0 1 -0.0125 -0.00125
-1 0 1 1 -0.0125 0.00125
-0 1 1 1 -0.0125 0.00125
-1 1 1 1 0.25 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.3333333333333333 0.0
 """
 
 [[observable.twosite]]
@@ -262,8 +264,18 @@ bonds = """
 3 1 1
 """
 elements = """
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.25 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.25 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[observable.twosite]]
@@ -277,14 +289,22 @@ bonds = """
 3 -1 1
 """
 elements = """
-0 0 0 0 0.125 0.0
+0 0 0 0 0.04166666666666667 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.125 0.0
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.125 0.0
 0 1 1 0 0.125 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.125 0.0
 0 1 0 1 -0.125 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.125 0.0
-1 1 1 1 0.125 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.20833333333333331 0.0
 """
 
 [[observable.twosite]]
@@ -298,8 +318,18 @@ bonds = """
 3 2 0
 """
 elements = """
+0 0 0 0 -0.08333333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 0 1 1 0 0.125 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.125 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.08333333333333333 0.0
 """
 
 [[observable.twosite]]
@@ -313,14 +343,22 @@ bonds = """
 3 0 2
 """
 elements = """
-0 0 0 0 0.0625 0.0
+0 0 0 0 -0.02083333333333333 0.0
+1 0 0 0 -0.004166666666666667 0.0004166666666666667
+0 1 0 0 -0.004166666666666667 0.0004166666666666667
 1 1 0 0 -0.0625 0.0
+0 0 1 0 -0.004166666666666667 -0.0004166666666666667
 1 0 1 0 -0.0625 0.0
 0 1 1 0 0.0625 0.0
+1 1 1 0 -0.004166666666666667 0.0004166666666666667
+0 0 0 1 -0.004166666666666667 -0.0004166666666666667
 1 0 0 1 0.0625 0.0
 0 1 0 1 -0.0625 0.0
+1 1 0 1 -0.004166666666666667 0.0004166666666666667
 0 0 1 1 -0.0625 0.0
-1 1 1 1 0.0625 0.0
+1 0 1 1 -0.004166666666666667 -0.0004166666666666667
+0 1 1 1 -0.004166666666666667 -0.0004166666666666667
+1 1 1 1 0.14583333333333331 0.0
 """
 
 [[observable.twosite]]

diff --git a/test/data/output_simple_mode_2.toml b/test/data/output_simple_mode_2.toml
@@ -171,8 +171,8 @@ group = 2
 sites = []
 dim = 2
 elements = """
-1 0 0.0 0.5
-0 1 0.0 -0.5
+1 0 0.0 -0.5
+0 1 0.0 0.5
 """
 
 [[observable.twosite]]

diff --git a/tool/tenes_simple.py b/tool/tenes_simple.py
@@ -635,7 +635,7 @@ def Splus(S: float) -> np.ndarray:
     ret = np.zeros((N, N))
     for i in range(1, N):
         m = S - i
-        ret[i - 1, i] = np.sqrt((S - m) * (S + m + 1.0))
+        ret[i, i - 1] = np.sqrt((S - m) * (S + m + 1.0))
     return ret
 
 
@@ -644,7 +644,7 @@ def Sminus(S: float) -> np.ndarray:
     ret = np.zeros((N, N))
     for i in range(N - 1):
         m = S - i
-        ret[i + 1, i] = np.sqrt((S + m) * (S - m + 1.0))
+        ret[i, i + 1] = np.sqrt((S + m) * (S - m + 1.0))
     return ret