ahrs_int_cmpl_euler_wrapper.c

Go to the documentation of this file.

/*
 * Copyright (C) 2015 Felix Ruess <felix.ruess@gmail.com>
 *
 * This file is part of paparazzi.
 *
 * paparazzi is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2, or (at your option)
 * any later version.
 *
 * paparazzi is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with Paparazzi; see the file COPYING.  If not, see
 * <http://www.gnu.org/licenses/>.
 */

#include "subsystems/ahrs/ahrs_int_cmpl_euler_wrapper.h"
#include "subsystems/ahrs.h"
#include "subsystems/abi.h"
#include "state.h"

#ifndef AHRS_ICE_OUTPUT_ENABLED
#define AHRS_ICE_OUTPUT_ENABLED TRUE
#endif
PRINT_CONFIG_VAR(AHRS_ICE_OUTPUT_ENABLED)


static bool ahrs_ice_output_enabled;
static uint32_t ahrs_ice_last_stamp;
static uint8_t ahrs_ice_id = AHRS_COMP_ID_ICE;

static void set_body_state_from_euler(void);

#if PERIODIC_TELEMETRY
#include "subsystems/datalink/telemetry.h"
#include "mcu_periph/sys_time.h"
#include "state.h"

static void send_filter(struct transport_tx *trans, struct link_device *dev)
{
  pprz_msg_send_FILTER(trans, dev, AC_ID,
                       &ahrs_ice.ltp_to_imu_euler.phi,
                       &ahrs_ice.ltp_to_imu_euler.theta,
                       &ahrs_ice.ltp_to_imu_euler.psi,
                       &ahrs_ice.measure.phi,
                       &ahrs_ice.measure.theta,
                       &ahrs_ice.measure.psi,
                       &ahrs_ice.hi_res_euler.phi,
                       &ahrs_ice.hi_res_euler.theta,
                       &ahrs_ice.hi_res_euler.psi,
                       &ahrs_ice.residual.phi,
                       &ahrs_ice.residual.theta,
                       &ahrs_ice.residual.psi,
                       &ahrs_ice.gyro_bias.p,
                       &ahrs_ice.gyro_bias.q,
                       &ahrs_ice.gyro_bias.r,
                       &ahrs_ice_id);
}

static void send_euler(struct transport_tx *trans, struct link_device *dev)
{
  struct Int32Eulers *eulers = stateGetNedToBodyEulers_i();
  pprz_msg_send_AHRS_EULER_INT(trans, dev, AC_ID,
                               &ahrs_ice.ltp_to_imu_euler.phi,
                               &ahrs_ice.ltp_to_imu_euler.theta,
                               &ahrs_ice.ltp_to_imu_euler.psi,
                               &(eulers->phi),
                               &(eulers->theta),
                               &(eulers->psi),
                               &ahrs_ice_id);
}

static void send_bias(struct transport_tx *trans, struct link_device *dev)
{
  pprz_msg_send_AHRS_GYRO_BIAS_INT(trans, dev, AC_ID,
                                   &ahrs_ice.gyro_bias.p, &ahrs_ice.gyro_bias.q,
                                   &ahrs_ice.gyro_bias.r, &ahrs_ice_id);
}

static void send_filter_status(struct transport_tx *trans, struct link_device *dev)
{
  uint8_t mde = 3;
  uint16_t val = 0;
  if (!ahrs_ice.is_aligned) { mde = 2; }
  uint32_t t_diff = get_sys_time_usec() - ahrs_ice_last_stamp;
  /* set lost if no new gyro measurements for 50ms */
  if (t_diff > 50000) { mde = 5; }
  pprz_msg_send_STATE_FILTER_STATUS(trans, dev, AC_ID, &ahrs_ice_id, &mde, &val);
}
#endif

#ifndef AHRS_ICE_IMU_ID
#define AHRS_ICE_IMU_ID ABI_BROADCAST
#endif
PRINT_CONFIG_VAR(AHRS_ICE_IMU_ID)
#ifndef AHRS_ICE_MAG_ID
#define AHRS_ICE_MAG_ID AHRS_ICE_IMU_ID
#endif
PRINT_CONFIG_VAR(AHRS_ICE_MAG_ID)
static abi_event gyro_ev;
static abi_event accel_ev;
static abi_event mag_ev;
static abi_event aligner_ev;
static abi_event body_to_imu_ev;


static void gyro_cb(uint8_t sender_id __attribute__((unused)),
                    uint32_t stamp, struct Int32Rates *gyro)
{
  ahrs_ice_last_stamp = stamp;
  if (ahrs_ice.is_aligned) {
    ahrs_ice_propagate(gyro);
    set_body_state_from_euler();
  }
}

static void accel_cb(uint8_t sender_id __attribute__((unused)),
                     uint32_t stamp __attribute__((unused)),
                     struct Int32Vect3 *accel)
{
  if (ahrs_ice.is_aligned) {
    ahrs_ice_update_accel(accel);
  }
}

static void mag_cb(uint8_t sender_id __attribute__((unused)),
                   uint32_t stamp __attribute__((unused)),
                   struct Int32Vect3 *mag)
{
  if (ahrs_ice.is_aligned) {
    ahrs_ice_update_mag(mag);
  }
}

static void aligner_cb(uint8_t __attribute__((unused)) sender_id,
                       uint32_t stamp __attribute__((unused)),
                       struct Int32Rates *lp_gyro, struct Int32Vect3 *lp_accel,
                       struct Int32Vect3 *lp_mag)
{
  if (!ahrs_ice.is_aligned) {
    if (ahrs_ice_align(lp_gyro, lp_accel, lp_mag)) {
      set_body_state_from_euler();
    }
  }
}

static void body_to_imu_cb(uint8_t sender_id __attribute__((unused)),
                           struct FloatQuat *q_b2i_f)
{
  ahrs_ice_set_body_to_imu_quat(q_b2i_f);
}

static bool ahrs_ice_enable_output(bool enable)
{
  ahrs_ice_output_enabled = enable;
  return ahrs_ice_output_enabled;
}

/* Rotate angles and rates from imu to body frame and set state */
static void set_body_state_from_euler(void)
{
  if (ahrs_ice_output_enabled) {
    struct Int32RMat *body_to_imu_rmat = orientationGetRMat_i(&ahrs_ice.body_to_imu);
    struct Int32RMat ltp_to_imu_rmat, ltp_to_body_rmat;
    /* Compute LTP to IMU rotation matrix */
    int32_rmat_of_eulers(&ltp_to_imu_rmat, &ahrs_ice.ltp_to_imu_euler);
    /* Compute LTP to BODY rotation matrix */
    int32_rmat_comp_inv(&ltp_to_body_rmat, &ltp_to_imu_rmat, body_to_imu_rmat);
    /* Set state */
    stateSetNedToBodyRMat_i(&ltp_to_body_rmat);

    struct Int32Rates body_rate;
    /* compute body rates */
    int32_rmat_transp_ratemult(&body_rate, body_to_imu_rmat, &ahrs_ice.imu_rate);
    /* Set state */
    stateSetBodyRates_i(&body_rate);
  }
}

void ahrs_ice_register(void)
{
  ahrs_ice_output_enabled = AHRS_ICE_OUTPUT_ENABLED;
  ahrs_ice_init();
  ahrs_register_impl(ahrs_ice_enable_output);

  /*
   * Subscribe to scaled IMU measurements and attach callbacks
   */
  AbiBindMsgIMU_GYRO_INT32(AHRS_ICE_IMU_ID, &gyro_ev, gyro_cb);
  AbiBindMsgIMU_ACCEL_INT32(AHRS_ICE_IMU_ID, &accel_ev, accel_cb);
  AbiBindMsgIMU_MAG_INT32(AHRS_ICE_MAG_ID, &mag_ev, mag_cb);
  AbiBindMsgIMU_LOWPASSED(ABI_BROADCAST, &aligner_ev, aligner_cb);
  AbiBindMsgBODY_TO_IMU_QUAT(ABI_BROADCAST, &body_to_imu_ev, body_to_imu_cb);

#if PERIODIC_TELEMETRY
  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_FILTER, send_filter);
  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_AHRS_EULER_INT, send_euler);
  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_AHRS_GYRO_BIAS_INT, send_bias);
  register_periodic_telemetry(DefaultPeriodic, PPRZ_MSG_ID_STATE_FILTER_STATUS, send_filter_status);
#endif
}