/*
 * Driver for Navcom receivers using proprietary NCT messages,
 * a binary protocol.
 *
 * Vendor website: http://www.navcomtech.com/
 * Technical references: Technical Reference Manual P/N 96-3120001-3001
 *
 * Tested with two SF-2040G models
 *
 * At this stage, this driver implements the following commands:
 *
 * 0x20: Data Request (tell the unit which responses you want)
 * 0x3f: LED Configuration (controls the front panel LEDs -- for testing)
 * 0x1c: Test Support Block (again, blinks the front panel lights)
 *
 * and it understands the following responses:
 *
 * 0x06: Acknowledgement (without error)
 * 0x15: Negative Acknowledge
 * 0x86: Channel Status
 * 0xae: Identification Block
 * 0xb0: Raw Meas. Data Block
 * 0xb1: PVT Block
 * 0xb5: Pseudorange Noise Statistics
 * 0xd3: LBM DSP Status Block
 * 0xef: Clock Drift and Offset
 *
 * By Diego Berge. Contact via web form at http://www.navlost.eu/contact
 *
 * Week counters are not limited to 10 bits. It's unknown what
 * the firmware is doing to disambiguate them, if anything; it might just
 * be adding a fixed offset based on a hidden epoch value, in which case
 * unhappy things will occur on the next rollover.
 *
 * This file is Copyright 2010 by the GPSD project
 * SPDX-License-Identifier: BSD-2-clause
 */

#include "../include/gpsd_config.h"  /* must be before all includes */

#include <math.h>
#include <stdbool.h>
#include <stdio.h>
#include <string.h>

#include "../include/gpsd.h"

#if defined(NAVCOM_ENABLE)
#include "../include/bits.h"
#include "../include/timespec.h"

/* Have data which is 24 bits long */
#define getles3224(buf,off)  (int32_t)(((uint32_t)getub((buf), (off)+2)<<24 | (uint32_t)getub((buf), (off)+1)<<16 | (uint32_t)getub((buf), (off))<<8)>>8)
#define getleu3224(buf,off) (uint32_t)(((uint32_t)getub((buf), (off)+2)<<24 | (uint32_t)getub((buf), (off)+1)<<16 | (uint32_t)getub((buf), (off))<<8)>>8)

/* And just to be difficult, Navcom is little endian but the GPS data stream
   is big endian.  Some messages contain raw GPS data */
#define getles16_be(buf, off)   (int16_t)((((uint16_t)getub(buf, (off)) << 8) \
                                    | (uint16_t)getub(buf, (off)+1)))
#define getleu16_be(buf, off)   (uint16_t)((((uint16_t)getub(buf, (off)) << 8) \
                                    | (uint16_t)getub(buf, (off)+1)))
#define getles32_be(buf, off)   (int32_t)((((uint16_t)getleu16_be(buf, (off)) << 16) \
                                    | getleu16_be(buf, (off)+2)))
#define getleu32_be(buf, off)   (uint32_t)((((uint16_t)getleu16_be(buf, (off)) << 16) \
                                    | getleu16_be(buf, (off)+2)))
#define getles64_be(buf, off)   (int64_t)((((uint64_t)getleu32_be(buf, (off)) << 32) \
                                    | getleu32_be(buf, (off)+4)))
#define getleu64_be(buf, off)   (uint64_t)((((uint64_t)getleu32_be(buf, (off)) << 32) \
                                    | getleu32_be(buf, (off)+4)))
#define getles3224_be(buf,off)     (int32_t)(((uint32_t)getub((buf), (off))<<24 \
                                    | (uint32_t)getub((buf), (off)+1)<<16 \
                                    | (uint32_t)getub((buf), (off)+2)<<8)>>8)

#define NAVCOM_CHANNELS 12

static uint8_t checksum(unsigned char *buf, size_t len)
{
    size_t n;
    uint8_t csum = (uint8_t) 0x00;
    for (n = 0; n < len; n++)
        csum ^= buf[n];
    return csum;
}

static bool navcom_send_cmd(struct gps_device_t *session, unsigned char *cmd,
                            size_t len)
{
    return (gpsd_write(session, (const char *)cmd, len) == (ssize_t) len);
}

/* Data Request */
static void navcom_cmd_0x20(struct gps_device_t *session, uint8_t block_id,
                            uint16_t rate)
{
    unsigned char msg[18];
    putbyte(msg, 0, 0x02);
    putbyte(msg, 1, 0x99);
    putbyte(msg, 2, 0x66);
    putbyte(msg, 3, 0x20);      /* Cmd ID */
    putle16(msg, 4, 0x000e);    /* Length */
    putbyte(msg, 6, 0x00);      /* Action */
    putbyte(msg, 7, 0x01);      /* Count of blocks */
    putbyte(msg, 8, block_id);  /* Data Block ID */
    putbyte(msg, 9, 0x02);      /* Logical Ports */
    putle16(msg, 10, rate);     /* Data rate */
    putbyte(msg, 12, 0x71);
    putbyte(msg, 13, 0x00);
    putle16(msg, 14, 0x0000);
    putbyte(msg, 16, checksum(msg + 3, 13));
    putbyte(msg, 17, 0x03);
    (void)navcom_send_cmd(session, msg, 18);
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: sent command 0x20 (Data Request) "
             "- data block id = %02x at rate %02x\n", block_id, rate);
}

// cppcheck-suppress unusedFunction
static void UNUSED navcom_cmd_0x3f(struct gps_device_t *session)
/* Changes the LED settings in the receiver */
{
    unsigned char msg[12];
    putbyte(msg, 0, 0x02);
    putbyte(msg, 1, 0x99);
    putbyte(msg, 2, 0x66);
    putbyte(msg, 3, 0x3f);      /* Cmd ID */
    putle16(msg, 4, 0x0008);
    putbyte(msg, 6, 0x01);      /* Action */
    putbyte(msg, 7, 0x00);      /* Reserved */
    putbyte(msg, 8, 0x02);      /* Link LED setting */
    putbyte(msg, 9, 0x0a);      /* Battery LED setting */
    putbyte(msg, 10, checksum(msg + 3, 7));
    putbyte(msg, 11, 0x03);
    (void)navcom_send_cmd(session, msg, 12);
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: sent command 0x3f (LED Configuration Block)\n");
}

/* Test Support Block - Blinks the LEDs */
static void navcom_cmd_0x1c(struct gps_device_t *session, uint8_t mode,
                            uint8_t length)
{
    unsigned char msg[12];
    putbyte(msg, 0, 0x02);
    putbyte(msg, 1, 0x99);
    putbyte(msg, 2, 0x66);
    putbyte(msg, 3, 0x1c);      /* Cmd ID */
    putle16(msg, 4, 0x0008);
    putbyte(msg, 6, 0x04);      /* Use ACK/NAK */
    putbyte(msg, 7, mode);      /* 0x01 or 0x02 */
    putbyte(msg, 8, length);    /* Only if mode == 0x01 */
    putbyte(msg, 9, 0x00);
    putbyte(msg, 10, checksum(msg + 3, 7));
    putbyte(msg, 11, 0x03);
    (void)navcom_send_cmd(session, msg, 12);
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: sent command 0x1c (Test Support Block)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: command 0x1c mode = %02x, length = %u\n",
             mode, length);
}

/* Serial Port Configuration */
static void navcom_cmd_0x11(struct gps_device_t *session,
                            uint8_t port_selection)
{
    /* NOTE - We only allow changing one port at a time,
     * although the message supports doing both at once. */
    unsigned char msg[12];
    putbyte(msg, 0, 0x02);
    putbyte(msg, 1, 0x99);
    putbyte(msg, 2, 0x66);
    putbyte(msg, 3, 0x11);      /* Cmd ID */
    putle16(msg, 4, 0x0008);    /* Length */
    putbyte(msg, 6, 0x04);      /* Action - Use ACK/NAK) */
    putbyte(msg, 7, port_selection);
    putbyte(msg, 8, 0x00);      /* Reserved */
    putbyte(msg, 9, 0x00);      /* Reserved */
    putbyte(msg, 10, checksum(msg + 3, 7));
    putbyte(msg, 11, 0x03);
    (void)navcom_send_cmd(session, msg, 12);
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: sent command 0x11 (Serial Port Configuration)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: serial port selection: 0x%02x\n", port_selection);
}

static void navcom_event_hook(struct gps_device_t *session, event_t event)
{
    if (session->context->readonly)
        return;

    /* Request the following messages: */
    if (event == event_identified) {
        /* NOTE - Channel Status allows us to know into which of the
         * unit's various serial ports we are connected.
         * Its value gets updated every time we receive a 0x06 (Ack)
         * message.  Note that if commands are being fed into the
         * unit from more than one port (which is entirely possible
         * although not necessarily a bright idea), there is a good
         * chance that we might misidentify our port */
        navcom_cmd_0x1c(session, 0x02, 0);      /* Test Support Block */
        navcom_cmd_0x20(session, 0xae, 0x0000); /* Identification Block */
        navcom_cmd_0x20(session, 0x86, 0x000a); /* Channel Status */
        navcom_cmd_0x1c(session, 0x01, 5);      /* Blink LEDs on receiver */
        navcom_cmd_0x20(session, 0xae, 0x1770); /* Identification Block - send every 10 min */
        navcom_cmd_0x20(session, 0xb1, 0x4000); /* PVT Block */
        navcom_cmd_0x20(session, 0xb5, 0x00c8); /* Pseudorange Noise Statistics - send every 20s */
        navcom_cmd_0x20(session, 0xb0, 0x4000); /* Raw Meas Data Block */
        navcom_cmd_0x20(session, 0x81, 0x0000); /* Packed Ephemeris Data - send once */
        navcom_cmd_0x20(session, 0x81, 0x4000); /* Packed Ephemeris Data */
        navcom_cmd_0x20(session, 0x86, 0x4000); /* Channel Status */
        navcom_cmd_0x20(session, 0x83, 0x4000); /* Ionosphere and UTC Data */
        navcom_cmd_0x20(session, 0xef, 0x0bb8); /* Clock Drift - send every 5 min */
    }
}

/* Ionosphere and UTC Data */
static gps_mask_t handle_0x83(struct gps_device_t *session)
{
    /* NOTE - At the present moment this is only being used
     * for determining the GPS-UTC time difference,
     * for which the iono data is not needed as far
     * as we are concerned.  However, I am still
     * reporting it (if debuglevel >= LOG_DATA) as a
     * matter of interest */
/* 2^-30 */
#define SF_A0 (0.000000000931322574615478515625)
/* 2^-50 */
#define SF_A1 (0.000000000000000888178419700125)
/* 2^12 */
#define SF_TOT (4096)
/* 2^-30 */
#define SF_ALPHA0 (0.000000000931322574615478515625)
/* 2^-27 */
#define SF_ALPHA1 (0.000000007450580596923828125)
/* 2^-24 */
#define SF_ALPHA2 (0.000000059604644775390625)
/* 2^-24 */
#define SF_ALPHA3 (0.000000059604644775390625)
/* 2^11 */
#define SF_BETA0 (2048)
/* 2^14 */
#define SF_BETA1 (16384)
/* 2^16 */
#define SF_BETA2 (65536)
/* 2^16 */
#define SF_BETA3 (65536)
    unsigned char *buf = session->lexer.outbuffer + 3;
    uint16_t week = getleu16(buf, 3);
    uint32_t tow = getleu32(buf, 5);
    int8_t alpha0 = getsb(buf, 9);
    int8_t alpha1 = getsb(buf, 10);
    int8_t alpha2 = getsb(buf, 11);
    int8_t alpha3 = getsb(buf, 12);
    int8_t beta0 = getsb(buf, 13);
    int8_t beta1 = getsb(buf, 14);
    int8_t beta2 = getsb(buf, 15);
    int8_t beta3 = getsb(buf, 16);
    int32_t a1 = getles32(buf, 17);
    int32_t a0 = getles32(buf, 21);
    uint8_t tot = getub(buf, 25);
    uint8_t wnt = getub(buf, 26);
    int8_t dtls = getsb(buf, 27);
    uint8_t wnlsf = getub(buf, 28);
    uint8_t dn = getub(buf, 29);
    int8_t dtlsf = getsb(buf, 30);

    /* Ref.: ICD-GPS-200C 20.3.3.5.2.4 */
    if ((week % 256) * 604800 + tow / 1000.0 < wnlsf * 604800 + dn * 86400) {
        /* Effectivity time is in the future, use dtls */
        session->context->leap_seconds = (int)dtls;
    } else {
        /* Effectivity time is not in the future, use dtlsf */
        session->context->leap_seconds = (int)dtlsf;
    }

    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0x83 (Ionosphere and UTC Data)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: Scaled parameters follow:\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: GPS Week: %u, GPS Time of Week: %u (GPS Time: %f)\n",
             week, tow, week * 604800 + tow / 1000.0);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: a0: %12.4E, a1: %12.4E, a2: %12.4E, a3: %12.4E, "
             "b0: %12.4E, b1: %12.4E, b2: %12.4E, b3: %12.4E\n",
             (double)alpha0 * SF_ALPHA0, (double)alpha1 * SF_ALPHA1,
             (double)alpha2 * SF_ALPHA2, (double)alpha3 * SF_ALPHA3,
             (double)beta0 * SF_BETA0, (double)beta1 * SF_BETA1,
             (double)beta2 * SF_BETA2, (double)beta3 * SF_BETA3);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: A0: %19.12E, A1: %19.12E\n", (double)a0 * SF_A0,
             (double)a1 * SF_A1);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: UTC Ref. Time: %lu, UTC Ref. Week: %u, dTls: %d\n",
             (unsigned long)tot * SF_TOT, wnt, dtls);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: Week of leap seconds: %u, Day number of leap seconds: %u, dTlsf: %d\n",
             wnlsf, dn, dtlsf);

    return 0;                   /* No flag for update of leap seconds (Not part of a fix) */

#undef SF_A0
#undef SF_A1
#undef SF_TOT
#undef SF_ALPHA0
#undef SF_ALPHA1
#undef SF_ALPHA2
#undef SF_ALPHA3
#undef SF_BETA0
#undef SF_BETA1
#undef SF_BETA2
#undef SF_BETA3
}

/* Acknowledgement (without error) */
static gps_mask_t handle_0x06(struct gps_device_t *session)
{
    unsigned char *buf = session->lexer.outbuffer + 3;
    uint8_t cmd_id = getub(buf, 3);
    uint8_t port = getub(buf, 4);
    session->driver.navcom.physical_port = port;        /* This tells us which serial port was used last */
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0x06 (Acknowledgement (without error))\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: acknowledged command id 0x%02x on port %c\n",
             cmd_id, (port == 0 ? 'A' : (port == 1 ? 'B' : '?')));
    return 0;                   /* Nothing updated */
}

/* Negative Acknowledge */
static gps_mask_t handle_0x15(struct gps_device_t *session)
{
    unsigned n;
    unsigned char *buf = session->lexer.outbuffer + 3;
    unsigned msg_len = getleu16(buf, 1);
    uint8_t port, cmd_id = getub(buf, 3);

    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0x15 (Negative Acknowledge)\n");

    if (18 < msg_len) {
        GPSD_LOG(LOG_PROG, &session->context->errout,
                 "Navcom: 0x15 too long %u\n", msg_len);
        return 0;
    }
    for (n = 4; n < (msg_len - 2); n += 2) {
        uint8_t err_id = getub(buf, n);
        uint8_t err_desc = getub(buf, n + 1);
        GPSD_LOG(LOG_DATA, &session->context->errout,
                 "Navcom: error id = 0x%02x, error description = 0x%02x\n",
                 err_id, err_desc);
    }
    port = getub(buf, n);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: negative acknowledge was for command id 0x%02x "
             "on port %c\n",
             cmd_id, (port == 0 ? 'A' : (port == 1 ? 'B' : '?')));
    return 0;                   /* Nothing updated */
}

/* PVT Block */
static gps_mask_t handle_0xb1(struct gps_device_t *session)
{
    gps_mask_t mask = 0;
    unsigned char *buf = session->lexer.outbuffer + 3;
    uint16_t week;
    uint32_t tow;
    timespec_t ts_tow = {0, 0};      // pacify codacy
    int32_t lat, lon;
    char ts_buf[TIMESPEC_LEN];
    /* Resolution of lat/lon values (2^-11) */
#define LL_RES (0.00048828125)
    uint8_t lat_fraction, lon_fraction;
    /* Resolution of lat/lon fractions (2^-15) */
#define LL_FRAC_RES (0.000030517578125)
    uint8_t nav_mode;
    int32_t ellips_height, altitude;
    /* Resolution of height and altitude values (2.0^-10) */
#define EL_RES (0.0009765625)
    double vel_north, vel_east, vel_up;
    uint8_t gdop, pdop, hdop, vdop, tdop;
    /* This value means "undefined" */
#define DOP_UNDEFINED (255)

    int16_t ant_height_adj;
    int32_t set_delta_up;
    /* Resolution of delta north, east, and up,
     * and ant. height adjustment values (1mm) */
#define D_RES (0.001)

#ifdef __UNUSED__
    /* Other values provided by the PVT block which we
     * may want to provide in the future.
     */
    uint8_t dgps_conf;
    uint16_t max_dgps_age;
    uint8_t ext_nav_mode;
    int32_t set_delta_north, set_delta_east;
    uint8_t nav_failure_code;
#endif /* __UNUSED__ */

    /* Timestamp */
    week = (uint16_t) getleu16(buf, 3);
    tow = (uint32_t) getleu32(buf, 5);        /* tow in ms */
    MSTOTS(&ts_tow, tow);
    session->newdata.time = gpsd_gpstime_resolv(session, week, ts_tow);

    /* Get latitude, longitude */
    lat = getles32(buf, 13);
    lon = getles32(buf, 17);
    lat_fraction = (uint8_t) (getub(buf, 21) >> 4);
    lon_fraction = (uint8_t) (getub(buf, 21) & 0x0f);

    session->newdata.latitude =
        (double)(lat * LL_RES + lat_fraction * LL_FRAC_RES) / 3600;
    session->newdata.longitude =
        (double)(lon * LL_RES + lon_fraction * LL_FRAC_RES) / 3600;

    /* Nav mode */
    nav_mode = (uint8_t) getub(buf, 22);
    if (-nav_mode & 0x80) {
        session->newdata.status = STATUS_UNK;
        session->newdata.mode = MODE_NO_FIX;
    } else {
        session->newdata.mode = ((nav_mode & 0x40)!=0 ? MODE_3D : MODE_2D);
        session->newdata.status =
            ((nav_mode & 0x03)!=0 ? STATUS_DGPS : STATUS_GPS);
    }

    /* altHAE */
    ellips_height = getles32(buf, 23);
    /* altMSL */
    altitude = getles32(buf, 27);

    ant_height_adj = getles16(buf, 51);
    set_delta_up = getles32(buf, 79);

    session->newdata.altMSL = (double)(altitude * EL_RES)
        + (ant_height_adj * D_RES) + (set_delta_up * D_RES);
    session->newdata.altHAE = (double)(ellips_height) * EL_RES
        + (ant_height_adj * D_RES) + (set_delta_up * D_RES);
    /* Let gpsd_error_model() deal with geoid_sep */

    /* Speed Data */
    vel_north = (double)getles3224(buf, 31);
    vel_east = (double)getles3224(buf, 34);
    vel_up = (double)getles3224(buf, 37);

    session->newdata.NED.velN = vel_north * 0.1;
    session->newdata.NED.velE = vel_east * 0.1;
    session->newdata.NED.velD = -vel_up * 0.1;

    /* Quality indicators */
    /* UNUSED fom = getub(buf, 40);     * FOM is DRMS */
    gdop = getub(buf, 41);
    pdop = getub(buf, 42);
    hdop = getub(buf, 43);
    vdop = getub(buf, 44);
    tdop = getub(buf, 45);
    /* UNUSED tfom = getub(buf, 46);    * tfom == 10 * TDOP */

    /* let gpsd_error_model() do the error estimates */

    if (gdop != DOP_UNDEFINED) {
        session->gpsdata.dop.gdop = gdop / 10.0;
        mask |= DOP_SET;
    }
    if (pdop != DOP_UNDEFINED) {
        session->gpsdata.dop.pdop = pdop / 10.0;
        mask |= DOP_SET;
    }
    if (hdop != DOP_UNDEFINED) {
        session->gpsdata.dop.hdop = hdop / 10.0;
        mask |= DOP_SET;
    }
    if (vdop != DOP_UNDEFINED) {
        session->gpsdata.dop.vdop = vdop / 10.0;
        mask |= DOP_SET;
    }
    if (tdop != DOP_UNDEFINED) {
        session->gpsdata.dop.tdop = tdop / 10.0;
        mask |= DOP_SET;
    }

    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0xb1 (PVT Report)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: navigation mode %s (0x%02x) - %s - %s\n",
             ((-nav_mode & 0x80)!='\0' ? "invalid" : "valid"), nav_mode,
             ((nav_mode & 0x40)!='\0' ? "3D" : "2D"),
             ((nav_mode & 0x03)!='\0' ? "DGPS" : "GPS"));
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: velocities: north = %f east = %f up = %f\n",
             session->newdata.NED.velN,
             session->newdata.NED.velE,
             -session->newdata.NED.velD);
#undef D_RES
#undef LL_RES
#undef LL_FRAC_RES
#undef EL_RES
#undef VEL_RES
#undef DOP_UNDEFINED

    mask |= LATLON_SET | ALTITUDE_SET | STATUS_SET | MODE_SET | USED_IS |
           HERR_SET | TIMERR_SET | VNED_SET | TIME_SET | NTPTIME_IS;
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "PVT 0xb1: time=%s, lat=%.2f lon=%.2f altHAE=%.2f "
             "altMSL %.2f mode=%d status=%d gdop=%.2f pdop=%.2f hdop=%.2f "
             "vdop=%.2f tdop=%.2f mask={%s}\n",
             timespec_str(&session->newdata.time, ts_buf, sizeof(ts_buf)),
             session->newdata.latitude,
             session->newdata.longitude,
             session->newdata.altHAE,
             session->newdata.altMSL,
             session->newdata.mode,
             session->newdata.status,
             session->gpsdata.dop.gdop,
             session->gpsdata.dop.pdop,
             session->gpsdata.dop.hdop,
             session->gpsdata.dop.vdop, session->gpsdata.dop.tdop,
             gps_maskdump(mask));
    return mask;
}

/* Packed Ephemeris Data */
static gps_mask_t handle_0x81(struct gps_device_t *session)
{
    /* Scale factors for everything */
    /* 2^-31 */
#define SF_TGD       (.000000000465661287307739257812)
    /* 2^4 */
#define SF_TOC     (16)
    /* 2^-55 */
#define SF_AF2       (.000000000000000027755575615628)
    /* 2^-43 */
#define SF_AF1       (.000000000000113686837721616029)
    /* 2^-31 */
#define SF_AF0       (.000000000465661287307739257812)
    /* 2^-5 */
#define SF_CRS       (.031250000000000000000000000000)
    /* 2^-43 */
#define SF_DELTA_N   (.000000000000113686837721616029)
    /* 2^-31 */
#define SF_M0   (.000000000465661287307739257812)
    /* 2^-29 */
#define SF_CUC       (.000000001862645149230957031250)
    /* 2^-33 */
#define SF_E     (.000000000116415321826934814453)
    /* 2^-29 */
#define SF_CUS       (.000000001862645149230957031250)
    /* 2^-19 */
#define SF_SQRT_A    (.000001907348632812500000000000)
    /* 2^4 */
#define SF_TOE     (16)
    /* 2^-29 */
#define SF_CIC       (.000000001862645149230957031250)
    /* 2^-31 */
#define SF_OMEGA0    (.000000000465661287307739257812)
    /* 2^-29 */
#define SF_CIS       (.000000001862645149230957031250)
    /* 2^-31 */
#define SF_I0   (.000000000465661287307739257812)
    /* 2^-5 */
#define SF_CRC       (.031250000000000000000000000000)
    /* 2^-31 */
#define SF_OMEGA     (.000000000465661287307739257812)
    /* 2^-43 */
#define SF_OMEGADOT  (.000000000000113686837721616029)
    /* 2^-43 */
#define SF_IDOT      (.000000000000113686837721616029)

    char ts_buf[TIMESPEC_LEN];
    unsigned char *buf = session->lexer.outbuffer + 3;
    uint8_t prn = getub(buf, 3);
    uint16_t week = getleu16(buf, 4);
    uint32_t tow = getleu32(buf, 6);
    uint16_t iodc = getleu16(buf, 10);
    /* And now the fun starts... everything that follows is
     * raw GPS data minus parity */
    /* Subframe 1, words 3 to 10 minus parity */
    uint16_t wn = (getleu16_be(buf, 12) & 0xffc0) >> 6;
    uint8_t cl2 = (getub(buf, 13) & 0x30) >> 4;
    uint8_t ura = getub(buf, 13) & 0x0f;
    uint8_t svh = (getub(buf, 14) & 0xfc) >> 2;
    /* We already have IODC from earlier in the message, so
     * we do not decode again */
/*    uint16_t iodc = (getub(buf, 14)&0x03)<<8;*/
    uint8_t l2pd = (getub(buf, 15) & 0x80) >> 7;
    int8_t tgd = getsb(buf, 26);
/*    iodc |= getub(buf, 27);*/
    uint16_t toc = getleu16_be(buf, 28);
    int8_t af2 = getsb(buf, 30);
    int16_t af1 = getles16_be(buf, 31);
    int32_t af0 = getles3224_be(buf, 33) >> 2;
    /* Subframe 2, words 3 to 10 minus parity */
    uint8_t iode = getub(buf, 36);
    int16_t crs = getles16_be(buf, 37);
    int16_t delta_n = getles16_be(buf, 39);
    int32_t m0 = getles32_be(buf, 41);
    int16_t cuc = getles16_be(buf, 45);
    uint32_t e = getleu32_be(buf, 47);
    int16_t cus = getles16_be(buf, 51);
    uint32_t sqrt_a = getleu32_be(buf, 53);
    uint16_t toe = getleu16_be(buf, 57);
    /* NOTE - Fit interval & AODO not collected */
    /* Subframe 3, words 3 to 10 minus parity */
    int16_t cic = getles16_be(buf, 60);
    int32_t Omega0 = getles32_be(buf, 62);
    int16_t cis = getles16_be(buf, 66);
    int32_t i0 = getles32_be(buf, 68);
    int16_t crc = getles16_be(buf, 72);
    int32_t omega = getles32_be(buf, 74);
    int32_t Omegadot = getles3224_be(buf, 78);
    /* Question: What is the proper way of shifting a signed int 2 bits to
     * the right, preserving sign? Answer: integer division by 4. */
    int16_t idot = (int16_t) (((getles16_be(buf, 82) & 0xfffc) / 4) |
                              ((getub(buf, 82) & 80) ? 0xc000 : 0x0000));
    session->context->gps_week = (unsigned short)wn;
    DTOTS(&session->context->gps_tow, (double)(toc * SF_TOC));
    /* leap second? */
    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0x81 (Packed Ephemeris Data)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: PRN: %u, Week: %u, TOW: %s "
             "SV clock bias/drift/drift rate: %#19.12E/%#19.12E/%#19.12E\n",
             prn,
             session->context->gps_week,
             timespec_str(&session->context->gps_tow, ts_buf, sizeof(ts_buf)),
             ((double)af0) * SF_AF0,
             ((double)af1) * SF_AF1, ((double)af2) * SF_AF2);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: IODE (!AODE): %u Crs: %19.12e, Delta n: %19.12e, M0: %19.12e\n",
             iode, (double)crs * SF_CRS,
             (double)delta_n * SF_DELTA_N * GPS_PI,
             (double)m0 * SF_M0 * GPS_PI);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: Cuc: %19.12e, Eccentricity: %19.12e, Cus: %19.12e, A^1/2: %19.12e\n",
             (double)cuc * SF_CUC, (double)e * SF_E, (double)cus * SF_CUS,
             (double)sqrt_a * SF_SQRT_A);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: TOE: %u, Cic: %19.12e, Omega %19.12e, Cis: %19.12e\n",
             toe * SF_TOE, (double)cic * SF_CIC,
             (double)Omega0 * SF_OMEGA0 * GPS_PI, (double)cis * SF_CIS);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: i0: %19.12e, Crc: %19.12e, omega: %19.12e, Omega dot: %19.12e\n",
             (double)i0 * SF_I0 * GPS_PI, (double)crc * SF_CRC,
             (double)omega * SF_OMEGA * GPS_PI,
             (double)Omegadot * SF_OMEGADOT * GPS_PI);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: IDOT: %19.12e, Codes on L2: 0x%x, GPS Week: %u, L2 P data flag: %x\n",
             (double)idot * SF_IDOT * GPS_PI, cl2,
             week - (week % 1024) + wn, l2pd);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: SV accuracy: 0x%x, SV health: 0x%x, TGD: %f, IODC (!AODC): %u\n",
             ura, svh, (double)tgd * SF_TGD, iodc);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: Transmission time: %u\n", tow);

#undef SF_TGD
#undef SF_TOC
#undef SF_AF2
#undef SF_AF1
#undef SF_AF0
#undef SF_CRS
#undef SF_DELTA_N
#undef SF_M0
#undef SF_CUC
#undef SF_E
#undef SF_CUS
#undef SF_SQRT_A
#undef SF_TOE
#undef SF_CIC
#undef SF_OMEGA0
#undef SF_CIS
#undef SF_I0
#undef SF_CRC
#undef SF_OMEGA
#undef SF_OMEGADOT
#undef SF_IDOT

    return 0;
}

/* Channel Status */
static gps_mask_t handle_0x86(struct gps_device_t *session)
{
    size_t n, i, nsu;
    unsigned char *buf = session->lexer.outbuffer + 3;
    unsigned msg_len = getleu16(buf, 1);
    unsigned short week = getleu16(buf, 3);
    uint32_t tow = getleu32(buf, 5);
    uint8_t eng_status = getub(buf, 9);
    uint16_t sol_status = getleu16(buf, 10);
    uint8_t sats_visible = getub(buf, 12);
    //uint8_t sats_tracked = getub(buf, 13);
    //uint8_t used_sats = getub(buf, 14);
    //uint8_t pdop = getub(buf, 15);
    timespec_t ts_tow;

    MSTOTS(&ts_tow, tow);

    /* Timestamp */
    session->gpsdata.skyview_time = gpsd_gpstime_resolv(session, week, ts_tow);

    /* Give this driver a single point of truth about DOPs */
    //session->gpsdata.dop.pdop = (int)pdop / 10.0;

    /* Satellite count */
    session->gpsdata.satellites_visible = (int)sats_visible;

    /* Fix mode */
    switch (sol_status & 0x05) {
    case 0x05:
        session->newdata.status = STATUS_DGPS;
        break;
    case 0x01:
        session->newdata.status = STATUS_GPS;
        break;
    default:
        session->newdata.status = STATUS_UNK;
    }

    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: engine status 0x%x almanac %s time 0x%x pos 0x%x\n",
             eng_status & 0x07, ((eng_status & 0x08) ? "valid" : "invalid"),
             eng_status & 0x30 >> 4, eng_status & 0xc0 >> 6);

    /* Satellite details */
    i = nsu = 0;
    if (298 < msg_len) {
        /* pasify coverity
         * msg_len = 18 + (14 * nsat)
         * assume 20 sats max */
        msg_len = 298;
    }
    for (n = 17; n < msg_len; n += 14) {
        uint8_t prn, ele, ca_snr, p2_snr, log_channel, hw_channel, s, stat;
        uint16_t azm, dgps_age;
        if (i >= MAXCHANNELS) {
            GPSD_LOG(LOG_ERROR, &session->context->errout,
                     "Navcom: packet type 0x86: too many satellites!\n");
            gpsd_zero_satellites(&session->gpsdata);
            return 0;
        }
        prn = getub(buf, n);
        /*
         * This field is described in the Technical Reference as follows:
         *
         * Channel Tracking Status:
         * B0-B1: C/A tracking status
         * B2-B3: P1 tracking status
         * B4-B5: P2 tracking status
         *    00 Acquisition or reacquisition
         *    01 Code loop locked
         *    02 Costas loop locked
         *    11 Full tracking with aiding and active
         *       multipath reduction - all data is valid
         * B6=1: C/A Bit sync
         * B7=1: C/A Frame sync
         *
         * By observation, the satellite is in use if this status is 0xff.
         * But errors here are not very serious, all they can affect is
         * the coverance-matrix calculation for error modeling.
         */
        stat = getub(buf, n + 1);
        log_channel = getub(buf, n + 2);
        ele = getub(buf, n + 5);
        azm = getleu16(buf, n + 6);
        ca_snr = getub(buf, n + 8);
        p2_snr = getub(buf, n + 10);
        dgps_age = getleu16(buf, n + 11);
        hw_channel = getub(buf, n + 13);
        s = (unsigned char)0;
        /* NOTE - In theory, I think one would check for hw channel number to
         * see if one is dealing with a GPS or other satellite, but the
         * channel numbers reported bear no resemblance to what the spec
         * says should be.  So I check for the fact that if all three
         * values below are zero, one is not interested in this satellite */
        if (!(ele == 0 && azm == 0 && dgps_age == 0)) {
            session->gpsdata.skyview[i].PRN = (short)prn;
            session->gpsdata.skyview[i].elevation = (double)ele;
            session->gpsdata.skyview[i].azimuth = (double)azm;
            s = session->gpsdata.skyview[i].ss =
                (p2_snr ? p2_snr : ca_snr) / 4.0;
            session->gpsdata.skyview[i++].used = (stat == 0xff);
            if (stat == 0xff)
                nsu++;
        }
        session->gpsdata.satellites_used = (int)nsu;
        GPSD_LOG(LOG_DATA, &session->context->errout,
                 "Navcom: prn = %3u, ele = %02u, azm = %03u, snr = %d (%s), "
                 "dgps age = %.1fs, log ch = %d, hw ch = 0x%02x\n",
                 prn, ele, azm, s, (p2_snr ? "P2" : "C/A"),
                 (double)dgps_age * 0.1, log_channel & 0x3f, hw_channel);
        GPSD_LOG(LOG_DATA, &session->context->errout,
                 "Navcom:           sol. valid = %c, clock = %s, pos. = %s, "
                 "height = %s, err. code = 0x%x\n",
                 ((sol_status & 0x01) ? 'Y' : 'N'),
                 ((sol_status & 0x02) ? "stable" : "unstable"),
                 ((sol_status & 0x04) ? "dgps" : "unaided"),
                 ((sol_status & 0x08) ? "solved" : "constrained"),
                 ((sol_status & 0x01) ? 0x00 : sol_status & 0x0f00 >> 8));
    }

    GPSD_LOG(LOG_DATA, &session->context->errout,
             "CS 0x86: visible=%d, used=%d, mask={SATELLITE|STATUS}\n",
             session->gpsdata.satellites_visible,
             session->gpsdata.satellites_used);
    return SATELLITE_SET | STATUS_SET;
}

/* Raw Meas. Data Block
 * Size 4 + 8 + (16 * numSat) = 524
 */
static gps_mask_t handle_0xb0(struct gps_device_t *session)
{
    char ts_buf[TIMESPEC_LEN];
    /* L1 wavelength (299792458m/s / 1575420000Hz)
     * from their Technical reference Manual */
#define LAMBDA_L1 (299792458.0 / 1575420000.0)
    unsigned n;
    unsigned char *buf = session->lexer.outbuffer + 3;
    unsigned msg_len = getleu16(buf, 1);
    uint16_t week = getleu16(buf, 3);
    uint32_t tow = getleu32(buf, 5);
    uint8_t tm_slew_acc = getub(buf, 9);
    uint8_t status = getub(buf, 10);

    session->context->gps_week = (unsigned short)week;
    MSTOTS(&session->context->gps_tow, tow);

    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0xb0 (Raw Meas. Data Block)\n");
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: week = %u, tow = %s "
             "time slew accumulator = %u (1/1023mS), status = 0x%02x "
             "(%sclock %s - %u blocks follow)\n",
             session->context->gps_week,
             timespec_str(&session->context->gps_tow, ts_buf, sizeof(ts_buf)),
             tm_slew_acc, status,
             ((status & 0x80) ? "channel time set - " : ""),
             ((status & 0x40) ? "stable" : "not stable"), status & 0x0f);

    if (530 < msg_len) {
        GPSD_LOG(LOG_PROG, &session->context->errout,
                 "Navcom: received packet type 0xb0, length %u too long\n",
                 msg_len);
        return 0;
    }

    for (n = 11; n < (msg_len - 1); n += 16) {
        uint8_t sv_status = getub(buf, n);
        uint8_t ch_status = getub(buf, n + 1);
        uint32_t ca_pseudorange = getleu32(buf, n + 2);
        /* integer division by 16 is a sign-preserving right shift of 4 bits */
        int32_t l1_phase = getles3224(buf, n + 6) / 16;
        uint8_t l1_slips = (uint8_t) (getles3224(buf, n + 6) & 0x0f);
        int16_t p1_ca_pseudorange = getles16(buf, n + 9);
        int16_t p2_ca_pseudorange = getles16(buf, n + 11);
        int32_t l2_phase = getles3224(buf, n + 13) / 16;
        uint8_t l2_slips = (uint8_t) (getles3224(buf, n + 13) & 0x0f);
        double c1 =
            ((sv_status & 0x80) ? (double)ca_pseudorange / 16.0 *
             LAMBDA_L1 : NAN);
        double l1 =
            ((sv_status & 0x80) ? (double)ca_pseudorange / 16.0 +
             (double)l1_phase / 256.0 : NAN);
        double l2 =
            ((sv_status & 0x20)
             ? ((double)ca_pseudorange / 16.0 +
                (double)p2_ca_pseudorange / 16.0) * (120.0 / 154.0)
             + (double)l2_phase / 256.0 : NAN);
        double p1 =
            ((sv_status & 0x40) ? c1 +
             (double)p1_ca_pseudorange / 16.0 * LAMBDA_L1 : NAN);
        double p2 =
            ((sv_status & 0x20) ? c1 +
             (double)p2_ca_pseudorange / 16.0 * LAMBDA_L1 : NAN);
        GPSD_LOG(LOG_SPIN, &session->context->errout,
                 "Navcom: >> sv status = 0x%02x (PRN %u - C/A & L1 %s - P1 %s - P2 & L2 %s)\n",
                 sv_status, (sv_status & 0x1f),
                 ((sv_status & 0x80) ? "valid" : "invalid"),
                 ((sv_status & 0x40) ? "valid" : "invalid"),
                 ((sv_status & 0x20) ? "valid" : "invalid"));
        GPSD_LOG(LOG_SPIN, &session->context->errout,
                 "Navcom: >>> ch status = 0x%02x "
                 "(Logical channel: %u - CA C/No: %u dBHz) "
                 "sL1: %u, sL2: %u\n", ch_status, ch_status & 0x0f,
                 ((ch_status & 0xf0) >> 4) + 35, l1_slips, l2_slips);
        GPSD_LOG(LOG_SPIN, &session->context->errout,
                 "Navcom: >>> C1: %14.3f, L1: %14.3f, L2: %14.3f, P1: %14.3f, P2: %14.3f\n",
                 c1, l1, l2, p1, p2);
    }
#undef LAMBDA_L1
    return 0;                   /* Raw measurements not yet implemented in gpsd */
}

/* Pseudorange Noise Statistics */
static gps_mask_t handle_0xb5(struct gps_device_t *session)
{
    if (sizeof(double) == 8) {
        gps_mask_t mask = TIME_SET;
        char *buf = (char *)session->lexer.outbuffer + 3;
        uint16_t week = getleu16(buf, 3);
        uint32_t tow = getleu32(buf, 5);
        timespec_t ts_tow;
#ifdef __UNUSED__
        double rms = getled64(buf, 9);
        /* Reason why it's unused is these figures do not agree
         * with those obtained from the PVT report (handle_0xb1).
         * The figures from 0xb1 do agree with the values reported
         * by Navcom's PC utility */
        /* let gpsd_error_model() handle this */
        //double ellips_maj = getled64(buf, 17);
        //double ellips_min = getled64(buf, 25);
        //double ellips_azm = getled64(buf, 33);
        double lat_sd = getled64(buf, 41);
        double lon_sd = getled64(buf, 49);
        double alt_sd = getled64(buf, 57);
        double hrms = sqrt(pow(lat_sd, 2) + pow(lon_sd, 2));
        /* Navcom doc unclear, this is likely sep? */
        session->newdata.sep = rms * 1.96;
        session->newdata.eph = hrms * 1.96;
        session->newdata.epv = alt_sd * 1.96;
        mask |= HERR_SET;
#endif /*  __UNUSED__ */
        MSTOTS(&ts_tow, tow);
        session->newdata.time = gpsd_gpstime_resolv(session,
                                                  (unsigned short)week,
                                                  ts_tow);
        GPSD_LOG(LOG_PROG, &session->context->errout,
                 "Navcom: received packet type 0xb5 (Pseudorange Noise Statistics)\n");
        GPSD_LOG(LOG_DATA, &session->context->errout,
                 "Navcom: sep = %f\n", session->newdata.sep);
        return mask;
    } else {
        /* Ignore this message block */
        if (!session->driver.navcom.warned) {
            GPSD_LOG(LOG_WARN, &session->context->errout,
                     "Navcom: received packet type 0xb5 (Pseudorange "
                     "Noise Statistics) ignored "
                     " - sizeof(double) == 64 bits required\n");
            session->driver.navcom.warned = true;
        }
        return 0;               /* Block ignored - wrong sizeof(double) */
    }
}

/* LBM DSP Status Block */
static gps_mask_t handle_0xd3(struct gps_device_t *session UNUSED)
{
    /* This block contains status information about the
     * unit's L-band (Inmarsat) module.  There is nothing
     * interesting in it for our purposes so we do not deal
     * with it.  This callback is purely to a) stop
     * "unrecognised packet" messages appearing in the log
     * and b) explain what it is for the curious */
    return 0;                   /* Nothing updated */
}

/* Identification Block */
static gps_mask_t handle_0xae(struct gps_device_t *session)
{
    char *engconfstr, *asicstr;
    unsigned char *buf = session->lexer.outbuffer + 3;
    size_t msg_len = (size_t) getleu16(buf, 1);
    uint8_t engconf = getub(buf, 3);
    uint8_t asic = getub(buf, 4);
    uint8_t swvermaj = getub(buf, 5);
    uint8_t swvermin = getub(buf, 6);
    uint16_t dcser = getleu16(buf, 7);
    uint8_t dcclass = getub(buf, 9);
    uint16_t rfcser = getleu16(buf, 10);
    uint8_t rfcclass = getub(buf, 12);
    uint8_t softtm[17] = "\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0";
    uint8_t bootstr[17] = "\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0";
    uint8_t ioptm[17] = "\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0\0";
    uint8_t iopvermaj = (uint8_t) 0x00;
    uint8_t iopvermin = (uint8_t) 0x00;
    uint8_t picver = (uint8_t) 0x00;
    uint8_t slsbn = (uint8_t) 0x00;
    uint8_t iopsbn = (uint8_t) 0x00;
    memcpy(softtm, &buf[13], 16);
    memcpy(bootstr, &buf[29], 16);
    if (msg_len == 0x0037) {    /* No IOP */
        slsbn = getub(buf, 53);
    } else {                    /* IOP Present */
        iopvermaj = getub(buf, 53);
        iopvermin = getub(buf, 54);
        memcpy(ioptm, &buf[55], 16);
        picver = getub(buf, 71);
        slsbn = getub(buf, 72);
        iopsbn = getub(buf, 73);
    }

    switch (engconf) {
    case 0x00:
        engconfstr = "Unknown/Undefined";
        break;
    case 0x01:
        engconfstr = "NCT 2000 S";
        break;
    case 0x02:
        engconfstr = "NCT 2000 D";
        break;
    case 0x03:
        engconfstr = "Startfire Single";
        break;
    case 0x04:
        engconfstr = "Starfire Dual";
        break;
    case 0x05:
        engconfstr = "Pole Mount RTK (Internal Radio)";
        break;
    case 0x06:
        engconfstr = "Pole Mount GIS (LBM)";
        break;
    case 0x07:
        engconfstr = "Black Box RTK (Internal Radio)";
        break;
    case 0x08:
        engconfstr = "Black Box GIS (LBM)";
        break;
    case 0x80:
        engconfstr = "R100";
        break;
    case 0x81:
        engconfstr = "R200";
        break;
    case 0x82:
        engconfstr = "R210";
        break;
    case 0x83:
        engconfstr = "R300";
        break;
    case 0x84:
        engconfstr = "R310";
        break;
    default:
        engconfstr = "?";
    }

    switch (asic) {
    case 0x01:
        asicstr = "A-ASIC";
        break;
    case 0x02:
        asicstr = "B-ASIC";
        break;
    case 0x03:
        asicstr = "C-ASIC";
        break;
    case 0x04:
        asicstr = "M-ASIC";
        break;
    default:
        asicstr = "?";
    }

    GPSD_LOG(LOG_PROG, &session->context->errout,
             "Navcom: received packet type 0xae (Identification Block)\n");
    if (msg_len == 0x0037) {
        GPSD_LOG(LOG_INF, &session->context->errout, "Navcom: ID Data: "
                 "%s %s Ver. %u.%u.%u, DC S/N: %u.%u, RF S/N: %u.%u, "
                 "Build ID: %s, Boot software: %s\n",
                 engconfstr, asicstr, swvermaj, swvermin, slsbn, dcser,
                 dcclass, rfcser, rfcclass, softtm, bootstr);
    } else {
        GPSD_LOG(LOG_INF, &session->context->errout, "Navcom: ID Data: "
                 "%s %s Ver. %u.%u.%u, DC S/N: %u.%u, RF S/N: %u.%u, "
                 "Build ID: %s, Boot software: %s, "
                 "IOP Ver.: %u.%u.%u, PIC: %u, IOP Build ID: %s\n",
                 engconfstr, asicstr, swvermaj, swvermin, slsbn, dcser,
                 dcclass, rfcser, rfcclass, softtm, bootstr, iopvermaj,
                 iopvermin, iopsbn, picver, ioptm);
    }

    (void)snprintf(session->subtype, sizeof(session->subtype),
                   "%s %s Ver. %u.%u.%u S/N %u.%u %u.%u",
                   engconfstr, asicstr, swvermaj, swvermin, slsbn, dcser,
                   dcclass, rfcser, rfcclass);
    return DEVICEID_SET;
}

/* Clock Drift and Offset */
static gps_mask_t handle_0xef(struct gps_device_t *session)
{
    unsigned char *buf = session->lexer.outbuffer + 3;
    //uint16_t week = getleu16(buf, 3);
    //uint32_t tow = getleu32(buf, 5);
    int8_t osc_temp = getsb(buf, 9);
    uint8_t nav_status = getub(buf, 10);
    double nav_clock_offset;
    float nav_clock_drift;
    float osc_filter_drift_est;
    char ts_buf[TIMESPEC_LEN];
    int32_t time_slew = (int32_t) getles32(buf, 27);
    if (sizeof(double) == 8) {
        nav_clock_offset = getled64((char *)buf, 11);
    } else {
        nav_clock_offset = NAN;
    }
    if (sizeof(float) == 4) {
        nav_clock_drift = getlef32((char *)buf, 19);
        osc_filter_drift_est = getlef32((char *)buf, 23);
    } else {
        nav_clock_drift = NAN;
        osc_filter_drift_est = NAN;
    }

    GPSD_LOG(LOG_DATA, &session->context->errout,
             "Navcom: oscillator temp. = %d, nav. status = 0x%02x, "
             "nav. clock offset = %f, nav. clock drift = %f, "
             "osc. filter drift est. = %f, acc.time slew value = %d\n",
             osc_temp, nav_status, nav_clock_offset, nav_clock_drift,
             osc_filter_drift_est, time_slew);
    GPSD_LOG(LOG_DATA, &session->context->errout,
             "CDO 0xef: time=%s mask={TIME}\n",
             timespec_str(&session->newdata.time, ts_buf, sizeof(ts_buf)));
    return 0;
}


gps_mask_t navcom_parse(struct gps_device_t * session, unsigned char *buf,
                        size_t len)
{
    unsigned char cmd_id;
    unsigned int msg_len;

    if (len == 0)
        return 0;

    cmd_id = (unsigned char)getub(buf, 3);
    //payload = &buf[6];
    msg_len = (unsigned int) getleu16(buf, 4);

    GPSD_LOG(LOG_RAW, &session->context->errout,
             "Navcom: packet type 0x%02x\n", cmd_id);

    session->cycle_end_reliable = true;

    switch (cmd_id) {
    case 0x06:
        return handle_0x06(session);
    case 0x15:
        return handle_0x15(session);
    case 0x81:
        return handle_0x81(session);
    case 0x83:
        return handle_0x83(session);
    case 0x86:
        return handle_0x86(session);
    case 0xae:
        return handle_0xae(session);
    case 0xb0:
        return handle_0xb0(session);
    case 0xb1:
        return handle_0xb1(session) | (CLEAR_IS | REPORT_IS);
    case 0xb5:
        return handle_0xb5(session);
    case 0xd3:
        return handle_0xd3(session);
    case 0xef:
        return handle_0xef(session);
    default:
        GPSD_LOG(LOG_PROG, &session->context->errout,
                 "Navcom: received packet type 0x%02x, length %d - "
                 "unknown or unimplemented\n",
                 cmd_id, msg_len);
        return 0;
    }
}


static gps_mask_t navcom_parse_input(struct gps_device_t *session)
{
    if (NAVCOM_PACKET == session->lexer.type) {
        return navcom_parse(session, session->lexer.outbuffer,
                          session->lexer.outbuflen);
    }
    if (NMEA_PACKET == session->lexer.type) {
        return nmea_parse((char *)session->lexer.outbuffer, session);;
    }
    return 0;
}

static ssize_t navcom_control_send(struct gps_device_t *session,
                                   char *buf, size_t len)
{
    putbyte(session->msgbuf, 0, 0x02);
    putbyte(session->msgbuf, 1, 0x99);
    putbyte(session->msgbuf, 2, 0x66);
    putbyte(session->msgbuf, 3, buf[0]);        /* Cmd ID */
    putle16(session->msgbuf, 4, len + 4);       /* Length */
    memcpy(session->msgbuf, buf + 6, len - 1);
    putbyte(session->msgbuf, 6 + len,
            checksum((unsigned char *)session->msgbuf + 3, len + 5));
    putbyte(session->msgbuf, 7 + len, 0x03);
    session->msgbuflen = len + 9;
    return gpsd_write(session, session->msgbuf, session->msgbuflen);
}

static bool navcom_speed(struct gps_device_t *session,
                         speed_t speed, char parity, int stopbits)
{
    /* parity and stopbit switching aren't implemented */
    if (parity != session->gpsdata.dev.parity
        || stopbits != (int)session->gpsdata.dev.parity) {
        return false;
    } else {
        uint8_t port, port_selection;
        uint8_t baud;
        if (session->driver.navcom.physical_port == (uint8_t) 0xFF) {
            /* We still don't know which port we're connected to */
            return false;
        }
        switch (speed) {
            /* NOTE - The spec says that certain baud combinations
             * on ports A and B are not allowed, those are
             * 1200/115200, 2400/57600, and 2400/115200.
             * To try and minimise the possibility of those
             * occurring, we do not allow baud rates below
             * 4800.  We could also disallow 57600 and 115200
             * to totally prevent this, but I do not consider
             * that reasonable.  Finding which baud speed the
             * other port is set at would also be too much
             * trouble, so we do not do it. */
        case 4800:
            baud = 0x04;
            break;
        case 9600:
            baud = 0x06;
            break;
        case 19200:
            baud = 0x08;
            break;
        case 38400:
            baud = 0x0a;
            break;
        case 57600:
            baud = 0x0c;
            break;
        case 115200:
            baud = 0x0e;
            break;
        default:
            /* Unsupported speed */
            return false;
        }

        /* Proceed to construct our message */
        port = session->driver.navcom.physical_port;
        port_selection = (port ? port : (uint8_t) 0xff) | baud;

        /* Send it off */
        navcom_cmd_0x11(session, port_selection);

        /* And cheekily return true, even though we have
         * no way to know if the speed change succeeded
         * until and if we receive an ACK (message 0x06),
         * which will be at the new baud speed if the
         * command was successful.  Bottom line, the client
         * should requery gpsd to see if the new speed is
         * different than the old one */
        return true;
    }
}

/* this is everything we export */
/* *INDENT-OFF* */
const struct gps_type_t driver_navcom =
{
    .type_name      = "Navcom NCT",             /* full name of type */
    .packet_type    = NAVCOM_PACKET,            /* lexer packet type */
    .flags          = DRIVER_STICKY,            /* remember this */
    .trigger        = NULL,                     /* none */
    .channels       = NAVCOM_CHANNELS,     // 12 L1 + 12 L2 + 2 Inmarsat L-Band
    .probe_detect   = NULL,                     /* no probe */
    .get_packet     = generic_get,              /* use generic one */
    .parse_packet   = navcom_parse_input,       /* parse message packets */
    .rtcm_writer    = gpsd_write,               /* send RTCM data straight */
    .init_query     = NULL,                     /* non-perturbing query */
    .event_hook     = navcom_event_hook,        /* lifetime event handler */
    .speed_switcher = navcom_speed,             /* we do change baud rates */
    .mode_switcher  = NULL,                     // there is not a mode switcher
    .rate_switcher  = NULL,                     /* no sample-rate switcher */
    .min_cycle.tv_sec  = 1,             /* not relevant, no rate switch */
    .min_cycle.tv_nsec = 0,             /* not relevant, no rate switch */
    .control_send   = navcom_control_send,      // how to send a control string
    .time_offset     = NULL,            /* no method for NTP fudge factor */
};
/* *INDENT-ON* */

#endif /* defined(NAVCOM_ENABLE) */
// vim: set expandtab shiftwidth=4
