Mirrors_Framework/PX4-Autopilot
1
0
Fork 0
mirror of https://github.com/PX4/PX4-Autopilot.git synced 2026-06-01 02:55:07 +08:00
Code Issues Actions 19 Packages Projects Releases Wiki Activity

testing float/double with Type template

Browse Source
This commit is contained in:
Bart Slinger
2018-09-14 13:15:52 +02:00
committed by Beat Küng
parent 707e288019
commit 30d5a79432
1 changed files with 24 additions and 48 deletions
Download Patch File Download Diff File Expand all files Collapse all files
test/least_squares.cpp
+24 -48
View File
@@ -4,7 +4,7 @@
using namespace matrix;
int test_4x3(void);
int test_4x4(void);
template<typename Type> int test_4x4(void);
int test_4x4_type_double(void);
int test_div_zero(void);
@@ -12,14 +12,15 @@ int main()
{
int ret;
ret = test_4x4();
ret = test_4x4<float>();
if (ret != 0) return ret;
ret = test_4x4<double>();
if (ret != 0) return ret;
ret = test_4x3();
if (ret != 0) return ret;
ret = test_4x4_type_double();
ret = test_div_zero();
if (ret != 0) return ret;
@@ -35,7 +36,7 @@ int test_4x3() {
};
Matrix<float, 4, 3> A(data);
float b_data[4] = {2.0, 3.0, 4.0, 5.0};
float b_data[4] = {2.0f, 3.0f, 4.0f, 5.0f};
Vector<float, 4> b(b_data);
float x_check_data[3] = {-0.69168233f,
@@ -51,54 +52,29 @@ int test_4x3() {
return 0;
}
template<typename Type>
int test_4x4() {
// Start with an (m x n) A matrix
float data[16] = { 20.f, -10.f, -13.f, 21.f,
17.f, 16.f, -18.f, -14.f,
0.7f, -0.8f, 0.9f, -0.5f,
-1.f, -1.1f, -1.2f, -1.3f
};
Matrix<float, 4, 4> A(data);
const Type data[16] = { 20.f, -10.f, -13.f, 21.f,
17.f, 16.f, -18.f, -14.f,
0.7f, -0.8f, 0.9f, -0.5f,
-1.f, -1.1f, -1.2f, -1.3f
};
Matrix<Type, 4, 4> A(data);
float b_data[4] = {2.0, 3.0, 4.0, 5.0};
Vector<float, 4> b(b_data);
Type b_data[4] = {2.0f, 3.0f, 4.0f, 5.0f};
Vector<Type, 4> b(b_data);
float x_check_data[4] = { 0.97893433f,
-2.80798701f,
-0.03175765f,
-2.19387649f
};
Vector<float, 4> x_check(x_check_data);
Type x_check_data[4] = { 0.97893433f,
-2.80798701f,
-0.03175765f,
-2.19387649f
};
Vector<Type, 4> x_check(x_check_data);
LeastSquaresSolver<float, 4, 4> qrd = LeastSquaresSolver<float, 4, 4>(A);
LeastSquaresSolver<Type, 4, 4> qrd = LeastSquaresSolver<Type, 4, 4>(A);
Vector<float, 4> x = qrd.solve(b);
TEST(isEqual(x, x_check));
return 0;
}
int test_4x4_type_double() {
// Start with an (m x n) A matrix
double data[16] = { 20., -10., -13., 21.,
17., 16., -18., -14.,
0.7, -0.8, 0.9, -0.5,
-1., -1.1, -1.2, -1.3
};
Matrix<double, 4, 4> A(data);
double b_data[4] = {2.0, 3.0, 4.0, 5.0};
Vector<double, 4> b(b_data);
double x_check_data[4] = { 0.97893433,
-2.80798701,
-0.03175765,
-2.19387649
};
Vector<double, 4> x_check(x_check_data);
LeastSquaresSolver<double, 4, 4> qrd = LeastSquaresSolver<double, 4, 4>(A);
Vector<double, 4> x = qrd.solve(b);
Vector<Type, 4> x = qrd.solve(b);
TEST(isEqual(x, x_check));
return 0;
}
@@ -107,7 +83,7 @@ int test_div_zero() {
float data[4] = {0.0f, 0.0f, 0.0f, 0.0f};
Matrix<float, 2, 2> A(data);
float b_data[2] = {1.0, 1.0};
float b_data[2] = {1.0f, 1.0f};
Vector<float, 2> b(b_data);
// Implement such that x returns zeros if it reaches div by zero