mirror of
https://gitee.com/mirrors_PX4/PX4-Autopilot.git
synced 2026-05-17 06:57:34 +08:00
Expand auto-format coverage and tiny style changes
This commit is contained in:
committed by
Mathieu Bresciani
parent
fbdd75da2e
commit
f20726d47f
+3
-3
@@ -76,7 +76,7 @@ IndentWidth: 8 # Modified
|
|||||||
IndentWrappedFunctionNames: false
|
IndentWrappedFunctionNames: false
|
||||||
JavaScriptQuotes: Leave
|
JavaScriptQuotes: Leave
|
||||||
JavaScriptWrapImports: true
|
JavaScriptWrapImports: true
|
||||||
KeepEmptyLinesAtTheStartOfBlocks: false
|
KeepEmptyLinesAtTheStartOfBlocks: true # Modified
|
||||||
MacroBlockBegin: ''
|
MacroBlockBegin: ''
|
||||||
MacroBlockEnd: ''
|
MacroBlockEnd: ''
|
||||||
MaxEmptyLinesToKeep: 1
|
MaxEmptyLinesToKeep: 1
|
||||||
@@ -84,8 +84,8 @@ NamespaceIndentation: None
|
|||||||
ObjCBlockIndentWidth: 2
|
ObjCBlockIndentWidth: 2
|
||||||
ObjCSpaceAfterProperty: false
|
ObjCSpaceAfterProperty: false
|
||||||
ObjCSpaceBeforeProtocolList: false
|
ObjCSpaceBeforeProtocolList: false
|
||||||
PenaltyBreakAssignment: 2
|
PenaltyBreakAssignment: 200 # Modified
|
||||||
PenaltyBreakBeforeFirstCallParameter: 1
|
PenaltyBreakBeforeFirstCallParameter: 20 # Modified
|
||||||
PenaltyBreakComment: 300
|
PenaltyBreakComment: 300
|
||||||
PenaltyBreakFirstLessLess: 120
|
PenaltyBreakFirstLessLess: 120
|
||||||
PenaltyBreakString: 1000
|
PenaltyBreakString: 1000
|
||||||
|
|||||||
@@ -60,6 +60,7 @@ public:
|
|||||||
RingBuffer &operator=(RingBuffer &&) = delete;
|
RingBuffer &operator=(RingBuffer &&) = delete;
|
||||||
|
|
||||||
bool allocate(uint8_t size) {
|
bool allocate(uint8_t size) {
|
||||||
|
|
||||||
if (_buffer != nullptr) {
|
if (_buffer != nullptr) {
|
||||||
delete[] _buffer;
|
delete[] _buffer;
|
||||||
}
|
}
|
||||||
@@ -92,6 +93,7 @@ public:
|
|||||||
}
|
}
|
||||||
|
|
||||||
void push(const data_type &sample) {
|
void push(const data_type &sample) {
|
||||||
|
|
||||||
uint8_t head_new = _head;
|
uint8_t head_new = _head;
|
||||||
|
|
||||||
if (!_first_write) {
|
if (!_first_write) {
|
||||||
|
|||||||
@@ -1,24 +1,18 @@
|
|||||||
#include "imu_down_sampler.hpp"
|
#include "imu_down_sampler.hpp"
|
||||||
|
|
||||||
ImuDownSampler::ImuDownSampler(float target_dt_sec):
|
ImuDownSampler::ImuDownSampler(float target_dt_sec) : _target_dt{target_dt_sec}, _imu_collection_time_adj{0.0f} {
|
||||||
_target_dt{target_dt_sec},
|
|
||||||
_imu_collection_time_adj{0.0f}
|
|
||||||
{
|
|
||||||
reset();
|
reset();
|
||||||
_imu_down_sampled.time_us = 0.0f;
|
_imu_down_sampled.time_us = 0.0f;
|
||||||
}
|
}
|
||||||
|
|
||||||
ImuDownSampler::~ImuDownSampler()
|
ImuDownSampler::~ImuDownSampler() {}
|
||||||
{
|
|
||||||
}
|
|
||||||
|
|
||||||
// integrate imu samples until target dt reached
|
// integrate imu samples until target dt reached
|
||||||
// assumes that dt of the gyroscope is close to the dt of the accelerometer
|
// assumes that dt of the gyroscope is close to the dt of the accelerometer
|
||||||
// returns true if target dt is reached
|
// returns true if target dt is reached
|
||||||
bool ImuDownSampler::update(imuSample imu_sample_new)
|
bool ImuDownSampler::update(imuSample imu_sample_new) {
|
||||||
{
|
|
||||||
if(_do_reset)
|
if (_do_reset) {
|
||||||
{
|
|
||||||
reset();
|
reset();
|
||||||
}
|
}
|
||||||
// accumulate time deltas
|
// accumulate time deltas
|
||||||
@@ -42,11 +36,11 @@ bool ImuDownSampler::update(imuSample imu_sample_new)
|
|||||||
|
|
||||||
// check if the target time delta between filter prediction steps has been exceeded
|
// check if the target time delta between filter prediction steps has been exceeded
|
||||||
if (_imu_down_sampled.delta_ang_dt >= _target_dt - _imu_collection_time_adj) {
|
if (_imu_down_sampled.delta_ang_dt >= _target_dt - _imu_collection_time_adj) {
|
||||||
|
|
||||||
// accumulate the amount of time to advance the IMU collection time so that we meet the
|
// accumulate the amount of time to advance the IMU collection time so that we meet the
|
||||||
// average EKF update rate requirement
|
// average EKF update rate requirement
|
||||||
_imu_collection_time_adj += 0.01f * (_imu_down_sampled.delta_ang_dt - _target_dt);
|
_imu_collection_time_adj += 0.01f * (_imu_down_sampled.delta_ang_dt - _target_dt);
|
||||||
_imu_collection_time_adj = math::constrain(_imu_collection_time_adj, -0.5f * _target_dt, 0.5f * _target_dt);
|
_imu_collection_time_adj = math::constrain(_imu_collection_time_adj, -0.5f * _target_dt,
|
||||||
|
0.5f * _target_dt);
|
||||||
|
|
||||||
_imu_down_sampled.delta_ang = _delta_angle_accumulated.to_axis_angle();
|
_imu_down_sampled.delta_ang = _delta_angle_accumulated.to_axis_angle();
|
||||||
|
|
||||||
@@ -57,14 +51,12 @@ bool ImuDownSampler::update(imuSample imu_sample_new)
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
imuSample ImuDownSampler::getDownSampledImuAndTriggerReset()
|
imuSample ImuDownSampler::getDownSampledImuAndTriggerReset() {
|
||||||
{
|
|
||||||
_do_reset = true;
|
_do_reset = true;
|
||||||
return _imu_down_sampled;
|
return _imu_down_sampled;
|
||||||
}
|
}
|
||||||
|
|
||||||
void ImuDownSampler::reset()
|
void ImuDownSampler::reset() {
|
||||||
{
|
|
||||||
_imu_down_sampled.delta_ang.setZero();
|
_imu_down_sampled.delta_ang.setZero();
|
||||||
_imu_down_sampled.delta_vel.setZero();
|
_imu_down_sampled.delta_vel.setZero();
|
||||||
_imu_down_sampled.delta_ang_dt = 0.0f;
|
_imu_down_sampled.delta_ang_dt = 0.0f;
|
||||||
|
|||||||
@@ -37,15 +37,14 @@
|
|||||||
*/
|
*/
|
||||||
#pragma once
|
#pragma once
|
||||||
|
|
||||||
#include <matrix/math.hpp>
|
|
||||||
#include <mathlib/mathlib.h>
|
#include <mathlib/mathlib.h>
|
||||||
|
#include <matrix/math.hpp>
|
||||||
|
|
||||||
#include "common.h"
|
#include "common.h"
|
||||||
|
|
||||||
using namespace estimator;
|
using namespace estimator;
|
||||||
|
|
||||||
class ImuDownSampler
|
class ImuDownSampler {
|
||||||
{
|
|
||||||
public:
|
public:
|
||||||
ImuDownSampler(float target_dt_sec);
|
ImuDownSampler(float target_dt_sec);
|
||||||
~ImuDownSampler();
|
~ImuDownSampler();
|
||||||
@@ -56,11 +55,9 @@ public:
|
|||||||
private:
|
private:
|
||||||
imuSample _imu_down_sampled;
|
imuSample _imu_down_sampled;
|
||||||
Quatf _delta_angle_accumulated;
|
Quatf _delta_angle_accumulated;
|
||||||
const float _target_dt; // [sec]
|
const float _target_dt; // [sec]
|
||||||
float _imu_collection_time_adj;
|
float _imu_collection_time_adj;
|
||||||
bool _do_reset;
|
bool _do_reset;
|
||||||
|
|
||||||
void reset();
|
void reset();
|
||||||
|
|
||||||
};
|
};
|
||||||
|
|
||||||
|
|||||||
+38
-43
@@ -41,39 +41,37 @@
|
|||||||
*
|
*
|
||||||
*/
|
*/
|
||||||
|
|
||||||
#include "ekf.h"
|
|
||||||
#include <ecl.h>
|
#include <ecl.h>
|
||||||
#include <mathlib/mathlib.h>
|
#include <mathlib/mathlib.h>
|
||||||
|
#include "ekf.h"
|
||||||
|
|
||||||
|
bool Ekf::fuseHorizontalVelocity(const Vector3f &innov, const Vector2f &innov_gate, const Vector3f &obs_var,
|
||||||
|
Vector3f &innov_var, Vector2f &test_ratio) {
|
||||||
|
|
||||||
bool Ekf::fuseHorizontalVelocity(const Vector3f &innov, const Vector2f &innov_gate,
|
innov_var(0) = P(4, 4) + obs_var(0);
|
||||||
const Vector3f &obs_var, Vector3f &innov_var, Vector2f &test_ratio)
|
innov_var(1) = P(5, 5) + obs_var(1);
|
||||||
{
|
|
||||||
innov_var(0) = P(4,4) + obs_var(0);
|
|
||||||
innov_var(1) = P(5,5) + obs_var(1);
|
|
||||||
test_ratio(0) = fmaxf(sq(innov(0)) / (sq(innov_gate(0)) * innov_var(0)),
|
test_ratio(0) = fmaxf(sq(innov(0)) / (sq(innov_gate(0)) * innov_var(0)),
|
||||||
sq(innov(1)) / (sq(innov_gate(0)) * innov_var(1)));
|
sq(innov(1)) / (sq(innov_gate(0)) * innov_var(1)));
|
||||||
|
|
||||||
const bool innov_check_pass = (test_ratio(0) <= 1.0f);
|
const bool innov_check_pass = (test_ratio(0) <= 1.0f);
|
||||||
if (innov_check_pass)
|
if (innov_check_pass) {
|
||||||
{
|
|
||||||
_time_last_hor_vel_fuse = _time_last_imu;
|
_time_last_hor_vel_fuse = _time_last_imu;
|
||||||
_innov_check_fail_status.flags.reject_hor_vel = false;
|
_innov_check_fail_status.flags.reject_hor_vel = false;
|
||||||
|
|
||||||
fuseVelPosHeight(innov(0),innov_var(0),0);
|
fuseVelPosHeight(innov(0), innov_var(0), 0);
|
||||||
fuseVelPosHeight(innov(1),innov_var(1),1);
|
fuseVelPosHeight(innov(1), innov_var(1), 1);
|
||||||
|
|
||||||
return true;
|
return true;
|
||||||
}else{
|
} else {
|
||||||
_innov_check_fail_status.flags.reject_hor_vel = true;
|
_innov_check_fail_status.flags.reject_hor_vel = true;
|
||||||
return false;
|
return false;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
bool Ekf::fuseVerticalVelocity(const Vector3f &innov, const Vector2f &innov_gate,
|
bool Ekf::fuseVerticalVelocity(const Vector3f &innov, const Vector2f &innov_gate, const Vector3f &obs_var,
|
||||||
const Vector3f &obs_var, Vector3f &innov_var, Vector2f &test_ratio)
|
Vector3f &innov_var, Vector2f &test_ratio) {
|
||||||
{
|
|
||||||
innov_var(2) = P(6,6) + obs_var(2);
|
innov_var(2) = P(6, 6) + obs_var(2);
|
||||||
test_ratio(1) = sq(innov(2)) / (sq(innov_gate(1)) * innov_var(2));
|
test_ratio(1) = sq(innov(2)) / (sq(innov_gate(1)) * innov_var(2));
|
||||||
|
|
||||||
const bool innov_check_pass = (test_ratio(1) <= 1.0f);
|
const bool innov_check_pass = (test_ratio(1) <= 1.0f);
|
||||||
@@ -81,20 +79,20 @@ bool Ekf::fuseVerticalVelocity(const Vector3f &innov, const Vector2f &innov_gate
|
|||||||
_time_last_ver_vel_fuse = _time_last_imu;
|
_time_last_ver_vel_fuse = _time_last_imu;
|
||||||
_innov_check_fail_status.flags.reject_ver_vel = false;
|
_innov_check_fail_status.flags.reject_ver_vel = false;
|
||||||
|
|
||||||
fuseVelPosHeight(innov(2),innov_var(2),2);
|
fuseVelPosHeight(innov(2), innov_var(2), 2);
|
||||||
|
|
||||||
return true;
|
return true;
|
||||||
}else{
|
} else {
|
||||||
_innov_check_fail_status.flags.reject_ver_vel = true;
|
_innov_check_fail_status.flags.reject_ver_vel = true;
|
||||||
return false;
|
return false;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
bool Ekf::fuseHorizontalPosition(const Vector3f &innov, const Vector2f &innov_gate,
|
bool Ekf::fuseHorizontalPosition(const Vector3f &innov, const Vector2f &innov_gate, const Vector3f &obs_var,
|
||||||
const Vector3f &obs_var, Vector3f &innov_var, Vector2f &test_ratio)
|
Vector3f &innov_var, Vector2f &test_ratio) {
|
||||||
{
|
|
||||||
innov_var(0) = P(7,7) + obs_var(0);
|
innov_var(0) = P(7, 7) + obs_var(0);
|
||||||
innov_var(1) = P(8,8) + obs_var(1);
|
innov_var(1) = P(8, 8) + obs_var(1);
|
||||||
test_ratio(0) = fmaxf(sq(innov(0)) / (sq(innov_gate(0)) * innov_var(0)),
|
test_ratio(0) = fmaxf(sq(innov(0)) / (sq(innov_gate(0)) * innov_var(0)),
|
||||||
sq(innov(1)) / (sq(innov_gate(0)) * innov_var(1)));
|
sq(innov(1)) / (sq(innov_gate(0)) * innov_var(1)));
|
||||||
|
|
||||||
@@ -108,20 +106,20 @@ bool Ekf::fuseHorizontalPosition(const Vector3f &innov, const Vector2f &innov_ga
|
|||||||
}
|
}
|
||||||
_innov_check_fail_status.flags.reject_hor_pos = false;
|
_innov_check_fail_status.flags.reject_hor_pos = false;
|
||||||
|
|
||||||
fuseVelPosHeight(innov(0),innov_var(0),3);
|
fuseVelPosHeight(innov(0), innov_var(0), 3);
|
||||||
fuseVelPosHeight(innov(1),innov_var(1),4);
|
fuseVelPosHeight(innov(1), innov_var(1), 4);
|
||||||
|
|
||||||
return true;
|
return true;
|
||||||
}else{
|
} else {
|
||||||
_innov_check_fail_status.flags.reject_hor_pos = true;
|
_innov_check_fail_status.flags.reject_hor_pos = true;
|
||||||
return false;
|
return false;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
bool Ekf::fuseVerticalPosition(const Vector3f &innov, const Vector2f &innov_gate,
|
bool Ekf::fuseVerticalPosition(const Vector3f &innov, const Vector2f &innov_gate, const Vector3f &obs_var,
|
||||||
const Vector3f &obs_var, Vector3f &innov_var, Vector2f &test_ratio)
|
Vector3f &innov_var, Vector2f &test_ratio) {
|
||||||
{
|
|
||||||
innov_var(2) = P(9,9) + obs_var(2);
|
innov_var(2) = P(9, 9) + obs_var(2);
|
||||||
test_ratio(1) = sq(innov(2)) / (sq(innov_gate(1)) * innov_var(2));
|
test_ratio(1) = sq(innov(2)) / (sq(innov_gate(1)) * innov_var(2));
|
||||||
|
|
||||||
const bool innov_check_pass = (test_ratio(1) <= 1.0f) || !_control_status.flags.tilt_align;
|
const bool innov_check_pass = (test_ratio(1) <= 1.0f) || !_control_status.flags.tilt_align;
|
||||||
@@ -129,31 +127,31 @@ bool Ekf::fuseVerticalPosition(const Vector3f &innov, const Vector2f &innov_gate
|
|||||||
_time_last_hgt_fuse = _time_last_imu;
|
_time_last_hgt_fuse = _time_last_imu;
|
||||||
_innov_check_fail_status.flags.reject_ver_pos = false;
|
_innov_check_fail_status.flags.reject_ver_pos = false;
|
||||||
|
|
||||||
fuseVelPosHeight(innov(2),innov_var(2),5);
|
fuseVelPosHeight(innov(2), innov_var(2), 5);
|
||||||
|
|
||||||
return true;
|
return true;
|
||||||
}else{
|
} else {
|
||||||
_innov_check_fail_status.flags.reject_ver_pos = true;
|
_innov_check_fail_status.flags.reject_ver_pos = true;
|
||||||
return false;
|
return false;
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
// Helper function that fuses a single velocity or position measurement
|
// Helper function that fuses a single velocity or position measurement
|
||||||
void Ekf::fuseVelPosHeight(const float innov, const float innov_var, const int obs_index)
|
void Ekf::fuseVelPosHeight(const float innov, const float innov_var, const int obs_index) {
|
||||||
{
|
|
||||||
float Kfusion[24] = {}; // Kalman gain vector for any single observation - sequential fusion is used.
|
float Kfusion[24] = {}; // Kalman gain vector for any single observation - sequential fusion is used.
|
||||||
const unsigned state_index = obs_index + 4; // we start with vx and this is the 4. state
|
const unsigned state_index = obs_index + 4; // we start with vx and this is the 4. state
|
||||||
|
|
||||||
// calculate kalman gain K = PHS, where S = 1/innovation variance
|
// calculate kalman gain K = PHS, where S = 1/innovation variance
|
||||||
for (int row = 0; row < _k_num_states; row++) {
|
for (int row = 0; row < _k_num_states; row++) {
|
||||||
Kfusion[row] = P(row,state_index) / innov_var;
|
Kfusion[row] = P(row, state_index) / innov_var;
|
||||||
}
|
}
|
||||||
|
|
||||||
matrix::SquareMatrix<float, _k_num_states> KHP;
|
matrix::SquareMatrix<float, _k_num_states> KHP;
|
||||||
|
|
||||||
for (unsigned row = 0; row < _k_num_states; row++) {
|
for (unsigned row = 0; row < _k_num_states; row++) {
|
||||||
for (unsigned column = 0; column < _k_num_states; column++) {
|
for (unsigned column = 0; column < _k_num_states; column++) {
|
||||||
KHP(row,column) = Kfusion[row] * P(state_index,column);
|
KHP(row, column) = Kfusion[row] * P(state_index, column);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
@@ -162,7 +160,7 @@ void Ekf::fuseVelPosHeight(const float innov, const float innov_var, const int o
|
|||||||
bool healthy = true;
|
bool healthy = true;
|
||||||
|
|
||||||
for (int i = 0; i < _k_num_states; i++) {
|
for (int i = 0; i < _k_num_states; i++) {
|
||||||
if (P(i,i) < KHP(i,i)) {
|
if (P(i, i) < KHP(i, i)) {
|
||||||
// zero rows and columns
|
// zero rows and columns
|
||||||
P.uncorrelateCovarianceSetVariance<1>(i, 0.0f);
|
P.uncorrelateCovarianceSetVariance<1>(i, 0.0f);
|
||||||
|
|
||||||
@@ -175,13 +173,12 @@ void Ekf::fuseVelPosHeight(const float innov, const float innov_var, const int o
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
|
||||||
// only apply covariance and state corrections if healthy
|
// only apply covariance and state corrections if healthy
|
||||||
if (healthy) {
|
if (healthy) {
|
||||||
// apply the covariance corrections
|
// apply the covariance corrections
|
||||||
for (unsigned row = 0; row < _k_num_states; row++) {
|
for (unsigned row = 0; row < _k_num_states; row++) {
|
||||||
for (unsigned column = 0; column < _k_num_states; column++) {
|
for (unsigned column = 0; column < _k_num_states; column++) {
|
||||||
P(row,column) = P(row,column) - KHP(row,column);
|
P(row, column) = P(row, column) - KHP(row, column);
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
@@ -190,12 +187,10 @@ void Ekf::fuseVelPosHeight(const float innov, const float innov_var, const int o
|
|||||||
|
|
||||||
// apply the state corrections
|
// apply the state corrections
|
||||||
fuse(Kfusion, innov);
|
fuse(Kfusion, innov);
|
||||||
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
void Ekf::setVelPosFaultStatus(const int index, const bool status)
|
void Ekf::setVelPosFaultStatus(const int index, const bool status) {
|
||||||
{
|
|
||||||
if (index == 0) {
|
if (index == 0) {
|
||||||
_fault_status.flags.bad_vel_N = status;
|
_fault_status.flags.bad_vel_N = status;
|
||||||
|
|
||||||
|
|||||||
+7
-23
@@ -42,34 +42,18 @@
|
|||||||
|
|
||||||
#ifdef ECL_STANDALONE
|
#ifdef ECL_STANDALONE
|
||||||
|
|
||||||
namespace math
|
namespace math {
|
||||||
{
|
|
||||||
|
|
||||||
float min(float val1, float val2)
|
float min(float val1, float val2) { return (val1 < val2) ? val1 : val2; }
|
||||||
{
|
|
||||||
return (val1 < val2) ? val1 : val2;
|
|
||||||
}
|
|
||||||
|
|
||||||
float max(float val1, float val2)
|
float max(float val1, float val2) { return (val1 > val2) ? val1 : val2; }
|
||||||
{
|
|
||||||
return (val1 > val2) ? val1 : val2;
|
|
||||||
}
|
|
||||||
|
|
||||||
float constrain(float val, float min, float max)
|
float constrain(float val, float min, float max) { return (val < min) ? min : ((val > max) ? max : val); }
|
||||||
{
|
|
||||||
return (val < min) ? min : ((val > max) ? max : val);
|
|
||||||
}
|
|
||||||
|
|
||||||
float radians(float degrees)
|
float radians(float degrees) { return (degrees / 180.0f) * M_PI_F; }
|
||||||
{
|
|
||||||
return (degrees / 180.0f) * M_PI_F;
|
|
||||||
}
|
|
||||||
|
|
||||||
float degrees(float radians)
|
float degrees(float radians) { return (radians * 180.0f) / M_PI_F; }
|
||||||
{
|
|
||||||
return (radians * 180.0f) / M_PI_F;
|
|
||||||
}
|
|
||||||
|
|
||||||
} // namespace math
|
} // namespace math
|
||||||
|
|
||||||
#endif /* ECL_STANDALONE */
|
#endif /* ECL_STANDALONE */
|
||||||
|
|||||||
+4
-5
@@ -51,8 +51,7 @@
|
|||||||
#define M_PI_2_F (M_PI / 2.0f)
|
#define M_PI_2_F (M_PI / 2.0f)
|
||||||
#endif
|
#endif
|
||||||
|
|
||||||
namespace math
|
namespace math {
|
||||||
{
|
|
||||||
// using namespace Eigen;
|
// using namespace Eigen;
|
||||||
|
|
||||||
float min(float val1, float val2);
|
float min(float val1, float val2);
|
||||||
@@ -61,10 +60,10 @@ float constrain(float val, float min, float max);
|
|||||||
float radians(float degrees);
|
float radians(float degrees);
|
||||||
float degrees(float radians);
|
float degrees(float radians);
|
||||||
|
|
||||||
}
|
} // namespace math
|
||||||
#else
|
#else
|
||||||
|
|
||||||
#include <mathlib/mathlib.h>
|
#include <mathlib/mathlib.h>
|
||||||
|
|
||||||
#endif //ECL_STANDALONE
|
#endif // ECL_STANDALONE
|
||||||
#endif //MATHLIB_H
|
#endif // MATHLIB_H
|
||||||
|
|||||||
@@ -3,6 +3,12 @@ do_format=$1
|
|||||||
files_to_format="""
|
files_to_format="""
|
||||||
EKF/AlphaFilter.hpp
|
EKF/AlphaFilter.hpp
|
||||||
EKF/RingBuffer.h
|
EKF/RingBuffer.h
|
||||||
|
EKF/vel_pos_fusion.cpp
|
||||||
|
EKF/imu_down_sampler.*pp
|
||||||
|
mathlib/*.cpp
|
||||||
|
mathlib/*.h
|
||||||
|
validation/*.cpp
|
||||||
|
validation/*.h
|
||||||
"""
|
"""
|
||||||
RED='\033[0;31m'
|
RED='\033[0;31m'
|
||||||
GREEN='\033[0;32m'
|
GREEN='\033[0;32m'
|
||||||
|
|||||||
@@ -43,17 +43,13 @@
|
|||||||
|
|
||||||
#include <ecl.h>
|
#include <ecl.h>
|
||||||
|
|
||||||
void
|
void DataValidator::put(uint64_t timestamp, float val, uint64_t error_count_in, int priority_in) {
|
||||||
DataValidator::put(uint64_t timestamp, float val, uint64_t error_count_in, int priority_in)
|
float data[dimensions] = {val}; // sets the first value and all others to 0
|
||||||
{
|
|
||||||
float data[dimensions] = { val }; //sets the first value and all others to 0
|
|
||||||
|
|
||||||
put(timestamp, data, error_count_in, priority_in);
|
put(timestamp, data, error_count_in, priority_in);
|
||||||
}
|
}
|
||||||
|
|
||||||
void
|
void DataValidator::put(uint64_t timestamp, const float val[dimensions], uint64_t error_count_in, int priority_in) {
|
||||||
DataValidator::put(uint64_t timestamp, const float val[dimensions], uint64_t error_count_in, int priority_in)
|
|
||||||
{
|
|
||||||
_event_count++;
|
_event_count++;
|
||||||
|
|
||||||
if (error_count_in > _error_count) {
|
if (error_count_in > _error_count) {
|
||||||
@@ -99,9 +95,8 @@ DataValidator::put(uint64_t timestamp, const float val[dimensions], uint64_t err
|
|||||||
_time_last = timestamp;
|
_time_last = timestamp;
|
||||||
}
|
}
|
||||||
|
|
||||||
float
|
float DataValidator::confidence(uint64_t timestamp) {
|
||||||
DataValidator::confidence(uint64_t timestamp)
|
|
||||||
{
|
|
||||||
float ret = 1.0f;
|
float ret = 1.0f;
|
||||||
|
|
||||||
/* check if we have any data */
|
/* check if we have any data */
|
||||||
@@ -128,7 +123,6 @@ DataValidator::confidence(uint64_t timestamp)
|
|||||||
/* cap error density counter at window size */
|
/* cap error density counter at window size */
|
||||||
_error_mask |= ERROR_FLAG_HIGH_ERRDENSITY;
|
_error_mask |= ERROR_FLAG_HIGH_ERRDENSITY;
|
||||||
_error_density = ERROR_DENSITY_WINDOW;
|
_error_density = ERROR_DENSITY_WINDOW;
|
||||||
|
|
||||||
}
|
}
|
||||||
|
|
||||||
/* no critical errors */
|
/* no critical errors */
|
||||||
@@ -144,17 +138,14 @@ DataValidator::confidence(uint64_t timestamp)
|
|||||||
return ret;
|
return ret;
|
||||||
}
|
}
|
||||||
|
|
||||||
void
|
void DataValidator::print() {
|
||||||
DataValidator::print()
|
|
||||||
{
|
|
||||||
if (_time_last == 0) {
|
if (_time_last == 0) {
|
||||||
ECL_INFO("\tno data");
|
ECL_INFO("\tno data");
|
||||||
return;
|
return;
|
||||||
}
|
}
|
||||||
|
|
||||||
for (unsigned i = 0; i < dimensions; i++) {
|
for (unsigned i = 0; i < dimensions; i++) {
|
||||||
ECL_INFO("\tval: %8.4f, lp: %8.4f mean dev: %8.4f RMS: %8.4f conf: %8.4f",
|
ECL_INFO("\tval: %8.4f, lp: %8.4f mean dev: %8.4f RMS: %8.4f conf: %8.4f", (double)_value[i],
|
||||||
(double) _value[i], (double)_lp[i], (double)_mean[i],
|
(double)_lp[i], (double)_mean[i], (double)_rms[i], (double)confidence(ecl_absolute_time()));
|
||||||
(double)_rms[i], (double)confidence(ecl_absolute_time()));
|
|
||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|||||||
+46
-44
@@ -44,8 +44,7 @@
|
|||||||
#include <math.h>
|
#include <math.h>
|
||||||
#include <stdint.h>
|
#include <stdint.h>
|
||||||
|
|
||||||
class DataValidator
|
class DataValidator {
|
||||||
{
|
|
||||||
public:
|
public:
|
||||||
static const unsigned dimensions = 3;
|
static const unsigned dimensions = 3;
|
||||||
|
|
||||||
@@ -57,146 +56,149 @@ public:
|
|||||||
*
|
*
|
||||||
* @param val Item to put
|
* @param val Item to put
|
||||||
*/
|
*/
|
||||||
void put(uint64_t timestamp, float val, uint64_t error_count, int priority);
|
void put(uint64_t timestamp, float val, uint64_t error_count, int priority);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Put a 3D item into the validator.
|
* Put a 3D item into the validator.
|
||||||
*
|
*
|
||||||
* @param val Item to put
|
* @param val Item to put
|
||||||
*/
|
*/
|
||||||
void put(uint64_t timestamp, const float val[dimensions], uint64_t error_count, int priority);
|
void put(uint64_t timestamp, const float val[dimensions], uint64_t error_count, int priority);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the next sibling in the group
|
* Get the next sibling in the group
|
||||||
*
|
*
|
||||||
* @return the next sibling
|
* @return the next sibling
|
||||||
*/
|
*/
|
||||||
DataValidator *sibling() { return _sibling; }
|
DataValidator *sibling() { return _sibling; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Set the sibling to the next node in the group
|
* Set the sibling to the next node in the group
|
||||||
*
|
*
|
||||||
*/
|
*/
|
||||||
void setSibling(DataValidator *new_sibling) { _sibling = new_sibling; }
|
void setSibling(DataValidator *new_sibling) { _sibling = new_sibling; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the confidence of this validator
|
* Get the confidence of this validator
|
||||||
* @return the confidence between 0 and 1
|
* @return the confidence between 0 and 1
|
||||||
*/
|
*/
|
||||||
float confidence(uint64_t timestamp);
|
float confidence(uint64_t timestamp);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the error count of this validator
|
* Get the error count of this validator
|
||||||
* @return the error count
|
* @return the error count
|
||||||
*/
|
*/
|
||||||
uint64_t error_count() const { return _error_count; }
|
uint64_t error_count() const { return _error_count; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the values of this validator
|
* Get the values of this validator
|
||||||
* @return the stored value
|
* @return the stored value
|
||||||
*/
|
*/
|
||||||
float *value() { return _value; }
|
float *value() { return _value; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the used status of this validator
|
* Get the used status of this validator
|
||||||
* @return true if this validator ever saw data
|
* @return true if this validator ever saw data
|
||||||
*/
|
*/
|
||||||
bool used() const { return (_time_last > 0); }
|
bool used() const { return (_time_last > 0); }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the priority of this validator
|
* Get the priority of this validator
|
||||||
* @return the stored priority
|
* @return the stored priority
|
||||||
*/
|
*/
|
||||||
int priority() const { return _priority; }
|
int priority() const { return _priority; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the error state of this validator
|
* Get the error state of this validator
|
||||||
* @return the bitmask with the error status
|
* @return the bitmask with the error status
|
||||||
*/
|
*/
|
||||||
uint32_t state() const { return _error_mask; }
|
uint32_t state() const { return _error_mask; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Reset the error state of this validator
|
* Reset the error state of this validator
|
||||||
*/
|
*/
|
||||||
void reset_state() { _error_mask = ERROR_FLAG_NO_ERROR; }
|
void reset_state() { _error_mask = ERROR_FLAG_NO_ERROR; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the RMS values of this validator
|
* Get the RMS values of this validator
|
||||||
* @return the stored RMS
|
* @return the stored RMS
|
||||||
*/
|
*/
|
||||||
float *rms() { return _rms; }
|
float *rms() { return _rms; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the vibration offset
|
* Get the vibration offset
|
||||||
* @return the stored vibration offset
|
* @return the stored vibration offset
|
||||||
*/
|
*/
|
||||||
float *vibration_offset() { return _vibe; }
|
float *vibration_offset() { return _vibe; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Print the validator value
|
* Print the validator value
|
||||||
*
|
*
|
||||||
*/
|
*/
|
||||||
void print();
|
void print();
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Set the timeout value
|
* Set the timeout value
|
||||||
*
|
*
|
||||||
* @param timeout_interval_us The timeout interval in microseconds
|
* @param timeout_interval_us The timeout interval in microseconds
|
||||||
*/
|
*/
|
||||||
void set_timeout(uint32_t timeout_interval_us) { _timeout_interval = timeout_interval_us; }
|
void set_timeout(uint32_t timeout_interval_us) { _timeout_interval = timeout_interval_us; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Set the equal count threshold
|
* Set the equal count threshold
|
||||||
*
|
*
|
||||||
* @param threshold The number of equal values before considering the sensor stale
|
* @param threshold The number of equal values before considering the sensor stale
|
||||||
*/
|
*/
|
||||||
void set_equal_value_threshold(uint32_t threshold) { _value_equal_count_threshold = threshold; }
|
void set_equal_value_threshold(uint32_t threshold) { _value_equal_count_threshold = threshold; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the timeout value
|
* Get the timeout value
|
||||||
*
|
*
|
||||||
* @return The timeout interval in microseconds
|
* @return The timeout interval in microseconds
|
||||||
*/
|
*/
|
||||||
uint32_t get_timeout() const { return _timeout_interval; }
|
uint32_t get_timeout() const { return _timeout_interval; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Data validator error states
|
* Data validator error states
|
||||||
*/
|
*/
|
||||||
static constexpr uint32_t ERROR_FLAG_NO_ERROR = (0x00000000U);
|
static constexpr uint32_t ERROR_FLAG_NO_ERROR = (0x00000000U);
|
||||||
static constexpr uint32_t ERROR_FLAG_NO_DATA = (0x00000001U);
|
static constexpr uint32_t ERROR_FLAG_NO_DATA = (0x00000001U);
|
||||||
static constexpr uint32_t ERROR_FLAG_STALE_DATA = (0x00000001U << 1);
|
static constexpr uint32_t ERROR_FLAG_STALE_DATA = (0x00000001U << 1);
|
||||||
static constexpr uint32_t ERROR_FLAG_TIMEOUT = (0x00000001U << 2);
|
static constexpr uint32_t ERROR_FLAG_TIMEOUT = (0x00000001U << 2);
|
||||||
static constexpr uint32_t ERROR_FLAG_HIGH_ERRCOUNT = (0x00000001U << 3);
|
static constexpr uint32_t ERROR_FLAG_HIGH_ERRCOUNT = (0x00000001U << 3);
|
||||||
static constexpr uint32_t ERROR_FLAG_HIGH_ERRDENSITY = (0x00000001U << 4);
|
static constexpr uint32_t ERROR_FLAG_HIGH_ERRDENSITY = (0x00000001U << 4);
|
||||||
|
|
||||||
private:
|
private:
|
||||||
uint32_t _error_mask{ERROR_FLAG_NO_ERROR}; /**< sensor error state */
|
uint32_t _error_mask{ERROR_FLAG_NO_ERROR}; /**< sensor error state */
|
||||||
|
|
||||||
uint32_t _timeout_interval{20000}; /**< interval in which the datastream times out in us */
|
uint32_t _timeout_interval{20000}; /**< interval in which the datastream times out in us */
|
||||||
|
|
||||||
uint64_t _time_last{0}; /**< last timestamp */
|
uint64_t _time_last{0}; /**< last timestamp */
|
||||||
uint64_t _event_count{0}; /**< total data counter */
|
uint64_t _event_count{0}; /**< total data counter */
|
||||||
uint64_t _error_count{0}; /**< error count */
|
uint64_t _error_count{0}; /**< error count */
|
||||||
|
|
||||||
int _error_density{0}; /**< ratio between successful reads and errors */
|
int _error_density{0}; /**< ratio between successful reads and errors */
|
||||||
|
|
||||||
int _priority{0}; /**< sensor nominal priority */
|
int _priority{0}; /**< sensor nominal priority */
|
||||||
|
|
||||||
float _mean[dimensions] {}; /**< mean of value */
|
float _mean[dimensions]{}; /**< mean of value */
|
||||||
float _lp[dimensions] {}; /**< low pass value */
|
float _lp[dimensions]{}; /**< low pass value */
|
||||||
float _M2[dimensions] {}; /**< RMS component value */
|
float _M2[dimensions]{}; /**< RMS component value */
|
||||||
float _rms[dimensions] {}; /**< root mean square error */
|
float _rms[dimensions]{}; /**< root mean square error */
|
||||||
float _value[dimensions] {}; /**< last value */
|
float _value[dimensions]{}; /**< last value */
|
||||||
float _vibe[dimensions] {}; /**< vibration level, in sensor unit */
|
float _vibe[dimensions]{}; /**< vibration level, in sensor unit */
|
||||||
|
|
||||||
unsigned _value_equal_count{0}; /**< equal values in a row */
|
unsigned _value_equal_count{0}; /**< equal values in a row */
|
||||||
unsigned _value_equal_count_threshold{VALUE_EQUAL_COUNT_DEFAULT}; /**< when to consider an equal count as a problem */
|
unsigned _value_equal_count_threshold{
|
||||||
|
VALUE_EQUAL_COUNT_DEFAULT}; /**< when to consider an equal count as a problem */
|
||||||
|
|
||||||
DataValidator *_sibling{nullptr}; /**< sibling in the group */
|
DataValidator *_sibling{nullptr}; /**< sibling in the group */
|
||||||
|
|
||||||
static const constexpr unsigned NORETURN_ERRCOUNT = 10000; /**< if the error count reaches this value, return sensor as invalid */
|
static const constexpr unsigned NORETURN_ERRCOUNT =
|
||||||
static const constexpr float ERROR_DENSITY_WINDOW = 100.0f; /**< window in measurement counts for errors */
|
10000; /**< if the error count reaches this value, return sensor as invalid */
|
||||||
static const constexpr unsigned VALUE_EQUAL_COUNT_DEFAULT = 100; /**< if the sensor value is the same (accumulated also between axes) this many times, flag it */
|
static const constexpr float ERROR_DENSITY_WINDOW = 100.0f; /**< window in measurement counts for errors */
|
||||||
|
static const constexpr unsigned VALUE_EQUAL_COUNT_DEFAULT =
|
||||||
|
100; /**< if the sensor value is the same (accumulated also between axes) this many times, flag it */
|
||||||
|
|
||||||
/* we don't want this class to be copied */
|
/* we don't want this class to be copied */
|
||||||
DataValidator(const DataValidator &) = delete;
|
DataValidator(const DataValidator &) = delete;
|
||||||
|
|||||||
@@ -43,8 +43,8 @@
|
|||||||
#include <ecl.h>
|
#include <ecl.h>
|
||||||
#include <float.h>
|
#include <float.h>
|
||||||
|
|
||||||
DataValidatorGroup::DataValidatorGroup(unsigned siblings)
|
DataValidatorGroup::DataValidatorGroup(unsigned siblings) {
|
||||||
{
|
|
||||||
DataValidator *next = nullptr;
|
DataValidator *next = nullptr;
|
||||||
DataValidator *prev = nullptr;
|
DataValidator *prev = nullptr;
|
||||||
|
|
||||||
@@ -68,8 +68,7 @@ DataValidatorGroup::DataValidatorGroup(unsigned siblings)
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
DataValidatorGroup::~DataValidatorGroup()
|
DataValidatorGroup::~DataValidatorGroup() {
|
||||||
{
|
|
||||||
while (_first) {
|
while (_first) {
|
||||||
DataValidator *next = _first->sibling();
|
DataValidator *next = _first->sibling();
|
||||||
delete (_first);
|
delete (_first);
|
||||||
@@ -77,8 +76,8 @@ DataValidatorGroup::~DataValidatorGroup()
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
DataValidator *DataValidatorGroup::add_new_validator()
|
DataValidator *DataValidatorGroup::add_new_validator() {
|
||||||
{
|
|
||||||
DataValidator *validator = new DataValidator();
|
DataValidator *validator = new DataValidator();
|
||||||
|
|
||||||
if (!validator) {
|
if (!validator) {
|
||||||
@@ -91,9 +90,8 @@ DataValidator *DataValidatorGroup::add_new_validator()
|
|||||||
return _last;
|
return _last;
|
||||||
}
|
}
|
||||||
|
|
||||||
void
|
void DataValidatorGroup::set_timeout(uint32_t timeout_interval_us) {
|
||||||
DataValidatorGroup::set_timeout(uint32_t timeout_interval_us)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
|
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
@@ -104,9 +102,8 @@ DataValidatorGroup::set_timeout(uint32_t timeout_interval_us)
|
|||||||
_timeout_interval_us = timeout_interval_us;
|
_timeout_interval_us = timeout_interval_us;
|
||||||
}
|
}
|
||||||
|
|
||||||
void
|
void DataValidatorGroup::set_equal_value_threshold(uint32_t threshold) {
|
||||||
DataValidatorGroup::set_equal_value_threshold(uint32_t threshold)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
|
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
@@ -115,10 +112,9 @@ DataValidatorGroup::set_equal_value_threshold(uint32_t threshold)
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
|
void DataValidatorGroup::put(unsigned index, uint64_t timestamp, const float val[3], uint64_t error_count,
|
||||||
|
int priority) {
|
||||||
|
|
||||||
void
|
|
||||||
DataValidatorGroup::put(unsigned index, uint64_t timestamp, const float val[3], uint64_t error_count, int priority)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
unsigned i = 0;
|
unsigned i = 0;
|
||||||
|
|
||||||
@@ -133,9 +129,8 @@ DataValidatorGroup::put(unsigned index, uint64_t timestamp, const float val[3],
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
float *
|
float *DataValidatorGroup::get_best(uint64_t timestamp, int *index) {
|
||||||
DataValidatorGroup::get_best(uint64_t timestamp, int *index)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
|
|
||||||
// XXX This should eventually also include voting
|
// XXX This should eventually also include voting
|
||||||
@@ -164,9 +159,8 @@ DataValidatorGroup::get_best(uint64_t timestamp, int *index)
|
|||||||
*/
|
*/
|
||||||
if ((((max_confidence < MIN_REGULAR_CONFIDENCE) && (confidence >= MIN_REGULAR_CONFIDENCE)) ||
|
if ((((max_confidence < MIN_REGULAR_CONFIDENCE) && (confidence >= MIN_REGULAR_CONFIDENCE)) ||
|
||||||
(confidence > max_confidence && (next->priority() >= max_priority)) ||
|
(confidence > max_confidence && (next->priority() >= max_priority)) ||
|
||||||
(fabsf(confidence - max_confidence) < 0.01f && (next->priority() > max_priority))
|
(fabsf(confidence - max_confidence) < 0.01f && (next->priority() > max_priority))) &&
|
||||||
) && (confidence > 0.0f)) {
|
(confidence > 0.0f)) {
|
||||||
|
|
||||||
max_index = i;
|
max_index = i;
|
||||||
max_confidence = confidence;
|
max_confidence = confidence;
|
||||||
max_priority = next->priority();
|
max_priority = next->priority();
|
||||||
@@ -180,13 +174,11 @@ DataValidatorGroup::get_best(uint64_t timestamp, int *index)
|
|||||||
/* the current best sensor is not matching the previous best sensor,
|
/* the current best sensor is not matching the previous best sensor,
|
||||||
* or the only sensor went bad */
|
* or the only sensor went bad */
|
||||||
if (max_index != _curr_best || ((max_confidence < FLT_EPSILON) && (_curr_best >= 0))) {
|
if (max_index != _curr_best || ((max_confidence < FLT_EPSILON) && (_curr_best >= 0))) {
|
||||||
|
|
||||||
bool true_failsafe = true;
|
bool true_failsafe = true;
|
||||||
|
|
||||||
/* check whether the switch was a failsafe or preferring a higher priority sensor */
|
/* check whether the switch was a failsafe or preferring a higher priority sensor */
|
||||||
if (pre_check_prio != -1 && pre_check_prio < max_priority &&
|
if (pre_check_prio != -1 && pre_check_prio < max_priority &&
|
||||||
fabsf(pre_check_confidence - max_confidence) < 0.1f) {
|
fabsf(pre_check_confidence - max_confidence) < 0.1f) {
|
||||||
|
|
||||||
/* this is not a failover */
|
/* this is not a failover */
|
||||||
true_failsafe = false;
|
true_failsafe = false;
|
||||||
|
|
||||||
@@ -222,16 +214,14 @@ DataValidatorGroup::get_best(uint64_t timestamp, int *index)
|
|||||||
return (best) ? best->value() : nullptr;
|
return (best) ? best->value() : nullptr;
|
||||||
}
|
}
|
||||||
|
|
||||||
float
|
float DataValidatorGroup::get_vibration_factor(uint64_t timestamp) {
|
||||||
DataValidatorGroup::get_vibration_factor(uint64_t timestamp)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
|
|
||||||
float vibe = 0.0f;
|
float vibe = 0.0f;
|
||||||
|
|
||||||
/* find the best RMS value of a non-timed out sensor */
|
/* find the best RMS value of a non-timed out sensor */
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
|
|
||||||
if (next->confidence(timestamp) > 0.5f) {
|
if (next->confidence(timestamp) > 0.5f) {
|
||||||
float *rms = next->rms();
|
float *rms = next->rms();
|
||||||
|
|
||||||
@@ -248,16 +238,14 @@ DataValidatorGroup::get_vibration_factor(uint64_t timestamp)
|
|||||||
return vibe;
|
return vibe;
|
||||||
}
|
}
|
||||||
|
|
||||||
float
|
float DataValidatorGroup::get_vibration_offset(uint64_t timestamp, int axis) {
|
||||||
DataValidatorGroup::get_vibration_offset(uint64_t timestamp, int axis)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
|
|
||||||
float vibe = -1.0f;
|
float vibe = -1.0f;
|
||||||
|
|
||||||
/* find the best vibration value of a non-timed out sensor */
|
/* find the best vibration value of a non-timed out sensor */
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
|
|
||||||
if (next->confidence(timestamp) > 0.5f) {
|
if (next->confidence(timestamp) > 0.5f) {
|
||||||
float *vibration_offset = next->vibration_offset();
|
float *vibration_offset = next->vibration_offset();
|
||||||
|
|
||||||
@@ -272,13 +260,10 @@ DataValidatorGroup::get_vibration_offset(uint64_t timestamp, int axis)
|
|||||||
return vibe;
|
return vibe;
|
||||||
}
|
}
|
||||||
|
|
||||||
void
|
void DataValidatorGroup::print() {
|
||||||
DataValidatorGroup::print()
|
|
||||||
{
|
ECL_INFO("validator: best: %d, prev best: %d, failsafe: %s (%u events)", _curr_best, _prev_best,
|
||||||
/* print the group's state */
|
(_toggle_count > 0) ? "YES" : "NO", _toggle_count);
|
||||||
ECL_INFO("validator: best: %d, prev best: %d, failsafe: %s (%u events)",
|
|
||||||
_curr_best, _prev_best, (_toggle_count > 0) ? "YES" : "NO",
|
|
||||||
_toggle_count);
|
|
||||||
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
unsigned i = 0;
|
unsigned i = 0;
|
||||||
@@ -303,14 +288,14 @@ DataValidatorGroup::print()
|
|||||||
}
|
}
|
||||||
}
|
}
|
||||||
|
|
||||||
int
|
int DataValidatorGroup::failover_index() {
|
||||||
DataValidatorGroup::failover_index()
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
unsigned i = 0;
|
unsigned i = 0;
|
||||||
|
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
if (next->used() && (next->state() != DataValidator::ERROR_FLAG_NO_ERROR) && (i == (unsigned)_prev_best)) {
|
if (next->used() && (next->state() != DataValidator::ERROR_FLAG_NO_ERROR) &&
|
||||||
|
(i == (unsigned)_prev_best)) {
|
||||||
return i;
|
return i;
|
||||||
}
|
}
|
||||||
|
|
||||||
@@ -321,14 +306,14 @@ DataValidatorGroup::failover_index()
|
|||||||
return -1;
|
return -1;
|
||||||
}
|
}
|
||||||
|
|
||||||
uint32_t
|
uint32_t DataValidatorGroup::failover_state() {
|
||||||
DataValidatorGroup::failover_state()
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
unsigned i = 0;
|
unsigned i = 0;
|
||||||
|
|
||||||
while (next != nullptr) {
|
while (next != nullptr) {
|
||||||
if (next->used() && (next->state() != DataValidator::ERROR_FLAG_NO_ERROR) && (i == (unsigned)_prev_best)) {
|
if (next->used() && (next->state() != DataValidator::ERROR_FLAG_NO_ERROR) &&
|
||||||
|
(i == (unsigned)_prev_best)) {
|
||||||
return next->state();
|
return next->state();
|
||||||
}
|
}
|
||||||
|
|
||||||
@@ -339,9 +324,8 @@ DataValidatorGroup::failover_state()
|
|||||||
return DataValidator::ERROR_FLAG_NO_ERROR;
|
return DataValidator::ERROR_FLAG_NO_ERROR;
|
||||||
}
|
}
|
||||||
|
|
||||||
uint32_t
|
uint32_t DataValidatorGroup::get_sensor_state(unsigned index) {
|
||||||
DataValidatorGroup::get_sensor_state(unsigned index)
|
|
||||||
{
|
|
||||||
DataValidator *next = _first;
|
DataValidator *next = _first;
|
||||||
unsigned i = 0;
|
unsigned i = 0;
|
||||||
|
|
||||||
|
|||||||
@@ -43,8 +43,7 @@
|
|||||||
|
|
||||||
#include "data_validator.h"
|
#include "data_validator.h"
|
||||||
|
|
||||||
class DataValidatorGroup
|
class DataValidatorGroup {
|
||||||
{
|
|
||||||
public:
|
public:
|
||||||
/**
|
/**
|
||||||
* @param siblings initial number of DataValidator's. Must be > 0.
|
* @param siblings initial number of DataValidator's. Must be > 0.
|
||||||
@@ -67,90 +66,89 @@ public:
|
|||||||
* @param error_count The current error count of the sensor
|
* @param error_count The current error count of the sensor
|
||||||
* @param priority The priority of the sensor
|
* @param priority The priority of the sensor
|
||||||
*/
|
*/
|
||||||
void put(unsigned index, uint64_t timestamp, const float val[3], uint64_t error_count, int priority);
|
void put(unsigned index, uint64_t timestamp, const float val[3], uint64_t error_count, int priority);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the best data triplet of the group
|
* Get the best data triplet of the group
|
||||||
*
|
*
|
||||||
* @return pointer to the array of best values
|
* @return pointer to the array of best values
|
||||||
*/
|
*/
|
||||||
float *get_best(uint64_t timestamp, int *index);
|
float *get_best(uint64_t timestamp, int *index);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the RMS / vibration factor
|
* Get the RMS / vibration factor
|
||||||
*
|
*
|
||||||
* @return float value representing the RMS, which a valid indicator for vibration
|
* @return float value representing the RMS, which a valid indicator for vibration
|
||||||
*/
|
*/
|
||||||
float get_vibration_factor(uint64_t timestamp);
|
float get_vibration_factor(uint64_t timestamp);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the vibration offset in the sensor unit
|
* Get the vibration offset in the sensor unit
|
||||||
*
|
*
|
||||||
* @return float value representing the vibration offset
|
* @return float value representing the vibration offset
|
||||||
*/
|
*/
|
||||||
float get_vibration_offset(uint64_t timestamp, int axis);
|
float get_vibration_offset(uint64_t timestamp, int axis);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the number of failover events
|
* Get the number of failover events
|
||||||
*
|
*
|
||||||
* @return the number of failovers
|
* @return the number of failovers
|
||||||
*/
|
*/
|
||||||
unsigned failover_count() const { return _toggle_count; }
|
unsigned failover_count() const { return _toggle_count; }
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the index of the failed sensor in the group
|
* Get the index of the failed sensor in the group
|
||||||
*
|
*
|
||||||
* @return index of the failed sensor
|
* @return index of the failed sensor
|
||||||
*/
|
*/
|
||||||
int failover_index();
|
int failover_index();
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the error state of the failed sensor in the group
|
* Get the error state of the failed sensor in the group
|
||||||
*
|
*
|
||||||
* @return bitmask with erro states of the failed sensor
|
* @return bitmask with erro states of the failed sensor
|
||||||
*/
|
*/
|
||||||
uint32_t failover_state();
|
uint32_t failover_state();
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Get the error state of the sensor with the specified index
|
* Get the error state of the sensor with the specified index
|
||||||
*
|
*
|
||||||
* @return bitmask with error states of the sensor
|
* @return bitmask with error states of the sensor
|
||||||
*/
|
*/
|
||||||
uint32_t get_sensor_state(unsigned index);
|
uint32_t get_sensor_state(unsigned index);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Print the validator value
|
* Print the validator value
|
||||||
*
|
*
|
||||||
*/
|
*/
|
||||||
void print();
|
void print();
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Set the timeout value for the whole group
|
* Set the timeout value for the whole group
|
||||||
*
|
*
|
||||||
* @param timeout_interval_us The timeout interval in microseconds
|
* @param timeout_interval_us The timeout interval in microseconds
|
||||||
*/
|
*/
|
||||||
void set_timeout(uint32_t timeout_interval_us);
|
void set_timeout(uint32_t timeout_interval_us);
|
||||||
|
|
||||||
/**
|
/**
|
||||||
* Set the equal count threshold for the whole group
|
* Set the equal count threshold for the whole group
|
||||||
*
|
*
|
||||||
* @param threshold The number of equal values before considering the sensor stale
|
* @param threshold The number of equal values before considering the sensor stale
|
||||||
*/
|
*/
|
||||||
void set_equal_value_threshold(uint32_t threshold);
|
void set_equal_value_threshold(uint32_t threshold);
|
||||||
|
|
||||||
|
|
||||||
private:
|
private:
|
||||||
DataValidator *_first{nullptr}; /**< first node in the group */
|
DataValidator *_first{nullptr}; /**< first node in the group */
|
||||||
DataValidator *_last{nullptr}; /**< last node in the group */
|
DataValidator *_last{nullptr}; /**< last node in the group */
|
||||||
|
|
||||||
uint32_t _timeout_interval_us{0}; /**< currently set timeout */
|
uint32_t _timeout_interval_us{0}; /**< currently set timeout */
|
||||||
|
|
||||||
int _curr_best{-1}; /**< currently best index */
|
int _curr_best{-1}; /**< currently best index */
|
||||||
int _prev_best{-1}; /**< the previous best index */
|
int _prev_best{-1}; /**< the previous best index */
|
||||||
|
|
||||||
uint64_t _first_failover_time{0}; /**< timestamp where the first failover occured or zero if none occured */
|
uint64_t _first_failover_time{0}; /**< timestamp where the first failover occured or zero if none occured */
|
||||||
|
|
||||||
unsigned _toggle_count{0}; /**< number of back and forth switches between two sensors */
|
unsigned _toggle_count{0}; /**< number of back and forth switches between two sensors */
|
||||||
|
|
||||||
static constexpr float MIN_REGULAR_CONFIDENCE = 0.9f;
|
static constexpr float MIN_REGULAR_CONFIDENCE = 0.9f;
|
||||||
|
|
||||||
|
|||||||
Reference in New Issue
Block a user