/*

 * 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

        // Identification Block - send every 10 min

        navcom_cmd_0x20(session, 0xae, 0x1770);

        navcom_cmd_0x20(session, 0xb1, 0x4000); // PVT Block

        // Pseudorange Noise Statistics - send every 20s

        navcom_cmd_0x20(session, 0xb5, 0x00c8);

        navcom_cmd_0x20(session, 0xb0, 0x4000); // Raw Meas Data Block

        // Packed Ephemeris Data - send once

        navcom_cmd_0x20(session, 0x81, 0x0000);

        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

        // Clock Drift - send every 5 min

        navcom_cmd_0x20(session, 0xef, 0x0bb8);

    }

}

// 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);

    // This tells us which serial port was used last

    session->driver.navcom.physical_port = port;

    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 +