mirror of
https://github.com/subsurface/subsurface.git
synced 2024-11-28 05:00:20 +00:00
56535e4b3b
In add dive mode simply bring the diver safely back to the surface (currently with a fixed ascent rate of 30ft/min (or 9m/min)). We should make that rate configurable (for the planner as well as the dive add function). Also, the dive add function should offer to automatically include a safety stop. Signed-off-by: Dirk Hohndel <dirk@hohndel.org>
1033 lines
26 KiB
C
1033 lines
26 KiB
C
/* planner.c
|
|
*
|
|
* code that allows us to plan future dives
|
|
*
|
|
* (c) Dirk Hohndel 2013
|
|
*/
|
|
#include <libintl.h>
|
|
#include <glib/gi18n.h>
|
|
#include <unistd.h>
|
|
#include <ctype.h>
|
|
#include "dive.h"
|
|
#include "divelist.h"
|
|
#include "planner.h"
|
|
|
|
int decostoplevels[] = { 0, 3000, 6000, 9000, 12000, 15000, 18000, 21000, 24000, 27000,
|
|
30000, 33000, 36000, 39000, 42000, 45000, 48000, 51000, 54000, 57000,
|
|
60000, 63000, 66000, 69000, 72000, 75000, 78000, 81000, 84000, 87000,
|
|
90000, 100000, 110000, 120000, 130000, 140000, 150000, 160000, 170000,
|
|
180000, 190000, 200000, 220000, 240000, 260000, 280000, 300000,
|
|
320000, 340000, 360000, 380000
|
|
};
|
|
double plangflow, plangfhigh;
|
|
char *disclaimer;
|
|
|
|
#if DEBUG_PLAN
|
|
void dump_plan(struct diveplan *diveplan)
|
|
{
|
|
struct divedatapoint *dp;
|
|
struct tm tm;
|
|
|
|
if (!diveplan) {
|
|
printf ("Diveplan NULL\n");
|
|
return;
|
|
}
|
|
utc_mkdate(diveplan->when, &tm);
|
|
|
|
printf("\nDiveplan @ %04d-%02d-%02d %02d:%02d:%02d (surfpres %dmbar):\n",
|
|
tm.tm_year+1900, tm.tm_mon+1, tm.tm_mday,
|
|
tm.tm_hour, tm.tm_min, tm.tm_sec,
|
|
diveplan->surface_pressure);
|
|
dp = diveplan->dp;
|
|
while (dp) {
|
|
printf("\t%3u:%02u: %dmm gas: %d o2 %d h2\n", FRACTION(dp->time, 60), dp->depth, dp->o2, dp->he);
|
|
dp = dp->next;
|
|
}
|
|
|
|
}
|
|
#endif
|
|
|
|
void set_last_stop(gboolean last_stop_6m)
|
|
{
|
|
if (last_stop_6m == TRUE)
|
|
decostoplevels[1] = 6000;
|
|
else
|
|
decostoplevels[1] = 3000;
|
|
}
|
|
|
|
void get_gas_from_events(struct divecomputer *dc, int time, int *o2, int *he)
|
|
{
|
|
struct event *event = dc->events;
|
|
while (event && event->time.seconds <= time) {
|
|
if (!strcmp(event->name, "gaschange")) {
|
|
*o2 = 10 * event->value & 0xffff;
|
|
*he = 10 * event->value >> 16;
|
|
}
|
|
event = event->next;
|
|
}
|
|
}
|
|
|
|
/* simple helper function to compare two permille values with
|
|
* (rounded) percent granularity */
|
|
static inline gboolean match_percent(int a, int b)
|
|
{
|
|
return (a + 5) / 10 == (b + 5) / 10;
|
|
}
|
|
|
|
static int get_gasidx(struct dive *dive, int o2, int he)
|
|
{
|
|
int gasidx = -1;
|
|
|
|
/* we treat air as 0/0 because it is special */
|
|
if (is_air(o2, he))
|
|
o2 = 0;
|
|
while (++gasidx < MAX_CYLINDERS)
|
|
if (match_percent(dive->cylinder[gasidx].gasmix.o2.permille, o2) &&
|
|
match_percent(dive->cylinder[gasidx].gasmix.he.permille, he))
|
|
return gasidx;
|
|
return -1;
|
|
}
|
|
|
|
void get_gas_string(int o2, int he, char *text, int len)
|
|
{
|
|
if (is_air(o2, he))
|
|
snprintf(text, len, _("air"));
|
|
else if (he == 0)
|
|
snprintf(text, len, _("EAN%d"), (o2 + 5) / 10);
|
|
else
|
|
snprintf(text, len, "(%d/%d)", (o2 + 5) / 10, (he + 5) / 10);
|
|
}
|
|
|
|
/* returns the tissue tolerance at the end of this (partial) dive */
|
|
double tissue_at_end(struct dive *dive, char **cached_datap, char **error_string_p)
|
|
{
|
|
struct divecomputer *dc;
|
|
struct sample *sample, *psample;
|
|
int i, j, t0, t1, gasidx, lastdepth;
|
|
int o2, he;
|
|
double tissue_tolerance;
|
|
static char buf[200];
|
|
|
|
if (!dive)
|
|
return 0.0;
|
|
if (*cached_datap) {
|
|
tissue_tolerance = restore_deco_state(*cached_datap);
|
|
} else {
|
|
tissue_tolerance = init_decompression(dive);
|
|
cache_deco_state(tissue_tolerance, cached_datap);
|
|
}
|
|
dc = &dive->dc;
|
|
if (!dc->samples)
|
|
return tissue_tolerance;
|
|
psample = sample = dc->sample;
|
|
lastdepth = t0 = 0;
|
|
/* we always start with gas 0 (unless an event tells us otherwise) */
|
|
o2 = dive->cylinder[0].gasmix.o2.permille;
|
|
he = dive->cylinder[0].gasmix.he.permille;
|
|
for (i = 0; i < dc->samples; i++, sample++) {
|
|
t1 = sample->time.seconds;
|
|
get_gas_from_events(&dive->dc, t0, &o2, &he);
|
|
if ((gasidx = get_gasidx(dive, o2, he)) == -1) {
|
|
snprintf(buf, sizeof(buf),_("Can't find gas %d/%d"), (o2 + 5) / 10, (he + 5) / 10);
|
|
*error_string_p = buf;
|
|
gasidx = 0;
|
|
}
|
|
if (i > 0)
|
|
lastdepth = psample->depth.mm;
|
|
for (j = t0; j < t1; j++) {
|
|
int depth = interpolate(lastdepth, sample->depth.mm, j - t0, t1 - t0);
|
|
tissue_tolerance = add_segment(depth_to_mbar(depth, dive) / 1000.0,
|
|
&dive->cylinder[gasidx].gasmix, 1, sample->po2, dive);
|
|
}
|
|
psample = sample;
|
|
t0 = t1;
|
|
}
|
|
return tissue_tolerance;
|
|
}
|
|
|
|
/* how many seconds until we can ascend to the next stop? */
|
|
int time_at_last_depth(struct dive *dive, int o2, int he, int next_stop, char **cached_data_p, char **error_string_p)
|
|
{
|
|
int depth, gasidx;
|
|
double surface_pressure, tissue_tolerance;
|
|
int wait = 0;
|
|
struct sample *sample;
|
|
|
|
if (!dive)
|
|
return 0;
|
|
surface_pressure = dive->dc.surface_pressure.mbar / 1000.0;
|
|
tissue_tolerance = tissue_at_end(dive, cached_data_p, error_string_p);
|
|
sample = &dive->dc.sample[dive->dc.samples - 1];
|
|
depth = sample->depth.mm;
|
|
gasidx = get_gasidx(dive, o2, he);
|
|
while (deco_allowed_depth(tissue_tolerance, surface_pressure, dive, 1) > next_stop) {
|
|
wait++;
|
|
tissue_tolerance = add_segment(depth_to_mbar(depth, dive) / 1000.0,
|
|
&dive->cylinder[gasidx].gasmix, 1, sample->po2, dive);
|
|
}
|
|
return wait;
|
|
}
|
|
|
|
int add_gas(struct dive *dive, int o2, int he)
|
|
{
|
|
int i;
|
|
struct gasmix *mix;
|
|
cylinder_t *cyl;
|
|
|
|
for (i = 0; i < MAX_CYLINDERS; i++) {
|
|
cyl = dive->cylinder + i;
|
|
mix = &cyl->gasmix;
|
|
if (cylinder_nodata(cyl))
|
|
break;
|
|
if (match_percent(o2, mix->o2.permille) && match_percent(he, mix->he.permille))
|
|
return i;
|
|
}
|
|
if (i == MAX_CYLINDERS) {
|
|
return -1;
|
|
}
|
|
mix->o2.permille = o2;
|
|
mix->he.permille = he;
|
|
/* since air is stored as 0/0 we need to set a name or an air cylinder
|
|
* would be seen as unset (by cylinder_nodata()) */
|
|
cyl->type.description = strdup(_("Cylinder for planning"));
|
|
return i;
|
|
}
|
|
|
|
struct dive *create_dive_from_plan(struct diveplan *diveplan, char **error_string)
|
|
{
|
|
struct dive *dive;
|
|
struct divedatapoint *dp;
|
|
struct divecomputer *dc;
|
|
struct sample *sample;
|
|
int oldo2 = O2_IN_AIR, oldhe = 0;
|
|
int oldpo2 = 0;
|
|
int lasttime = 0;
|
|
|
|
*error_string = NULL;
|
|
if (!diveplan || !diveplan->dp)
|
|
return NULL;
|
|
#if DEBUG_PLAN & 4
|
|
printf("in create_dive_from_plan\n");
|
|
dump_plan(diveplan);
|
|
#endif
|
|
dive = alloc_dive();
|
|
dive->when = diveplan->when;
|
|
dive->dc.surface_pressure.mbar = diveplan->surface_pressure;
|
|
dc = &dive->dc;
|
|
dc->model = strdup(_("Simulated Dive"));
|
|
dp = diveplan->dp;
|
|
|
|
/* let's start with the gas given on the first segment */
|
|
if (dp->o2 || dp->he) {
|
|
oldo2 = dp->o2;
|
|
oldhe = dp->he;
|
|
}
|
|
sample = prepare_sample(dc);
|
|
sample->po2 = dp->po2;
|
|
finish_sample(dc);
|
|
if (add_gas(dive, oldo2, oldhe) < 0)
|
|
goto gas_error_exit;
|
|
while (dp) {
|
|
int o2 = dp->o2, he = dp->he;
|
|
int po2 = dp->po2;
|
|
int time = dp->time;
|
|
int depth = dp->depth;
|
|
|
|
if (time == 0) {
|
|
/* special entries that just inform the algorithm about
|
|
* additional gases that are available */
|
|
if (add_gas(dive, o2, he) < 0)
|
|
goto gas_error_exit;
|
|
dp = dp->next;
|
|
continue;
|
|
}
|
|
if (!o2 && !he) {
|
|
o2 = oldo2;
|
|
he = oldhe;
|
|
}
|
|
|
|
/* Check for SetPoint change */
|
|
if (oldpo2 != po2) {
|
|
if (lasttime)
|
|
add_event(dc, lasttime, 20, 0, po2, "SP change"); // SAMPLE_EVENT_PO2
|
|
oldpo2 = po2;
|
|
}
|
|
|
|
/* Create new gas, and gas change event if necessary;
|
|
* Sadly, we inherited our gaschange event from libdivecomputer which only
|
|
* support percentage values, so round the entries */
|
|
if (o2 != oldo2 || he != oldhe) {
|
|
int plano2 = (o2 + 5) / 10 * 10;
|
|
int planhe = (he + 5) / 10 * 10;
|
|
int value;
|
|
if (add_gas(dive, plano2, planhe) < 0)
|
|
goto gas_error_exit;
|
|
value = (plano2 / 10) | ((planhe / 10) << 16);
|
|
add_event(dc, lasttime, 25, 0, value, "gaschange"); // SAMPLE_EVENT_GASCHANGE2
|
|
oldo2 = o2; oldhe = he;
|
|
}
|
|
/* Create sample */
|
|
sample = prepare_sample(dc);
|
|
/* set po2 at beginning of this segment */
|
|
/* and keep it valid for last sample - where it likely doesn't matter */
|
|
sample[-1].po2 = po2;
|
|
sample->po2 = po2;
|
|
sample->time.seconds = time;
|
|
sample->depth.mm = depth;
|
|
finish_sample(dc);
|
|
lasttime = time;
|
|
dp = dp->next;
|
|
}
|
|
if (dc->samples <= 1) {
|
|
/* not enough there yet to create a dive - most likely the first time is missing */
|
|
free(dive);
|
|
dive = NULL;
|
|
}
|
|
#if DEBUG_PLAN & 32
|
|
if (dive)
|
|
save_dive(stdout, dive);
|
|
#endif
|
|
return dive;
|
|
|
|
gas_error_exit:
|
|
free(dive);
|
|
*error_string = _("Too many gas mixes");
|
|
return NULL;
|
|
}
|
|
|
|
void free_dps(struct divedatapoint *dp)
|
|
{
|
|
while (dp) {
|
|
struct divedatapoint *ndp = dp->next;
|
|
free(dp);
|
|
dp = ndp;
|
|
}
|
|
}
|
|
|
|
struct divedatapoint *create_dp(int time_incr, int depth, int o2, int he, int po2)
|
|
{
|
|
struct divedatapoint *dp;
|
|
|
|
dp = malloc(sizeof(struct divedatapoint));
|
|
dp->time = time_incr;
|
|
dp->depth = depth;
|
|
dp->o2 = o2;
|
|
dp->he = he;
|
|
dp->po2 = po2;
|
|
dp->entered = FALSE;
|
|
dp->next = NULL;
|
|
return dp;
|
|
}
|
|
|
|
struct divedatapoint *get_nth_dp(struct diveplan *diveplan, int idx)
|
|
{
|
|
struct divedatapoint **ldpp, *dp = diveplan->dp;
|
|
int i = 0;
|
|
ldpp = &diveplan->dp;
|
|
|
|
while (dp && i++ < idx) {
|
|
ldpp = &dp->next;
|
|
dp = dp->next;
|
|
}
|
|
while (i++ <= idx) {
|
|
*ldpp = dp = create_dp(0, 0, 0, 0, 0);
|
|
ldpp = &((*ldpp)->next);
|
|
}
|
|
return dp;
|
|
}
|
|
|
|
/* return -1 to warn about potentially very long calculation */
|
|
int add_duration_to_nth_dp(struct diveplan *diveplan, int idx, int duration, gboolean is_rel)
|
|
{
|
|
struct divedatapoint *pdp, *dp = get_nth_dp(diveplan, idx);
|
|
if (idx > 0) {
|
|
pdp = get_nth_dp(diveplan, idx - 1);
|
|
if (duration && (is_rel || duration <= pdp->time))
|
|
duration += pdp->time;
|
|
}
|
|
dp->time = duration;
|
|
if (duration > 180 * 60)
|
|
return -1;
|
|
return 0;
|
|
}
|
|
|
|
/* this function is ONLY called from the dialog callback - so it
|
|
* marks this entry as 'entered'.
|
|
* Do NOT call from other parts of the planning code without changing
|
|
* that logic */
|
|
void add_depth_to_nth_dp(struct diveplan *diveplan, int idx, int depth)
|
|
{
|
|
struct divedatapoint *dp = get_nth_dp(diveplan, idx);
|
|
dp->depth = depth;
|
|
dp->entered = TRUE;
|
|
}
|
|
|
|
void add_gas_to_nth_dp(struct diveplan *diveplan, int idx, int o2, int he)
|
|
{
|
|
struct divedatapoint *dp = get_nth_dp(diveplan, idx);
|
|
dp->o2 = o2;
|
|
dp->he = he;
|
|
}
|
|
|
|
void add_po2_to_nth_dp(struct diveplan *diveplan, int idx, int po2)
|
|
{
|
|
struct divedatapoint *dp = get_nth_dp(diveplan, idx);
|
|
dp->po2 = po2;
|
|
}
|
|
|
|
void add_to_end_of_diveplan(struct diveplan *diveplan, struct divedatapoint *dp)
|
|
{
|
|
struct divedatapoint **lastdp = &diveplan->dp;
|
|
struct divedatapoint *ldp = *lastdp;
|
|
int lasttime = 0;
|
|
while (*lastdp) {
|
|
ldp = *lastdp;
|
|
if (ldp->time > lasttime)
|
|
lasttime = ldp->time;
|
|
lastdp = &(*lastdp)->next;
|
|
}
|
|
*lastdp = dp;
|
|
if (ldp)
|
|
dp->time += lasttime;
|
|
}
|
|
|
|
struct divedatapoint * plan_add_segment(struct diveplan *diveplan, int duration, int depth, int o2, int he, int po2)
|
|
{
|
|
struct divedatapoint *dp = create_dp(duration, depth, o2, he, po2);
|
|
add_to_end_of_diveplan(diveplan, dp);
|
|
return(dp);
|
|
}
|
|
|
|
struct gaschanges {
|
|
int depth;
|
|
int gasidx;
|
|
};
|
|
|
|
static struct gaschanges *analyze_gaslist(struct diveplan *diveplan, struct dive *dive, int *gaschangenr)
|
|
{
|
|
int nr = 0;
|
|
struct gaschanges *gaschanges = NULL;
|
|
struct divedatapoint *dp = diveplan->dp;
|
|
|
|
while (dp) {
|
|
if (dp->time == 0) {
|
|
int i = 0, j = 0;
|
|
nr++;
|
|
gaschanges = realloc(gaschanges, nr * sizeof(struct gaschanges));
|
|
while (i < nr - 1) {
|
|
if (dp->depth < gaschanges[i].depth) {
|
|
memmove(gaschanges + i + 1, gaschanges + i, (nr - i - 1) * sizeof(struct gaschanges));
|
|
break;
|
|
}
|
|
i++;
|
|
}
|
|
gaschanges[i].depth = dp->depth;
|
|
do {
|
|
if (dive->cylinder[j].gasmix.o2.permille == dp->o2 &&
|
|
dive->cylinder[j].gasmix.he.permille == dp->he) {
|
|
gaschanges[i].gasidx = j;
|
|
break;
|
|
}
|
|
j++;
|
|
} while (j < MAX_CYLINDERS);
|
|
}
|
|
dp = dp->next;
|
|
}
|
|
*gaschangenr = nr;
|
|
#if DEBUG_PLAN & 16
|
|
for (nr = 0; nr < *gaschangenr; nr++)
|
|
printf("gaschange nr %d: @ %5.2lfm gasidx %d (%d/%d)\n", nr, gaschanges[nr].depth / 1000.0,
|
|
gaschanges[nr].gasidx, (dive->cylinder[gaschanges[nr].gasidx].gasmix.o2.permille + 5) / 10,
|
|
(dive->cylinder[gaschanges[nr].gasidx].gasmix.he.permille + 5) / 10);
|
|
#endif
|
|
return gaschanges;
|
|
}
|
|
|
|
/* sort all the stops into one ordered list */
|
|
static int *sort_stops(int *dstops, int dnr, struct gaschanges *gstops, int gnr)
|
|
{
|
|
int i, gi, di;
|
|
int total = dnr + gnr;
|
|
int *stoplevels = malloc(total * sizeof(int));
|
|
|
|
/* no gaschanges */
|
|
if (gnr == 0) {
|
|
memcpy(stoplevels, dstops, dnr * sizeof(int));
|
|
return stoplevels;
|
|
}
|
|
i = total - 1;
|
|
gi = gnr - 1;
|
|
di = dnr - 1;
|
|
while (i >= 0) {
|
|
if (dstops[di] > gstops[gi].depth) {
|
|
stoplevels[i] = dstops[di];
|
|
di--;
|
|
} else if (dstops[di] == gstops[gi].depth) {
|
|
stoplevels[i] = dstops[di];
|
|
di--;
|
|
gi--;
|
|
} else {
|
|
stoplevels[i] = gstops[gi].depth;
|
|
gi--;
|
|
}
|
|
i--;
|
|
if (di < 0) {
|
|
while (gi >= 0)
|
|
stoplevels[i--] = gstops[gi--].depth;
|
|
break;
|
|
}
|
|
if (gi < 0) {
|
|
while (di >= 0)
|
|
stoplevels[i--] = dstops[di--];
|
|
break;
|
|
}
|
|
}
|
|
while (i >= 0)
|
|
stoplevels[i--] = 0;
|
|
|
|
#if DEBUG_PLAN & 16
|
|
int k;
|
|
for (k = gnr + dnr -1; k >= 0; k--) {
|
|
printf("stoplevel[%d]: %5.2lfm\n", k, stoplevels[k]/1000.0);
|
|
if (stoplevels[k] == 0)
|
|
break;
|
|
}
|
|
#endif
|
|
return stoplevels;
|
|
}
|
|
|
|
#if USE_GTK_UI
|
|
static void add_plan_to_notes(struct diveplan *diveplan, struct dive *dive)
|
|
{
|
|
char buffer[20000];
|
|
int consumption[MAX_CYLINDERS] = { 0, };
|
|
int len, gasidx, lastdepth = 0, lasttime = 0;
|
|
struct divedatapoint *dp = diveplan->dp;
|
|
int o2, he;
|
|
|
|
if (!dp)
|
|
return;
|
|
|
|
snprintf(buffer, sizeof(buffer), _("%s\nSubsurface dive plan\nbased on GFlow = %.0f and GFhigh = %.0f\n\n"),
|
|
disclaimer, plangflow * 100, plangfhigh * 100);
|
|
/* we start with gas 0, then check if that was changed */
|
|
o2 = dive->cylinder[0].gasmix.o2.permille;
|
|
he = dive->cylinder[0].gasmix.he.permille;
|
|
do {
|
|
const char *depth_unit;
|
|
char gas[64];
|
|
double depthvalue;
|
|
int decimals;
|
|
double used;
|
|
int newo2 = o2, newhe = he;
|
|
struct divedatapoint *nextdp;
|
|
|
|
if (dp->time == 0)
|
|
continue;
|
|
depthvalue = get_depth_units(dp->depth, &decimals, &depth_unit);
|
|
/* do we change gas after this segment? We need to look at the gas
|
|
* for the next segment (that isn't just a record of available gas !!)
|
|
* to find out */
|
|
nextdp = dp->next;
|
|
while (nextdp && nextdp->time == 0)
|
|
nextdp = nextdp->next;
|
|
if (nextdp) {
|
|
newo2 = nextdp->o2;
|
|
newhe = nextdp->he;
|
|
if (newhe == 0 && newo2 == 0) {
|
|
/* same as last segment */
|
|
newo2 = o2;
|
|
newhe = he;
|
|
}
|
|
}
|
|
/* do we want to skip this leg as it is devoid of anything useful? */
|
|
if (!dp->entered && o2 == newo2 && he == newhe && nextdp && dp->depth != lastdepth && nextdp->depth != dp->depth)
|
|
continue;
|
|
get_gas_string(o2, he, gas, sizeof(gas));
|
|
gasidx = get_gasidx(dive, o2, he);
|
|
len = strlen(buffer);
|
|
if (dp->depth != lastdepth) {
|
|
used = diveplan->bottomsac / 1000.0 * depth_to_mbar((dp->depth + lastdepth) / 2, dive) *
|
|
(dp->time - lasttime) / 60;
|
|
snprintf(buffer + len, sizeof(buffer) - len, _("Transition to %.*f %s in %d:%02d min - runtime %d:%02u on %s\n"),
|
|
decimals, depthvalue, depth_unit,
|
|
FRACTION(dp->time - lasttime, 60),
|
|
FRACTION(dp->time, 60),
|
|
gas);
|
|
} else {
|
|
/* we use deco SAC rate during the calculated deco stops, bottom SAC rate everywhere else */
|
|
int sac = dp->entered ? diveplan->bottomsac : diveplan->decosac;
|
|
used = sac / 1000.0 * depth_to_mbar(dp->depth, dive) * (dp->time - lasttime) / 60;
|
|
snprintf(buffer + len, sizeof(buffer) - len, _("Stay at %.*f %s for %d:%02d min - runtime %d:%02u on %s\n"),
|
|
decimals, depthvalue, depth_unit,
|
|
FRACTION(dp->time - lasttime, 60),
|
|
FRACTION(dp->time, 60),
|
|
gas);
|
|
}
|
|
if (gasidx != -1)
|
|
consumption[gasidx] += used;
|
|
get_gas_string(newo2, newhe, gas, sizeof(gas));
|
|
if (o2 != newo2 || he != newhe) {
|
|
len = strlen(buffer);
|
|
snprintf(buffer + len, sizeof(buffer) - len, _("Switch gas to %s\n"), gas);
|
|
}
|
|
o2 = newo2;
|
|
he = newhe;
|
|
lasttime = dp->time;
|
|
lastdepth = dp->depth;
|
|
} while ((dp = dp->next) != NULL);
|
|
len = strlen(buffer);
|
|
snprintf(buffer + len, sizeof(buffer) - len, _("Gas consumption:\n"));
|
|
for (gasidx = 0; gasidx < MAX_CYLINDERS; gasidx++) {
|
|
double volume;
|
|
const char *unit;
|
|
char gas[64];
|
|
if (consumption[gasidx] == 0)
|
|
continue;
|
|
len = strlen(buffer);
|
|
volume = get_volume_units(consumption[gasidx], NULL, &unit);
|
|
get_gas_string(dive->cylinder[gasidx].gasmix.o2.permille,
|
|
dive->cylinder[gasidx].gasmix.he.permille, gas, sizeof(gas));
|
|
snprintf(buffer + len, sizeof(buffer) - len, _("%.0f%s of %s\n"), volume, unit, gas);
|
|
}
|
|
dive->notes = strdup(buffer);
|
|
}
|
|
#endif
|
|
|
|
void plan(struct diveplan *diveplan, char **cached_datap, struct dive **divep, bool add_deco, char **error_string_p)
|
|
{
|
|
struct dive *dive;
|
|
struct sample *sample;
|
|
int wait_time, o2, he, po2;
|
|
int ceiling, depth, transitiontime;
|
|
int stopidx, gi;
|
|
double tissue_tolerance;
|
|
struct gaschanges *gaschanges;
|
|
int gaschangenr;
|
|
int *stoplevels;
|
|
|
|
set_gf(diveplan->gflow, diveplan->gfhigh);
|
|
if (!diveplan->surface_pressure)
|
|
diveplan->surface_pressure = SURFACE_PRESSURE;
|
|
if (*divep)
|
|
delete_single_dive(dive_table.nr - 1);
|
|
*divep = dive = create_dive_from_plan(diveplan, error_string_p);
|
|
if (!dive)
|
|
return;
|
|
record_dive(dive);
|
|
|
|
sample = &dive->dc.sample[dive->dc.samples - 1];
|
|
/* we start with gas 0, then check if that was changed */
|
|
o2 = dive->cylinder[0].gasmix.o2.permille;
|
|
he = dive->cylinder[0].gasmix.he.permille;
|
|
get_gas_from_events(&dive->dc, sample->time.seconds, &o2, &he);
|
|
po2 = dive->dc.sample[dive->dc.samples - 1].po2;
|
|
depth = dive->dc.sample[dive->dc.samples - 1].depth.mm;
|
|
|
|
/* if all we wanted was the dive just get us back to the surface */
|
|
if (!add_deco) {
|
|
transitiontime = depth / 150; /* this still needs to be made configurable */
|
|
plan_add_segment(diveplan, transitiontime, 0, o2, he, po2);
|
|
/* re-create the dive */
|
|
delete_single_dive(dive_table.nr - 1);
|
|
*divep = dive = create_dive_from_plan(diveplan, error_string_p);
|
|
if (dive)
|
|
record_dive(dive);
|
|
return;
|
|
}
|
|
|
|
tissue_tolerance = tissue_at_end(dive, cached_datap, error_string_p);
|
|
ceiling = deco_allowed_depth(tissue_tolerance, diveplan->surface_pressure / 1000.0, dive, 1);
|
|
#if DEBUG_PLAN & 4
|
|
printf("gas %d/%d\n", o2, he);
|
|
printf("depth %5.2lfm ceiling %5.2lfm\n", depth / 1000.0, ceiling / 1000.0);
|
|
#endif
|
|
if (depth < ceiling) /* that's not good... */
|
|
depth = ceiling;
|
|
for (stopidx = 0; stopidx < sizeof(decostoplevels) / sizeof(int); stopidx++)
|
|
if (decostoplevels[stopidx] >= depth)
|
|
break;
|
|
stopidx--;
|
|
|
|
/* so now we know the first decostop level above us
|
|
* NOTE, this could be the surface or a long list of potential stops
|
|
* we do NOT start only at the ceiling, as the ceiling may come down
|
|
* further during the ascent.
|
|
* Next we need to figure out if there are better gases available
|
|
* and at which depths we are supposed to switch to them */
|
|
gaschanges = analyze_gaslist(diveplan, dive, &gaschangenr);
|
|
stoplevels = sort_stops(decostoplevels, stopidx + 1, gaschanges, gaschangenr);
|
|
|
|
gi = gaschangenr - 1;
|
|
stopidx += gaschangenr;
|
|
if (depth > stoplevels[stopidx]) {
|
|
transitiontime = (depth - stoplevels[stopidx]) / 150;
|
|
#if DEBUG_PLAN & 2
|
|
printf("transitiontime %d:%02d to depth %5.2lfm\n", FRACTION(transitiontime, 60), stoplevels[stopidx] / 1000.0);
|
|
#endif
|
|
plan_add_segment(diveplan, transitiontime, stoplevels[stopidx], o2, he, po2);
|
|
/* re-create the dive */
|
|
delete_single_dive(dive_table.nr - 1);
|
|
*divep = dive = create_dive_from_plan(diveplan, error_string_p);
|
|
if (!dive)
|
|
goto error_exit;
|
|
record_dive(dive);
|
|
}
|
|
while (stopidx > 0) { /* this indicates that we aren't surfacing directly */
|
|
/* if we are in a double-step, eg, when 3m/10ft stop is disabled,
|
|
* just skip the first stop at that depth */
|
|
if (stoplevels[stopidx] == stoplevels[stopidx - 1]) {
|
|
stopidx--;
|
|
continue;
|
|
}
|
|
if (gi >= 0 && stoplevels[stopidx] == gaschanges[gi].depth) {
|
|
o2 = dive->cylinder[gaschanges[gi].gasidx].gasmix.o2.permille;
|
|
he = dive->cylinder[gaschanges[gi].gasidx].gasmix.he.permille;
|
|
#if DEBUG_PLAN & 16
|
|
printf("switch to gas %d (%d/%d) @ %5.2lfm\n", gaschanges[gi].gasidx,
|
|
(o2 + 5) / 10, (he + 5) / 10, gaschanges[gi].depth / 1000.0);
|
|
#endif
|
|
gi--;
|
|
}
|
|
wait_time = time_at_last_depth(dive, o2, he, stoplevels[stopidx - 1], cached_datap, error_string_p);
|
|
/* typically deco plans are done in one minute increments; we may want to
|
|
* make this configurable at some point */
|
|
wait_time = ((wait_time + 59) / 60) * 60;
|
|
#if DEBUG_PLAN & 2
|
|
tissue_tolerance = tissue_at_end(dive, cached_datap, error_string_p);
|
|
ceiling = deco_allowed_depth(tissue_tolerance, diveplan->surface_pressure / 1000.0, dive, 1);
|
|
printf("waittime %d:%02d at depth %5.2lfm; ceiling %5.2lfm\n", FRACTION(wait_time, 60),
|
|
stoplevels[stopidx] / 1000.0, ceiling / 1000.0);
|
|
#endif
|
|
if (wait_time)
|
|
plan_add_segment(diveplan, wait_time, stoplevels[stopidx], o2, he, po2);
|
|
transitiontime = (stoplevels[stopidx] - stoplevels[stopidx - 1]) / 150;
|
|
#if DEBUG_PLAN & 2
|
|
printf("transitiontime %d:%02d to depth %5.2lfm\n", FRACTION(transitiontime, 60), stoplevels[stopidx - 1] / 1000.0);
|
|
#endif
|
|
plan_add_segment(diveplan, transitiontime, stoplevels[stopidx - 1], o2, he, po2);
|
|
/* re-create the dive */
|
|
delete_single_dive(dive_table.nr - 1);
|
|
*divep = dive = create_dive_from_plan(diveplan, error_string_p);
|
|
if (!dive)
|
|
goto error_exit;
|
|
record_dive(dive);
|
|
stopidx--;
|
|
}
|
|
#if USE_GTK_UI
|
|
add_plan_to_notes(diveplan, dive);
|
|
/* now make the dive visible in the dive list */
|
|
report_dives(FALSE, FALSE);
|
|
show_and_select_dive(dive);
|
|
#endif
|
|
error_exit:
|
|
free(stoplevels);
|
|
free(gaschanges);
|
|
}
|
|
|
|
/*
|
|
* Get a value in tenths (so "10.2" == 102, "9" = 90)
|
|
*
|
|
* Return negative for errors.
|
|
*/
|
|
static int get_tenths(const char *begin, const char **endp)
|
|
{
|
|
char *end;
|
|
int value = strtol(begin, &end, 10);
|
|
|
|
if (begin == end)
|
|
return -1;
|
|
value *= 10;
|
|
|
|
/* Fraction? We only look at the first digit */
|
|
if (*end == '.') {
|
|
end++;
|
|
if (!isdigit(*end))
|
|
return -1;
|
|
value += *end - '0';
|
|
do {
|
|
end++;
|
|
} while (isdigit(*end));
|
|
}
|
|
*endp = end;
|
|
return value;
|
|
}
|
|
|
|
static int get_thousandths(const char *begin, const char **endp)
|
|
{
|
|
char *end;
|
|
int value = strtol(begin, &end, 10);
|
|
|
|
if (begin == end)
|
|
return -1;
|
|
value *= 1000;
|
|
|
|
/* now deal with up to three more digits after decimal point */
|
|
if (*end == '.') {
|
|
int factor = 100;
|
|
do {
|
|
++end;
|
|
if (!isdigit(*end))
|
|
break;
|
|
value += (*end - '0') * factor;
|
|
factor /= 10;
|
|
} while (factor);
|
|
}
|
|
*endp = end;
|
|
return value;
|
|
}
|
|
|
|
static int get_permille(const char *begin, const char **end)
|
|
{
|
|
int value = get_tenths(begin, end);
|
|
if (value >= 0) {
|
|
/* Allow a percentage sign */
|
|
if (**end == '%')
|
|
++*end;
|
|
}
|
|
return value;
|
|
}
|
|
|
|
int validate_gas(const char *text, int *o2_p, int *he_p)
|
|
{
|
|
int o2, he;
|
|
|
|
if (!text)
|
|
return 0;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
|
|
if (!*text)
|
|
return 0;
|
|
|
|
if (!strcasecmp(text, _("air"))) {
|
|
o2 = O2_IN_AIR; he = 0; text += strlen(_("air"));
|
|
} else if (!strncasecmp(text, _("ean"), 3)) {
|
|
o2 = get_permille(text+3, &text); he = 0;
|
|
} else {
|
|
o2 = get_permille(text, &text); he = 0;
|
|
if (*text == '/')
|
|
he = get_permille(text+1, &text);
|
|
}
|
|
|
|
/* We don't want any extra crud */
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
if (*text)
|
|
return 0;
|
|
|
|
/* Validate the gas mix */
|
|
if (*text || o2 < 1 || o2 > 1000 || he < 0 || o2+he > 1000)
|
|
return 0;
|
|
|
|
/* Let it rip */
|
|
*o2_p = o2;
|
|
*he_p = he;
|
|
return 1;
|
|
}
|
|
|
|
int validate_time(const char *text, int *sec_p, int *rel_p)
|
|
{
|
|
int min, sec, rel;
|
|
char *end;
|
|
|
|
if (!text)
|
|
return 0;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
|
|
rel = 1;
|
|
if (*text == '+') {
|
|
rel = 1;
|
|
text++;
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
} else if (*text == '@') {
|
|
rel = 0;
|
|
text++;
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
}
|
|
|
|
min = strtol(text, &end, 10);
|
|
if (text == end)
|
|
return 0;
|
|
|
|
if (min < 0 || min > 1000)
|
|
return 0;
|
|
|
|
/* Ok, minutes look ok */
|
|
text = end;
|
|
sec = 0;
|
|
if (*text == ':') {
|
|
text++;
|
|
sec = strtol(text, &end, 10);
|
|
if (end != text+2)
|
|
return 0;
|
|
if (sec < 0)
|
|
return 0;
|
|
text = end;
|
|
if (*text == ':') {
|
|
if (sec >= 60)
|
|
return 0;
|
|
min = min*60 + sec;
|
|
text++;
|
|
sec = strtol(text, &end, 10);
|
|
if (end != text+2)
|
|
return 0;
|
|
if (sec < 0)
|
|
return 0;
|
|
text = end;
|
|
}
|
|
}
|
|
|
|
/* Maybe we should accept 'min' at the end? */
|
|
if (g_ascii_isspace(*text))
|
|
text++;
|
|
if (*text)
|
|
return 0;
|
|
|
|
*sec_p = min*60 + sec;
|
|
*rel_p = rel;
|
|
return 1;
|
|
}
|
|
|
|
int validate_depth(const char *text, int *mm_p)
|
|
{
|
|
int depth, imperial;
|
|
|
|
if (!text)
|
|
return 0;
|
|
|
|
depth = get_tenths(text, &text);
|
|
if (depth < 0)
|
|
return 0;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
|
|
imperial = get_units()->length == FEET;
|
|
if (*text == 'm') {
|
|
imperial = 0;
|
|
text++;
|
|
} else if (!strcasecmp(text, _("ft"))) {
|
|
imperial = 1;
|
|
text += 2;
|
|
}
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
if (*text)
|
|
return 0;
|
|
|
|
if (imperial) {
|
|
depth = feet_to_mm(depth / 10.0);
|
|
} else {
|
|
depth *= 100;
|
|
}
|
|
*mm_p = depth;
|
|
/* we don't support extreme depths */
|
|
if (depth > 400000)
|
|
return 0;
|
|
return 1;
|
|
}
|
|
|
|
int validate_po2(const char *text, int *mbar_po2)
|
|
{
|
|
int po2;
|
|
|
|
if (!text)
|
|
return 0;
|
|
|
|
po2 = get_tenths(text, &text);
|
|
if (po2 < 0)
|
|
return 0;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
if (*text)
|
|
return 0;
|
|
|
|
*mbar_po2 = po2 * 100;
|
|
return 1;
|
|
}
|
|
|
|
int validate_volume(const char *text, int *sac)
|
|
{
|
|
int volume, imperial;
|
|
|
|
if (!text)
|
|
return 0;
|
|
|
|
volume = get_thousandths(text, &text);
|
|
if (volume < 0)
|
|
return 0;
|
|
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
|
|
imperial = get_units()->volume == CUFT;
|
|
if (*text == 'l') {
|
|
imperial = 0;
|
|
text++;
|
|
} else if (!strncasecmp(text, _("cuft"), 4)) {
|
|
imperial = 1;
|
|
text += 4;
|
|
}
|
|
while (g_ascii_isspace(*text) || *text == '/')
|
|
text++;
|
|
if (!strncasecmp(text, _("min"), 3))
|
|
text += 3;
|
|
while (g_ascii_isspace(*text))
|
|
text++;
|
|
if (*text)
|
|
return 0;
|
|
if (imperial)
|
|
volume = cuft_to_l(volume) + 0.5; /* correct for mcuft -> ml */
|
|
*sac = volume;
|
|
return 1;
|
|
}
|
|
|
|
#if USE_GTK_UI
|
|
struct diveplan diveplan = {};
|
|
char *cache_data = NULL;
|
|
struct dive *planned_dive = NULL;
|
|
|
|
/* make a copy of the diveplan so far and display the corresponding dive */
|
|
void show_planned_dive(char **error_string_p)
|
|
{
|
|
struct diveplan tempplan;
|
|
struct divedatapoint *dp, **dpp;
|
|
|
|
memcpy(&tempplan, &diveplan, sizeof(struct diveplan));
|
|
dpp = &tempplan.dp;
|
|
dp = diveplan.dp;
|
|
while (dp && *dpp) {
|
|
*dpp = malloc(sizeof(struct divedatapoint));
|
|
memcpy(*dpp, dp, sizeof(struct divedatapoint));
|
|
dp = dp->next;
|
|
dpp = &(*dpp)->next;
|
|
}
|
|
#if DEBUG_PLAN & 1
|
|
printf("in show_planned_dive:\n");
|
|
dump_plan(&tempplan);
|
|
#endif
|
|
plan(&tempplan, &cache_data, &planned_dive, error_string_p);
|
|
free_dps(tempplan.dp);
|
|
}
|
|
|
|
/* Subsurface follows the lead of most divecomputers to use times
|
|
* without timezone - so all times are implicitly assumed to be
|
|
* local time of the dive location; so in order to give the current
|
|
* time in that way we actually need to add the timezone offset */
|
|
timestamp_t current_time_notz(void)
|
|
{
|
|
time_t now = time(NULL);
|
|
struct tm *local = localtime(&now);
|
|
return utc_mktime(local);
|
|
}
|
|
#endif
|