}
typedef struct {
- double r,g,b;
- enum {CENTER,LEFT} allign;
+ double r,g,b;
+ enum {CENTER,LEFT} allign;
} text_render_options_t;
static void plot_text(cairo_t *cr, text_render_options_t *tro,
cairo_text_extents(cr, buffer, &extents);
if (tro->allign == CENTER)
- x -= extents.width/2 + extents.x_bearing;
+ x -= extents.width/2 + extents.x_bearing;
y += extents.height * 1.2;
cairo_move_to(cr, x, y);
}
/* Scale to 0,0 -> maxx,maxy */
-#define SCALE(x,y) (x)*maxx/scalex+topx,(y)*maxy/scaley+topy
+#define SCALE(x,y) (x)*maxx/scalex,(y)*maxy/scaley
static void plot_depth_text(struct dive *dive, cairo_t *cr,
- double topx, double topy, double maxx, double maxy)
+ double maxx, double maxy)
{
double scalex, scaley;
int maxtime, maxdepth;
}
static void plot_depth_profile(struct dive *dive, cairo_t *cr,
- double topx, double topy, double maxx, double maxy)
+ double maxx, double maxy)
{
double scalex, scaley;
int begins, sec, depth;
for (i = 1; i < dive->samples; i++) {
sample++;
sec = sample->time.seconds;
- depth = to_feet(sample->depth);
- cairo_line_to(cr, SCALE(sec, depth));
+ if (sec <= maxtime) {
+ depth = to_feet(sample->depth);
+ cairo_line_to(cr, SCALE(sec, depth));
+ }
}
scaley = 1.0;
- cairo_line_to(cr, SCALE(sec, 0));
+ cairo_line_to(cr, SCALE(MIN(sec,maxtime), 0));
cairo_line_to(cr, SCALE(begins, 0));
cairo_close_path(cr);
cairo_set_source_rgba(cr, 1, 0.2, 0.2, 0.20);
/* gets both the actual start and end pressure as well as the scaling factors */
static int get_cylinder_pressure_range(struct dive *dive, double *scalex, double *scaley,
- double *startp, double *endp)
+ pressure_t *startp, pressure_t *endp)
{
int i;
- double min, max;
- double bar;
+ int min, max, mbar;
*scalex = round_seconds_up(dive->duration.seconds);
max = 0;
- min = 5000;
+ min = 5000000;
if (startp)
- *startp = *endp = 0.0;
+ startp->mbar = endp->mbar = 0;
for (i = 0; i < dive->samples; i++) {
struct sample *sample = dive->sample + i;
/* FIXME! We only track cylinder 0 right now */
if (sample->cylinderindex)
continue;
- if (!sample->cylinderpressure.mbar)
+ mbar = sample->cylinderpressure.mbar;
+ if (!mbar)
continue;
- bar = sample->cylinderpressure.mbar;
- if (bar != 0.0 && startp && *startp == 0.0)
- *startp = bar;
- if (bar < min)
- min = bar;
- if (bar > max)
- max = bar;
+ if (mbar && startp && !startp->mbar)
+ startp->mbar = mbar;
+ if (mbar < min)
+ min = mbar;
+ if (mbar > max)
+ max = mbar;
}
if (endp)
- *endp = bar;
+ endp->mbar = mbar;
if (!max)
return 0;
*scaley = max * 1.5;
}
static void plot_cylinder_pressure(struct dive *dive, cairo_t *cr,
- double topx, double topy, double maxx, double maxy)
+ double maxx, double maxy)
{
int i, sec = -1;
double scalex, scaley;
if (!mbar)
continue;
sec = sample->time.seconds;
- cairo_line_to(cr, SCALE(sec, mbar));
+ if (sec <= dive->duration.seconds)
+ cairo_line_to(cr, SCALE(sec, mbar));
}
/*
* We may have "surface time" events, in which case we don't go
cairo_stroke(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, cairo_t *cr,
+ double maxx, double maxy)
+{
+ text_render_options_t tro = {0.2, 1.0, 0.2, LEFT};
+ const double liters_per_cuft = 28.317;
+ double airuse;
+
+ airuse = calculate_airuse(dive);
+ if (!airuse)
+ return;
+
+ /* I really need to start addign some unit setting thing */
+ airuse /= liters_per_cuft;
+ plot_text(cr, &tro, maxx*0.8, maxy*0.8, "cuft: %4.2f", airuse);
+ if (dive->duration.seconds) {
+ double pressure = 1 + (dive->meandepth.mm / 10000.0);
+ double sac = airuse / pressure * 60 / dive->duration.seconds;
+ plot_text(cr, &tro, maxx*0.8, maxy*0.85, "SAC: %4.2f", sac);
+ }
+}
static void plot_cylinder_pressure_text(struct dive *dive, cairo_t *cr,
- double topx, double topy, double maxx, double maxy)
+ double maxx, double maxy)
{
double scalex, scaley;
- double startp,endp;
+ pressure_t startp, endp;
cairo_set_font_size(cr, 10);
if (get_cylinder_pressure_range(dive, &scalex, &scaley,
&startp, &endp)) {
text_render_options_t tro = {0.2, 1.0, 0.2, LEFT};
- plot_text(cr, &tro, SCALE(0, startp), "%3.0f bar", startp/1000.0);
- plot_text(cr, &tro, SCALE(dive->duration.seconds, endp),
- "%3.0f bar", endp/1000.0);
+ plot_text(cr, &tro, SCALE(0, startp.mbar), "%3.0f bar", startp.mbar/1000.0);
+ plot_text(cr, &tro, SCALE(dive->duration.seconds, endp.mbar),
+ "%3.0f bar", endp.mbar/1000.0);
}
}
topy = h / 20.0;
maxx = (w - 2*topx);
maxy = (h - 2*topy);
+ cairo_translate(cr, topx, topy);
/* Cylinder pressure plot */
- plot_cylinder_pressure(dive, cr, topx, topy, maxx, maxy);
+ plot_cylinder_pressure(dive, cr, maxx, maxy);
/* Depth profile */
- plot_depth_profile(dive, cr, topx, topy, maxx, maxy);
+ plot_depth_profile(dive, cr, maxx, maxy);
/* Text on top of all graphs.. */
- plot_depth_text(dive, cr, topx, topy, maxx, maxy);
- plot_cylinder_pressure_text(dive, cr, topx, topy, maxx, maxy);
+ plot_depth_text(dive, cr, maxx, maxy);
+ plot_cylinder_pressure_text(dive, cr, maxx, maxy);
+
+ /* And info box in the lower right corner.. */
+ plot_info(dive, cr, maxx, maxy);
/* Bounding box last */
scalex = scaley = 1.0;
projects / ext / subsurface.git / blobdiff