2019-06-12 21:35:47 +08:00
|
|
|
// include first, needed by safety policies
|
|
|
|
#include "safety_declarations.h"
|
2017-07-15 12:17:32 +08:00
|
|
|
// Include the actual safety policies.
|
2017-07-18 01:20:08 +08:00
|
|
|
#include "safety/safety_defaults.h"
|
|
|
|
#include "safety/safety_honda.h"
|
2017-08-25 13:31:34 +08:00
|
|
|
#include "safety/safety_toyota.h"
|
2018-10-05 01:22:10 +08:00
|
|
|
#include "safety/safety_tesla.h"
|
2019-02-06 10:37:02 +08:00
|
|
|
#include "safety/safety_gm_ascm.h"
|
2018-02-23 16:19:22 +08:00
|
|
|
#include "safety/safety_gm.h"
|
2018-05-03 13:45:23 +08:00
|
|
|
#include "safety/safety_ford.h"
|
2018-08-18 12:31:00 +08:00
|
|
|
#include "safety/safety_hyundai.h"
|
2018-11-07 04:28:33 +08:00
|
|
|
#include "safety/safety_chrysler.h"
|
2019-01-29 13:42:06 +08:00
|
|
|
#include "safety/safety_subaru.h"
|
2019-09-10 06:58:41 +08:00
|
|
|
#include "safety/safety_mazda.h"
|
2020-02-27 03:11:56 +08:00
|
|
|
#include "safety/safety_nissan.h"
|
2019-10-10 04:54:22 +08:00
|
|
|
#include "safety/safety_volkswagen.h"
|
2017-08-12 07:17:43 +08:00
|
|
|
#include "safety/safety_elm327.h"
|
2017-07-15 12:17:32 +08:00
|
|
|
|
2019-10-05 04:30:00 +08:00
|
|
|
// from cereal.car.CarParams.SafetyModel
|
2019-11-21 03:56:26 +08:00
|
|
|
#define SAFETY_SILENT 0U
|
2019-12-21 07:51:32 +08:00
|
|
|
#define SAFETY_HONDA_NIDEC 1U
|
2019-10-05 04:30:00 +08:00
|
|
|
#define SAFETY_TOYOTA 2U
|
|
|
|
#define SAFETY_ELM327 3U
|
|
|
|
#define SAFETY_GM 4U
|
2019-12-21 07:51:32 +08:00
|
|
|
#define SAFETY_HONDA_BOSCH_GIRAFFE 5U
|
2019-10-05 04:30:00 +08:00
|
|
|
#define SAFETY_FORD 6U
|
|
|
|
#define SAFETY_HYUNDAI 8U
|
|
|
|
#define SAFETY_CHRYSLER 9U
|
|
|
|
#define SAFETY_TESLA 10U
|
|
|
|
#define SAFETY_SUBARU 11U
|
|
|
|
#define SAFETY_MAZDA 13U
|
2020-02-27 03:11:56 +08:00
|
|
|
#define SAFETY_NISSAN 14U
|
2020-02-20 04:37:38 +08:00
|
|
|
#define SAFETY_VOLKSWAGEN_MQB 15U
|
2019-10-05 04:30:00 +08:00
|
|
|
#define SAFETY_ALLOUTPUT 17U
|
|
|
|
#define SAFETY_GM_ASCM 18U
|
2019-11-21 03:56:26 +08:00
|
|
|
#define SAFETY_NOOUTPUT 19U
|
2019-12-21 07:51:32 +08:00
|
|
|
#define SAFETY_HONDA_BOSCH_HARNESS 20U
|
2020-04-02 03:38:27 +08:00
|
|
|
#define SAFETY_VOLKSWAGEN_PQ 21U
|
2020-02-28 07:28:24 +08:00
|
|
|
#define SAFETY_SUBARU_LEGACY 22U
|
2019-10-05 04:30:00 +08:00
|
|
|
|
2019-11-21 03:56:26 +08:00
|
|
|
uint16_t current_safety_mode = SAFETY_SILENT;
|
2017-07-15 12:17:32 +08:00
|
|
|
const safety_hooks *current_hooks = &nooutput_hooks;
|
|
|
|
|
2019-12-21 17:25:54 +08:00
|
|
|
int safety_rx_hook(CAN_FIFOMailBox_TypeDef *to_push){
|
|
|
|
return current_hooks->rx(to_push);
|
2017-07-15 12:17:32 +08:00
|
|
|
}
|
|
|
|
|
2017-07-21 14:36:06 +08:00
|
|
|
int safety_tx_hook(CAN_FIFOMailBox_TypeDef *to_send) {
|
|
|
|
return current_hooks->tx(to_send);
|
2017-07-15 12:17:32 +08:00
|
|
|
}
|
|
|
|
|
2017-07-21 14:36:06 +08:00
|
|
|
int safety_tx_lin_hook(int lin_num, uint8_t *data, int len){
|
|
|
|
return current_hooks->tx_lin(lin_num, data, len);
|
2017-07-15 12:17:32 +08:00
|
|
|
}
|
|
|
|
|
2018-02-28 15:27:47 +08:00
|
|
|
int safety_fwd_hook(int bus_num, CAN_FIFOMailBox_TypeDef *to_fwd) {
|
|
|
|
return current_hooks->fwd(bus_num, to_fwd);
|
|
|
|
}
|
|
|
|
|
2020-02-21 05:57:07 +08:00
|
|
|
// Given a CRC-8 poly, generate a static lookup table to use with a fast CRC-8
|
|
|
|
// algorithm. Called at init time for safety modes using CRC-8.
|
|
|
|
void gen_crc_lookup_table(uint8_t poly, uint8_t crc_lut[]) {
|
|
|
|
for (int i = 0; i < 256; i++) {
|
|
|
|
uint8_t crc = i;
|
|
|
|
for (int j = 0; j < 8; j++) {
|
|
|
|
if ((crc & 0x80U) != 0U)
|
|
|
|
crc = (uint8_t)((crc << 1) ^ poly);
|
|
|
|
else
|
|
|
|
crc <<= 1;
|
|
|
|
}
|
|
|
|
crc_lut[i] = crc;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2019-12-18 16:17:25 +08:00
|
|
|
bool msg_allowed(int addr, int bus, const AddrBus addr_list[], int len) {
|
2019-11-17 16:24:19 +08:00
|
|
|
bool allowed = false;
|
|
|
|
for (int i = 0; i < len; i++) {
|
|
|
|
if ((addr == addr_list[i].addr) && (bus == addr_list[i].bus)) {
|
|
|
|
allowed = true;
|
|
|
|
break;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
return allowed;
|
|
|
|
}
|
|
|
|
|
2019-12-21 17:25:54 +08:00
|
|
|
// compute the time elapsed (in microseconds) from 2 counter samples
|
|
|
|
// case where ts < ts_last is ok: overflow is properly re-casted into uint32_t
|
|
|
|
uint32_t get_ts_elapsed(uint32_t ts, uint32_t ts_last) {
|
|
|
|
return ts - ts_last;
|
|
|
|
}
|
|
|
|
|
|
|
|
int get_addr_check_index(CAN_FIFOMailBox_TypeDef *to_push, AddrCheckStruct addr_list[], const int len) {
|
|
|
|
int bus = GET_BUS(to_push);
|
|
|
|
int addr = GET_ADDR(to_push);
|
|
|
|
|
|
|
|
int index = -1;
|
|
|
|
for (int i = 0; i < len; i++) {
|
|
|
|
for (uint8_t j = 0U; addr_list[i].addr[j] != 0; j++) {
|
|
|
|
if ((addr == addr_list[i].addr[j]) && (bus == addr_list[i].bus)) {
|
|
|
|
index = i;
|
|
|
|
goto Return;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
Return:
|
|
|
|
return index;
|
|
|
|
}
|
|
|
|
|
|
|
|
// 1Hz safety function called by main. Now just a check for lagging safety messages
|
|
|
|
void safety_tick(const safety_hooks *hooks) {
|
|
|
|
uint32_t ts = TIM2->CNT;
|
|
|
|
if (hooks->addr_check != NULL) {
|
|
|
|
for (int i=0; i < hooks->addr_check_len; i++) {
|
|
|
|
uint32_t elapsed_time = get_ts_elapsed(ts, hooks->addr_check[i].last_timestamp);
|
|
|
|
// lag threshold is max of: 1s and MAX_MISSED_MSGS * expected timestep.
|
|
|
|
// Quite conservative to not risk false triggers.
|
|
|
|
// 2s of lag is worse case, since the function is called at 1Hz
|
|
|
|
bool lagging = elapsed_time > MAX(hooks->addr_check[i].expected_timestep * MAX_MISSED_MSGS, 1e6);
|
|
|
|
hooks->addr_check[i].lagging = lagging;
|
|
|
|
if (lagging) {
|
|
|
|
controls_allowed = 0;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
void update_counter(AddrCheckStruct addr_list[], int index, uint8_t counter) {
|
|
|
|
if (index != -1) {
|
|
|
|
uint8_t expected_counter = (addr_list[index].last_counter + 1U) % (addr_list[index].max_counter + 1U);
|
|
|
|
addr_list[index].wrong_counters += (expected_counter == counter) ? -1 : 1;
|
|
|
|
addr_list[index].wrong_counters = MAX(MIN(addr_list[index].wrong_counters, MAX_WRONG_COUNTERS), 0);
|
|
|
|
addr_list[index].last_counter = counter;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
bool is_msg_valid(AddrCheckStruct addr_list[], int index) {
|
|
|
|
bool valid = true;
|
|
|
|
if (index != -1) {
|
|
|
|
if ((!addr_list[index].valid_checksum) || (addr_list[index].wrong_counters >= MAX_WRONG_COUNTERS)) {
|
|
|
|
valid = false;
|
|
|
|
controls_allowed = 0;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
return valid;
|
|
|
|
}
|
|
|
|
|
|
|
|
void update_addr_timestamp(AddrCheckStruct addr_list[], int index) {
|
|
|
|
if (index != -1) {
|
|
|
|
uint32_t ts = TIM2->CNT;
|
|
|
|
addr_list[index].last_timestamp = ts;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
bool addr_safety_check(CAN_FIFOMailBox_TypeDef *to_push,
|
|
|
|
AddrCheckStruct *rx_checks,
|
|
|
|
const int rx_checks_len,
|
|
|
|
uint8_t (*get_checksum)(CAN_FIFOMailBox_TypeDef *to_push),
|
|
|
|
uint8_t (*compute_checksum)(CAN_FIFOMailBox_TypeDef *to_push),
|
|
|
|
uint8_t (*get_counter)(CAN_FIFOMailBox_TypeDef *to_push)) {
|
|
|
|
|
|
|
|
int index = get_addr_check_index(to_push, rx_checks, rx_checks_len);
|
|
|
|
update_addr_timestamp(rx_checks, index);
|
|
|
|
|
|
|
|
if (index != -1) {
|
|
|
|
// checksum check
|
2020-01-19 05:26:34 +08:00
|
|
|
if ((get_checksum != NULL) && (compute_checksum != NULL) && rx_checks[index].check_checksum) {
|
|
|
|
uint8_t checksum = get_checksum(to_push);
|
|
|
|
uint8_t checksum_comp = compute_checksum(to_push);
|
|
|
|
rx_checks[index].valid_checksum = checksum_comp == checksum;
|
|
|
|
} else {
|
|
|
|
rx_checks[index].valid_checksum = true;
|
2019-12-21 17:25:54 +08:00
|
|
|
}
|
|
|
|
|
2020-01-19 05:26:34 +08:00
|
|
|
// counter check (max_counter == 0 means skip check)
|
|
|
|
if ((get_counter != NULL) && (rx_checks[index].max_counter > 0U)) {
|
|
|
|
uint8_t counter = get_counter(to_push);
|
|
|
|
update_counter(rx_checks, index, counter);
|
|
|
|
} else {
|
|
|
|
rx_checks[index].wrong_counters = 0U;
|
2019-12-21 17:25:54 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
return is_msg_valid(rx_checks, index);
|
|
|
|
}
|
|
|
|
|
2020-04-10 06:53:46 +08:00
|
|
|
void relay_malfunction_set(void) {
|
|
|
|
relay_malfunction = true;
|
|
|
|
fault_occurred(FAULT_RELAY_MALFUNCTION);
|
|
|
|
}
|
|
|
|
|
|
|
|
void relay_malfunction_reset(void) {
|
|
|
|
relay_malfunction = false;
|
|
|
|
fault_recovered(FAULT_RELAY_MALFUNCTION);
|
|
|
|
}
|
2019-12-21 17:25:54 +08:00
|
|
|
|
2017-07-15 12:17:32 +08:00
|
|
|
typedef struct {
|
|
|
|
uint16_t id;
|
|
|
|
const safety_hooks *hooks;
|
|
|
|
} safety_hook_config;
|
|
|
|
|
|
|
|
const safety_hook_config safety_hook_registry[] = {
|
2019-11-21 03:56:26 +08:00
|
|
|
{SAFETY_SILENT, &nooutput_hooks},
|
2019-12-21 07:51:32 +08:00
|
|
|
{SAFETY_HONDA_NIDEC, &honda_nidec_hooks},
|
2017-08-25 13:31:34 +08:00
|
|
|
{SAFETY_TOYOTA, &toyota_hooks},
|
2019-10-03 09:20:32 +08:00
|
|
|
{SAFETY_ELM327, &elm327_hooks},
|
2018-05-03 13:45:23 +08:00
|
|
|
{SAFETY_GM, &gm_hooks},
|
2019-12-21 15:18:20 +08:00
|
|
|
{SAFETY_HONDA_BOSCH_GIRAFFE, &honda_bosch_giraffe_hooks},
|
|
|
|
{SAFETY_HONDA_BOSCH_HARNESS, &honda_bosch_harness_hooks},
|
2018-08-18 12:31:00 +08:00
|
|
|
{SAFETY_HYUNDAI, &hyundai_hooks},
|
2018-11-07 04:28:33 +08:00
|
|
|
{SAFETY_CHRYSLER, &chrysler_hooks},
|
2019-01-29 13:42:06 +08:00
|
|
|
{SAFETY_SUBARU, &subaru_hooks},
|
2020-02-28 07:28:24 +08:00
|
|
|
{SAFETY_SUBARU_LEGACY, &subaru_legacy_hooks},
|
2019-09-10 06:58:41 +08:00
|
|
|
{SAFETY_MAZDA, &mazda_hooks},
|
2020-02-20 04:37:38 +08:00
|
|
|
{SAFETY_VOLKSWAGEN_MQB, &volkswagen_mqb_hooks},
|
2020-04-02 03:38:27 +08:00
|
|
|
{SAFETY_VOLKSWAGEN_PQ, &volkswagen_pq_hooks},
|
2019-11-21 03:56:26 +08:00
|
|
|
{SAFETY_NOOUTPUT, &nooutput_hooks},
|
2019-11-16 04:20:15 +08:00
|
|
|
#ifdef ALLOW_DEBUG
|
|
|
|
{SAFETY_TESLA, &tesla_hooks},
|
2020-02-27 03:11:56 +08:00
|
|
|
{SAFETY_NISSAN, &nissan_hooks},
|
2017-07-18 01:48:16 +08:00
|
|
|
{SAFETY_ALLOUTPUT, &alloutput_hooks},
|
2019-10-03 09:20:32 +08:00
|
|
|
{SAFETY_GM_ASCM, &gm_ascm_hooks},
|
2019-11-16 04:20:15 +08:00
|
|
|
{SAFETY_FORD, &ford_hooks},
|
|
|
|
#endif
|
2017-07-15 12:17:32 +08:00
|
|
|
};
|
|
|
|
|
2019-11-07 06:30:22 +08:00
|
|
|
int set_safety_hooks(uint16_t mode, int16_t param) {
|
2019-11-27 13:19:54 +08:00
|
|
|
safety_mode_cnt = 0U; // reset safety mode timer
|
|
|
|
int set_status = -1; // not set
|
2019-06-15 01:57:14 +08:00
|
|
|
int hook_config_count = sizeof(safety_hook_registry) / sizeof(safety_hook_config);
|
|
|
|
for (int i = 0; i < hook_config_count; i++) {
|
2017-07-18 01:27:34 +08:00
|
|
|
if (safety_hook_registry[i].id == mode) {
|
2017-07-15 12:17:32 +08:00
|
|
|
current_hooks = safety_hook_registry[i].hooks;
|
2019-10-05 04:30:00 +08:00
|
|
|
current_safety_mode = safety_hook_registry[i].id;
|
2019-11-27 13:19:54 +08:00
|
|
|
set_status = 0; // set
|
2019-06-13 09:18:07 +08:00
|
|
|
break;
|
2017-07-15 12:17:32 +08:00
|
|
|
}
|
|
|
|
}
|
2019-06-13 09:18:07 +08:00
|
|
|
if ((set_status == 0) && (current_hooks->init != NULL)) {
|
|
|
|
current_hooks->init(param);
|
|
|
|
}
|
|
|
|
return set_status;
|
2017-07-15 12:17:32 +08:00
|
|
|
}
|
2017-08-24 03:49:56 +08:00
|
|
|
|
2018-05-26 08:52:37 +08:00
|
|
|
// convert a trimmed integer to signed 32 bit int
|
|
|
|
int to_signed(int d, int bits) {
|
2019-06-11 15:58:22 +08:00
|
|
|
int d_signed = d;
|
2019-06-18 04:48:14 +08:00
|
|
|
if (d >= (1 << MAX((bits - 1), 0))) {
|
|
|
|
d_signed = d - (1 << MAX(bits, 0));
|
2018-05-26 08:52:37 +08:00
|
|
|
}
|
2019-06-11 15:58:22 +08:00
|
|
|
return d_signed;
|
2018-05-26 08:52:37 +08:00
|
|
|
}
|
2018-06-01 09:49:12 +08:00
|
|
|
|
|
|
|
// given a new sample, update the smaple_t struct
|
|
|
|
void update_sample(struct sample_t *sample, int sample_new) {
|
2019-06-15 01:57:14 +08:00
|
|
|
int sample_size = sizeof(sample->values) / sizeof(sample->values[0]);
|
|
|
|
for (int i = sample_size - 1; i > 0; i--) {
|
2018-06-01 09:49:12 +08:00
|
|
|
sample->values[i] = sample->values[i-1];
|
|
|
|
}
|
|
|
|
sample->values[0] = sample_new;
|
|
|
|
|
|
|
|
// get the minimum and maximum measured samples
|
2019-06-11 16:10:45 +08:00
|
|
|
sample->min = sample->values[0];
|
|
|
|
sample->max = sample->values[0];
|
2019-06-15 01:57:14 +08:00
|
|
|
for (int i = 1; i < sample_size; i++) {
|
2019-06-11 16:07:41 +08:00
|
|
|
if (sample->values[i] < sample->min) {
|
|
|
|
sample->min = sample->values[i];
|
|
|
|
}
|
|
|
|
if (sample->values[i] > sample->max) {
|
|
|
|
sample->max = sample->values[i];
|
|
|
|
}
|
2018-06-01 09:49:12 +08:00
|
|
|
}
|
|
|
|
}
|
2018-06-13 15:56:29 +08:00
|
|
|
|
2019-06-18 09:48:38 +08:00
|
|
|
bool max_limit_check(int val, const int MAX_VAL, const int MIN_VAL) {
|
|
|
|
return (val > MAX_VAL) || (val < MIN_VAL);
|
2018-06-13 16:37:36 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
// check that commanded value isn't too far from measured
|
2019-06-13 11:12:48 +08:00
|
|
|
bool dist_to_meas_check(int val, int val_last, struct sample_t *val_meas,
|
2018-06-13 16:37:36 +08:00
|
|
|
const int MAX_RATE_UP, const int MAX_RATE_DOWN, const int MAX_ERROR) {
|
|
|
|
|
|
|
|
// *** val rate limit check ***
|
2019-06-18 06:29:04 +08:00
|
|
|
int highest_allowed_rl = MAX(val_last, 0) + MAX_RATE_UP;
|
|
|
|
int lowest_allowed_rl = MIN(val_last, 0) - MAX_RATE_UP;
|
2018-06-13 16:37:36 +08:00
|
|
|
|
|
|
|
// if we've exceeded the meas val, we must start moving toward 0
|
2019-06-18 06:29:04 +08:00
|
|
|
int highest_allowed = MIN(highest_allowed_rl, MAX(val_last - MAX_RATE_DOWN, MAX(val_meas->max, 0) + MAX_ERROR));
|
|
|
|
int lowest_allowed = MAX(lowest_allowed_rl, MIN(val_last + MAX_RATE_DOWN, MIN(val_meas->min, 0) - MAX_ERROR));
|
2018-06-13 16:37:36 +08:00
|
|
|
|
|
|
|
// check for violation
|
2019-06-18 06:29:04 +08:00
|
|
|
return (val < lowest_allowed) || (val > highest_allowed);
|
2018-06-13 16:37:36 +08:00
|
|
|
}
|
|
|
|
|
2018-06-14 06:23:56 +08:00
|
|
|
// check that commanded value isn't fighting against driver
|
2019-06-13 11:12:48 +08:00
|
|
|
bool driver_limit_check(int val, int val_last, struct sample_t *val_driver,
|
2019-06-18 09:48:38 +08:00
|
|
|
const int MAX_VAL, const int MAX_RATE_UP, const int MAX_RATE_DOWN,
|
2018-06-14 06:23:56 +08:00
|
|
|
const int MAX_ALLOWANCE, const int DRIVER_FACTOR) {
|
|
|
|
|
2019-06-18 06:29:04 +08:00
|
|
|
int highest_allowed_rl = MAX(val_last, 0) + MAX_RATE_UP;
|
|
|
|
int lowest_allowed_rl = MIN(val_last, 0) - MAX_RATE_UP;
|
2018-06-14 06:23:56 +08:00
|
|
|
|
2019-06-18 09:48:38 +08:00
|
|
|
int driver_max_limit = MAX_VAL + (MAX_ALLOWANCE + val_driver->max) * DRIVER_FACTOR;
|
|
|
|
int driver_min_limit = -MAX_VAL + (-MAX_ALLOWANCE + val_driver->min) * DRIVER_FACTOR;
|
2018-06-14 06:23:56 +08:00
|
|
|
|
|
|
|
// if we've exceeded the applied torque, we must start moving toward 0
|
2019-06-18 06:29:04 +08:00
|
|
|
int highest_allowed = MIN(highest_allowed_rl, MAX(val_last - MAX_RATE_DOWN,
|
2019-06-18 04:48:14 +08:00
|
|
|
MAX(driver_max_limit, 0)));
|
2019-06-18 06:29:04 +08:00
|
|
|
int lowest_allowed = MAX(lowest_allowed_rl, MIN(val_last + MAX_RATE_DOWN,
|
2019-06-18 04:48:14 +08:00
|
|
|
MIN(driver_min_limit, 0)));
|
2018-06-14 06:23:56 +08:00
|
|
|
|
|
|
|
// check for violation
|
|
|
|
return (val < lowest_allowed) || (val > highest_allowed);
|
|
|
|
}
|
|
|
|
|
|
|
|
|
2018-06-13 15:56:29 +08:00
|
|
|
// real time check, mainly used for steer torque rate limiter
|
2019-06-13 11:12:48 +08:00
|
|
|
bool rt_rate_limit_check(int val, int val_last, const int MAX_RT_DELTA) {
|
2018-06-13 15:56:29 +08:00
|
|
|
|
2018-06-13 15:57:36 +08:00
|
|
|
// *** torque real time rate limit check ***
|
2019-06-18 04:48:14 +08:00
|
|
|
int highest_val = MAX(val_last, 0) + MAX_RT_DELTA;
|
|
|
|
int lowest_val = MIN(val_last, 0) - MAX_RT_DELTA;
|
2018-06-13 15:56:29 +08:00
|
|
|
|
2018-06-13 16:37:36 +08:00
|
|
|
// check for violation
|
2018-06-13 15:57:36 +08:00
|
|
|
return (val < lowest_val) || (val > highest_val);
|
2018-06-13 15:56:29 +08:00
|
|
|
}
|
2018-08-14 13:52:31 +08:00
|
|
|
|
|
|
|
|
|
|
|
// interp function that holds extreme values
|
|
|
|
float interpolate(struct lookup_t xy, float x) {
|
2019-06-13 09:18:07 +08:00
|
|
|
|
2018-08-14 13:52:31 +08:00
|
|
|
int size = sizeof(xy.x) / sizeof(xy.x[0]);
|
2019-06-13 09:18:07 +08:00
|
|
|
float ret = xy.y[size - 1]; // default output is last point
|
|
|
|
|
2018-08-14 13:52:31 +08:00
|
|
|
// x is lower than the first point in the x array. Return the first point
|
|
|
|
if (x <= xy.x[0]) {
|
2019-06-13 09:18:07 +08:00
|
|
|
ret = xy.y[0];
|
2018-08-14 13:52:31 +08:00
|
|
|
|
|
|
|
} else {
|
|
|
|
// find the index such that (xy.x[i] <= x < xy.x[i+1]) and linearly interp
|
2019-06-12 12:23:48 +08:00
|
|
|
for (int i=0; i < (size - 1); i++) {
|
2018-08-14 13:52:31 +08:00
|
|
|
if (x < xy.x[i+1]) {
|
|
|
|
float x0 = xy.x[i];
|
|
|
|
float y0 = xy.y[i];
|
|
|
|
float dx = xy.x[i+1] - x0;
|
|
|
|
float dy = xy.y[i+1] - y0;
|
2019-12-13 10:23:03 +08:00
|
|
|
// dx should not be zero as xy.x is supposed to be monotonic
|
2019-06-11 16:07:41 +08:00
|
|
|
if (dx <= 0.) {
|
|
|
|
dx = 0.0001;
|
|
|
|
}
|
2019-06-13 09:18:07 +08:00
|
|
|
ret = (dy * (x - x0) / dx) + y0;
|
|
|
|
break;
|
2018-08-14 13:52:31 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
2019-06-13 09:18:07 +08:00
|
|
|
return ret;
|
2018-08-14 13:52:31 +08:00
|
|
|
}
|