X-Git-Url: http://git.tdb.fi/?a=blobdiff_plain;f=profile.c;h=20317b93717704e922000528cf1b5060df1d7bd1;hb=682135838ff313594c7f67fabd9be8f88a33883b;hp=0add26d0844d98fb824f669c4c0593d360d4e8b7;hpb=1937df188a2928478a9012a382ba1360d5e309b7;p=ext%2Fsubsurface.git diff --git a/profile.c b/profile.c index 0add26d..20317b9 100644 --- a/profile.c +++ b/profile.c @@ -1,3 +1,7 @@ +/* profile.c */ +/* creates all the necessary data for drawing the dive profile + * uses cairo to draw it + */ #include #include #include @@ -10,6 +14,7 @@ int selected_dive = 0; +typedef enum { STABLE, SLOW, MODERATE, FAST, CRAZY } velocity_t; /* Plot info with smoothing, velocity indication * and one-, two- and three-minute minimums and maximums */ struct plot_info { @@ -24,7 +29,7 @@ struct plot_info { /* Depth info */ int val; int smoothed; - enum { STABLE, SLOW, MODERATE, FAST, CRAZY } velocity; + velocity_t velocity; struct plot_data *min[3]; struct plot_data *max[3]; int avg[3]; @@ -71,7 +76,7 @@ static void set_source_rgba(struct graphics_context *gc, double r, double g, dou cairo_set_source_rgba(gc->cr, r, g, b, a); } -static void set_source_rgb(struct graphics_context *gc, double r, double g, double b) +void set_source_rgb(struct graphics_context *gc, double r, double g, double b) { set_source_rgba(gc, r, g, b, 1); } @@ -441,7 +446,7 @@ static int get_cylinder_pressure_range(struct graphics_context *gc, struct plot_ gc->leftx = 0; gc->rightx = get_maxtime(pi); - gc->topy = 0; gc->bottomy = pi->maxpressure * 1.5; + gc->bottomy = 0; gc->topy = pi->maxpressure * 1.5; return pi->maxpressure != 0; } @@ -468,83 +473,6 @@ static void plot_cylinder_pressure(struct graphics_context *gc, struct plot_info cairo_stroke(gc->cr); } -/* - * Return air usage (in liters). - */ -static double calculate_airuse(struct dive *dive) -{ - double airuse = 0; - int i; - - for (i = 0; i < MAX_CYLINDERS; i++) { - cylinder_t *cyl = dive->cylinder + i; - int size = cyl->type.size.mliter; - double kilo_atm; - - if (!size) - continue; - - kilo_atm = (cyl->start.mbar - cyl->end.mbar) / 1013250.0; - - /* Liters of air at 1 atm == milliliters at 1k atm*/ - airuse += kilo_atm * size; - } - return airuse; -} - -static void plot_info(struct dive *dive, struct graphics_context *gc) -{ - text_render_options_t tro = {10, 0.2, 1.0, 0.2, RIGHT, BOTTOM}; - const double liters_per_cuft = 28.317; - const char *unit, *format, *desc; - double airuse; - char buffer1[80]; - char buffer2[80]; - int len; - - airuse = calculate_airuse(dive); - if (!airuse) { - update_air_info(" \n "); - return; - } - switch (output_units.volume) { - case LITER: - unit = "l"; - format = "vol: %4.0f %s"; - break; - case CUFT: - unit = "cuft"; - format = "vol: %4.2f %s"; - airuse /= liters_per_cuft; - break; - } - tro.vpos = -1.0; - plot_text(gc, &tro, 0.98, 0.98, format, airuse, unit); - len = snprintf(buffer1, sizeof(buffer1), format, airuse, unit); - tro.vpos = -2.2; - if (dive->duration.seconds) { - double pressure = 1 + (dive->meandepth.mm / 10000.0); - double sac = airuse / pressure * 60 / dive->duration.seconds; - plot_text(gc, &tro, 0.98, 0.98, "SAC: %4.2f %s/min", sac, unit); - snprintf(buffer1+len, sizeof(buffer1)-len, - "\nSAC: %4.2f %s/min", sac, unit); - } - len = 0; - tro.vpos = -3.4; - desc = dive->cylinder[0].type.description; - if (desc || dive->cylinder[0].gasmix.o2.permille) { - int o2 = dive->cylinder[0].gasmix.o2.permille / 10; - if (!desc) - desc = ""; - if (!o2) - o2 = 21; - plot_text(gc, &tro, 0.98, 0.98, "%s (%d%%)", desc, o2); - len = snprintf(buffer2, sizeof(buffer2), "%s (%d%%): used ", desc, o2); - } - snprintf(buffer2+len, sizeof(buffer2)-len, buffer1); - update_air_info(buffer2); -} - static int mbar_to_PSI(int mbar) { pressure_t p = {mbar}; @@ -624,6 +552,33 @@ static void analyze_plot_info_minmax(struct plot_data *entry, struct plot_data * analyze_plot_info_minmax_minute(entry, first, last, 2); } +static velocity_t velocity(int speed) +{ + velocity_t v; + + if (speed < -304) /* ascent faster than -60ft/min */ + v = CRAZY; + else if (speed < -152) /* above -30ft/min */ + v = FAST; + else if (speed < -76) /* -15ft/min */ + v = MODERATE; + else if (speed < -25) /* -5ft/min */ + v = SLOW; + else if (speed < 25) /* very hard to find data, but it appears that the recommendations + for descent are usually about 2x ascent rate; still, we want + stable to mean stable */ + v = STABLE; + else if (speed < 152) /* between 5 and 30ft/min is considered slow */ + v = SLOW; + else if (speed < 304) /* up to 60ft/min is moderate */ + v = MODERATE; + else if (speed < 507) /* up to 100ft/min is fast */ + v = FAST; + else /* more than that is just crazy - you'll blow your ears out */ + v = CRAZY; + + return v; +} static struct plot_info *analyze_plot_info(struct plot_info *pi) { int i; @@ -660,28 +615,17 @@ static struct plot_info *analyze_plot_info(struct plot_info *pi) entry->smoothed = (val+4) / 9; } /* vertical velocity in mm/sec */ + /* Linus wants to smooth this - let's at least look at the samples that aren't FAST or CRAZY */ if (entry[0].sec - entry[-1].sec) { - val = (entry[0].val - entry[-1].val) / (entry[0].sec - entry[-1].sec); - if (val < -304) /* ascent faster than -60ft/min */ - entry->velocity = CRAZY; - else if (val < -152) /* above -30ft/min */ - entry->velocity = FAST; - else if (val < -76) /* -15ft/min */ - entry->velocity = MODERATE; - else if (val < -25) /* -5ft/min */ - entry->velocity = SLOW; - else if (val < 25) /* very hard to find data, but it appears that the recommendations - for descent are usually about 2x ascent rate; still, we want - stable to mean stable */ - entry->velocity = STABLE; - else if (val < 152) /* between 5 and 30ft/min is considered slow */ - entry->velocity = SLOW; - else if (val < 304) /* up to 60ft/min is moderate */ - entry->velocity = MODERATE; - else if (val < 507) /* up to 100ft/min is fast */ - entry->velocity = FAST; - else /* more than that is just crazy - you'll blow your ears out */ - entry->velocity = CRAZY; + entry->velocity = velocity((entry[0].val - entry[-1].val) / (entry[0].sec - entry[-1].sec)); + /* if our samples are short and we aren't too FAST*/ + if (entry[0].sec - entry[-1].sec < 30 && entry->velocity < FAST) { + int past = -2; + while (pi->entry <= entry-past && entry[0].sec - entry[past].sec < 30) + past--; + entry->velocity = velocity((entry[0].val - entry[past].val) / + (entry[0].sec - entry[past].sec)); + } } else entry->velocity = STABLE; } @@ -788,12 +732,10 @@ void plot(struct graphics_context *gc, int w, int h, struct dive *dive) plot_depth_text(gc, pi); plot_cylinder_pressure_text(gc, pi); - /* And info box in the lower right corner.. */ + /* Bounding box last */ gc->leftx = 0; gc->rightx = 1.0; gc->topy = 0; gc->bottomy = 1.0; - plot_info(dive, gc); - /* Bounding box last */ set_source_rgb(gc, 1, 1, 1); move_to(gc, 0, 0); line_to(gc, 0, 1); @@ -803,35 +745,3 @@ void plot(struct graphics_context *gc, int w, int h, struct dive *dive) cairo_stroke(gc->cr); } - -static gboolean expose_event(GtkWidget *widget, GdkEventExpose *event, gpointer data) -{ - struct dive *dive = current_dive; - struct graphics_context gc = { .printer = 0 }; - int w,h; - - w = widget->allocation.width; - h = widget->allocation.height; - - gc.cr = gdk_cairo_create(widget->window); - set_source_rgb(&gc, 0, 0, 0); - cairo_paint(gc.cr); - - if (dive) - plot(&gc, w, h, dive); - - cairo_destroy(gc.cr); - - return FALSE; -} - -GtkWidget *dive_profile_widget(void) -{ - GtkWidget *da; - - da = gtk_drawing_area_new(); - gtk_widget_set_size_request(da, 350, 250); - g_signal_connect(da, "expose_event", G_CALLBACK(expose_event), NULL); - - return da; -}