From 2954c0ee3a41d3fbcd83984ea93964adaffc5a14 Mon Sep 17 00:00:00 2001 From: Nico Date: Fri, 6 Feb 2026 14:43:43 +0100 Subject: [PATCH] Punktliste 2 --- GHA_triaxial/gha1_ES.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/GHA_triaxial/gha1_ES.py b/GHA_triaxial/gha1_ES.py index adeda08..fc9a1e8 100644 --- a/GHA_triaxial/gha1_ES.py +++ b/GHA_triaxial/gha1_ES.py @@ -227,7 +227,7 @@ def gha1_ES(ell: EllipsoidTriaxial, beta0: float, omega0: float, alpha0: float, gamma0 = jacobi_konstante(beta, omega, alpha, ell) # Referenz-γ0 - P_all: NDArray[NDArray] = np.array([ell.ell2cart_karney(beta, omega)]) + P_all: List[NDArray] = [ell.ell2cart_karney(beta, omega)] alpha_end: List[float] = [alpha] s_acc = 0.0 @@ -248,12 +248,12 @@ def gha1_ES(ell: EllipsoidTriaxial, beta0: float, omega0: float, alpha0: float, Pk = P_all[-1] alpha1 = float(alpha_end[-1]) - return Pk, alpha1, P_all + return Pk, alpha1, np.array(P_all) if __name__ == "__main__": ell = EllipsoidTriaxial.init_name("BursaSima1980round") - s = 18000 + s = 180000 #alpha0 = 3 alpha0 = wu.gms2rad([5 ,0 ,0]) beta = 0 @@ -267,7 +267,7 @@ if __name__ == "__main__": print(point1) print(res) print(alpha) - #print(points) + print(points) #print("alpha1 (am Endpunkt):", res.alpha1) print(res - point1) print(point1app - point1, "approx") \ No newline at end of file