7 #include <libxml/parser.h>
8 #include <libxml/tree.h>
13 * Some silly typedefs to make our units very explicit.
15 * Also, the units are chosen so that values can be expressible as
16 * integers, so that we never have FP rounding issues. And they
17 * are small enough that converting to/from imperial units doesn't
20 * We also strive to make '0' a meaningless number saying "not
21 * initialized", since many values are things that may not have
22 * been reported (eg tank pressure or temperature from dive
23 * computers that don't support them). But sometimes -1 is an even
24 * more explicit way of saying "not there".
26 * Thus "millibar" for pressure, for example, or "millikelvin" for
27 * temperatures. Doing temperatures in celsius or fahrenheit would
28 * make for loss of precision when converting from one to the other,
29 * and using millikelvin is SI-like but also means that a temperature
30 * of '0' is clearly just a missing temperature or tank pressure.
32 * Also strive to use units that can not possibly be mistaken for a
33 * valid value in a "normal" system without conversion. If the max
34 * depth of a dive is '20000', you probably didn't convert from mm on
35 * output, or if the max depth gets reported as "0.2ft" it was either
36 * a really boring dive, or there was some missing input conversion,
37 * and a 60-ft dive got recorded as 60mm.
39 * Doing these as "structs containing value" means that we always
40 * have to explicitly write out those units in order to get at the
41 * actual value. So there is hopefully little fear of using a value
42 * in millikelvin as Fahrenheit by mistake.
44 * We don't actually use these all yet, so maybe they'll change, but
45 * I made a number of types as guidelines.
86 static int to_feet(depth_t depth)
88 return depth.mm * 0.00328084 + 0.5;
91 static int to_C(temperature_t temp)
95 return (temp.mkelvin - 273150) / 1000;
98 static int to_PSI(pressure_t pressure)
100 return pressure.mbar * 0.0145037738 + 0.5;
106 temperature_t temperature;
107 pressure_t tankpressure;
113 depth_t maxdepth, meandepth;
114 duration_t duration, surfacetime;
116 temperature_t airtemp, watertemp;
117 pressure_t beginning_pressure, end_pressure;
119 struct sample sample[];
122 static struct dive **dive_table;
123 static int nr_dives, nr_allocated;
125 static void record_dive(struct dive *dive)
127 if (nr_dives >= nr_allocated) {
128 nr_allocated = (nr_dives + 32) * 3 / 2;
129 dive_table = realloc(dive_table, nr_allocated * sizeof(struct dive *));
133 dive_table[nr_dives++] = dive;
136 static void show_dive(int nr, struct dive *dive)
141 tm = gmtime(&dive->when);
143 printf("Dive %d with %d samples at %02d:%02d:%02d %04d-%02d-%02d\n",
145 tm->tm_hour, tm->tm_min, tm->tm_sec,
146 tm->tm_year+1900, tm->tm_mon+1, tm->tm_mday);
151 for (i = 0; i < dive->samples; i++) {
152 struct sample *s = dive->sample + i;
154 printf("%4d:%02d: %3d ft, %2d C, %4d PSI\n",
155 s->time.seconds / 60,
156 s->time.seconds % 60,
158 to_C(s->temperature),
159 to_PSI(s->tankpressure));
163 static int sortfn(const void *_a, const void *_b)
165 const struct dive *a = *(void **)_a;
166 const struct dive *b = *(void **)_b;
168 if (a->when < b->when)
170 if (a->when > b->when)
175 static void report_dives(void)
178 qsort(dive_table, nr_dives, sizeof(struct dive *), sortfn);
180 for (i = 0; i < nr_dives; i++)
181 show_dive(i+1, dive_table[i]);
184 static void nonmatch(const char *type, const char *fullname, const char *name, char *buffer)
187 printf("Unable to match %s '(%.*s)%s' (%s)\n", type,
188 (int) (name - fullname), fullname, name,
193 static const char *last_part(const char *name)
195 const char *p = strrchr(name, '.');
196 return p ? p+1 : name;
199 typedef void (*matchfn_t)(char *buffer, void *);
201 static int match(const char *pattern, const char *name, matchfn_t fn, char *buf, void *data)
203 if (strcasecmp(pattern, name))
210 * Dive info as it is being built up..
212 static int alloc_samples;
213 static struct dive *dive;
214 static struct sample *sample;
217 static time_t utc_mktime(struct tm *tm)
219 static const int mdays[] = {
220 0, 31, 59, 90, 120, 151, 181, 212, 243, 273, 304, 334
222 int year = tm->tm_year;
223 int month = tm->tm_mon;
224 int day = tm->tm_mday;
226 /* First normalize relative to 1900 */
229 else if (year > 1900)
232 /* Normalized to Jan 1, 1970: unix time */
235 if (year < 0 || year > 129) /* algo only works for 1970-2099 */
237 if (month < 0 || month > 11) /* array bounds */
239 if (month < 2 || (year + 2) % 4)
241 if (tm->tm_hour < 0 || tm->tm_min < 0 || tm->tm_sec < 0)
243 return (year * 365 + (year + 1) / 4 + mdays[month] + day) * 24*60*60UL +
244 tm->tm_hour * 60*60 + tm->tm_min * 60 + tm->tm_sec;
247 static void divedate(char *buffer, void *_when)
250 time_t *when = _when;
252 if (sscanf(buffer, "%d.%d.%d", &d, &m, &y) == 3) {
256 if (tm.tm_sec | tm.tm_min | tm.tm_hour)
257 *when = utc_mktime(&tm);
262 static void divetime(char *buffer, void *_when)
265 time_t *when = _when;
267 if (sscanf(buffer, "%d:%d:%d", &h, &m, &s) >= 2) {
272 *when = utc_mktime(&tm);
277 /* Libdivecomputer: "2011-03-20 10:22:38" */
278 static void divedatetime(char *buffer, void *_when)
282 time_t *when = _when;
284 if (sscanf(buffer, "%d-%d-%d %d:%d:%d",
285 &y, &m, &d, &hr, &min, &sec) == 6) {
292 *when = utc_mktime(&tm);
308 static enum number_type integer_or_float(char *buffer, union int_or_float *res)
314 /* Integer or floating point? */
315 val = strtol(buffer, &end, 10);
316 if (val < 0 || end == buffer)
319 /* Looks like it might be floating point? */
322 fp = strtod(buffer, &end);
333 static void pressure(char *buffer, void *_press)
335 pressure_t *pressure = _press;
336 union int_or_float val;
338 switch (integer_or_float(buffer, &val)) {
340 /* Maybe it's in Bar? */
341 if (val.fp < 500.0) {
342 pressure->mbar = val.fp * 1000;
345 printf("Unknown fractional pressure reading %s\n", buffer);
350 * Random integer? Maybe in PSI? Or millibar already?
352 * We assume that 5 bar is a ridiculous tank pressure,
353 * so if it's smaller than 5000, it's in PSI..
356 pressure->mbar = val.i * 68.95;
359 pressure->mbar = val.i;
362 printf("Strange pressure reading %s\n", buffer);
367 static void depth(char *buffer, void *_depth)
369 depth_t *depth = _depth;
370 union int_or_float val;
372 switch (integer_or_float(buffer, &val)) {
373 /* Integer values are probably in feet */
375 depth->mm = 304.8 * val.i;
377 /* Float? Probably meters.. */
379 depth->mm = val.fp * 1000;
382 printf("Strange depth reading %s\n", buffer);
387 static void temperature(char *buffer, void *_temperature)
389 temperature_t *temperature = _temperature;
390 union int_or_float val;
392 switch (integer_or_float(buffer, &val)) {
393 /* C or F? Who knows? Let's default to Celsius */
398 /* Ignore zero. It means "none" */
403 temperature->mkelvin = (val.fp + 273.16) * 1000;
407 if (val.fp < 212.0) {
408 temperature->mkelvin = (val.fp + 459.67) * 5000/9;
411 /* Kelvin or already millikelvin */
414 temperature->mkelvin = val.fp;
417 printf("Strange temperature reading %s\n", buffer);
422 static void sampletime(char *buffer, void *_time)
424 duration_t *time = _time;
425 union int_or_float val;
427 switch (integer_or_float(buffer, &val)) {
429 time->seconds = val.i;
432 printf("Strange sample time reading %s\n", buffer);
437 /* We're in samples - try to convert the random xml value to something useful */
438 static void try_to_fill_sample(struct sample *sample, const char *name, char *buf)
440 const char *last = last_part(name);
442 if (match("pressure", last, pressure, buf, &sample->tankpressure))
444 if (match("cylpress", last, pressure, buf, &sample->tankpressure))
446 if (match("depth", last, depth, buf, &sample->depth))
448 if (match("temperature", last, temperature, buf, &sample->temperature))
450 if (match("sampletime", last, sampletime, buf, &sample->time))
452 if (match("time", last, sampletime, buf, &sample->time))
455 nonmatch("sample", name, last, buf);
458 /* We're in the top-level dive xml. Try to convert whatever value to a dive value */
459 static void try_to_fill_dive(struct dive *dive, const char *name, char *buf)
461 const char *last = last_part(name);
463 if (match("date", last, divedate, buf, &dive->when))
465 if (match("time", last, divetime, buf, &dive->when))
467 if (match("datetime", last, divedatetime, buf, &dive->when))
469 nonmatch("dive", name, last, buf);
472 static unsigned int dive_size(int samples)
474 return sizeof(struct dive) + samples*sizeof(struct sample);
478 * File boundaries are dive boundaries. But sometimes there are
479 * multiple dives per file, so there can be other events too that
480 * trigger a "new dive" marker and you may get some nesting due
481 * to that. Just ignore nesting levels.
483 static void dive_start(void)
491 size = dive_size(alloc_samples);
495 memset(dive, 0, size);
496 memset(&tm, 0, sizeof(tm));
499 static void dive_end(void)
507 static void sample_start(void)
514 if (nr >= alloc_samples) {
517 alloc_samples = (alloc_samples * 3)/2 + 10;
518 size = dive_size(alloc_samples);
519 dive = realloc(dive, size);
523 sample = dive->sample + nr;
524 memset(sample, 0, sizeof(*sample));
527 static void sample_end(void)
535 static void entry(const char *name, int size, const char *raw)
537 char *buf = malloc(size+1);
541 memcpy(buf, raw, size);
544 try_to_fill_sample(sample, name, buf);
548 try_to_fill_dive(dive, name, buf);
553 static const char *nodename(xmlNode *node, char *buf, int len)
556 if (!node || !node->name)
564 const char *name = node->name;
565 int i = strlen(name);
567 unsigned char c = name[i];
573 if (!node || !node->name)
583 static void visit_one_node(xmlNode *node)
586 const unsigned char *content;
587 char buffer[MAXNAME];
590 content = node->content;
594 /* Trim whitespace at beginning */
595 while (isspace(*content))
598 /* Trim whitespace at end */
599 len = strlen(content);
600 while (len && isspace(content[len-1]))
606 /* Don't print out the node name if it is "text" */
607 if (!strcmp(node->name, "text"))
610 name = nodename(node, buffer, sizeof(buffer));
612 entry(name, len, content);
615 static void traverse(xmlNode *node)
619 for (n = node; n; n = n->next) {
620 /* XML from libdivecomputer: 'dive' per new dive */
621 if (!strcmp(n->name, "dive")) {
623 traverse(n->children);
629 * At least both libdivecomputer and Suunto
632 * Well - almost. Ignore case.
634 if (!strcasecmp(n->name, "sample")) {
636 traverse(n->children);
641 /* Anything else - just visit it and recurse */
643 traverse(n->children);
647 static void parse(const char *filename)
651 doc = xmlReadFile(filename, NULL, 0);
653 fprintf(stderr, "Failed to parse '%s'.\n", filename);
658 traverse(xmlDocGetRootElement(doc));
664 static void parse_argument(const char *arg)
666 const char *p = arg+1;
674 fprintf(stderr, "Bad argument '%s'\n", arg);
680 int main(int argc, char **argv)
686 for (i = 1; i < argc; i++) {
687 const char *a = argv[i];