admin/omath
1
0
Fork 0
mirror of https://github.com/orange-cpp/omath.git synced 2026-02-13 07:03:25 +00:00
Code Issues Packages Projects Releases Wiki Activity

added additional methods

Browse Source
This commit is contained in:
Orange
2025-04-29 20:10:17 +03:00
parent 4ef674f7b4
commit 4615769682
2 changed files with 18 additions and 5 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
include/omath/mat.hpp
View File

@@ -257,9 +257,9 @@ namespace omath
if constexpr (Rows > 2)
{
Type det = 0;
for (size_t i = 0; i < Columns; ++i)
for (size_t column = 0; column < Columns; ++column)
{
const Type cofactor = (i % 2 == 0 ? 1 : -1) * At(0, i) * Minor(0, i).Determinant();
const Type cofactor = At(0, column) * AlgComplement(0, column);
det += cofactor;
}
return det;
@@ -268,7 +268,7 @@ namespace omath
}
[[nodiscard]]
constexpr Mat<Rows - 1, Columns - 1, Type, StoreType> Minor(const size_t row, const size_t column) const
constexpr Mat<Rows - 1, Columns - 1, Type, StoreType> Strip(const size_t row, const size_t column) const
{
static_assert(Rows-1 > 0 && Columns-1 > 0);
Mat<Rows - 1, Columns - 1, Type, StoreType> result;
@@ -288,6 +288,19 @@ namespace omath
return result;
}
[[nodiscard]]
constexpr Type Minor(const size_t row, const size_t column) const
{
return Strip(row, column).Determinant();
}
[[nodiscard]]
constexpr Type AlgComplement(const size_t row, const size_t column) const
{
const auto minor = Minor(row, column);
return (row + column + 2) % 2 == 0 ? minor: -minor;
}
[[nodiscard]]
constexpr const std::array<Type, Rows * Columns>& RawArray() const
{

4
tests/general/unit_test_mat.cpp
View File

@@ -180,10 +180,10 @@ TEST(UnitTestMatStandalone, Determinant_3x3)
}
// Test Minor for 3x3 matrix
TEST(UnitTestMatStandalone, Minor_3x3)
TEST(UnitTestMatStandalone, Strip_3x3)
{
constexpr Mat<3, 3> m{{3, 0, 2}, {2, 0, -2}, {0, 1, 1}};
auto minor = m.Minor(0, 0);
auto minor = m.Strip(0, 0);
EXPECT_EQ(minor.RowCount(), 2);
EXPECT_EQ(minor.ColumnsCount(), 2);
EXPECT_FLOAT_EQ(minor.At(0, 0), 0.0f);