statistics.c

   1 /* statistics.c */
   2 /* creates the UI for the Info & Stats page -
   3  * controlled through the following interfaces:
   4  *
   5  * void show_dive_stats(struct dive *dive)
   6  * void flush_dive_stats_changes(struct dive *dive)
   7  *
   8  * called from gtk-ui:
   9  * GtkWidget *stats_widget(void)
  10  */
  11 #include <stdio.h>
  12 #include <string.h>
  13 #include <stdlib.h>
  14 #include <stdarg.h>
  15 #include <time.h>
  16
  17 #include "dive.h"
  18 #include "display.h"
  19 #include "display-gtk.h"
  20 #include "divelist.h"
  21
  22 typedef struct {
  23         GtkWidget *date,
  24                 *dive_time,
  25                 *surf_intv,
  26                 *max_depth,
  27                 *avg_depth,
  28                 *water_temp,
  29                 *sac,
  30                 *otu,
  31                 *o2he,
  32                 *gas_used,
  33                 *total_time,
  34                 *avg_time,
  35                 *max_overall_depth,
  36                 *avg_overall_depth,
  37                 *min_sac,
  38                 *avg_sac,
  39                 *max_sac;
  40 } info_stat_widget_t;
  41
  42 static info_stat_widget_t info_stat_w;
  43
  44 typedef struct {
  45         duration_t total_time;
  46         /* avg_time is simply total_time / nr -- let's not keep this */
  47         depth_t max_depth;
  48         depth_t avg_depth;
  49         volume_t max_sac;
  50         volume_t min_sac;
  51         volume_t avg_sac;
  52 } info_stat_t;
  53
  54 static info_stat_t info_stat;
  55
  56 static void process_all_dives(struct dive *dive, struct dive **prev_dive)
  57 {
  58         int idx;
  59         struct dive *dp;
  60         int old_tt, sac_time = 0;
  61
  62         *prev_dive = NULL;
  63         memset(&info_stat, 0, sizeof(info_stat));
  64         /* this relies on the fact that the dives in the dive_table
  65          * are in chronological order */
  66         for (idx = 0; idx < dive_table.nr; idx++) {
  67                 dp = dive_table.dives[idx];
  68                 if (dp->when == dive->when) {
  69                         /* that's the one we are showing */
  70                         if (idx > 0)
  71                                 *prev_dive = dive_table.dives[idx-1];
  72                 }
  73                 old_tt = info_stat.total_time.seconds;
  74                 info_stat.total_time.seconds += dp->duration.seconds;
  75                 if (dp->maxdepth.mm > info_stat.max_depth.mm)
  76                         info_stat.max_depth.mm = dp->maxdepth.mm;
  77                 info_stat.avg_depth.mm = (1.0 * old_tt * info_stat.avg_depth.mm +
  78                                 dp->duration.seconds * dp->meandepth.mm) / info_stat.total_time.seconds;
  79                 if (dp->sac > 2800) { /* less than .1 cuft/min (2800ml/min) is bogus */
  80                         int old_sac_time = sac_time;
  81                         sac_time += dp->duration.seconds;
  82                         info_stat.avg_sac.mliter = (1.0 * old_sac_time * info_stat.avg_sac.mliter +
  83                                                 dp->duration.seconds * dp->sac) / sac_time ;
  84                         if (dp->sac > info_stat.max_sac.mliter)
  85                                 info_stat.max_sac.mliter = dp->sac;
  86                         if (info_stat.min_sac.mliter == 0 || dp->sac < info_stat.min_sac.mliter)
  87                                 info_stat.min_sac.mliter = dp->sac;
  88                 }
  89         }
  90 }
  91
  92 static void set_label(GtkWidget *w, const char *fmt, ...)
  93 {
  94         char buf[80];
  95         va_list args;
  96
  97         va_start(args, fmt);
  98         vsnprintf(buf, sizeof(buf), fmt, args);
  99         va_end(args);
 100         gtk_label_set_text(GTK_LABEL(w), buf);
 101 }
 102
 103 static char * get_time_string(int seconds, int maxdays)
 104 {
 105         static char buf[80];
 106         if (maxdays && seconds > 3600 * 24 * maxdays)
 107                 snprintf(buf, sizeof(buf), "more than %d days", maxdays);
 108         else {
 109                 int days = seconds / 3600 / 24;
 110                 int hours = (seconds - days * 3600 * 24) / 3600;
 111                 int minutes = (seconds - days * 3600 * 24 - hours * 3600) / 60;
 112                 if (days > 0)
 113                         snprintf(buf, sizeof(buf), "%dd %dh %dmin", days, hours, minutes);
 114                 else
 115                         snprintf(buf, sizeof(buf), "%dh %dmin", hours, minutes);
 116         }
 117         return buf;
 118 }
 119
 120 void show_dive_stats(struct dive *dive)
 121 {
 122         char buf[80];
 123         double value;
 124         int decimals;
 125         const char *unit;
 126         int idx, offset, gas_used;
 127         struct dive *prev_dive;
 128         struct tm *tm;
 129
 130         process_all_dives(dive, &prev_dive);
 131
 132         tm = gmtime(&dive->when);
 133         snprintf(buf, sizeof(buf),
 134                 "%s, %s %d, %d %2d:%02d",
 135                 weekday(tm->tm_wday),
 136                 monthname(tm->tm_mon),
 137                 tm->tm_mday, tm->tm_year + 1900,
 138                 tm->tm_hour, tm->tm_min);
 139
 140         set_label(info_stat_w.date, buf);
 141         set_label(info_stat_w.dive_time, "%d min", (dive->duration.seconds + 30) / 60);
 142         if (prev_dive)
 143                 set_label(info_stat_w.surf_intv,
 144                         get_time_string(dive->when - (prev_dive->when + prev_dive->duration.seconds), 4));
 145         else
 146                 set_label(info_stat_w.surf_intv, "unknown");
 147         value = get_depth_units(dive->maxdepth.mm, &decimals, &unit);
 148         set_label(info_stat_w.max_depth, "%.*f %s", decimals, value, unit);
 149         value = get_depth_units(dive->meandepth.mm, &decimals, &unit);
 150         set_label(info_stat_w.avg_depth, "%.*f %s", decimals, value, unit);
 151         if (dive->watertemp.mkelvin) {
 152                 value = get_temp_units(dive->watertemp.mkelvin, &unit);
 153                 set_label(info_stat_w.water_temp, "%.1f %s", value, unit);
 154         } else
 155                 set_label(info_stat_w.water_temp, "");
 156         value = get_volume_units(dive->sac, &decimals, &unit);
 157         if (value > 0) {
 158                 set_label(info_stat_w.sac, "%.*f %s/min", decimals, value, unit);
 159         } else
 160                 set_label(info_stat_w.sac, "");
 161         set_label(info_stat_w.otu, "%d", dive->otu);
 162         offset = 0;
 163         gas_used = 0;
 164         buf[0] = '\0';
 165         /* for the O2/He readings just create a list of them */
 166         for (idx = 0; idx < MAX_CYLINDERS; idx++) {
 167                 cylinder_t *cyl = &dive->cylinder[idx];
 168                 /* we assume that every valid cylinder has either a working pressure
 169                  * or a size; but for good measure let's also accept cylinders with
 170                  * a starting or ending pressure*/
 171                 if (cyl->type.workingpressure.mbar || cyl->type.size.mliter ||
 172                         cyl->start.mbar || cyl->end.mbar) {
 173                         /* 0% O2 strangely means air, so 21% - I don't like that at all */
 174                         int o2 = cyl->gasmix.o2.permille ? : 209;
 175                         if (offset > 0) {
 176                                 snprintf(buf+offset, 80-offset, ", ");
 177                                 offset += 2;
 178                         }
 179                         snprintf(buf+offset, 80-offset, "%d/%d", (o2 + 5) / 10,
 180                                 (cyl->gasmix.he.permille + 5) / 10);
 181                         offset = strlen(buf);
 182                 }
 183                 /* and if we have size, start and end pressure, we can
 184                  * calculate the total gas used */
 185                 if (cyl->type.size.mliter && cyl->start.mbar && cyl->end.mbar)
 186                         gas_used += cyl->type.size.mliter / 1000.0 *
 187                                 (cyl->start.mbar - cyl->end.mbar);
 188         }
 189         set_label(info_stat_w.o2he, buf);
 190         if (gas_used) {
 191                 value = get_volume_units(gas_used, &decimals, &unit);
 192                 set_label(info_stat_w.gas_used, "%.*f %s", decimals, value, unit);
 193         } else
 194                 set_label(info_stat_w.gas_used, "");
 195         /* and now do the statistics */
 196         set_label(info_stat_w.total_time, get_time_string(info_stat.total_time.seconds, 0));
 197         set_label(info_stat_w.avg_time, get_time_string(info_stat.total_time.seconds / dive_table.nr, 0));
 198         value = get_depth_units(info_stat.max_depth.mm, &decimals, &unit);
 199         set_label(info_stat_w.max_overall_depth, "%.*f %s", decimals, value, unit);
 200         value = get_depth_units(info_stat.avg_depth.mm, &decimals, &unit);
 201         set_label(info_stat_w.avg_overall_depth, "%.*f %s", decimals, value, unit);
 202         value = get_volume_units(info_stat.max_sac.mliter, &decimals, &unit);
 203         set_label(info_stat_w.max_sac, "%.*f %s/min", decimals, value, unit);
 204         value = get_volume_units(info_stat.min_sac.mliter, &decimals, &unit);
 205         set_label(info_stat_w.min_sac, "%.*f %s/min", decimals, value, unit);
 206         value = get_volume_units(info_stat.avg_sac.mliter, &decimals, &unit);
 207         set_label(info_stat_w.avg_sac, "%.*f %s/min", decimals, value, unit);
 208 }
 209
 210 void flush_dive_stats_changes(struct dive *dive)
 211 {
 212         /* We do nothing: we require the "Ok" button press */
 213 }
 214
 215 static GtkWidget *new_info_label_in_frame(GtkWidget *box, const char *label)
 216 {
 217         GtkWidget *label_widget;
 218         GtkWidget *frame;
 219
 220         frame = gtk_frame_new(label);
 221         label_widget = gtk_label_new(NULL);
 222         gtk_box_pack_start(GTK_BOX(box), frame, TRUE, TRUE, 3);
 223         gtk_container_add(GTK_CONTAINER(frame), label_widget);
 224
 225         return label_widget;
 226 }
 227
 228 GtkWidget *stats_widget(void)
 229 {
 230
 231         GtkWidget *vbox, *hbox, *infoframe, *statsframe, *framebox;
 232
 233         vbox = gtk_vbox_new(FALSE, 3);
 234
 235         infoframe = gtk_frame_new("Dive Info");
 236         gtk_box_pack_start(GTK_BOX(vbox), infoframe, TRUE, FALSE, 3);
 237         framebox = gtk_vbox_new(FALSE, 3);
 238         gtk_container_add(GTK_CONTAINER(infoframe), framebox);
 239
 240         /* first row */
 241         hbox = gtk_hbox_new(FALSE, 3);
 242         gtk_box_pack_start(GTK_BOX(framebox), hbox, TRUE, FALSE, 3);
 243
 244         info_stat_w.date = new_info_label_in_frame(hbox, "Date");
 245         info_stat_w.dive_time = new_info_label_in_frame(hbox, "Dive Time");
 246         info_stat_w.surf_intv = new_info_label_in_frame(hbox, "Surf Intv");
 247
 248         /* second row */
 249         hbox = gtk_hbox_new(FALSE, 3);
 250         gtk_box_pack_start(GTK_BOX(framebox), hbox, TRUE, FALSE, 3);
 251
 252         info_stat_w.max_depth = new_info_label_in_frame(hbox, "Max Depth");
 253         info_stat_w.avg_depth = new_info_label_in_frame(hbox, "Avg Depth");
 254         info_stat_w.water_temp = new_info_label_in_frame(hbox, "Water Temp");
 255
 256         /* third row */
 257         hbox = gtk_hbox_new(FALSE, 3);
 258         gtk_box_pack_start(GTK_BOX(framebox), hbox, TRUE, FALSE, 3);
 259
 260         info_stat_w.sac = new_info_label_in_frame(hbox, "SAC");
 261         info_stat_w.otu = new_info_label_in_frame(hbox, "OTU");
 262         info_stat_w.o2he = new_info_label_in_frame(hbox, "O" UTF8_SUBSCRIPT_2 " / He");
 263         info_stat_w.gas_used = new_info_label_in_frame(hbox, "Gas Used");
 264
 265         statsframe = gtk_frame_new("Statistics");
 266         gtk_box_pack_start(GTK_BOX(vbox), statsframe, TRUE, FALSE, 3);
 267         framebox = gtk_vbox_new(FALSE, 3);
 268         gtk_container_add(GTK_CONTAINER(statsframe), framebox);
 269
 270         /* first row */
 271         hbox = gtk_hbox_new(FALSE, 3);
 272         gtk_box_pack_start(GTK_BOX(framebox), hbox, TRUE, FALSE, 3);
 273
 274         info_stat_w.total_time = new_info_label_in_frame(hbox, "Total Time");
 275         info_stat_w.avg_time = new_info_label_in_frame(hbox, "Avg Time");
 276         info_stat_w.max_overall_depth = new_info_label_in_frame(hbox, "Max Depth");
 277         info_stat_w.avg_overall_depth = new_info_label_in_frame(hbox, "Avg Depth");
 278
 279         /* second row */
 280         hbox = gtk_hbox_new(FALSE, 3);
 281         gtk_box_pack_start(GTK_BOX(framebox), hbox, TRUE, FALSE, 3);
 282
 283         info_stat_w.max_sac = new_info_label_in_frame(hbox, "Max SAC");
 284         info_stat_w.min_sac = new_info_label_in_frame(hbox, "Min SAC");
 285         info_stat_w.avg_sac = new_info_label_in_frame(hbox, "Avg SAC");
 286
 287         return vbox;
 288 }