git.gag.com Git - fw/altos/blob - src/test/ao_flight_test.c

   1 /*
   2  * Copyright © 2009 Keith Packard <keithp@keithp.com>
   3  *
   4  * This program is free software; you can redistribute it and/or modify
   5  * it under the terms of the GNU General Public License as published by
   6  * the Free Software Foundation; version 2 of the License.
   7  *
   8  * This program is distributed in the hope that it will be useful, but
   9  * WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  11  * General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU General Public License along
  14  * with this program; if not, write to the Free Software Foundation, Inc.,
  15  * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
  16  */
  17
  18 #define _GNU_SOURCE
  19
  20 #include <stdint.h>
  21 #include <stdio.h>
  22 #include <stdlib.h>
  23 #include <string.h>
  24 #include <getopt.h>
  25 #include <math.h>
  26
  27 #define AO_HERTZ        100
  28
  29 #define HAS_ADC 1
  30 #define AO_DATA_RING    64
  31 #define ao_data_ring_next(n)    (((n) + 1) & (AO_DATA_RING - 1))
  32 #define ao_data_ring_prev(n)    (((n) - 1) & (AO_DATA_RING - 1))
  33
  34 #define AO_M_TO_HEIGHT(m)       ((int16_t) (m))
  35 #define AO_MS_TO_SPEED(ms)      ((int16_t) ((ms) * 16))
  36 #define AO_MSS_TO_ACCEL(mss)    ((int16_t) ((mss) * 16))
  37
  38 #if MEGAMETRUM
  39 #define AO_ADC_NUM_SENSE        6
  40 #define HAS_MS5607              1
  41 #define HAS_MPU6000             1
  42 #define HAS_MMA655X             0
  43
  44 struct ao_adc {
  45         int16_t                 sense[AO_ADC_NUM_SENSE];
  46         int16_t                 v_batt;
  47         int16_t                 v_pbatt;
  48         int16_t                 accel_ref;
  49         int16_t                 accel;
  50         int16_t                 temp;
  51 };
  52 #else
  53 /*
  54  * One set of samples read from the A/D converter
  55  */
  56 struct ao_adc {
  57         int16_t         accel;          /* accelerometer */
  58         int16_t         pres;           /* pressure sensor */
  59         int16_t         pres_real;      /* unclipped */
  60         int16_t         temp;           /* temperature sensor */
  61         int16_t         v_batt;         /* battery voltage */
  62         int16_t         sense_d;        /* drogue continuity sense */
  63         int16_t         sense_m;        /* main continuity sense */
  64 };
  65
  66 #ifndef HAS_ACCEL
  67 #define HAS_ACCEL 1
  68 #define HAS_ACCEL_REF 0
  69 #endif
  70
  71 #endif
  72
  73 #define __pdata
  74 #define __data
  75 #define __xdata
  76 #define __code
  77 #define __reentrant
  78
  79 #define HAS_FLIGHT 1
  80 #define HAS_IGNITE 1
  81 #define HAS_USB 1
  82 #define HAS_GPS 1
  83
  84 #include <ao_data.h>
  85 #include <ao_log.h>
  86
  87 #define to_fix16(x) ((int16_t) ((x) * 65536.0 + 0.5))
  88 #define to_fix32(x) ((int32_t) ((x) * 65536.0 + 0.5))
  89 #define from_fix(x)     ((x) >> 16)
  90
  91 /*
  92  * Above this height, the baro sensor doesn't work
  93  */
  94 #define AO_BARO_SATURATE        13000
  95 #define AO_MIN_BARO_VALUE       ao_altitude_to_pres(AO_BARO_SATURATE)
  96
  97 /*
  98  * Above this speed, baro measurements are unreliable
  99  */
 100 #define AO_MAX_BARO_SPEED       200
 101
 102 #define ACCEL_NOSE_UP   (ao_accel_2g >> 2)
 103
 104 extern enum ao_flight_state ao_flight_state;
 105
 106 #define FALSE 0
 107 #define TRUE 1
 108
 109 volatile struct ao_data ao_data_ring[AO_DATA_RING];
 110 volatile uint8_t ao_data_head;
 111 int     ao_summary = 0;
 112
 113 #define ao_led_on(l)
 114 #define ao_led_off(l)
 115 #define ao_timer_set_adc_interval(i)
 116 #define ao_wakeup(wchan) ao_dump_state()
 117 #define ao_cmd_register(c)
 118 #define ao_usb_disable()
 119 #define ao_telemetry_set_interval(x)
 120 #define ao_rdf_set(rdf)
 121 #define ao_packet_slave_start()
 122 #define ao_packet_slave_stop()
 123
 124 enum ao_igniter {
 125         ao_igniter_drogue = 0,
 126         ao_igniter_main = 1
 127 };
 128
 129 struct ao_data ao_data_static;
 130
 131 int     drogue_height;
 132 double  drogue_time;
 133 int     main_height;
 134 double  main_time;
 135
 136 int     tick_offset;
 137
 138 static int32_t  ao_k_height;
 139
 140 void
 141 ao_ignite(enum ao_igniter igniter)
 142 {
 143         double time = (double) (ao_data_static.tick + tick_offset) / 100;
 144
 145         if (igniter == ao_igniter_drogue) {
 146                 drogue_time = time;
 147                 drogue_height = ao_k_height >> 16;
 148         } else {
 149                 main_time = time;
 150                 main_height = ao_k_height >> 16;
 151         }
 152 }
 153
 154 struct ao_task {
 155         int dummy;
 156 };
 157
 158 #define ao_add_task(t,f,n) ((void) (t))
 159
 160 #define ao_log_start()
 161 #define ao_log_stop()
 162
 163 #define AO_MS_TO_TICKS(ms)      ((ms) / 10)
 164 #define AO_SEC_TO_TICKS(s)      ((s) * 100)
 165
 166 #define AO_FLIGHT_TEST
 167
 168 int     ao_flight_debug;
 169
 170 FILE *emulator_in;
 171 char *emulator_app;
 172 char *emulator_name;
 173 char *emulator_info;
 174 double emulator_error_max = 4;
 175 double emulator_height_error_max = 20;  /* noise in the baro sensor */
 176
 177 void
 178 ao_dump_state(void);
 179
 180 void
 181 ao_sleep(void *wchan);
 182
 183 const char const * const ao_state_names[] = {
 184         "startup", "idle", "pad", "boost", "fast",
 185         "coast", "drogue", "main", "landed", "invalid"
 186 };
 187
 188 struct ao_cmds {
 189         void            (*func)(void);
 190         const char      *help;
 191 };
 192
 193 #define ao_xmemcpy(d,s,c) memcpy(d,s,c)
 194 #define ao_xmemset(d,v,c) memset(d,v,c)
 195 #define ao_xmemcmp(d,s,c) memcmp(d,s,c)
 196
 197 #define AO_NEED_ALTITUDE_TO_PRES 1
 198 #if MEGAMETRUM
 199 #include "ao_convert_pa.c"
 200 #include <ao_ms5607.h>
 201 struct ao_ms5607_prom   ms5607_prom;
 202 #include "ao_ms5607_convert.c"
 203 #else
 204 #include "ao_convert.c"
 205 #endif
 206
 207 struct ao_config {
 208         uint16_t        main_deploy;
 209         int16_t         accel_plus_g;
 210         int16_t         accel_minus_g;
 211         uint8_t         pad_orientation;
 212         uint16_t        apogee_lockout;
 213 };
 214
 215 #define AO_PAD_ORIENTATION_ANTENNA_UP   0
 216 #define AO_PAD_ORIENTATION_ANTENNA_DOWN 1
 217
 218 #define ao_config_get()
 219
 220 struct ao_config ao_config;
 221
 222 #define DATA_TO_XDATA(x) (x)
 223
 224
 225 #define GRAVITY 9.80665
 226 extern int16_t ao_ground_accel, ao_flight_accel;
 227 extern int16_t ao_accel_2g;
 228
 229 typedef int16_t accel_t;
 230
 231 extern uint16_t ao_sample_tick;
 232
 233 extern alt_t    ao_sample_height;
 234 extern accel_t  ao_sample_accel;
 235 extern int32_t  ao_accel_scale;
 236 extern alt_t    ao_ground_height;
 237 extern alt_t    ao_sample_alt;
 238
 239 int ao_sample_prev_tick;
 240 uint16_t        prev_tick;
 241
 242 #include "ao_kalman.c"
 243 #include "ao_sample.c"
 244 #include "ao_flight.c"
 245
 246 #define to_double(f)    ((f) / 65536.0)
 247
 248 static int      ao_records_read = 0;
 249 static int      ao_eof_read = 0;
 250 static int      ao_flight_ground_accel;
 251 static int      ao_flight_started = 0;
 252 static int      ao_test_max_height;
 253 static double   ao_test_max_height_time;
 254 static int      ao_test_main_height;
 255 static double   ao_test_main_height_time;
 256 static double   ao_test_landed_time;
 257 static double   ao_test_landed_height;
 258 static double   ao_test_landed_time;
 259 static int      landed_set;
 260 static double   landed_time;
 261 static double   landed_height;
 262
 263 #if HAS_MPU6000
 264 static struct ao_mpu6000_sample ao_ground_mpu6000;
 265 #endif
 266
 267 void
 268 ao_test_exit(void)
 269 {
 270         double  drogue_error;
 271         double  main_error;
 272         double  landed_error;
 273         double  landed_time_error;
 274
 275         if (!ao_test_main_height_time) {
 276                 ao_test_main_height_time = ao_test_max_height_time;
 277                 ao_test_main_height = ao_test_max_height;
 278         }
 279         drogue_error = fabs(ao_test_max_height_time - drogue_time);
 280         main_error = fabs(ao_test_main_height_time - main_time);
 281         landed_error = fabs(ao_test_landed_height - landed_height);
 282         landed_time_error = ao_test_landed_time - landed_time;
 283         if (drogue_error > emulator_error_max || main_error > emulator_error_max) {
 284                 printf ("%s %s\n",
 285                         emulator_app, emulator_name);
 286                 if (emulator_info)
 287                         printf ("\t%s\n", emulator_info);
 288                 printf ("\tApogee error %g\n", drogue_error);
 289                 printf ("\tMain error %g\n", main_error);
 290                 printf ("\tLanded height error %g\n", landed_error);
 291                 printf ("\tLanded time error %g\n", landed_time_error);
 292                 printf ("\tActual: apogee: %d at %7.2f main: %d at %7.2f landed %7.2f at %7.2f\n",
 293                         ao_test_max_height, ao_test_max_height_time,
 294                         ao_test_main_height, ao_test_main_height_time,
 295                         ao_test_landed_height, ao_test_landed_time);
 296                 printf ("\tComputed: apogee: %d at %7.2f main: %d at %7.2f landed %7.2f at %7.2f\n",
 297                         drogue_height, drogue_time, main_height, main_time,
 298                         landed_height, landed_time);
 299                 exit (1);
 300         }
 301         exit(0);
 302 }
 303
 304 #if HAS_MPU6000
 305 static double
 306 ao_mpu6000_accel(int16_t sensor)
 307 {
 308         return sensor / 32767.0 * MPU6000_ACCEL_FULLSCALE * GRAVITY;
 309 }
 310
 311 static double
 312 ao_mpu6000_gyro(int16_t sensor)
 313 {
 314         return sensor / 32767.0 * MPU6000_GYRO_FULLSCALE;
 315 }
 316 #endif
 317
 318 void
 319 ao_insert(void)
 320 {
 321         double  time;
 322
 323         ao_data_ring[ao_data_head] = ao_data_static;
 324         ao_data_head = ao_data_ring_next(ao_data_head);
 325         if (ao_flight_state != ao_flight_startup) {
 326 #if HAS_ACCEL
 327                 double  accel = ((ao_flight_ground_accel - ao_data_accel_cook(&ao_data_static)) * GRAVITY * 2.0) /
 328                         (ao_config.accel_minus_g - ao_config.accel_plus_g);
 329 #else
 330                 double  accel = 0.0;
 331 #endif
 332 #if MEGAMETRUM
 333                 double  height;
 334
 335                 ao_ms5607_convert(&ao_data_static.ms5607_raw, &ao_data_static.ms5607_cooked);
 336                 height = ao_pa_to_altitude(ao_data_static.ms5607_cooked.pres) - ao_ground_height;
 337 #else
 338                 double  height = ao_pres_to_altitude(ao_data_static.adc.pres_real) - ao_ground_height;
 339 #endif
 340
 341                 if (!tick_offset)
 342                         tick_offset = -ao_data_static.tick;
 343                 if ((prev_tick - ao_data_static.tick) > 0x400)
 344                         tick_offset += 65536;
 345                 prev_tick = ao_data_static.tick;
 346                 time = (double) (ao_data_static.tick + tick_offset) / 100;
 347
 348                 if (ao_test_max_height < height) {
 349                         ao_test_max_height = height;
 350                         ao_test_max_height_time = time;
 351                         ao_test_landed_height = height;
 352                         ao_test_landed_time = time;
 353                 }
 354                 if (height > ao_config.main_deploy) {
 355                         ao_test_main_height_time = time;
 356                         ao_test_main_height = height;
 357                 }
 358
 359                 if (ao_test_landed_height > height) {
 360                         ao_test_landed_height = height;
 361                         ao_test_landed_time = time;
 362                 }
 363
 364                 if (ao_flight_state == ao_flight_landed && !landed_set) {
 365                         landed_set = 1;
 366                         landed_time = time;
 367                         landed_height = height;
 368                 }
 369
 370                 if (!ao_summary) {
 371                         printf("%7.2f height %8.2f accel %8.3f "
 372 #if MEGAMETRUM
 373                                "accel_x %8.3f accel_y %8.3f accel_z %8.3f gyro_x %8.3f gyro_y %8.3f gyro_z %8.3f "
 374 #endif
 375                                "state %-8.8s k_height %8.2f k_speed %8.3f k_accel %8.3f avg_height %5d drogue %4d main %4d error %5d\n",
 376                                time,
 377                                height,
 378                                accel,
 379 #if MEGAMETRUM
 380                                ao_mpu6000_accel(ao_data_static.mpu6000.accel_x),
 381                                ao_mpu6000_accel(ao_data_static.mpu6000.accel_y),
 382                                ao_mpu6000_accel(ao_data_static.mpu6000.accel_z),
 383                                ao_mpu6000_gyro(ao_data_static.mpu6000.gyro_x - ao_ground_mpu6000.gyro_x),
 384                                ao_mpu6000_gyro(ao_data_static.mpu6000.gyro_y - ao_ground_mpu6000.gyro_y),
 385                                ao_mpu6000_gyro(ao_data_static.mpu6000.gyro_z - ao_ground_mpu6000.gyro_z),
 386 #endif
 387                                ao_state_names[ao_flight_state],
 388                                ao_k_height / 65536.0,
 389                                ao_k_speed / 65536.0 / 16.0,
 390                                ao_k_accel / 65536.0 / 16.0,
 391                                ao_avg_height,
 392                                drogue_height,
 393                                main_height,
 394                                ao_error_h_sq_avg);
 395
 396 //                      if (ao_flight_state == ao_flight_landed)
 397 //                              ao_test_exit();
 398                 }
 399         }
 400 }
 401
 402 #define AO_MAX_CALLSIGN                 8
 403 #define AO_MAX_VERSION                  8
 404 #define AO_MAX_TELEMETRY                128
 405
 406 struct ao_telemetry_generic {
 407         uint16_t        serial;         /* 0 */
 408         uint16_t        tick;           /* 2 */
 409         uint8_t         type;           /* 4 */
 410         uint8_t         payload[27];    /* 5 */
 411         /* 32 */
 412 };
 413
 414 #define AO_TELEMETRY_SENSOR_TELEMETRUM  0x01
 415 #define AO_TELEMETRY_SENSOR_TELEMINI    0x02
 416 #define AO_TELEMETRY_SENSOR_TELENANO    0x03
 417
 418 struct ao_telemetry_sensor {
 419         uint16_t        serial;         /*  0 */
 420         uint16_t        tick;           /*  2 */
 421         uint8_t         type;           /*  4 */
 422
 423         uint8_t         state;          /*  5 flight state */
 424         int16_t         accel;          /*  6 accelerometer (TM only) */
 425         int16_t         pres;           /*  8 pressure sensor */
 426         int16_t         temp;           /* 10 temperature sensor */
 427         int16_t         v_batt;         /* 12 battery voltage */
 428         int16_t         sense_d;        /* 14 drogue continuity sense (TM/Tm) */
 429         int16_t         sense_m;        /* 16 main continuity sense (TM/Tm) */
 430
 431         int16_t         acceleration;   /* 18 m/s² * 16 */
 432         int16_t         speed;          /* 20 m/s * 16 */
 433         int16_t         height;         /* 22 m */
 434
 435         int16_t         ground_pres;    /* 24 average pres on pad */
 436         int16_t         ground_accel;   /* 26 average accel on pad */
 437         int16_t         accel_plus_g;   /* 28 accel calibration at +1g */
 438         int16_t         accel_minus_g;  /* 30 accel calibration at -1g */
 439         /* 32 */
 440 };
 441
 442 #define AO_TELEMETRY_CONFIGURATION      0x04
 443
 444 struct ao_telemetry_configuration {
 445         uint16_t        serial;                         /*  0 */
 446         uint16_t        tick;                           /*  2 */
 447         uint8_t         type;                           /*  4 */
 448
 449         uint8_t         device;                         /*  5 device type */
 450         uint16_t        flight;                         /*  6 flight number */
 451         uint8_t         config_major;                   /*  8 Config major version */
 452         uint8_t         config_minor;                   /*  9 Config minor version */
 453         uint16_t        apogee_delay;                   /* 10 Apogee deploy delay in seconds */
 454         uint16_t        main_deploy;                    /* 12 Main deploy alt in meters */
 455         uint16_t        flight_log_max;                 /* 14 Maximum flight log size in kB */
 456         char            callsign[AO_MAX_CALLSIGN];      /* 16 Radio operator identity */
 457         char            version[AO_MAX_VERSION];        /* 24 Software version */
 458         /* 32 */
 459 };
 460
 461 #define AO_TELEMETRY_LOCATION           0x05
 462
 463 #define AO_GPS_MODE_NOT_VALID           'N'
 464 #define AO_GPS_MODE_AUTONOMOUS          'A'
 465 #define AO_GPS_MODE_DIFFERENTIAL        'D'
 466 #define AO_GPS_MODE_ESTIMATED           'E'
 467 #define AO_GPS_MODE_MANUAL              'M'
 468 #define AO_GPS_MODE_SIMULATED           'S'
 469
 470 struct ao_telemetry_location {
 471         uint16_t        serial;         /*  0 */
 472         uint16_t        tick;           /*  2 */
 473         uint8_t         type;           /*  4 */
 474
 475         uint8_t         flags;          /*  5 Number of sats and other flags */
 476         int16_t         altitude;       /*  6 GPS reported altitude (m) */
 477         int32_t         latitude;       /*  8 latitude (degrees * 10⁷) */
 478         int32_t         longitude;      /* 12 longitude (degrees * 10⁷) */
 479         uint8_t         year;           /* 16 (- 2000) */
 480         uint8_t         month;          /* 17 (1-12) */
 481         uint8_t         day;            /* 18 (1-31) */
 482         uint8_t         hour;           /* 19 (0-23) */
 483         uint8_t         minute;         /* 20 (0-59) */
 484         uint8_t         second;         /* 21 (0-59) */
 485         uint8_t         pdop;           /* 22 (m * 5) */
 486         uint8_t         hdop;           /* 23 (m * 5) */
 487         uint8_t         vdop;           /* 24 (m * 5) */
 488         uint8_t         mode;           /* 25 */
 489         uint16_t        ground_speed;   /* 26 cm/s */
 490         int16_t         climb_rate;     /* 28 cm/s */
 491         uint8_t         course;         /* 30 degrees / 2 */
 492         uint8_t         unused[1];      /* 31 */
 493         /* 32 */
 494 };
 495
 496 #define AO_TELEMETRY_SATELLITE          0x06
 497
 498 struct ao_telemetry_satellite_info {
 499         uint8_t         svid;
 500         uint8_t         c_n_1;
 501 };
 502
 503 struct ao_telemetry_satellite {
 504         uint16_t                                serial;         /*  0 */
 505         uint16_t                                tick;           /*  2 */
 506         uint8_t                                 type;           /*  4 */
 507         uint8_t                                 channels;       /*  5 number of reported sats */
 508
 509         struct ao_telemetry_satellite_info      sats[12];       /* 6 */
 510         uint8_t                                 unused[2];      /* 30 */
 511         /* 32 */
 512 };
 513
 514 union ao_telemetry_all {
 515         struct ao_telemetry_generic             generic;
 516         struct ao_telemetry_sensor              sensor;
 517         struct ao_telemetry_configuration       configuration;
 518         struct ao_telemetry_location            location;
 519         struct ao_telemetry_satellite           satellite;
 520 };
 521
 522 uint16_t
 523 uint16(uint8_t *bytes, int off)
 524 {
 525         return (uint16_t) bytes[off] | (((uint16_t) bytes[off+1]) << 8);
 526 }
 527
 528 int16_t
 529 int16(uint8_t *bytes, int off)
 530 {
 531         return (int16_t) uint16(bytes, off);
 532 }
 533
 534 uint32_t
 535 uint32(uint8_t *bytes, int off)
 536 {
 537         return (uint32_t) bytes[off] | (((uint32_t) bytes[off+1]) << 8) |
 538                 (((uint32_t) bytes[off+2]) << 16) |
 539                 (((uint32_t) bytes[off+3]) << 24);
 540 }
 541
 542 int32_t
 543 int32(uint8_t *bytes, int off)
 544 {
 545         return (int32_t) uint32(bytes, off);
 546 }
 547
 548 static int log_format;
 549
 550 void
 551 ao_sleep(void *wchan)
 552 {
 553         if (wchan == &ao_data_head) {
 554                 char            type = 0;
 555                 uint16_t        tick = 0;
 556                 uint16_t        a = 0, b = 0;
 557                 uint8_t         bytes[1024];
 558                 union ao_telemetry_all  telem;
 559                 char            line[1024];
 560                 char            *saveptr;
 561                 char            *l;
 562                 char            *words[64];
 563                 int             nword;
 564
 565                 for (;;) {
 566                         if (ao_records_read > 2 && ao_flight_state == ao_flight_startup)
 567                         {
 568 #if MEGAMETRUM
 569                                 ao_data_static.mpu6000 = ao_ground_mpu6000;
 570 #else
 571                                 ao_data_static.adc.accel = ao_flight_ground_accel;
 572 #endif
 573                                 ao_insert();
 574                                 return;
 575                         }
 576
 577                         if (!fgets(line, sizeof (line), emulator_in)) {
 578                                 if (++ao_eof_read >= 1000) {
 579                                         if (!ao_summary)
 580                                                 printf ("no more data, exiting simulation\n");
 581                                         ao_test_exit();
 582                                 }
 583                                 ao_data_static.tick += 10;
 584                                 ao_insert();
 585                                 return;
 586                         }
 587                         l = line;
 588                         for (nword = 0; nword < 64; nword++) {
 589                                 words[nword] = strtok_r(l, " \t\n", &saveptr);
 590                                 l = NULL;
 591                                 if (words[nword] == NULL)
 592                                         break;
 593                         }
 594 #if MEGAMETRUM
 595                         if (log_format == AO_LOG_FORMAT_MEGAMETRUM && nword == 30 && strlen(words[0]) == 1) {
 596                                 int     i;
 597                                 struct ao_ms5607_value  value;
 598
 599                                 type = words[0][0];
 600                                 tick = strtoul(words[1], NULL, 16);
 601 //                              printf ("%c %04x", type, tick);
 602                                 for (i = 2; i < nword; i++) {
 603                                         bytes[i - 2] = strtoul(words[i], NULL, 16);
 604 //                                      printf(" %02x", bytes[i-2]);
 605                                 }
 606 //                              printf ("\n");
 607                                 switch (type) {
 608                                 case 'F':
 609                                         ao_flight_ground_accel = int16(bytes, 2);
 610                                         ao_flight_started = 1;
 611                                         ao_ground_pres = int32(bytes, 4);
 612                                         ao_ground_height = ao_pa_to_altitude(ao_ground_pres);
 613                                         break;
 614                                 case 'A':
 615                                         ao_data_static.tick = tick;
 616                                         ao_data_static.ms5607_raw.pres = int32(bytes, 0);
 617                                         ao_data_static.ms5607_raw.temp = int32(bytes, 4);
 618                                         ao_ms5607_convert(&ao_data_static.ms5607_raw, &value);
 619                                         ao_data_static.mpu6000.accel_x = int16(bytes, 8);
 620                                         ao_data_static.mpu6000.accel_y = -int16(bytes, 10);
 621                                         ao_data_static.mpu6000.accel_z = int16(bytes, 12);
 622                                         ao_data_static.mpu6000.gyro_x = int16(bytes, 14);
 623                                         ao_data_static.mpu6000.gyro_y = -int16(bytes, 16);
 624                                         ao_data_static.mpu6000.gyro_z = int16(bytes, 18);
 625                                         if (ao_records_read == 0)
 626                                                 ao_ground_mpu6000 = ao_data_static.mpu6000;
 627                                         else if (ao_records_read < 10) {
 628 #define f(f) ao_ground_mpu6000.f = ao_ground_mpu6000.f + ((ao_data_static.mpu6000.f - ao_ground_mpu6000.f) >> 2)
 629                                                 f(accel_x);
 630                                                 f(accel_y);
 631                                                 f(accel_z);
 632                                                 f(gyro_x);
 633                                                 f(gyro_y);
 634                                                 f(gyro_z);
 635                                         }
 636                                         ao_records_read++;
 637                                         ao_insert();
 638                                         return;
 639                                 }
 640                                 continue;
 641                         } else if (nword == 3 && strcmp(words[0], "ms5607") == 0) {
 642                                 if (strcmp(words[1], "reserved:") == 0)
 643                                         ms5607_prom.reserved = strtoul(words[2], NULL, 10);
 644                                 else if (strcmp(words[1], "sens:") == 0)
 645                                         ms5607_prom.sens = strtoul(words[2], NULL, 10);
 646                                 else if (strcmp(words[1], "off:") == 0)
 647                                         ms5607_prom.off = strtoul(words[2], NULL, 10);
 648                                 else if (strcmp(words[1], "tcs:") == 0)
 649                                         ms5607_prom.tcs = strtoul(words[2], NULL, 10);
 650                                 else if (strcmp(words[1], "tco:") == 0)
 651                                         ms5607_prom.tco = strtoul(words[2], NULL, 10);
 652                                 else if (strcmp(words[1], "tref:") == 0)
 653                                         ms5607_prom.tref = strtoul(words[2], NULL, 10);
 654                                 else if (strcmp(words[1], "tempsens:") == 0)
 655                                         ms5607_prom.tempsens = strtoul(words[2], NULL, 10);
 656                                 else if (strcmp(words[1], "crc:") == 0)
 657                                         ms5607_prom.crc = strtoul(words[2], NULL, 10);
 658                                 continue;
 659                         }
 660 #else
 661                         if (nword == 4 && log_format != AO_LOG_FORMAT_MEGAMETRUM) {
 662                                 type = words[0][0];
 663                                 tick = strtoul(words[1], NULL, 16);
 664                                 a = strtoul(words[2], NULL, 16);
 665                                 b = strtoul(words[3], NULL, 16);
 666                                 if (type == 'P')
 667                                         type = 'A';
 668                         }
 669 #endif
 670                         else if (nword == 2 && strcmp(words[0], "log-format") == 0) {
 671                                 log_format = strtoul(words[1], NULL, 10);
 672                         } else if (nword >= 6 && strcmp(words[0], "Accel") == 0) {
 673                                 ao_config.accel_plus_g = atoi(words[3]);
 674                                 ao_config.accel_minus_g = atoi(words[5]);
 675                         } else if (nword >= 4 && strcmp(words[0], "Main") == 0) {
 676                                 ao_config.main_deploy = atoi(words[2]);
 677                         } else if (nword >= 3 && strcmp(words[0], "Apogee") == 0 &&
 678                                    strcmp(words[1], "lockout:") == 0) {
 679                                 ao_config.apogee_lockout = atoi(words[2]);
 680                         } else if (nword >= 36 && strcmp(words[0], "CALL") == 0) {
 681                                 tick = atoi(words[10]);
 682                                 if (!ao_flight_started) {
 683                                         type = 'F';
 684                                         a = atoi(words[26]);
 685                                         ao_flight_started = 1;
 686                                 } else {
 687                                         type = 'A';
 688                                         a = atoi(words[12]);
 689                                         b = atoi(words[14]);
 690                                 }
 691                         } else if (nword == 3 && strcmp(words[0], "BARO") == 0) {
 692                                 tick = strtol(words[1], NULL, 16);
 693                                 a = 16384 - 328;
 694                                 b = strtol(words[2], NULL, 10);
 695                                 type = 'A';
 696                                 if (!ao_flight_started) {
 697                                         ao_flight_ground_accel = 16384 - 328;
 698                                         ao_config.accel_plus_g = 16384 - 328;
 699                                         ao_config.accel_minus_g = 16384 + 328;
 700                                         ao_flight_started = 1;
 701                                 }
 702                         } else if (nword == 2 && strcmp(words[0], "TELEM") == 0) {
 703                                 __xdata char    *hex = words[1];
 704                                 char    elt[3];
 705                                 int     i, len;
 706                                 uint8_t sum;
 707
 708                                 len = strlen(hex);
 709                                 if (len > sizeof (bytes) * 2) {
 710                                         len = sizeof (bytes)*2;
 711                                         hex[len] = '\0';
 712                                 }
 713                                 for (i = 0; i < len; i += 2) {
 714                                         elt[0] = hex[i];
 715                                         elt[1] = hex[i+1];
 716                                         elt[2] = '\0';
 717                                         bytes[i/2] = (uint8_t) strtol(elt, NULL, 16);
 718                                 }
 719                                 len = i/2;
 720                                 if (bytes[0] != len - 2) {
 721                                         printf ("bad length %d != %d\n", bytes[0], len - 2);
 722                                         continue;
 723                                 }
 724                                 sum = 0x5a;
 725                                 for (i = 1; i < len-1; i++)
 726                                         sum += bytes[i];
 727                                 if (sum != bytes[len-1]) {
 728                                         printf ("bad checksum\n");
 729                                         continue;
 730                                 }
 731                                 if ((bytes[len-2] & 0x80) == 0) {
 732                                         continue;
 733                                 }
 734                                 if (len == 36) {
 735                                         ao_xmemcpy(&telem, bytes + 1, 32);
 736                                         tick = telem.generic.tick;
 737                                         switch (telem.generic.type) {
 738                                         case AO_TELEMETRY_SENSOR_TELEMETRUM:
 739                                         case AO_TELEMETRY_SENSOR_TELEMINI:
 740                                         case AO_TELEMETRY_SENSOR_TELENANO:
 741                                                 if (!ao_flight_started) {
 742                                                         ao_flight_ground_accel = telem.sensor.ground_accel;
 743                                                         ao_config.accel_plus_g = telem.sensor.accel_plus_g;
 744                                                         ao_config.accel_minus_g = telem.sensor.accel_minus_g;
 745                                                         ao_flight_started = 1;
 746                                                 }
 747                                                 type = 'A';
 748                                                 a = telem.sensor.accel;
 749                                                 b = telem.sensor.pres;
 750                                                 break;
 751                                         }
 752                                 } else if (len == 99) {
 753                                         ao_flight_started = 1;
 754                                         tick = uint16(bytes+1, 21);
 755                                         ao_flight_ground_accel = int16(bytes+1, 7);
 756                                         ao_config.accel_plus_g = int16(bytes+1, 17);
 757                                         ao_config.accel_minus_g = int16(bytes+1, 19);
 758                                         type = 'A';
 759                                         a = int16(bytes+1, 23);
 760                                         b = int16(bytes+1, 25);
 761                                 } else if (len == 98) {
 762                                         ao_flight_started = 1;
 763                                         tick = uint16(bytes+1, 20);
 764                                         ao_flight_ground_accel = int16(bytes+1, 6);
 765                                         ao_config.accel_plus_g = int16(bytes+1, 16);
 766                                         ao_config.accel_minus_g = int16(bytes+1, 18);
 767                                         type = 'A';
 768                                         a = int16(bytes+1, 22);
 769                                         b = int16(bytes+1, 24);
 770                                 } else {
 771                                         printf("unknown len %d\n", len);
 772                                         continue;
 773                                 }
 774                         }
 775                         if (type != 'F' && !ao_flight_started)
 776                                 continue;
 777
 778 #if MEGAMETRUM
 779 #else
 780                         switch (type) {
 781                         case 'F':
 782                                 ao_flight_ground_accel = a;
 783                                 if (ao_config.accel_plus_g == 0) {
 784                                         ao_config.accel_plus_g = a;
 785                                         ao_config.accel_minus_g = a + 530;
 786                                 }
 787                                 if (ao_config.main_deploy == 0)
 788                                         ao_config.main_deploy = 250;
 789                                 ao_flight_started = 1;
 790                                 break;
 791                         case 'S':
 792                                 break;
 793                         case 'A':
 794                                 ao_data_static.tick = tick;
 795                                 ao_data_static.adc.accel = a;
 796                                 ao_data_static.adc.pres_real = b;
 797                                 if (b < AO_MIN_BARO_VALUE)
 798                                         b = AO_MIN_BARO_VALUE;
 799                                 ao_data_static.adc.pres = b;
 800                                 ao_records_read++;
 801                                 ao_insert();
 802                                 return;
 803                         case 'T':
 804                                 ao_data_static.tick = tick;
 805                                 ao_data_static.adc.temp = a;
 806                                 ao_data_static.adc.v_batt = b;
 807                                 break;
 808                         case 'D':
 809                         case 'G':
 810                         case 'N':
 811                         case 'W':
 812                         case 'H':
 813                                 break;
 814                         }
 815 #endif
 816                 }
 817
 818         }
 819 }
 820 #define COUNTS_PER_G 264.8
 821
 822 void
 823 ao_dump_state(void)
 824 {
 825 }
 826
 827 static const struct option options[] = {
 828         { .name = "summary", .has_arg = 0, .val = 's' },
 829         { .name = "debug", .has_arg = 0, .val = 'd' },
 830         { .name = "info", .has_arg = 1, .val = 'i' },
 831         { 0, 0, 0, 0},
 832 };
 833
 834 void run_flight_fixed(char *name, FILE *f, int summary, char *info)
 835 {
 836         emulator_name = name;
 837         emulator_in = f;
 838         emulator_info = info;
 839         ao_summary = summary;
 840         ao_flight_init();
 841         ao_flight();
 842 }
 843
 844 int
 845 main (int argc, char **argv)
 846 {
 847         int     summary = 0;
 848         int     c;
 849         int     i;
 850         char    *info = NULL;
 851
 852 #if HAS_ACCEL
 853         emulator_app="full";
 854 #else
 855         emulator_app="baro";
 856 #endif
 857         while ((c = getopt_long(argc, argv, "sdi:", options, NULL)) != -1) {
 858                 switch (c) {
 859                 case 's':
 860                         summary = 1;
 861                         break;
 862                 case 'd':
 863                         ao_flight_debug = 1;
 864                         break;
 865                 case 'i':
 866                         info = optarg;
 867                         break;
 868                 }
 869         }
 870
 871         if (optind == argc)
 872                 run_flight_fixed("<stdin>", stdin, summary, info);
 873         else
 874                 for (i = optind; i < argc; i++) {
 875                         FILE    *f = fopen(argv[i], "r");
 876                         if (!f) {
 877                                 perror(argv[i]);
 878                                 continue;
 879                         }
 880                         run_flight_fixed(argv[i], f, summary, info);
 881                         fclose(f);
 882                 }
 883         exit(0);
 884 }