Mercurial > hg > constant-q-cpp

annotate test/TestMathUtilities.cpp @ 196:da283326bcd3 tip master

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Update plugin versions in RDF
author Chris Cannam <cannam@all-day-breakfast.com>
date Fri, 28 Feb 2020 09:43:02 +0000
parents 6b13f9c694a8
children
rev   line source
c@131 1 /* -*- c-basic-offset: 4 indent-tabs-mode: nil -*- vi:set ts=8 sts=4 sw=4: */
c@131 2
c@131 3 #include "dsp/MathUtilities.h"
c@131 4
c@131 5 #include <cmath>
c@131 6
c@131 7 #define BOOST_TEST_DYN_LINK
c@131 8 #define BOOST_TEST_MAIN
c@131 9
c@131 10 #include <boost/test/unit_test.hpp>
c@131 11
c@131 12 BOOST_AUTO_TEST_SUITE(TestMathUtilities)
c@131 13
c@131 14 BOOST_AUTO_TEST_CASE(round)
c@131 15 {
c@131 16 BOOST_CHECK_EQUAL(MathUtilities::round(0.5), 1.0);
c@131 17 BOOST_CHECK_EQUAL(MathUtilities::round(0.49), 0.0);
c@131 18 BOOST_CHECK_EQUAL(MathUtilities::round(0.99), 1.0);
c@131 19 BOOST_CHECK_EQUAL(MathUtilities::round(0.01), 0.0);
c@131 20 BOOST_CHECK_EQUAL(MathUtilities::round(0.0), 0.0);
c@131 21 BOOST_CHECK_EQUAL(MathUtilities::round(100.0), 100.0);
c@131 22 BOOST_CHECK_EQUAL(MathUtilities::round(-0.2), 0.0);
c@131 23 BOOST_CHECK_EQUAL(MathUtilities::round(-0.5), -1.0);
c@131 24 BOOST_CHECK_EQUAL(MathUtilities::round(-0.99), -1.0);
c@131 25 BOOST_CHECK_EQUAL(MathUtilities::round(-1.0), -1.0);
c@131 26 BOOST_CHECK_EQUAL(MathUtilities::round(-1.1), -1.0);
c@131 27 BOOST_CHECK_EQUAL(MathUtilities::round(-1.5), -2.0);
c@131 28 }
c@131 29
c@131 30 BOOST_AUTO_TEST_CASE(mean)
c@131 31 {
c@131 32 BOOST_CHECK_EQUAL(MathUtilities::mean(0, 0), 0);
c@131 33 double d0[] = { 0, 4, 3, -1 };
c@131 34 BOOST_CHECK_EQUAL(MathUtilities::mean(d0, 4), 1.5);
c@131 35 double d1[] = { -2.6 };
c@131 36 BOOST_CHECK_EQUAL(MathUtilities::mean(d1, 1), -2.6);
c@131 37 std::vector<double> v;
c@131 38 v.push_back(0);
c@131 39 v.push_back(4);
c@131 40 v.push_back(3);
c@131 41 v.push_back(-1);
c@131 42 BOOST_CHECK_EQUAL(MathUtilities::mean(v, 0, 4), 1.5);
c@131 43 BOOST_CHECK_EQUAL(MathUtilities::mean(v, 1, 2), 3.5);
c@131 44 BOOST_CHECK_EQUAL(MathUtilities::mean(v, 3, 1), -1);
c@131 45 BOOST_CHECK_EQUAL(MathUtilities::mean(v, 3, 0), 0);
c@131 46 }
c@131 47
c@131 48 BOOST_AUTO_TEST_CASE(sum)
c@131 49 {
c@131 50 BOOST_CHECK_EQUAL(MathUtilities::sum(0, 0), 0);
c@131 51 double d0[] = { 0, 4, 3, -1 };
c@131 52 BOOST_CHECK_EQUAL(MathUtilities::sum(d0, 4), 6);
c@131 53 double d1[] = { -2.6 };
c@131 54 BOOST_CHECK_EQUAL(MathUtilities::sum(d1, 1), -2.6);
c@131 55 }
c@131 56
c@131 57 BOOST_AUTO_TEST_CASE(median)
c@131 58 {
c@131 59 BOOST_CHECK_EQUAL(MathUtilities::median(0, 0), 0);
c@131 60 double d0[] = { 0, 4, 3, -1 };
c@131 61 BOOST_CHECK_EQUAL(MathUtilities::median(d0, 4), 1.5);
c@131 62 double d1[] = { 0, 4, 3, -1, -1 };
c@131 63 BOOST_CHECK_EQUAL(MathUtilities::median(d1, 5), 0);
c@131 64 double d2[] = { 1.0, -2.0 };
c@131 65 BOOST_CHECK_EQUAL(MathUtilities::median(d2, 2), -0.5);
c@131 66 double d3[] = { -2.6 };
c@131 67 BOOST_CHECK_EQUAL(MathUtilities::median(d3, 1), -2.6);
c@131 68 }
c@131 69
c@131 70 BOOST_AUTO_TEST_CASE(princarg)
c@131 71 {
c@131 72 BOOST_CHECK_EQUAL(MathUtilities::princarg(M_PI), M_PI);
c@131 73 BOOST_CHECK_EQUAL(MathUtilities::princarg(-M_PI), M_PI);
c@131 74 BOOST_CHECK_EQUAL(MathUtilities::princarg(2 * M_PI), 0.0);
c@131 75 BOOST_CHECK_EQUAL(MathUtilities::princarg(5 * M_PI), M_PI);
c@131 76 BOOST_CHECK_EQUAL(MathUtilities::princarg(1.0), 1.0);
c@131 77 BOOST_CHECK_EQUAL(MathUtilities::princarg(-1.0), -1.0);
c@131 78 BOOST_CHECK_EQUAL(MathUtilities::princarg(-10.0), -10.0 + 4 * M_PI);
c@131 79 }
c@131 80
c@131 81 BOOST_AUTO_TEST_CASE(isPowerOfTwo)
c@131 82 {
c@131 83 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(0), false);
c@131 84 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(1), true);
c@131 85 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(-2), false);
c@131 86 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(2), true);
c@131 87 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(3), false);
c@131 88 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(12), false);
c@131 89 BOOST_CHECK_EQUAL(MathUtilities::isPowerOfTwo(16), true);
c@131 90 }
c@131 91
c@131 92 BOOST_AUTO_TEST_CASE(nextPowerOfTwo)
c@131 93 {
c@131 94 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(0), 1);
c@131 95 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(1), 1);
c@131 96 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(-2), 1);
c@131 97 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(2), 2);
c@131 98 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(3), 4);
c@131 99 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(12), 16);
c@131 100 BOOST_CHECK_EQUAL(MathUtilities::nextPowerOfTwo(16), 16);
c@131 101 }
c@131 102
c@131 103 BOOST_AUTO_TEST_CASE(previousPowerOfTwo)
c@131 104 {
c@131 105 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(0), 1);
c@131 106 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(1), 1);
c@131 107 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(-2), 1);
c@131 108 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(2), 2);
c@131 109 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(3), 2);
c@131 110 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(12), 8);
c@131 111 BOOST_CHECK_EQUAL(MathUtilities::previousPowerOfTwo(16), 16);
c@131 112 }
c@131 113
c@131 114 BOOST_AUTO_TEST_CASE(nearestPowerOfTwo)
c@131 115 {
c@131 116 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(0), 1);
c@131 117 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(1), 1);
c@131 118 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(-2), 1);
c@131 119 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(2), 2);
c@131 120 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(3), 4);
c@131 121 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(11), 8);
c@131 122 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(12), 16);
c@131 123 BOOST_CHECK_EQUAL(MathUtilities::nearestPowerOfTwo(16), 16);
c@131 124 }
c@131 125
c@131 126 BOOST_AUTO_TEST_CASE(factorial)
c@131 127 {
c@131 128 BOOST_CHECK_EQUAL(MathUtilities::factorial(-10), 0.0);
c@131 129 BOOST_CHECK_EQUAL(MathUtilities::factorial(0), 1.0);
c@131 130 BOOST_CHECK_EQUAL(MathUtilities::factorial(1), 1.0);
c@131 131 BOOST_CHECK_EQUAL(MathUtilities::factorial(2), 2.0);
c@131 132 BOOST_CHECK_EQUAL(MathUtilities::factorial(3), 6.0);
c@131 133 BOOST_CHECK_EQUAL(MathUtilities::factorial(4), 24.0);
c@131 134
c@131 135 // Too big for an int, hence double return value from factorial
c@131 136 BOOST_CHECK_EQUAL(MathUtilities::factorial(20), 2432902008176640000.0);
c@131 137 }
c@131 138
c@131 139 BOOST_AUTO_TEST_CASE(gcd)
c@131 140 {
c@131 141 BOOST_CHECK_EQUAL(MathUtilities::gcd(1, 1), 1);
c@131 142 BOOST_CHECK_EQUAL(MathUtilities::gcd(2, 1), 1);
c@131 143 BOOST_CHECK_EQUAL(MathUtilities::gcd(2, 3), 1);
c@131 144 BOOST_CHECK_EQUAL(MathUtilities::gcd(4, 2), 2);
c@131 145 BOOST_CHECK_EQUAL(MathUtilities::gcd(18, 24), 6);
c@131 146 BOOST_CHECK_EQUAL(MathUtilities::gcd(27, 18), 9);
c@131 147 BOOST_CHECK_EQUAL(MathUtilities::gcd(18, 36), 18);
c@131 148 BOOST_CHECK_EQUAL(MathUtilities::gcd(37, 18), 1);
c@131 149 }
c@131 150
c@131 151 BOOST_AUTO_TEST_SUITE_END()
c@131 152
c@131 153