From 1528d942934cdafeaac311270f5ed131859f0b25 Mon Sep 17 00:00:00 2001 From: Mateusz Pusz Date: Sat, 27 May 2023 13:16:16 +0200 Subject: [PATCH] chore: `linear_algebra` example removed --- example/linear_algebra.cpp | 450 ------------------------------------- 1 file changed, 450 deletions(-) delete mode 100644 example/linear_algebra.cpp diff --git a/example/linear_algebra.cpp b/example/linear_algebra.cpp deleted file mode 100644 index 9ec0a94b..00000000 --- a/example/linear_algebra.cpp +++ /dev/null @@ -1,450 +0,0 @@ -// The MIT License (MIT) -// -// Copyright (c) 2018 Mateusz Pusz -// -// Permission is hereby granted, free of charge, to any person obtaining a copy -// of this software and associated documentation files (the "Software"), to deal -// in the Software without restriction, including without limitation the rights -// to use, copy, modify, merge, publish, distribute, sublicense, and/or sell -// copies of the Software, and to permit persons to whom the Software is -// furnished to do so, subject to the following conditions: -// -// The above copyright notice and this permission notice shall be included in all -// copies or substantial portions of the Software. -// -// THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR -// IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, -// FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE -// AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER -// LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, -// OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE -// SOFTWARE. - -#include -// linear_algebra.hpp has to be included first otherwise the header will fail to compile! -#include -#include -#include -#include -#include -#include -#include -#include - -namespace STD_LA { - -template -std::ostream& operator<<(std::ostream& os, const vector& v) -{ - os << "|"; - for (auto i = 0U; i < v.size(); ++i) { - os << STD_FMT::format(" {:>9}", v(i)); - } - os << " |"; - return os; -} - -template -std::ostream& operator<<(std::ostream& os, const matrix& v) -{ - for (auto i = 0U; i < v.rows(); ++i) { - os << "|"; - for (auto j = 0U; j < v.columns(); ++j) { - os << STD_FMT::format(" {:>9}", v(i, j)); - } - os << (i != v.rows() - 1U ? " |\n" : " |"); - } - return os; -} - -} // namespace STD_LA - -template -using vector = std::math::fs_vector; - -template -inline constexpr bool mp_units::is_vector> = true; - -namespace { - -using namespace mp_units; -using namespace mp_units::si::unit_symbols; - -// template -// using matrix = std::math::fs_matrix; - - -// constexpr auto mm = isq::length[si::milli]; -// constexpr auto m = isq::length[si::metre]; -// constexpr auto km = isq::length[si::kilo]; -// constexpr auto N = isq::force[si::newton]; - -// template -// using length_v = quantity>; - -// template -// using force_v = quantity>; - -void quantity_of_vector_cast() -{ - std::cout << "\nquantity_of_vector_cast:\n"; - - auto v = vector{1001, 1002, 1003} * isq::position_vector[m]; - auto u = vector{3, 2, 1} * isq::position_vector[km]; - - std::cout << "v = " << v << "\n"; - std::cout << "u = " << u << "\n"; - - std::cout << "v[km] = " << value_cast(v) << "\n"; - std::cout << "u[m] = " << u[m] << "\n"; -} - -void quantity_of_vector_multiply_by_scalar_value() -{ - std::cout << "\nquantity_of_vector_multiply_scalar_value:\n"; - - auto v = vector{1, 2, 3} * isq::position_vector[m]; - - std::cout << "v = " << v << "\n"; - - std::cout << "2 * v = " << 2 * v << "\n"; -} - -void quantity_of_vector_divide_by_scalar_value() -{ - std::cout << "\nquantity_of_vector_divide_scalar_value:\n"; - - auto v = vector{2, 4, 6} * isq::position_vector[m]; - - std::cout << "v = " << v << "\n"; - - std::cout << "v / 2 = " << v / 2 << "\n"; -} - -void quantity_of_vector_add() -{ - std::cout << "\nquantity_of_vector_add:\n"; - - // TODO Does not compile (LA bug) - // auto v = vector{1, 2, 3} * isq::position_vector[m]; - // auto u = vector{3, 2, 1} * isq::position_vector[m]; - // auto t = vector{3, 2, 1} * isq::position_vector[km]; - - // quantity> v{vector{1, 2, 3}}; // should not compile - auto v = vector{1, 2, 3} * isq::position_vector[m]; - auto u = vector{3, 2, 1} * isq::position_vector[m]; - auto t = vector{3, 2, 1} * isq::position_vector[km]; - - std::cout << "v = " << v << "\n"; - std::cout << "u = " << u << "\n"; - std::cout << "t = " << t << "\n"; - - std::cout << "v + u = " << v + u << "\n"; - std::cout << "v + t = " << v + t << "\n"; -} - -void quantity_of_vector_subtract() -{ - std::cout << "\nquantity_of_vector_subtract:\n"; - - auto v = vector{1, 2, 3} * isq::position_vector[m]; - auto u = vector{3, 2, 1} * isq::position_vector[m]; - auto t = vector{3, 2, 1} * isq::position_vector[km]; - - std::cout << "v = " << v << "\n"; - std::cout << "u = " << u << "\n"; - std::cout << "t = " << t << "\n"; - - std::cout << "v - u = " << v - u << "\n"; - std::cout << "v - t = " << v - t << "\n"; -} - -void quantity_of_vector_divide_by_scalar_quantity() -{ - std::cout << "\nquantity_of_vector_divide_scalar:\n"; - - auto pos = vector{30, 20, 10} * isq::position_vector[km]; - auto dur = 2 * isq::duration[h]; - quantity> v = pos / dur; - - std::cout << "pos = " << pos << "\n"; - std::cout << "dur = " << dur << "\n"; - - std::cout << "pos / dur = " << pos / dur << "\n"; - std::cout << "v = " << v << "\n"; -} - - -// void quantity_of_vector_multiply_same() -// { -// std::cout << "\nquantity_of_vector_multiply_same:\n"; - -// length_v<> v(vector<>{1, 2, 3}); -// length_v<> u(vector<>{3, 2, 1}); - -// std::cout << "v = " << v << "\n"; -// std::cout << "u = " << u << "\n"; - -// std::cout << "v * u = " << v * u << "\n"; -// std::cout << "2 * m * v = " << 2. * m * v << "\n"; -// } - -// void quantity_of_vector_multiply_different() -// { -// std::cout << "\nquantity_of_vector_multiply_different:\n"; - -// force_v<> v(vector<>{1, 2, 3}); -// length_v<> u(vector<>{3, 2, 1}); - -// std::cout << "v = " << v << "\n"; -// std::cout << "u = " << u << "\n"; - -// std::cout << "v * u = " << v * u << "\n"; -// std::cout << "2 * N * u = " << 2. * N * u << "\n"; -// std::cout << "2 * u = " << 2 * u << "\n"; -// } - -// void quantity_of_vector_divide_by_scalar() -// { -// std::cout << "\nquantity_of_vector_divide_by_scalar:\n"; - -// length_v<> v(vector<>{4, 8, 12}); - -// std::cout << "v = " << v << "\n"; - -// // TODO Uncomment when bug in the LA is fixed -// // std::cout << "v / 2 * s = " << v / 2 * s << "\n"; -// // std::cout << "v / 2 = " << v / 2 << "\n"; -// } - -void quantity_of_vector_tests() -{ - quantity_of_vector_cast(); - quantity_of_vector_multiply_by_scalar_value(); - quantity_of_vector_divide_by_scalar_value(); - quantity_of_vector_add(); - quantity_of_vector_subtract(); - quantity_of_vector_divide_by_scalar_quantity(); - // quantity_of_vector_multiply_same(); - // quantity_of_vector_multiply_different(); - // quantity_of_vector_divide_by_scalar(); -} - -// template -// using length_m = quantity>; - -// void quantity_of_matrix_add() -// { -// std::cout << "\nquantity_of_matrix_add:\n"; - -// length_m<> v(matrix<>{{1, 2, 3}, {4, 5, 6}, {7, 8, 9}}); -// length_m<> u(matrix<>{{3, 2, 1}, {3, 2, 1}, {3, 2, 1}}); -// length_m> t(matrix<>{{3, 2, 1}, {3, 2, 1}, {3, 2, 1}}); - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; -// std::cout << "t =\n" << t << "\n"; - -// std::cout << "v + u =\n" << v + u << "\n"; -// std::cout << "v + t =\n" << v + t << "\n"; - -// // TODO Uncomment when fixed in the LA lib -// // std::cout << "v[mm] =\n" << matrix>(v) << "\n"; -// } - -// void quantity_of_matrix_multiply_same() -// { -// std::cout << "\nquantity_of_matrix_multiply_same:\n"; - -// length_m<> v(matrix<>{{1, 2, 3}, {4, 5, 6}, {7, 8, 9}}); -// length_v<> u(vector<>{3, 2, 1}); - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; - -// std::cout << "v * u =\n" << v * u << "\n"; -// std::cout << "2 * m * u =\n" << 2. * m * u << "\n"; -// } - -// void quantity_of_matrix_multiply_different() -// { -// std::cout << "\nquantity_of_matrix_multiply_different:\n"; - -// force_v<> v(vector<>{1, 2, 3}); -// length_m<> u(matrix<>{{1, 2, 3}, {4, 5, 6}, {7, 8, 9}}); - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; - -// std::cout << "v * u =\n" << v * u << "\n"; -// std::cout << "2 * N * u =\n" << 2. * N * u << "\n"; -// std::cout << "2 * u =\n" << 2 * u << "\n"; -// } - -// void quantity_of_matrix_divide_by_scalar() -// { -// std::cout << "\nquantity_of_matrix_divide_by_scalar:\n"; - -// length_m<> v(matrix<>{{2, 4, 6}, {4, 6, 8}, {8, 4, 2}}); - -// std::cout << "v =\n" << v << "\n"; - -// // TODO Uncomment when bug in the LA is fixed -// // std::cout << "v / (2 * s) =\n" << v / (2 * s) << "\n"; -// // std::cout << "v / 2 =\n" << v / 2 << "\n"; -// } - -// void quantity_of_matrix_tests() -// { -// quantity_of_matrix_add(); -// quantity_of_matrix_multiply_same(); -// quantity_of_matrix_multiply_different(); -// quantity_of_matrix_divide_by_scalar(); -// } - - -// void vector_of_quantity_add() -// { -// std::cout << "\nvector_of_quantity_add:\n"; - -// vector> v = {1 * m, 2 * m, 3 * m}; -// vector> u = {3 * m, 2 * m, 1 * m}; -// vector> t = {3 * km, 2 * km, 1 * km}; - -// std::cout << "v = " << v << "\n"; -// std::cout << "u = " << u << "\n"; -// std::cout << "t = " << t << "\n"; - -// std::cout << "v + u = " << v + u << "\n"; -// std::cout << "v + t = " << v + t << "\n"; -// std::cout << "t[m] = " << vector>(t) << "\n"; -// } - -// void vector_of_quantity_multiply_same() -// { -// std::cout << "\nvector_of_quantity_multiply_same:\n"; - -// vector> v = {1 * m, 2 * m, 3 * m}; -// vector> u = {3 * m, 2 * m, 1 * m}; - -// std::cout << "v = " << v << "\n"; -// std::cout << "u = " << u << "\n"; - -// std::cout << "v * u = " << v * u << "\n"; -// std::cout << "2 * m * v = " << 2. * m * v << "\n"; -// } - -// void vector_of_quantity_multiply_different() -// { -// std::cout << "\nvector_of_quantity_multiply_different:\n"; - -// vector> v = {1 * N, 2 * N, 3 * N}; -// vector> u = {3 * m, 2 * m, 1 * m}; - -// std::cout << "v = " << v << "\n"; -// std::cout << "u = " << u << "\n"; - -// std::cout << "v * u = " << v * u << "\n"; -// std::cout << "2 * N * u = " << 2. * N * u << "\n"; -// std::cout << "2 * u = " << 2 * u << "\n"; -// } - -// void vector_of_quantity_divide_by_scalar() -// { -// std::cout << "\nvector_of_quantity_divide_by_scalar:\n"; - -// vector> v = {4 * m, 8 * m, 12 * m}; - -// std::cout << "v = " << v << "\n"; - -// // TODO Uncomment when bug in the LA is fixed -// // std::cout << "v / (2 * s) = " << v / (2 * s) << "\n"; -// // std::cout << "v / 2 = " << v / 2 << "\n"; -// } - -// void vector_of_quantity_tests() -// { -// vector_of_quantity_add(); -// vector_of_quantity_multiply_same(); -// vector_of_quantity_multiply_different(); -// vector_of_quantity_divide_by_scalar(); -// } - -// void matrix_of_quantity_add() -// { -// std::cout << "\nmatrix_of_quantity_add:\n"; - -// matrix> v = {{1 * m, 2 * m, 3 * m}, {4 * m, 5 * m, 6 * m}, {7 * m, 8 * m, 9 * m}}; -// matrix> u = {{3 * m, 2 * m, 1 * m}, {3 * m, 2 * m, 1 * m}, {3 * m, 2 * m, 1 * m}}; -// matrix> t = {{3 * mm, 2 * mm, 1 * mm}, {3 * mm, 2 * mm, 1 * mm}, {3 * mm, 2 * mm, 1 * mm}}; - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; -// std::cout << "t =\n" << t << "\n"; - -// std::cout << "v + u =\n" << v + u << "\n"; -// std::cout << "v + t =\n" << v + t << "\n"; - -// // TODO Uncomment when fixed in the LA lib -// // std::cout << "v[mm] =\n" << matrix>(v) << "\n"; -// } - -// void matrix_of_quantity_multiply_same() -// { -// std::cout << "\nmatrix_of_quantity_multiply_same:\n"; - -// matrix> v = {{1 * m, 2 * m, 3 * m}, {4 * m, 5 * m, 6 * m}, {7 * m, 8 * m, 9 * m}}; -// vector> u = {3 * m, 2 * m, 1 * m}; - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; - -// std::cout << "v * u =\n" << v * u << "\n"; -// std::cout << "2 * m * u =\n" << 2. * m * u << "\n"; -// } - -// void matrix_of_quantity_multiply_different() -// { -// std::cout << "\nmatrix_of_quantity_multiply_different:\n"; - -// vector> v = {1 * N, 2 * N, 3 * N}; -// matrix> u = {{1 * m, 2 * m, 3 * m}, {4 * m, 5 * m, 6 * m}, {7 * m, 8 * m, 9 * m}}; - -// std::cout << "v =\n" << v << "\n"; -// std::cout << "u =\n" << u << "\n"; - -// std::cout << "v * u =\n" << v * u << "\n"; -// std::cout << "2 * N * u =\n" << 2. * N * u << "\n"; -// std::cout << "2 * u =\n" << 2 * u << "\n"; -// } - -// void matrix_of_quantity_divide_by_scalar() -// { -// std::cout << "\nmatrix_of_quantity_divide_by_scalar:\n"; - -// matrix> v = {{2 * m, 4 * m, 6 * m}, {4 * m, 6 * m, 8 * m}, {8 * m, 4 * m, 2 * m}}; - -// std::cout << "v =\n" << v << "\n"; - -// // TODO Uncomment when bug in the LA is fixed -// // std::cout << "v / (2 * s) =\n" << v / (2 * s) << "\n"; -// // std::cout << "v / 2 =\n" << v / 2 << "\n"; -// } - -// void matrix_of_quantity_tests() -// { -// matrix_of_quantity_add(); -// matrix_of_quantity_multiply_same(); -// matrix_of_quantity_multiply_different(); -// matrix_of_quantity_divide_by_scalar(); -// } - -} // namespace - -int main() -{ - quantity_of_vector_tests(); - // quantity_of_matrix_tests(); - // vector_of_quantity_tests(); - // matrix_of_quantity_tests(); -}