kiam-astro 6.0.7__cp311-cp311-win_amd64.whl → 6.0.8__cp311-cp311-win_amd64.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- kiam_astro/FKIAMToolbox.cp311-win_amd64.pyd +0 -0
- kiam_astro/optimal.py +2 -56
- {kiam_astro-6.0.7.dist-info → kiam_astro-6.0.8.dist-info}/METADATA +1 -1
- {kiam_astro-6.0.7.dist-info → kiam_astro-6.0.8.dist-info}/RECORD +8 -8
- kiam_astro.libs/{libMathAndP.KCBHYULMHMDE63WRVB7DHUMKW2SH6367.gfortran-win_amd64-c9c106af8ff67abea1d3a3b079643807.dll → libMathAndP.NUETR3I7GYNVL2EE46PJHCKVGI4LBW5G.gfortran-win_amd64-7c9d0f2ddc1fc49b42de02176ec5d33c.dll} +0 -0
- {kiam_astro-6.0.7.dist-info → kiam_astro-6.0.8.dist-info}/DELVEWHEEL +0 -0
- {kiam_astro-6.0.7.dist-info → kiam_astro-6.0.8.dist-info}/WHEEL +0 -0
- {kiam_astro-6.0.7.dist-info → kiam_astro-6.0.8.dist-info}/top_level.txt +0 -0
|
Binary file
|
kiam_astro/optimal.py
CHANGED
|
@@ -958,7 +958,7 @@ def r2bp_pontr_mopt_bnd_f_rv(t: float, y: numpy.ndarray, stmreq: bool, fmax: flo
|
|
|
958
958
|
Extended equations of motion:
|
|
959
959
|
Two-body problem, position-velocity variables, bounded thrust force, mass-optimal problem.
|
|
960
960
|
|
|
961
|
-
Parameters
|
|
961
|
+
Parameters:
|
|
962
962
|
-----------
|
|
963
963
|
|
|
964
964
|
`t` : float
|
|
@@ -2026,7 +2026,7 @@ def objective_r2bp_pontr_mass_bnd_f_ee(z: numpy.ndarray, x0: numpy.ndarray, x1:
|
|
|
2026
2026
|
```
|
|
2027
2027
|
"""
|
|
2028
2028
|
|
|
2029
|
-
f = residue_r2bp_pontr_mass_bnd_f_ee(z, x0, x1, fmax, vex, tof, switch_width,
|
|
2029
|
+
f = residue_r2bp_pontr_mass_bnd_f_ee(z, x0, x1, fmax, vex, tof, switch_width, mu, atol, rtol)
|
|
2030
2030
|
return norm(f)**2
|
|
2031
2031
|
def control_r2bp_pontr_mass_bnd_f_ee(y: numpy.ndarray, fmax: float, vex: float, switch_width: float):
|
|
2032
2032
|
"""
|
|
@@ -2337,57 +2337,3 @@ def estimate_tof(x0: numpy.ndarray, x1: numpy.ndarray, umax: float, mu: float =
|
|
|
2337
2337
|
p1 = oe1[0] * (1 - oe1[1] ** 2)
|
|
2338
2338
|
tof = (sqrt(mu / p0) - sqrt(mu / p1)) / umax
|
|
2339
2339
|
return tof
|
|
2340
|
-
|
|
2341
|
-
def test_rv():
|
|
2342
|
-
units = kiam.units('sun')
|
|
2343
|
-
x0 = numpy.array([1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 80.0])
|
|
2344
|
-
x1 = numpy.array([1.5, 0.0, 0.0, 0.0, 1 / numpy.sqrt(1.5), 0.0])
|
|
2345
|
-
fmax = 14.8e-03 / units['AccUnit']
|
|
2346
|
-
vex = 9.3 / units['VelUnit']
|
|
2347
|
-
tof = 2.8 * numpy.pi
|
|
2348
|
-
z0 = numpy.zeros(7, )
|
|
2349
|
-
z0[0:3] = numpy.array([1.41454267e+02, -1.69272674e+00, -4.26704189e-23])
|
|
2350
|
-
z0[3:6] = numpy.array([-1.69272675e+00, 1.41378451e+02, 2.47113618e-23])
|
|
2351
|
-
z0[6] = 5.55630093e-01
|
|
2352
|
-
switch_width = 2.5e-04
|
|
2353
|
-
zopt, res, _ = solve_r2bp_pontr_mass_bnd_f_rv(x0, x1, fmax, vex, tof, switch_width, z0, 1.0, method='least-squares', atol=1e-10, rtol=1e-10)
|
|
2354
|
-
print(f'Оптимальный вектор сопряженных переменных: {zopt}')
|
|
2355
|
-
T, Y = propagate_r2bp_pontr_mass_bnd_f_rv(numpy.linspace(0.0, tof, 10000), numpy.concatenate((x0, zopt)), fmax, vex, switch_width)
|
|
2356
|
-
fig = kiam.plot(Y[0, :], Y[1, :], None, 'x', 'y', axis_equal=True)
|
|
2357
|
-
fig.show()
|
|
2358
|
-
force_value, force_vector = control_r2bp_pontr_mass_bnd_f_rv(Y, fmax, vex, switch_width)
|
|
2359
|
-
fig = kiam.plot(T, force_value, None, 't', 'f')
|
|
2360
|
-
fig.show()
|
|
2361
|
-
h = hamiltonian_r2bp_pontr_mass_bnd_f_rv(Y, fmax, vex, switch_width, 1.0)
|
|
2362
|
-
fig = kiam.plot(T, h, None, 't', 'f')
|
|
2363
|
-
fig.show()
|
|
2364
|
-
print(f'Масса: {Y[6, 0]:.2f} ---> {Y[6, -1]:.2f}, затрачивается {(1 - Y[6, -1] / Y[6, 0]) * 100:.2f}%')
|
|
2365
|
-
def test_ee():
|
|
2366
|
-
units = kiam.units('sun')
|
|
2367
|
-
x0 = numpy.array([1.0, 0.0, 0.0, 0.0, 1.0, 0.0, 80.0])
|
|
2368
|
-
x1 = numpy.array([1.5, 0.0, 0.0, 0.0, 1 / numpy.sqrt(1.5), 0.0])
|
|
2369
|
-
x0[0:6] = kiam.rv2ee(x0[0:6], 1.0, False)
|
|
2370
|
-
x1 = kiam.rv2ee(x1, 1.0, False)
|
|
2371
|
-
fmax = 14.8e-03 / units['AccUnit']
|
|
2372
|
-
vex = 9.3 / units['VelUnit']
|
|
2373
|
-
tof = 2.8 * numpy.pi
|
|
2374
|
-
z0 = numpy.zeros(7, )
|
|
2375
|
-
z0[0:3] = numpy.array([1.41530246e+02, -7.69020160e-02, 1.69380406e+00])
|
|
2376
|
-
z0[5] = 8.25721397e-10
|
|
2377
|
-
z0[6] = 5.55628531e-01
|
|
2378
|
-
switch_width = 0.0
|
|
2379
|
-
zopt, res = solve_r2bp_pontr_mass_bnd_f_ee(x0, x1, fmax, vex, tof, switch_width, z0, 1.0, method='least-squares', atol=1e-10, rtol=1e-10)
|
|
2380
|
-
print(f'Оптимальный вектор сопряженных переменных: {zopt}')
|
|
2381
|
-
T, Y = propagate_r2bp_pontr_mass_bnd_f_ee(numpy.linspace(0.0, tof, 10000), numpy.concatenate((x0, zopt)), fmax, vex, switch_width)
|
|
2382
|
-
RV = kiam.ee2rv(Y[0:6, :], 1.0, False)
|
|
2383
|
-
fig = kiam.plot(RV[0, :], RV[1, :], None, 'x', 'y', axis_equal=True)
|
|
2384
|
-
fig.show()
|
|
2385
|
-
force_value, force_vector = control_r2bp_pontr_mass_bnd_f_ee(Y, fmax, vex, switch_width)
|
|
2386
|
-
fig = kiam.plot(T, force_value, None, 't', 'f')
|
|
2387
|
-
fig.show()
|
|
2388
|
-
h = hamiltonian_r2bp_pontr_mass_bnd_f_ee(Y, fmax, vex, switch_width, 1.0)
|
|
2389
|
-
fig = kiam.plot(T, h, None, 't', 'f')
|
|
2390
|
-
fig.show()
|
|
2391
|
-
print(f'Масса: {Y[6, 0]:.2f} ---> {Y[6, -1]:.2f}, затрачивается {(1 - Y[6, -1] / Y[6, 0]) * 100:.2f}%')
|
|
2392
|
-
|
|
2393
|
-
test_ee()
|
|
@@ -1,8 +1,8 @@
|
|
|
1
1
|
kiam_astro/dace.py,sha256=QhTK6ZNitEr_fIapBo84_CTQLi8ui1_IpsgI-rAkpHU,42721
|
|
2
|
-
kiam_astro/FKIAMToolbox.cp311-win_amd64.pyd,sha256
|
|
2
|
+
kiam_astro/FKIAMToolbox.cp311-win_amd64.pyd,sha256=-U1ki5j6bymG9rTFh5k89iTjDLCNOaNYJ0GaZ6xlpyU,321536
|
|
3
3
|
kiam_astro/JPLEPH,sha256=aB_-XyPHuwyVjtH-a_E0KRL9aQxG5IKdbTBV0yQ_s2g,74468736
|
|
4
4
|
kiam_astro/kiam.py,sha256=171bMkKI6D4HhdBc5lT1CncRe0THdiafiIGhX9I6rGg,210201
|
|
5
|
-
kiam_astro/optimal.py,sha256=
|
|
5
|
+
kiam_astro/optimal.py,sha256=fGGzG9qvdO9n9UqpYd93hKyWiBmkle7LOZ5wVFCbAV8,67625
|
|
6
6
|
kiam_astro/trajectory.py,sha256=b4MdrT-bRJvGgQN2Ln_uaUQchuoRMUKWmxMDNtLaP7o,94477
|
|
7
7
|
kiam_astro/__init__.py,sha256=rKfNulUOTqYZ6q8wwuPTTlkn29Rks_5hLtzdAmc5xFQ,330
|
|
8
8
|
kiam_astro/images/Earth.jpg,sha256=kpalQJ5CV7K6Ch5kOAisL-QdcXNxwwgyoQaDFutEdKs,109661
|
|
@@ -10,11 +10,11 @@ kiam_astro/images/Earth2.jpg,sha256=W1TMWGxsvysodi701AEfbPQieouTpje4GKDFQJDObCw,
|
|
|
10
10
|
kiam_astro/images/Moon - Plate Caree Projection.jpg,sha256=rJ347Lqza2Ya4dzJGGg0xqi89qbBhcIbllXp78zAczY,151314
|
|
11
11
|
kiam_astro/images/Moon1.jpg,sha256=eWyMqBjGonmhVYqU4weAd-n4Nsw3cywScLQxtCFSjcc,541900
|
|
12
12
|
kiam_astro/images/Moon2.jpg,sha256=fv03jz1mpcmlRK4Y7-AayyC6m-wtEIIpDUddORGhhyE,245713
|
|
13
|
-
kiam_astro-6.0.
|
|
14
|
-
kiam_astro-6.0.
|
|
15
|
-
kiam_astro-6.0.
|
|
16
|
-
kiam_astro-6.0.
|
|
17
|
-
kiam_astro-6.0.
|
|
18
|
-
kiam_astro.libs/libMathAndP.
|
|
13
|
+
kiam_astro-6.0.8.dist-info/DELVEWHEEL,sha256=LzS3S8-_bYTE7LnpOxMCES4-H8FVblPsKAzzqawlIOY,222
|
|
14
|
+
kiam_astro-6.0.8.dist-info/METADATA,sha256=Ec6Zil8bxoyZYO2gdoswc7nffRuA6tZrycczqR5bYq4,1321
|
|
15
|
+
kiam_astro-6.0.8.dist-info/RECORD,,
|
|
16
|
+
kiam_astro-6.0.8.dist-info/top_level.txt,sha256=qNlsAYWQr6J_l2w9AhaNb0r02w-JJF1ccO0uB5FVKpo,17
|
|
17
|
+
kiam_astro-6.0.8.dist-info/WHEEL,sha256=U4Q8B8nLmo6Sd27xJLdTpEbQfCLTtXIbLhvUDHobvIE,102
|
|
18
|
+
kiam_astro.libs/libMathAndP.NUETR3I7GYNVL2EE46PJHCKVGI4LBW5G.gfortran-win_amd64-7c9d0f2ddc1fc49b42de02176ec5d33c.dll,sha256=fJ0PLdwfxJtC3gIXbsXTPApt7llPM7nqDfSdEycOtY8,2362893
|
|
19
19
|
tests/test_kiam_basic.py,sha256=OrPukJOxtUecKzQNKKKXulaS8W6ocYHeIQoBTBp877s,26838
|
|
20
20
|
tests/__init__.py,sha256=rKfNulUOTqYZ6q8wwuPTTlkn29Rks_5hLtzdAmc5xFQ,330
|
|
index 4f728d1..4802359 100644
|
|
|
Binary file
|
|
File without changes
|
|
File without changes
|
|
File without changes
|