add tetra

.gitignore

@@ -52,3 +52,4 @@ docs/_build/
 .vscode/
 .vscode/
 
 
 *.out
 *.out
+*.full

genax.py

@@ -98,7 +98,9 @@ def angles_alap(anglestopick, plot = False):
     tCx = tCy / tgtA
     tCx = tCy / tgtA
 
 
     Cx = cp.concatenate((hCx, tCx), axis=None)
     Cx = cp.concatenate((hCx, tCx), axis=None)
+    Cx = cp.append(Cx, [0.5], axis=False)
     Cy = cp.concatenate((hCy, tCy), axis=None)
     Cy = cp.concatenate((hCy, tCy), axis=None)
+    Cy = cp.append(Cy, [cp.sqrt(3.0)/2.0], axis=False)
 
 
     return Cx, Cy
     return Cx, Cy
 
 
@@ -140,7 +142,10 @@ def angles_ratet(anglestopick, plot = False):
     sin = cp.sin(anglestopick)
     sin = cp.sin(anglestopick)
 
 
     Dx = cp.outer(cp.full(anglestopick.size, 1.0, dtype=cp.float64), Ex).flatten()
     Dx = cp.outer(cp.full(anglestopick.size, 1.0, dtype=cp.float64), Ex).flatten()
+    Dx = cp.append(Dx, [0.5], axis=False)
     Dy = cp.outer(cos, Ey).flatten()
     Dy = cp.outer(cos, Ey).flatten()
+    Dy = cp.append(Dy, [cp.sqrt(3.0)/6.0], axis=False)
     Dz = cp.outer(sin, Ey).flatten()
     Dz = cp.outer(sin, Ey).flatten()
+    Dz = cp.append(Dz, [cp.sqrt(2.0/3.0)], axis=False)
 
 
     return Dx, Dy, Dz
     return Dx, Dy, Dz
\ No newline at end of file