he2851/Masterprojekt
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Ausgabe der Zwischenpunkte GHA1_num

Browse Source
This commit is contained in:
Tammo.Weber
2026-01-06 15:12:21 +01:00
parent fcae02a0d9
commit cf756e3d9a
1 changed files with 5 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files

7
GHA_triaxial/panou.py
View File

@@ -65,7 +65,7 @@ def gha1_num(ell, point, alpha0, s, num):
werte = rk.rk_chat(F, v_init, s, num) werte = rk.rk_chat(F, v_init, s, num)
x1, _, y1, _, z1, _ = werte[-1] x1, _, y1, _, z1, _ = werte[-1]
return x1, y1, z1 return x1, y1, z1, werte
def checkLiouville(ell: ellipsoide.EllipsoidTriaxial, points): def checkLiouville(ell: ellipsoide.EllipsoidTriaxial, points):
@@ -181,11 +181,14 @@ if __name__ == "__main__":
re = ellipsoide.EllipsoidBiaxial.init_name("Bessel") re = ellipsoide.EllipsoidBiaxial.init_name("Bessel")
# Panou 2013, 7, Table 1, beta0=60° # Panou 2013, 7, Table 1, beta0=60°
beta0, lamb0, beta1, lamb1, c, alpha0, alpha1, s = nep.get_example(table=1, example=5) beta0, lamb0, beta1, lamb1, c, alpha0, alpha1, s = nep.get_example(table=1, example=3)
P0 = ell.ell2cart(beta0, lamb0) P0 = ell.ell2cart(beta0, lamb0)
P1 = ell.ell2cart(beta1, lamb1) P1 = ell.ell2cart(beta1, lamb1)
# P1_num = gha1_num(ell, P0, alpha0, s, 1000) # P1_num = gha1_num(ell, P0, alpha0, s, 1000)
P1_num = gha1_num(ell, P0, alpha0, s, 10000) P1_num = gha1_num(ell, P0, alpha0, s, 10000)
P1_ana = gha1_ana(ell, P0, alpha0, s, 30) P1_ana = gha1_ana(ell, P0, alpha0, s, 30)
beta, lamb = ellipsoide.EllipsoidTriaxial.cart2ell(ell, P1_num)
pass pass