X-Git-Url: http://git.tdb.fi/?a=blobdiff_plain;f=profile.c;h=20317b93717704e922000528cf1b5060df1d7bd1;hb=682135838ff313594c7f67fabd9be8f88a33883b;hp=b2679c5d57b30f3eee228d1cc901a24604f0f4c5;hpb=89fe2c723f10649d2d98a88b385bd4733e5acbc4;p=ext%2Fsubsurface.git diff --git a/profile.c b/profile.c index b2679c5..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,7 +14,9 @@ int selected_dive = 0; -/* Plot info with smoothing and one-, two- and three-minute minimums and maximums */ +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 { int nr; int maxtime; @@ -23,11 +29,23 @@ struct plot_info { /* Depth info */ int val; int smoothed; + velocity_t velocity; struct plot_data *min[3]; struct plot_data *max[3]; int avg[3]; } entry[]; }; + +/* convert velocity to colors */ +typedef struct { double r, g, b; } rgb_t; +static const rgb_t rgb[] = { + [STABLE] = {0.0, 0.4, 0.0}, + [SLOW] = {0.4, 0.8, 0.0}, + [MODERATE] = {0.8, 0.8, 0.0}, + [FAST] = {0.8, 0.5, 0.0}, + [CRAZY] = {1.0, 0.0, 0.0}, +}; + #define plot_info_size(nr) (sizeof(struct plot_info) + (nr)*sizeof(struct plot_data)) /* Scale to 0,0 -> maxx,maxy */ @@ -58,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); } @@ -70,6 +88,8 @@ static void set_source_rgb(struct graphics_context *gc, double r, double g, doub * current dive. However, we don't scale past less than * 30 minutes or 90 ft, just so that small dives show * up as such. + * we also need to add 180 seconds at the end so the min/max + * plots correctly */ static int get_maxtime(struct plot_info *pi) { @@ -234,7 +254,9 @@ static void plot_depth_profile(struct graphics_context *gc, struct plot_info *pi { int i; cairo_t *cr = gc->cr; - int begins, sec, depth; + int ends, sec, depth; + int *secs; + int *depths; struct plot_data *entry; int maxtime, maxdepth, marker; @@ -278,24 +300,47 @@ static void plot_depth_profile(struct graphics_context *gc, struct plot_info *pi entry = pi->entry; set_source_rgba(gc, 1, 0.2, 0.2, 0.80); - begins = entry->sec; - move_to(gc, entry->sec, entry->val); + secs = (int *) malloc(sizeof(int) * pi->nr); + depths = (int *) malloc(sizeof(int) * pi->nr); + secs[0] = entry->sec; + depths[0] = entry->val; for (i = 1; i < pi->nr; i++) { entry++; sec = entry->sec; - if (sec <= maxtime) { + if (sec <= maxtime || entry->val > 0) { + /* we want to draw the segments in different colors + * representing the vertical velocity, so we need to + * chop this into short segments */ + rgb_t color = rgb[entry->velocity]; depth = entry->val; + set_source_rgb(gc, color.r, color.g, color.b); + move_to(gc, secs[i-1], depths[i-1]); line_to(gc, sec, depth); + cairo_stroke(cr); + ends = i; } + secs[i] = sec; + depths[i] = depth; } + move_to(gc, secs[ends], depths[ends]); gc->topy = 0; gc->bottomy = 1.0; - line_to(gc, MIN(sec,maxtime), 0); - line_to(gc, begins, 0); + line_to(gc, secs[ends], 0); + line_to(gc, secs[0], 0); cairo_close_path(cr); - set_source_rgba(gc, 1, 0.2, 0.2, 0.20); - cairo_fill_preserve(cr); set_source_rgba(gc, 1, 0.2, 0.2, 0.80); cairo_stroke(cr); + /* now do it again for the neat fill */ + gc->topy = 0; gc->bottomy = maxdepth; + set_source_rgba(gc, 1, 0.2, 0.2, 0.20); + move_to(gc, secs[0], depths[0]); + for (i = 1; i <= ends; i++) { + line_to(gc, secs[i],depths[i]); + } + gc->topy = 0; gc->bottomy = 1.0; + line_to(gc, secs[ends], 0); + line_to(gc, secs[0], 0); + cairo_close_path(gc->cr); + cairo_fill(gc->cr); } static int setup_temperature_limits(struct graphics_context *gc, struct plot_info *pi) @@ -401,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; } @@ -428,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(""); - 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}; @@ -584,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; @@ -611,12 +606,28 @@ static struct plot_info *analyze_plot_info(struct plot_info *pi) } /* Smoothing function: 5-point triangular smooth */ - for (i = 2; i < nr-2; i++) { + for (i = 2; i < nr-1; i++) { struct plot_data *entry = pi->entry+i; int val; - val = entry[-2].val + 2*entry[-1].val + 3*entry[0].val + 2*entry[1].val + entry[2].val; - entry->smoothed = (val+4) / 9; + if (i < nr-2) { + val = entry[-2].val + 2*entry[-1].val + 3*entry[0].val + 2*entry[1].val + entry[2].val; + 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) { + 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; } /* One-, two- and three-minute minmax data */ @@ -721,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); @@ -736,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; -}