aboutsummaryrefslogtreecommitdiffstatshomepage
path: root/schedule.c
diff options
context:
space:
mode:
Diffstat (limited to 'schedule.c')
-rw-r--r--schedule.c270
1 files changed, 0 insertions, 270 deletions
diff --git a/schedule.c b/schedule.c
deleted file mode 100644
index d13a9c9..0000000
--- a/schedule.c
+++ /dev/null
@@ -1,270 +0,0 @@
-/* SPDX-License-Identifier: MIT-0 */
-
-#include <assert.h>
-#include <math.h>
-
-#include "schedule.h"
-
-#define SWITCH_INTERMEDIATE 1
-
-#define STOPLEN_ROUGH 10
-#define STOPLEN_FINE 1
-
-const gas_t *best_gas(const double depth, const gas_t *gasses, const int nof_gasses)
-{
- const gas_t *best = NULL;
- double mod_best = -1;
-
- for (int i = 0; i < nof_gasses; i++) {
- double mod = gas_mod(&gasses[i]);
-
- if (depth <= mod && (mod_best == -1 || mod < mod_best)) {
- best = &gasses[i];
- mod_best = mod;
- }
- }
-
- return best;
-}
-
-const gas_t *next_gas(const double depth, const gas_t *gasses, const int nof_gasses)
-{
- const gas_t *next = NULL;
- double mod_best = 0;
-
- for (int i = 0; i < nof_gasses; i++) {
- double mod = gas_mod(&gasses[i]);
-
- if (depth > mod && mod > mod_best) {
- next = &gasses[i];
- mod_best = mod;
- }
- }
-
- return next;
-}
-
-int direct_ascent(const decostate_t *ds, const double depth, const double time, const gas_t *gas)
-{
- decostate_t ds_ = *ds;
- assert(ds_.firststop == -1);
-
- add_segment_ascdec(&ds_, depth, abs_depth(0), time, gas);
-
- return gauge_depth(ceiling(&ds_, ds_.gfhi)) <= 0;
-}
-
-void simulate_dive(decostate_t *ds, waypoint_t *waypoints, const int nof_waypoints, waypoint_callback_t wp_cb)
-{
- double depth = abs_depth(0);
- double runtime = 0;
-
- for (int i = 0; i < nof_waypoints; i++) {
- double d = waypoints[i].depth;
- double t = waypoints[i].time;
- const gas_t *g = waypoints[i].gas;
-
- if (d != depth) {
- runtime += add_segment_ascdec(ds, depth, d, t, g);
- depth = d;
- } else {
- runtime += add_segment_const(ds, d, t, g);
- }
-
- if (wp_cb)
- wp_cb(ds, (waypoint_t){.depth = d, .time = t, .gas = g}, SEG_DIVE);
- }
-}
-
-double calc_ndl(decostate_t *ds, const double depth, const double ascrate, const gas_t *gas)
-{
- decostate_t ds_ = *ds;
- double ndl = 0;
-
- while (ndl < 360) {
- double tmp = add_segment_const(&ds_, depth, 1, gas);
-
- if (!direct_ascent(&ds_, depth, gauge_depth(depth) / ascrate, gas))
- break;
-
- ndl += tmp;
- }
-
- return ndl;
-}
-
-double calc_stoplen_rough(const decostate_t *ds, const double depth, const double current_gf, const gas_t *gas)
-{
- decostate_t ds_ = *ds;
- double stoplen = 0;
-
- for (;;) {
- double tmp = add_segment_const(&ds_, depth, STOPLEN_ROUGH, gas);
-
- if (ceiling(&ds_, current_gf) != depth)
- break;
-
- stoplen += tmp;
- }
-
- return stoplen;
-}
-
-double deco_stop(decostate_t *ds, const double depth, const double current_gf, const gas_t *gas)
-{
- double stoplen = 0;
-
- /* rough steps */
- double stoplen_rough = calc_stoplen_rough(ds, depth, current_gf, gas);
-
- if (stoplen_rough) {
- add_segment_const(ds, depth, stoplen_rough, gas);
- stoplen += stoplen_rough;
- }
-
- /* fine steps */
- while (ceiling(ds, current_gf) == depth)
- stoplen += add_segment_const(ds, depth, STOPLEN_FINE, gas);
-
- return stoplen;
-}
-
-int should_switch(const gas_t **next, double *switch_depth, const decostate_t *ds, const double depth,
- const double next_stop, const gas_t *deco_gasses, const int nof_gasses)
-{
- /* check if we switch at MOD or just at stops */
- if (!SWITCH_INTERMEDIATE)
- return 0;
-
- /* check if there is a gas to switch to */
- *next = next_gas(depth, deco_gasses, nof_gasses);
-
- if (*next == NULL)
- return 0;
-
- /* check if the switch happens before the current ceiling */
- *switch_depth = round_ceiling(ds, (*next)->mod) - ds->ceil_multiple;
- assert(*switch_depth <= (*next)->mod);
-
- if (*switch_depth <= next_stop)
- return 0;
-
- return 1;
-}
-
-decoinfo_t calc_deco(decostate_t *ds, const double start_depth, const gas_t *start_gas, const gas_t *deco_gasses,
- const int nof_gasses, waypoint_callback_t wp_cb)
-{
- decoinfo_t ret = {.tts = 0, .ndl = 0};
-
- /* setup start parameters */
- double depth = start_depth;
- const gas_t *gas = start_gas;
-
- const double asc_per_min = msw_to_bar(9);
-
- /* check if direct ascent is possible */
- if (direct_ascent(ds, depth, gauge_depth(depth) / asc_per_min, gas)) {
- ret.ndl = calc_ndl(ds, depth, asc_per_min, gas);
- return ret;
- }
-
- /* determine first stop */
- double current_gf = get_gf(ds, depth);
-
- if (ds->firststop == -1)
- ds->firststop = ceiling(ds, current_gf);
-
- /* switch to best deco gas if there is one available */
- const gas_t *best = best_gas(depth, deco_gasses, nof_gasses);
-
- if (best)
- gas = best;
-
- /* alternate between ascending and stopping */
- for (;;) {
- /* extra bookkeeping because waypoints and segments do not match 1:1 */
- double last_waypoint_depth = depth;
- double waypoint_time;
-
- /* ascend */
- for (;;) {
- /* determine next stop */
- double next_stop = ceiling(ds, current_gf);
-
- /* find out if we need to switch gas on the way */
- const gas_t *next;
- double switch_depth;
-
- if (should_switch(&next, &switch_depth, ds, depth, next_stop, deco_gasses, nof_gasses)) {
- /* ascend to gas switch */
- ret.tts += add_segment_ascdec(ds, depth, switch_depth, fabs(depth - switch_depth) / asc_per_min, gas);
- depth = switch_depth;
- current_gf = get_gf(ds, depth);
-
- /*
- * since we're stopping for a switch next, this is the last of
- * any number of consecutive travel segments and the waypoint
- * callback should be called.
- */
- waypoint_time = fabs(last_waypoint_depth - depth) / asc_per_min;
-
- if (wp_cb)
- wp_cb(ds, (waypoint_t){.depth = depth, .time = waypoint_time, .gas = gas}, SEG_TRAVEL);
-
- last_waypoint_depth = depth;
-
- /* switch gas */
- gas = next;
-
- ret.tts += add_segment_const(ds, switch_depth, 1, gas);
-
- if (wp_cb)
- wp_cb(ds, (waypoint_t){.depth = depth, .time = 1, .gas = gas}, SEG_GAS_SWITCH);
-
- continue;
- }
-
- /* ascend to current ceiling */
- ret.tts += add_segment_ascdec(ds, depth, next_stop, fabs(depth - next_stop) / asc_per_min, gas);
- depth = next_stop;
- current_gf = get_gf(ds, depth);
-
- /* if the ceiling moved while we ascended, keep ascending */
- if (depth > ceiling(ds, current_gf))
- continue;
-
- /*
- * since we've actually reached the ceiling, this is the last of
- * any number of consecutive travel segments and the waypoint
- * callback should be called.
- */
- waypoint_time = fabs(last_waypoint_depth - depth) / asc_per_min;
- enum segtype_t segtype = depth <= abs_depth(0) ? SEG_SURFACE : SEG_TRAVEL;
-
- if (wp_cb)
- wp_cb(ds, (waypoint_t){.depth = depth, .time = waypoint_time, .gas = gas}, segtype);
-
- break;
- }
-
- /* terminate if we surfaced */
- if (depth <= abs_depth(0))
- break;
-
- /* switch to better gas if available */
- const gas_t *best = best_gas(depth, deco_gasses, nof_gasses);
-
- if (best)
- gas = best;
-
- /* stop */
- double stoplen = deco_stop(ds, depth, current_gf, gas);
- ret.tts += stoplen;
-
- if (wp_cb)
- wp_cb(ds, (waypoint_t){.depth = depth, .time = stoplen, .gas = gas}, SEG_DECO_STOP);
- }
-
- return ret;
-}
generated by cgit v1.2.3-70-g09d2 (git 2.47.0) at 2025-01-12 21:59:09 +0000