2 * Copyright © 2013 Keith Packard <keithp@keithp.com>
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; either version 2 of the License, or
7 * (at your option) any later version.
9 * This program is distributed in the hope that it will be useful, but
10 * WITHOUT ANY WARRANTY; without even the implied warranty of
11 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 * General Public License for more details.
14 * You should have received a copy of the GNU General Public License along
15 * with this program; if not, write to the Free Software Foundation, Inc.,
16 * 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
19 #ifndef AO_FLIGHT_TEST
22 #include <ao_micropeak.h>
23 #include <ao_log_micro.h>
33 ao_microkalman_predict();
34 ao_microkalman_correct();
37 #define NUM_PA_HIST (GROUND_AVG)
39 #define SKIP_PA_HIST(i,j) (((i) + (j)) & (NUM_PA_HIST - 1))
41 static uint32_t pa_hist[NUM_PA_HIST];
49 uint32_t pa_interval_min, pa_interval_max;
53 uint8_t accel_lock = 0;
56 /* Wait for motion, averaging values to get ground pressure */
60 ao_microkalman_init();
67 if ((sample_count & 0x1f) == 0)
68 ao_led_on(AO_LED_REPORT);
73 if ((sample_count & 0x1f) == 0)
74 ao_led_off(AO_LED_REPORT);
77 h = SKIP_PA_HIST(h,1);
78 pa_diff = pa_ground - ao_pa;
81 /* Check for a significant pressure change */
82 if (pa_diff > BOOST_DETECT)
86 if (sample_count < GROUND_AVG * 2) {
87 if (sample_count < GROUND_AVG)
91 pa_ground = pa_sum >> GROUND_AVG_SHIFT;
100 /* Go back and find the last sample close to the ground */
101 pa_min = pa_ground - LAND_DETECT;
102 for (i = SKIP_PA_HIST(h,-2); i != SKIP_PA_HIST(h,2); i = SKIP_PA_HIST(i,-2)) {
103 if (pa_hist[i] >= pa_min)
107 /* Log the remaining samples so we get a complete history since leaving the ground */
108 #if LOG_INTERVAL < NUM_PA_HIST
109 for (; i != h; i = SKIP_PA_HIST(i,2))
116 /* Now sit around until the pressure is stable again and record the max */
121 pa_interval_min = ao_pa;
122 pa_interval_max = ao_pa;
124 time += SAMPLE_SLEEP;
125 ao_delay_until(time);
126 if ((sample_count & 3) == 0)
127 ao_led_on(AO_LED_REPORT);
129 if ((sample_count & 3) == 0)
130 ao_led_off(AO_LED_REPORT);
131 if (log_count == LOG_INTERVAL - 1) {
137 /* If accelerating upwards, don't look for min pressure */
138 if (ao_pa_accel < ACCEL_LOCK_PA)
139 accel_lock = ACCEL_LOCK_TIME;
142 else if (ao_pa < pa_min)
145 if (sample_count == (GROUND_AVG - 1)) {
146 pa_diff = pa_interval_max - pa_interval_min;
148 /* Check to see if the pressure is now stable */
149 if (pa_diff < LAND_DETECT)
152 pa_interval_min = ao_pa;
153 pa_interval_max = ao_pa;
155 if (ao_pa < pa_interval_min)
156 pa_interval_min = ao_pa;
157 if (ao_pa > pa_interval_max)
158 pa_interval_max = ao_pa;