10 int selected_dive = 0;
13 * Cairo scaling really is horribly horribly mis-designed.
15 * Which is sad, because I really like Cairo otherwise. But
16 * the fact that the line width is scaled with the same scale
17 * as the coordinate system is a f*&%ing disaster. So we
18 * can't use it, and instead have this butt-ugly wrapper thing..
20 struct graphics_context {
23 double scalex, scaley;
26 /* Scale to 0,0 -> maxx,maxy */
27 #define SCALE(gc,x,y) (x)*gc->maxx/gc->scalex,(y)*gc->maxy/gc->scaley
29 static void move_to(struct graphics_context *gc, double x, double y)
31 cairo_move_to(gc->cr, SCALE(gc, x, y));
34 static void line_to(struct graphics_context *gc, double x, double y)
36 cairo_line_to(gc->cr, SCALE(gc, x, y));
39 #define ROUND_UP(x,y) ((((x)+(y)-1)/(y))*(y))
42 * When showing dive profiles, we scale things to the
43 * current dive. However, we don't scale past less than
44 * 30 minutes or 90 ft, just so that small dives show
47 static int round_seconds_up(int seconds)
49 return MAX(30*60, ROUND_UP(seconds, 60*10));
52 static int round_depth_up(depth_t depth)
54 unsigned mm = depth.mm;
56 return MAX(30000, ROUND_UP(mm+3000, 10000));
61 enum {CENTER,LEFT} allign;
62 } text_render_options_t;
64 static void plot_text(struct graphics_context *gc, text_render_options_t *tro,
65 double x, double y, const char *fmt, ...)
68 cairo_text_extents_t extents;
74 vsnprintf(buffer, sizeof(buffer), fmt, args);
77 cairo_text_extents(cr, buffer, &extents);
79 if (tro->allign == CENTER)
80 dx = -(extents.width/2 + extents.x_bearing);
81 dy = extents.height * 1.2;
84 cairo_rel_move_to(cr, dx, dy);
86 cairo_text_path(cr, buffer);
87 cairo_set_source_rgb(cr, 0, 0, 0);
91 cairo_rel_move_to(cr, dx, dy);
93 cairo_set_source_rgb(cr, tro->r, tro->g, tro->b);
94 cairo_show_text(cr, buffer);
98 * Find the next maximum point in a 10-minute window.
100 * We exit early if we hit "enough" of a depth reversal,
101 * which is roughly 10 feet.
103 static struct sample *next_minmax(struct sample *sample, struct sample *end, int minmax)
105 const int enough = 3000;
106 struct sample *result;
107 int timelimit, depthlimit;
112 timelimit = 24*60*60;
113 depthlimit = sample->depth.mm;
122 time = sample->time.seconds;
123 depth = sample->depth.mm;
124 if (time > timelimit)
128 if (depth <= depthlimit) {
129 if (depthlimit - depth > enough)
134 if (depth >= depthlimit) {
135 if (depth - depthlimit > enough)
143 /* Look up to ten minutes into the future */
144 timelimit = time + 600;
149 static void render_depth_sample(struct graphics_context *gc, struct sample *sample)
151 text_render_options_t tro = {1.0, 0.2, 0.2, CENTER};
152 int sec = sample->time.seconds;
153 depth_t depth = sample->depth;
157 switch (output_units.length) {
159 d = depth.mm / 1000.0;
167 plot_text(gc, &tro, sec, depth.mm, fmt, d);
170 static void plot_text_samples(struct graphics_context *gc, struct sample *a, struct sample *b)
172 struct sample *max, *min;
176 if (b->time.seconds - a->time.seconds < 3*60)
179 max = next_minmax(a, b, 1);
181 render_depth_sample(gc, max);
182 min = next_minmax(max, b, 0);
184 plot_text_samples(gc, min, b);
190 static void plot_depth_text(struct dive *dive, struct graphics_context *gc)
192 struct sample *sample, *end;
193 int maxtime, maxdepth;
195 /* Get plot scaling limits */
196 maxtime = round_seconds_up(dive->duration.seconds);
197 maxdepth = round_depth_up(dive->maxdepth);
199 gc->scalex = maxtime;
200 gc->scaley = maxdepth;
202 cairo_set_font_size(gc->cr, 14);
205 * We never take the last sample into account.
206 * It should be a surface event anyway, although
207 * there are buggy cases where it isn't..
209 sample = dive->sample;
210 end = dive->sample + dive->samples - 1;
212 plot_text_samples(gc, sample, end);
215 static void plot_depth_profile(struct dive *dive, struct graphics_context *gc)
217 cairo_t *cr = gc->cr;
218 int begins, sec, depth;
220 struct sample *sample;
221 int maxtime, maxdepth, marker;
223 samples = dive->samples;
227 cairo_set_line_width(gc->cr, 2);
229 /* Get plot scaling limits */
230 maxtime = round_seconds_up(dive->duration.seconds);
231 maxdepth = round_depth_up(dive->maxdepth);
233 /* Time markers: every 5 min */
234 gc->scalex = maxtime;
236 for (i = 5*60; i < maxtime; i += 5*60) {
241 /* Depth markers: every 30 ft or 10 m*/
243 gc->scaley = maxdepth;
244 switch (output_units.length) {
245 case METERS: marker = 10000; break;
246 case FEET: marker = 9144; break; /* 30 ft */
249 cairo_set_source_rgba(cr, 1, 1, 1, 0.5);
250 for (i = marker; i < maxdepth; i += marker) {
256 /* Show mean depth */
257 cairo_set_source_rgba(cr, 1, 0.2, 0.2, 0.40);
258 move_to(gc, 0, dive->meandepth.mm);
259 line_to(gc, 1, dive->meandepth.mm);
262 gc->scalex = maxtime;
264 sample = dive->sample;
265 cairo_set_source_rgba(cr, 1, 0.2, 0.2, 0.80);
266 begins = sample->time.seconds;
267 move_to(gc, sample->time.seconds, sample->depth.mm);
268 for (i = 1; i < dive->samples; i++) {
270 sec = sample->time.seconds;
271 if (sec <= maxtime) {
272 depth = sample->depth.mm;
273 line_to(gc, sec, depth);
277 line_to(gc, MIN(sec,maxtime), 0);
278 line_to(gc, begins, 0);
279 cairo_close_path(cr);
280 cairo_set_source_rgba(cr, 1, 0.2, 0.2, 0.20);
281 cairo_fill_preserve(cr);
282 cairo_set_source_rgba(cr, 1, 0.2, 0.2, 0.80);
286 /* gets both the actual start and end pressure as well as the scaling factors */
287 static int get_cylinder_pressure_range(struct dive *dive, struct graphics_context *gc,
288 pressure_t *startp, pressure_t *endp)
293 gc->scalex = round_seconds_up(dive->duration.seconds);
298 startp->mbar = endp->mbar = 0;
300 for (i = 0; i < dive->samples; i++) {
302 struct sample *sample = dive->sample + i;
304 /* FIXME! We only track cylinder 0 right now */
305 if (sample->cylinderindex)
307 mbar = sample->cylinderpressure.mbar;
321 gc->scaley = max * 1.5;
325 static void plot_cylinder_pressure(struct dive *dive, struct graphics_context *gc)
329 if (!get_cylinder_pressure_range(dive, gc, NULL, NULL))
332 cairo_set_source_rgba(gc->cr, 0.2, 1.0, 0.2, 0.80);
334 move_to(gc, 0, dive->cylinder[0].start.mbar);
335 for (i = 1; i < dive->samples; i++) {
337 struct sample *sample = dive->sample + i;
339 mbar = sample->cylinderpressure.mbar;
342 sec = sample->time.seconds;
343 if (sec <= dive->duration.seconds)
344 line_to(gc, sec, mbar);
347 * We may have "surface time" events, in which case we don't go
348 * back to dive duration
350 if (sec < dive->duration.seconds)
351 line_to(gc, dive->duration.seconds, dive->cylinder[0].end.mbar);
352 cairo_stroke(gc->cr);
356 * Return air usage (in liters).
358 static double calculate_airuse(struct dive *dive)
363 for (i = 0; i < MAX_CYLINDERS; i++) {
364 cylinder_t *cyl = dive->cylinder + i;
365 int size = cyl->type.size.mliter;
371 kilo_atm = (cyl->start.mbar - cyl->end.mbar) / 1013250.0;
373 /* Liters of air at 1 atm == milliliters at 1k atm*/
374 airuse += kilo_atm * size;
379 static void plot_info(struct dive *dive, struct graphics_context *gc)
381 text_render_options_t tro = {0.2, 1.0, 0.2, LEFT};
382 const double liters_per_cuft = 28.317;
386 airuse = calculate_airuse(dive);
390 /* I really need to start addign some unit setting thing */
391 switch (output_units.volume) {
397 airuse /= liters_per_cuft;
400 plot_text(gc, &tro, 0.8, 0.8, "vol: %4.2f %s", airuse, unit);
401 if (dive->duration.seconds) {
402 double pressure = 1 + (dive->meandepth.mm / 10000.0);
403 double sac = airuse / pressure * 60 / dive->duration.seconds;
404 plot_text(gc, &tro, 0.8, 0.85, "SAC: %4.2f %s/min", sac, unit);
408 static void plot_cylinder_pressure_text(struct dive *dive, struct graphics_context *gc)
410 pressure_t startp, endp;
412 cairo_set_font_size(gc->cr, 10);
414 if (get_cylinder_pressure_range(dive, gc, &startp, &endp)) {
416 const char *unit = "bar";
418 switch (output_units.pressure) {
420 start = startp.mbar * 100;
421 end = startp.mbar * 100;
425 start = (startp.mbar + 500) / 1000;
426 end = (endp.mbar + 500) / 1000;
430 start = to_PSI(startp);
436 text_render_options_t tro = {0.2, 1.0, 0.2, LEFT};
437 plot_text(gc, &tro, 0, startp.mbar, "%d %s", start, unit);
438 plot_text(gc, &tro, dive->duration.seconds, endp.mbar,
443 static void plot(struct graphics_context *gc, int w, int h, struct dive *dive)
449 cairo_translate(gc->cr, topx, topy);
452 * We can use "cairo_translate()" because that doesn't
453 * scale line width etc. But the actual scaling we need
454 * do set up ourselves..
458 gc->maxx = (w - 2*topx);
459 gc->maxy = (h - 2*topy);
461 /* Cylinder pressure plot */
462 plot_cylinder_pressure(dive, gc);
465 plot_depth_profile(dive, gc);
467 /* Text on top of all graphs.. */
468 plot_depth_text(dive, gc);
469 plot_cylinder_pressure_text(dive, gc);
471 /* And info box in the lower right corner.. */
472 gc->scalex = gc->scaley = 1.0;
475 /* Bounding box last */
476 cairo_set_source_rgb(gc->cr, 1, 1, 1);
481 cairo_close_path(gc->cr);
482 cairo_stroke(gc->cr);
486 static gboolean expose_event(GtkWidget *widget, GdkEventExpose *event, gpointer data)
488 struct dive *dive = current_dive;
489 struct graphics_context gc;
492 w = widget->allocation.width;
493 h = widget->allocation.height;
495 gc.cr = gdk_cairo_create(widget->window);
496 cairo_set_source_rgb(gc.cr, 0, 0, 0);
500 plot(&gc, w, h, dive);
502 cairo_destroy(gc.cr);
507 GtkWidget *dive_profile_widget(void)
511 da = gtk_drawing_area_new();
512 gtk_widget_set_size_request(da, 450, 350);
513 g_signal_connect(da, "expose_event", G_CALLBACK(expose_event), NULL);