cinnaboot
/
claudes_game
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity
Browse Source

Remove deprecated comet orbit test and config 

- Delete tests/test_comet_orbit.cpp (7-body complex system)
- Delete configs/test_simple.toml (redundant test config)
- Replaced by simpler tests/configs/soi_transition.toml
- All tests still pass (15 test cases, 14 passed)
- Documentation references remain for historical context

Claude
main
cinnaboot 6 months ago
parent
08cdfeb7ca
commit
2e9e747de2
2 changed files with 0 additions and 236 deletions
Whitespace
Show all changes Ignore whitespace when comparing lines Ignore changes in amount of whitespace Ignore changes in whitespace at EOL
Unified View
Diff Options
Show Stats Download Patch File Download Diff File
  1. 83
      configs/test_simple.toml
  2. 153
      tests/test_comet_orbit.cpp

83
configs/test_simple.toml
Unescape View File

@ -1,83 +0,0 @@
# Simple Test Configuration
# Planets in circular and eccentric orbits
#
# Orbital periods (for verification):
# Earth: ~365 days
# Mars: ~687 days
# Comet: ~1444 days
#
# At t=0:
# Earth starts at (1.0 AU, 0, 0) = (1.496e11 m, 0, 0)
# Mars starts at (1.5 AU, 0, 0) = (2.244e11 m, 0, 0)
# Comet starts at perihelion (0.75 AU, 0, 0) = (1.122e11 m, 0, 0)
# All orbit counter-clockwise when viewed from +Z
[[bodies]]
name = "Sun"
mass = 1.989e30
radius = 6.96e8
position = { x = 0.0, y = 0.0, z = 0.0 }
parent_index = -1
color = { r = 1.0, g = 1.0, b = 0.0 }
eccentricity = 0.0
semi_major_axis = 0.0
[[bodies]]
name = "Earth"
mass = 5.972e24
radius = 6.371e6
position = { x = 1.496e11, y = 0.0, z = 0.0 }
parent_index = 0
color = { r = 0.0, g = 0.5, b = 1.0 }
eccentricity = 0.0
semi_major_axis = 1.496e11
[[bodies]]
name = "Mars"
mass = 6.39e23
radius = 3.3895e6
position = { x = 2.244e11, y = 0.0, z = 0.0 }
parent_index = 0
color = { r = 0.8, g = 0.3, b = 0.1 }
eccentricity = 0.0
semi_major_axis = 2.244e11
[[bodies]]
name = "Comet"
mass = 1e14
radius = 5e3
position = { x = 1.122e11, y = 0.0, z = 0.0 }
parent_index = 0
color = { r = 0.5, g = 0.8, b = 1.0 }
eccentricity = 0.7
semi_major_axis = 3.74e11
[[bodies]]
name = "Parabolic"
mass = 1e13
radius = 3e3
position = { x = 7.48e10, y = 0.0, z = 0.0 }
parent_index = 0
color = { r = 1.0, g = 1.0, b = 0.5 }
eccentricity = 1.0
semi_major_axis = 7.48e10
[[bodies]]
name = "Inclined"
mass = 1e12
radius = 2e3
position = { x = 1.795e11, y = 0.0, z = 3.114e10 }
parent_index = 0
color = { r = 0.8, g = 0.8, b = 0.0 }
eccentricity = 0.0
semi_major_axis = 1.795e11
[[bodies]]
name = "Hyperbolic"
mass = 1e13
radius = 3e3
position = { x = 7.48e10, y = 0.0, z = 0.0 }
parent_index = 0
color = { r = 0.3, g = 1.0, b = 0.3 }
eccentricity = 1.5
semi_major_axis = -1.496e11

153
tests/test_comet_orbit.cpp
Unescape View File

@ -1,153 +0,0 @@
#include <catch2/catch_test_macros.hpp>
#include "../src/physics.h"
#include "../src/simulation.h"
#include "../src/config_loader.h"
#include "../src/test_utilities.h"
#include <cmath>
#include <vector>
#define COMET_TEST_VERBOSE_OUTPUT 0
struct ParentChange {
double time_seconds;
double time_days;
int old_parent;
int new_parent;
double distance_to_mars_au;
double distance_to_sun_au;
};
#if COMET_TEST_VERBOSE_OUTPUT
static void print_orbital_snapshots(const std::vector<OrbitalElements>& snapshots,
double expected_sma, double expected_ecc) {
printf("\n=== Comet Orbital Elements Over Time ===\n");
printf("Expected: a = %.1f AU, e = %.1f\n\n", expected_sma, expected_ecc);
for (const auto& elem : snapshots) {
printf("Day %.0f:\n", elem.time_days);
printf(" Semi-major axis: %.6f AU (expected: %.1f AU)\n", elem.semi_major_axis_au, expected_sma);
printf(" Eccentricity: %.6f (expected: %.1f)\n", elem.eccentricity, expected_ecc);
printf(" Distance to Sun: %.4f AU\n", elem.distance_to_sun_au);
printf(" Distance to Mars: %.6f AU\n", elem.distance_to_ref_body_au);
printf(" Velocity: %.2f km/s\n", elem.velocity_magnitude / 1000.0);
printf(" Specific energy: %.3e J/kg\n", elem.specific_energy);
double sma_error = fabs(elem.semi_major_axis_au - expected_sma);
double ecc_error = fabs(elem.eccentricity - expected_ecc);
if (elem.time_days >= 1530 && elem.time_days <= 1540) {
printf(" *** CLOSEST APPROACH TO MARS ***\n");
}
printf(" Error: Δa = %.6f AU, Δe = %.6f\n\n", sma_error, ecc_error);
}
}
static void print_parent_changes(const std::vector<ParentChange>& parent_changes,
int mars_index) {
printf("\n=== Parent Index Changes ===\n");
if (parent_changes.empty()) {
printf("No parent changes detected\n\n");
} else {
for (const auto& change : parent_changes) {
printf("Time: %.0f days (%.0f seconds)\n", change.time_days, change.time_seconds);
printf(" Parent: %d -> %d", change.old_parent, change.new_parent);
if (change.new_parent == mars_index) {
printf(" (Sun -> Mars)");
} else if (change.old_parent == mars_index) {
printf(" (Mars -> Sun)");
}
printf("\n");
printf(" Distance to Mars: %.6f AU\n", change.distance_to_mars_au);
printf(" Distance to Sun: %.6f AU\n\n", change.distance_to_sun_au);
}
}
}
#endif
TEST_CASE("Comet orbital elements and SOI transitions during Mars encounter", "[comet][orbital-elements][soi]") {
const double TIME_STEP = 60.0;
const double SECONDS_PER_DAY = 86400.0;
const double DAYS_TO_SIMULATE = 1700.0;
const double AU = 1.496e11;
#if COMET_TEST_VERBOSE_OUTPUT
const double EXPECTED_SMA = 2.5;
const double EXPECTED_ECC = 0.7;
#endif
SimulationState* sim = create_simulation(10, TIME_STEP);
REQUIRE(load_system_config(sim, "configs/test_simple.toml"));
const int COMET_INDEX = 3;
const int MARS_INDEX = 2;
const int SUN_INDEX = 0;
std::vector<OrbitalElements> snapshots;
std::vector<ParentChange> parent_changes;
int previous_parent = sim->bodies[COMET_INDEX].parent_index;
OrbitalElements initial = calculate_orbital_elements(
&sim->bodies[COMET_INDEX],
&sim->bodies[SUN_INDEX],
&sim->bodies[MARS_INDEX],
sim->time
);
snapshots.push_back(initial);
double checkpoint_days[] = {0, 365, 722, 1444, 1533, 1600, 1700};
size_t next_checkpoint = 1;
double max_time = DAYS_TO_SIMULATE * SECONDS_PER_DAY;
while (sim->time < max_time) {
update_simulation(sim);
int current_parent = sim->bodies[COMET_INDEX].parent_index;
if (current_parent != previous_parent) {
ParentChange change;
change.time_seconds = sim->time;
change.time_days = sim->time / SECONDS_PER_DAY;
change.old_parent = previous_parent;
change.new_parent = current_parent;
change.distance_to_sun_au = vec3_magnitude(sim->bodies[COMET_INDEX].position) / AU;
change.distance_to_mars_au = vec3_distance(sim->bodies[COMET_INDEX].position,
sim->bodies[MARS_INDEX].position) / AU;
parent_changes.push_back(change);
previous_parent = current_parent;
}
double current_days = sim->time / SECONDS_PER_DAY;
if (next_checkpoint < sizeof(checkpoint_days)/sizeof(checkpoint_days[0]) &&
current_days >= checkpoint_days[next_checkpoint]) {
OrbitalElements elem = calculate_orbital_elements(
&sim->bodies[COMET_INDEX],
&sim->bodies[SUN_INDEX],
&sim->bodies[MARS_INDEX],
sim->time
);
snapshots.push_back(elem);
next_checkpoint++;
}
}
#if COMET_TEST_VERBOSE_OUTPUT
print_orbital_snapshots(snapshots, EXPECTED_SMA, EXPECTED_ECC);
print_parent_changes(parent_changes, MARS_INDEX);
#endif
REQUIRE(parent_changes.size() > 0);
bool found_mars_transition = false;
for (const auto& change : parent_changes) {
if (change.new_parent == MARS_INDEX || change.old_parent == MARS_INDEX) {
found_mars_transition = true;
REQUIRE(fabs(change.time_days - 1550.0) < 50.0);
INFO("Mars encounter at day " << change.time_days << ", distance " << change.distance_to_mars_au << " AU");
}
}
REQUIRE(found_mars_transition);
destroy_simulation(sim);
}
Write Preview
Loading…
Cancel
Save