forked from ATTron/astroz
-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathtransfer_propagation.zig
More file actions
104 lines (91 loc) · 4.15 KB
/
transfer_propagation.zig
File metadata and controls
104 lines (91 loc) · 4.15 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
const std = @import("std");
const astroz = @import("astroz");
const Mission = astroz.Mission;
const OrbitalMechanics = astroz.OrbitalMechanics;
const constants = astroz.constants;
pub fn main() !void {
var gpa = std.heap.GeneralPurposeAllocator(.{}){};
defer _ = gpa.deinit();
const allocator = gpa.allocator();
const orbitalMechanics = OrbitalMechanics.init(constants.sun.mu, constants.sun);
std.debug.print("=== Earth to Mars Transfer ===\n", .{});
const earthMarsParams = try Mission.MissionParameters.init(constants.earth, constants.mars, 0.0, null, "hohmann");
var earthMarsMission = Mission.init(allocator, earthMarsParams, orbitalMechanics);
defer earthMarsMission.deinit();
try earthMarsMission.propagateTransfer(520.0, 5.0);
std.debug.print("Generated {d} trajectory points\n", .{earthMarsMission.trajectoryPredictions.items.len});
std.debug.print("Waypoints:\n", .{});
for (earthMarsMission.trajectoryPredictions.items) |point| {
if (std.mem.eql(u8, point.label, "waypoint")) {
std.debug.print("Day {d:>6.1}: {s} at ({d:>8.0}, {d:>8.0}) km\n", .{
point.time,
point.body,
point.position.x(),
point.position.y(),
});
}
}
std.debug.print("\n=== Mars to Jupiter Transfer ===\n", .{});
const marsJupiterParams = try Mission.MissionParameters.init(constants.mars, constants.jupiter, 0.0, null, "hohmann");
var marsJupiterMission = Mission.init(allocator, marsJupiterParams, orbitalMechanics);
defer marsJupiterMission.deinit();
try marsJupiterMission.propagateTransfer(1200.0, 10.0);
std.debug.print("Generated {d} trajectory points\n", .{marsJupiterMission.trajectoryPredictions.items.len});
std.debug.print("Waypoints:\n", .{});
for (marsJupiterMission.trajectoryPredictions.items) |point| {
if (std.mem.eql(u8, point.label, "waypoint")) {
std.debug.print("Day {d:>6.1}: {s} at ({d:>8.0}, {d:>8.0}) km\n", .{
point.time,
point.body,
point.position.x(),
point.position.y(),
});
}
}
std.debug.print("\n=== Venus to Earth Transfer ===\n", .{});
const venusEarthParams = try Mission.MissionParameters.init(constants.venus, constants.earth, 0.0, null, "hohmann");
var venusEarthMission = Mission.init(allocator, venusEarthParams, orbitalMechanics);
defer venusEarthMission.deinit();
try venusEarthMission.propagateTransfer(400.0, 2.0);
std.debug.print("Generated {d} trajectory points\n", .{venusEarthMission.trajectoryPredictions.items.len});
std.debug.print("Waypoints:\n", .{});
for (venusEarthMission.trajectoryPredictions.items) |point| {
if (std.mem.eql(u8, point.label, "waypoint")) {
std.debug.print("Day {d:>6.1}: {s} at ({d:>8.0}, {d:>8.0}) km\n", .{
point.time,
point.body,
point.position.x(),
point.position.y(),
});
}
}
std.debug.print("\n=== Programmatic Access Example ===\n", .{});
std.debug.print("First 5 Earth-Mars transfer points:\n", .{});
std.debug.print("Time (days) | Body | X (km) | Y (km) | Label\n", .{});
std.debug.print("-----------|------|--------|--------|----------\n", .{});
var count: usize = 0;
for (earthMarsMission.trajectoryPredictions.items) |point| {
if (count < 5) {
std.debug.print("{d:>10.1} | {s:<8} | {d:>8.0} | {d:>8.0} | {s}\n", .{
point.time,
point.body,
point.position.x(),
point.position.y(),
point.label,
});
count += 1;
}
}
std.debug.print("\nTransfer trajectory points (first 5):\n", .{});
count = 0;
for (earthMarsMission.trajectoryPredictions.items) |point| {
if (std.mem.eql(u8, point.body, "Transfer") and count < 5) {
std.debug.print("Day {d:>6.1}: Transfer at ({d:>8.0}, {d:>8.0}) km\n", .{
point.time,
point.position.x(),
point.position.y(),
});
count += 1;
}
}
}