lisyarus/psemek
1
0
Fork 1
Code Issues Pull requests 1 Projects Releases Packages Wiki Activity Actions

Add geom::matrix::diagonal

Browse source
This commit is contained in:
Nikita Lisitsa 2021-12-30 17:21:10 +03:00
parent 64b3e876ad
commit cae538d7a5
1 changed files with 10 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

10
libs/geom/include/psemek/geom/matrix.hpp
View file

@ -48,6 +48,7 @@ namespace psemek::geom
static matrix zero(); static matrix zero();
static matrix identity(); static matrix identity();
static matrix scalar(T const & s); static matrix scalar(T const & s);
static matrix diagonal(vector<T, std::min(R, C)> const & d);
}; };
template <typename T, std::size_t R, std::size_t C> template <typename T, std::size_t R, std::size_t C>
@ -74,6 +75,15 @@ namespace psemek::geom
return m; return m;
} }
template <typename T, std::size_t R, std::size_t C>
matrix<T, R, C> matrix<T, R, C>::diagonal(vector<T, std::min(R, C)> const & d)
{
matrix<T, R, C> m = zero();
for (std::size_t i = 0; i < std::min(R, C); ++i)
m[i][i] = d[i];
return m;
}
template <typename T, std::size_t R, std::size_t C> template <typename T, std::size_t R, std::size_t C>
bool operator == (matrix<T, R, C> const & m1, matrix<T, R, C> const & m2) bool operator == (matrix<T, R, C> const & m1, matrix<T, R, C> const & m2)
{ {