]> git.tdb.fi Git - ext/subsurface.git/blobdiff - uemis.c
Fix profile and average depth for freedives
[ext/subsurface.git] / uemis.c
diff --git a/uemis.c b/uemis.c
index 259d30c53758f9d8ad3e8f2868ef47fa30534efa..06ef5b4231c6f85e6474d182b04a5cf7911d8f2b 100644 (file)
--- a/uemis.c
+++ b/uemis.c
@@ -14,7 +14,6 @@
 #include <stdlib.h>
 #define __USE_XOPEN
 #include <time.h>
-#include <regex.h>
 
 #include "dive.h"
 #include "uemis.h"
@@ -73,64 +72,6 @@ static void decode( uint8_t *inbuf, uint8_t *outbuf, int inbuf_len ) {
 }
 /* end code from Bob Trower */
 
-/* small helper functions */
-/* simpleregex allocates (and reallocates) the found buffer
- * don't forget to free it when you are done
- */
-static int simpleregex(char *buffer, char *regex, char **found) {
-       int status;
-       regex_t re;
-       regmatch_t match[5];
-
-       if (regcomp(&re, regex, 0) !=0) {
-               fprintf(stderr,"internal error, regex failed!\n");
-               exit(1);
-       }
-       status = regexec(&re,buffer,5,match,0);
-       if(status == 0) {
-               *found = realloc(*found,match[1].rm_eo-match[1].rm_so + 1);
-               strncpy(*found,buffer+match[1].rm_so,match[1].rm_eo-match[1].rm_so);
-               (*found)[match[1].rm_eo-match[1].rm_so] = '\0';
-       }
-       return(status == 0);
-}
-
-/* read in line of arbitrary length (important for SDA files that can
- * have lines that are tens of kB long
- * don't forget to free it when you are done
- */
-#define MYGETL_INCR 1024
-static char * mygetline(FILE * f) {
-       size_t size = 0;
-       size_t len  = 0;
-       char * buf  = NULL;
-
-       do {
-               size += MYGETL_INCR;
-               if ((buf = realloc(buf,size)) == NULL)
-                       break;
-               fgets(buf+len,MYGETL_INCR,f);
-               len = strlen(buf);
-       } while (!feof(f) && buf[len-1]!='\n');
-       return buf;
-}
-
-/* text matching, used to build very poor man's XML parser */
-int matchit(FILE *infd, char *regex, char *typeregex, char **found) {
-       char *buffer;
-
-       while (!feof(infd)) {
-               buffer = mygetline(infd);
-               if (buffer && simpleregex(buffer,regex,found)) {
-                       buffer = mygetline(infd);
-                       if (buffer && simpleregex(buffer,typeregex,found)) {
-                               return 1;
-                       }
-               }
-       }
-       return 0;
-}
-
 /*
  * pressure_to_depth: In centibar. And when converting to
  * depth, I'm just going to always use saltwater, because I
@@ -145,9 +86,9 @@ static int pressure_to_depth(uint16_t value)
 {
        double atm, cm;
 
-       atm = (value / 100.0) / 1.01325;
+       atm = bar_to_atm(value / 100.0);
        cm = 100 * atm + 0.5;
-       return( (cm > 0) ? 10 * (long)cm : 0);
+       return (cm > 0) ? 10 * (long)cm : 0;
 }
 
 /*
@@ -177,14 +118,69 @@ bail:
        return datalen;
 }
 
+/* Create events from the flag bits;
+ * These bits basically represent what is displayed on screen at sample time.
+ * Many of these 'warnings' are way hyper-active and seriously clutter the
+ * profile plot - so these are disabled by default
+ */
+void uemis_event(struct dive *dive, struct sample *sample, uemis_sample_t *u_sample)
+{
+       uint8_t *flags = u_sample->flags;
+
+       if (flags[1] & 0x01)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Safety Stop Violation");
+       if (flags[1] & 0x08)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Speed Alarm");
+#if WANT_CRAZY_WARNINGS
+       if (flags[1] & 0x06) /* both bits 1 and 2 are a warning */
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Speed Warning");
+       if (flags[1] & 0x10)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "PO2 Green Warning");
+#endif
+       if (flags[1] & 0x20)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "PO2 Ascend Warning");
+       if (flags[1] & 0x40)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "PO2 Ascend Alarm");
+       /* flags[2] reflects the deco / time bar
+        * flags[3] reflects more display details on deco and pO2 */
+       if (flags[4] & 0x01)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Tank Pressure Info");
+       if (flags[4] & 0x04)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "RGT Warning");
+       if (flags[4] & 0x08)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "RGT Alert");
+       if (flags[4] & 0x40)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Tank Change Suggested");
+       if (flags[4] & 0x80)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Depth Limit Exceeded");
+       if (flags[5] & 0x01)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Max Deco Time Warning");
+       if (flags[5] & 0x04)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Dive Time Info");
+       if (flags[5] & 0x08)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Dive Time Alert");
+       if (flags[5] & 0x10)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Marker");
+       if (flags[6] & 0x02)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "No Tank Data");
+       if (flags[6] & 0x04)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Low Battery Warning");
+       if (flags[6] & 0x08)
+               add_event(dive, sample->time.seconds, 0, 0, 0, "Low Battery Alert");
+       /* flags[7] reflects the little on screen icons that remind of previous
+        * warnings / alerts - not useful for events */
+}
+
 /*
  * parse uemis base64 data blob into struct dive
  */
-void uemis_parse_divelog_binary(char *base64, struct dive **divep) {
+void uemis_parse_divelog_binary(char *base64, void *datap) {
        int datalen;
        int i;
        uint8_t *data;
        struct sample *sample;
+       uemis_sample_t *u_sample;
+       struct dive **divep = datap;
        struct dive *dive = *divep;
        int template, gasoffset;
 
@@ -201,89 +197,44 @@ void uemis_parse_divelog_binary(char *base64, struct dive **divep) {
        gasoffset = template = *(uint8_t *)(data+115);
        if (template == 3)
                gasoffset = 4;
+       if (template == 0)
+               template = 1;
        for (i = 0; i < template; i++) {
                float volume = *(float *)(data+116+25*(gasoffset + i)) * 1000.0;
-               /* uemis always assumes a working pressure of 3000psi / 206bar - even though that's wrong */
-               /* I also think that the unit that it stores (cuft for me) might change with SDA settings */
-               // dive->cylinder[i].type.size.mliter = volume * 206.84 / 28.317;
-               dive->cylinder[i].type.size.mliter = volume * 200 / 28.317;
+               /* uemis always assumes a working pressure of 202.6bar (!?!?) - I first thought
+                * it was 3000psi, but testing against all my dives gets me that strange number.
+                * Still, that's of course completely bogus and shows they don't get how
+                * cylinders are named in non-metric parts of the world...
+                * we store the incorrect working pressure to get the SAC calculations "close"
+                * but the user will have to correct this manually
+                */
+               dive->cylinder[i].type.size.mliter = volume;
+               dive->cylinder[i].type.workingpressure.mbar = 202600;
                dive->cylinder[i].gasmix.o2.permille = *(uint8_t *)(data+120+25*(gasoffset + i)) * 10 + 0.5;
                dive->cylinder[i].gasmix.he.permille = 0;
        }
        /* first byte of divelog data is at offset 0x123 */
        i = 0x123;
-       while ((i < datalen) && (*(uint16_t *)(data+i))) {
+       u_sample = (uemis_sample_t *)(data + i);
+       while ((i < datalen) && (u_sample->dive_time)) {
                /* it seems that a dive_time of 0 indicates the end of the valid readings */
                /* the SDA usually records more samples after the end of the dive --
                 * we want to discard those, but not cut the dive short; sadly the dive
                 * duration in the header is a) in minutes and b) up to 3 minutes short */
-               if (*(uint16_t *)(data+i) > dive->duration.seconds + 180)
+               if (u_sample->dive_time > dive->duration.seconds + 180)
                        break;
                sample = prepare_sample(divep);
                dive = *divep; /* prepare_sample might realloc the dive */
-               sample->time.seconds = *(uint16_t *)(data+i);
-               sample->depth.mm = pressure_to_depth(*(uint16_t *)(data+i+2));
-               sample->temperature.mkelvin = (*(uint16_t *)(data+i+4) * 100) + 273150;
-               sample->cylinderpressure.mbar= *(uint16_t *)(data+i+23) * 10;
-               sample->cylinderindex = *(uint8_t *)(data+i+22);
-               finish_sample(dive, sample);
+               sample->time.seconds = u_sample->dive_time;
+               sample->depth.mm = pressure_to_depth(u_sample->water_pressure);
+               sample->temperature.mkelvin = (u_sample->dive_temperature * 100) + 273150;
+               sample->cylinderindex = u_sample->active_tank;
+               sample->cylinderpressure.mbar= u_sample->tank_pressure * 10;
+               uemis_event(dive, sample, u_sample);
+               finish_sample(dive);
                i += 0x25;
+               u_sample++;
        }
        dive->duration.seconds = sample->time.seconds - 1;
        return;
 }
-
-/* parse a single file
- * TODO: we don't report any errors when the parse fails - we simply don't add them to the list
- */
-void
-parse_uemis_file(char *divelogfilename) {
-       char *found=NULL;
-       struct tm tm;
-       struct dive *dive;
-
-       FILE *divelogfile = fopen(divelogfilename,"r");
-
-       dive = alloc_dive();
-
-       if (! matchit(divelogfile,"val key=\"date\"","<ts>\\([^<]*\\)</ts>",&found)) {
-               /* some error handling */
-               goto bail;
-       }
-       strptime(found, "%Y-%m-%dT%H:%M:%S", &tm);
-       dive->when = utc_mktime(&tm);
-       if (! matchit(divelogfile,"<val key=\"duration\">",
-                       "<float>\\([0-9.]*\\)</float>", &found)) {
-               /* some error handling */
-               goto bail;
-       }
-       dive->duration.seconds = 60.0 * atof(found);
-
-       if (! matchit(divelogfile,"<val key=\"depth\">",
-                       "<int>\\([0-9.]*\\)</int>", &found)) {
-               /* some error handling */
-               goto bail;
-       }
-       dive->maxdepth.mm = pressure_to_depth(atoi(found));
-
-       if (! matchit(divelogfile,"<val key=\"file_content\">",
-                       ">\\([a-zA-Z0-9+/=]*\\)<", &found)) {
-               /* some error handling */
-               goto bail;
-       }
-       uemis_parse_divelog_binary(found,&dive);
-       record_dive(dive);
-bail:
-       if (found)
-               free(found);
-}
-
-/*
- * parse the two files extracted from the SDA
- */
-void
-uemis_import() {
-       if (open_import_file_dialog("*.SDA","uemis Zurich SDA files",
-                                       &parse_uemis_file))
-               report_dives();
-}