]> git.tdb.fi Git - ext/subsurface.git/blobdiff - profile.c
Change plot_info to use depth (instead of val) for depth value
[ext/subsurface.git] / profile.c
index 0e5361704fc88cca0214919732b63894b4bb7fa5..eef68f20bbd9a12b9424e20590c548cc5a8a578f 100644 (file)
--- a/profile.c
+++ b/profile.c
@@ -22,12 +22,14 @@ struct plot_info {
        int maxtime;
        int meandepth, maxdepth;
        int minpressure, maxpressure;
+       int endpressure; /* start pressure better be max pressure */
        int mintemp, maxtemp;
        struct plot_data {
+               unsigned int same_cylinder:1;
                int sec;
                int pressure, temperature;
                /* Depth info */
-               int val;
+               int depth;
                int smoothed;
                velocity_t velocity;
                struct plot_data *min[3];
@@ -159,14 +161,31 @@ static void plot_text(struct graphics_context *gc, const text_render_options_t *
 static void plot_one_event(struct graphics_context *gc, struct plot_info *pi, struct event *event, const text_render_options_t *tro)
 {
        int i, depth = 0;
+       int x,y;
 
        for (i = 0; i < pi->nr; i++) {
                struct plot_data *data = pi->entry + i;
                if (event->time.seconds < data->sec)
                        break;
-               depth = data->val;
+               depth = data->depth;
        }
-       plot_text(gc, tro, event->time.seconds, depth, "%s", event->name);
+       /* draw a little tirangular marker and attach tooltip */
+       x = SCALEX(gc, event->time.seconds);
+       y = SCALEY(gc, depth);
+       set_source_rgba(gc, 1.0, 1.0, 0.1, 0.8);
+       cairo_move_to(gc->cr, x-15, y+6);
+       cairo_line_to(gc->cr, x-3  , y+6);
+       cairo_line_to(gc->cr, x-9, y-6);
+       cairo_line_to(gc->cr, x-15, y+6);
+       cairo_stroke_preserve(gc->cr);
+       cairo_fill(gc->cr);
+       set_source_rgba(gc, 0.0, 0.0, 0.0, 0.8);
+       cairo_move_to(gc->cr, x-9, y-3);
+       cairo_line_to(gc->cr, x-9, y+1);
+       cairo_move_to(gc->cr, x-9, y+4);
+       cairo_line_to(gc->cr, x-9, y+4);
+       cairo_stroke(gc->cr);
+       attach_tooltip(x-15, y-6, 12, 12, event->name);
 }
 
 static void plot_events(struct graphics_context *gc, struct plot_info *pi, struct dive *dive)
@@ -188,9 +207,9 @@ static void render_depth_sample(struct graphics_context *gc, struct plot_data *e
        int sec = entry->sec, decimals;
        double d;
 
-       d = get_depth_units(entry->val, &decimals, NULL);
+       d = get_depth_units(entry->depth, &decimals, NULL);
 
-       plot_text(gc, tro, sec, entry->val, "%.*f", decimals, d);
+       plot_text(gc, tro, sec, entry->depth, "%.*f", decimals, d);
 }
 
 static void plot_text_samples(struct graphics_context *gc, struct plot_info *pi)
@@ -202,7 +221,7 @@ static void plot_text_samples(struct graphics_context *gc, struct plot_info *pi)
        for (i = 0; i < pi->nr; i++) {
                struct plot_data *entry = pi->entry + i;
 
-               if (entry->val < 2000)
+               if (entry->depth < 2000)
                        continue;
 
                if (entry == entry->max[2])
@@ -248,13 +267,13 @@ static void plot_minmax_profile_minute(struct graphics_context *gc, struct plot_
        struct plot_data *entry = pi->entry;
 
        set_source_rgba(gc, 1, 0.2, 1, a);
-       move_to(gc, entry->sec, entry->min[index]->val);
+       move_to(gc, entry->sec, entry->min[index]->depth);
        for (i = 1; i < pi->nr; i++) {
                entry++;
-               line_to(gc, entry->sec, entry->min[index]->val);
+               line_to(gc, entry->sec, entry->min[index]->depth);
        }
        for (i = 1; i < pi->nr; i++) {
-               line_to(gc, entry->sec, entry->max[index]->val);
+               line_to(gc, entry->sec, entry->max[index]->depth);
                entry--;
        }
        cairo_close_path(gc->cr);
@@ -306,10 +325,12 @@ static void plot_depth_profile(struct graphics_context *gc, struct plot_info *pi
        cairo_stroke(cr);
 
        /* Show mean depth */
-       set_source_rgba(gc, 1, 0.2, 0.2, 0.40);
-       move_to(gc, 0, pi->meandepth);
-       line_to(gc, 1, pi->meandepth);
-       cairo_stroke(cr);
+       if (! gc->printer) {
+               set_source_rgba(gc, 1, 0.2, 0.2, 0.40);
+               move_to(gc, 0, pi->meandepth);
+               line_to(gc, 1, pi->meandepth);
+               cairo_stroke(cr);
+       }
 
        gc->leftx = 0; gc->rightx = maxtime;
 
@@ -331,7 +352,7 @@ static void plot_depth_profile(struct graphics_context *gc, struct plot_info *pi
        entry = pi->entry;
        move_to(gc, 0, 0);
        for (i = 0; i < pi->nr; i++, entry++)
-               line_to(gc, entry->sec, entry->val);
+               line_to(gc, entry->sec, entry->depth);
        cairo_close_path(gc->cr);
        if (gc->printer) {
                set_source_rgba(gc, 1, 1, 1, 0.2);
@@ -351,9 +372,9 @@ static void plot_depth_profile(struct graphics_context *gc, struct plot_info *pi
                 * representing the vertical velocity, so we need to
                 * chop this into short segments */
                rgb_t color = rgb[entry->velocity];
-               depth = entry->val;
+               depth = entry->depth;
                set_source_rgb(gc, color.r, color.g, color.b);
-               move_to(gc, entry[-1].sec, entry[-1].val);
+               move_to(gc, entry[-1].sec, entry[-1].depth);
                line_to(gc, sec, depth);
                cairo_stroke(cr);
        }
@@ -469,23 +490,32 @@ static int get_cylinder_pressure_range(struct graphics_context *gc, struct plot_
 static void plot_cylinder_pressure(struct graphics_context *gc, struct plot_info *pi)
 {
        int i;
+       int have_pressure = FALSE;
 
        if (!get_cylinder_pressure_range(gc, pi))
                return;
 
        set_source_rgba(gc, 0.2, 1.0, 0.2, 0.80);
 
-       move_to(gc, 0, pi->maxpressure);
-       for (i = 1; i < pi->nr; i++) {
+       for (i = 0; i < pi->nr; i++) {
                int mbar;
                struct plot_data *entry = pi->entry + i;
 
                mbar = entry->pressure;
                if (!mbar)
                        continue;
-               line_to(gc, entry->sec, mbar);
+               have_pressure = TRUE;
+               if (entry->same_cylinder)
+                       line_to(gc, entry->sec, mbar);
+               else
+                       move_to(gc, entry->sec, mbar);
        }
-       line_to(gc, pi->maxtime, pi->minpressure);
+       /* if we have valid samples, we don't want to draw a line to the minpressure
+        * but just end wherever the dive ended (think valve shutdowns during dive)
+        * but that doesn't work so well if we have only max and min
+        */
+       if (! have_pressure)
+               line_to(gc, pi->maxtime, pi->minpressure);
        cairo_stroke(gc->cr);
 }
 
@@ -504,24 +534,24 @@ static void plot_cylinder_pressure_text(struct graphics_context *gc, struct plot
                switch (output_units.pressure) {
                case PASCAL:
                        start = pi->maxpressure * 100;
-                       end = pi->minpressure * 100;
+                       end = pi->endpressure * 100;
                        unit = "pascal";
                        break;
                case BAR:
                        start = (pi->maxpressure + 500) / 1000;
-                       end = (pi->minpressure + 500) / 1000;
+                       end = (pi->endpressure + 500) / 1000;
                        unit = "bar";
                        break;
                case PSI:
                        start = mbar_to_PSI(pi->maxpressure);
-                       end = mbar_to_PSI(pi->minpressure);
+                       end = mbar_to_PSI(pi->endpressure);
                        unit = "psi";
                        break;
                }
 
                text_render_options_t tro = {10, 0.2, 1.0, 0.2, LEFT, TOP};
                plot_text(gc, &tro, 0, pi->maxpressure, "%d %s", start, unit);
-               plot_text(gc, &tro, pi->maxtime, pi->minpressure,
+               plot_text(gc, &tro, pi->maxtime, pi->endpressure,
                          "%d %s", end, unit);
        }
 }
@@ -543,17 +573,17 @@ static void analyze_plot_info_minmax_minute(struct plot_data *entry, struct plot
 
        /* Then go forward until we hit an entry past the time */
        min = max = p;
-       avg = p->val;
+       avg = p->depth;
        nr = 1;
        while (++p < last) {
-               int val = p->val;
+               int depth = p->depth;
                if (p->sec > time + seconds)
                        break;
-               avg += val;
+               avg += depth;
                nr ++;
-               if (val < min->val)
+               if (depth < min->depth)
                        min = p;
-               if (val > max->val)
+               if (depth > max->depth)
                        max = p;
        }
        entry->min[index] = min;
@@ -622,24 +652,24 @@ static struct plot_info *analyze_plot_info(struct plot_info *pi)
        }
 
        /* Smoothing function: 5-point triangular smooth */
-       for (i = 2; i < nr-1; i++) {
+       for (i = 2; i < nr; i++) {
                struct plot_data *entry = pi->entry+i;
-               int val;
+               int depth;
 
                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;
+                       depth = entry[-2].depth + 2*entry[-1].depth + 3*entry[0].depth + 2*entry[1].depth + entry[2].depth;
+                       entry->smoothed = (depth+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));
+                       entry->velocity = velocity((entry[0].depth - entry[-1].depth) / (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) { 
+                       if (entry[0].sec - entry[-1].sec < 15 && entry->velocity < FAST) {
                                int past = -2;
-                               while (i+past > 0 && entry[0].sec - entry[past].sec < 30)
+                               while (i+past > 0 && entry[0].sec - entry[past].sec < 15)
                                        past--;
-                               entry->velocity = velocity((entry[0].val - entry[past].val) / 
+                               entry->velocity = velocity((entry[0].depth - entry[past].depth) / 
                                                        (entry[0].sec - entry[past].sec));
                        }
                } else
@@ -664,6 +694,7 @@ static struct plot_info *analyze_plot_info(struct plot_info *pi)
  */
 static struct plot_info *create_plot_info(struct dive *dive)
 {
+       int cylinderindex = -1;
        int lastdepth, lastindex;
        int i, nr = dive->samples + 4, sec;
        size_t alloc_size = plot_info_size(nr);
@@ -683,8 +714,13 @@ static struct plot_info *create_plot_info(struct dive *dive)
                struct plot_data *entry = pi->entry + i + 2;
 
                sec = entry->sec = sample->time.seconds;
-               depth = entry->val = sample->depth.mm;
+               depth = entry->depth = sample->depth.mm;
+
+               entry->same_cylinder = sample->cylinderindex == cylinderindex;
+               cylinderindex = sample->cylinderindex;
                entry->pressure = sample->cylinderpressure.mbar;
+               if (!entry->same_cylinder && !entry->pressure)
+                       entry->pressure = dive->cylinder[cylinderindex].start.mbar;
                entry->temperature = sample->temperature.mkelvin;
 
                if (depth || lastdepth)
@@ -700,11 +736,15 @@ static struct plot_info *create_plot_info(struct dive *dive)
        i = dive->samples + 2;
        pi->entry[i].sec = sec + 20;
        pi->entry[i+1].sec = sec + 40;
+       if (cylinderindex >= 0) {
+               pi->entry[i].pressure = dive->cylinder[cylinderindex].end.mbar;
+               pi->entry[i].same_cylinder = 1;
+       }
 
        pi->nr = lastindex+1;
        pi->maxtime = pi->entry[lastindex].sec;
 
-       pi->minpressure = dive->cylinder[0].end.mbar;
+       pi->endpressure = pi->minpressure = dive->cylinder[0].end.mbar;
        pi->maxpressure = dive->cylinder[0].start.mbar;
 
        pi->meandepth = dive->meandepth.mm;
@@ -712,14 +752,11 @@ static struct plot_info *create_plot_info(struct dive *dive)
        return analyze_plot_info(pi);
 }
 
-void plot(struct graphics_context *gc, int w, int h, struct dive *dive)
+void plot(struct graphics_context *gc, cairo_rectangle_int_t *drawing_area, struct dive *dive)
 {
-       double topx, topy;
        struct plot_info *pi = create_plot_info(dive);
 
-       topx = w / 20.0;
-       topy = h / 20.0;
-       cairo_translate(gc->cr, topx, topy);
+       cairo_translate(gc->cr, drawing_area->x, drawing_area->y);
        cairo_set_line_width(gc->cr, 2);
        cairo_set_line_cap(gc->cr, CAIRO_LINE_CAP_ROUND);
        cairo_set_line_join(gc->cr, CAIRO_LINE_JOIN_ROUND);
@@ -731,8 +768,8 @@ void plot(struct graphics_context *gc, int w, int h, struct dive *dive)
         *
         * Snif. What a pity.
         */
-       gc->maxx = (w - 2*topx);
-       gc->maxy = (h - 2*topy);
+       gc->maxx = (drawing_area->width - 2*drawing_area->x);
+       gc->maxy = (drawing_area->height - 2*drawing_area->y);
 
        /* Temperature profile */
        plot_temperature_profile(gc, pi);