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

Implement matrix Frobenius (L2) norm

Browse source
This commit is contained in:
Nikita Lisitsa 2021-01-16 22:43:22 +03:00
parent 7698e54b18
commit 71bafff64d
1 changed files with 17 additions and 1 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

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

@ -2,6 +2,7 @@
#include <psemek/geom/detail/array.hpp> #include <psemek/geom/detail/array.hpp>
#include <psemek/geom/vector.hpp> #include <psemek/geom/vector.hpp>
#include <psemek/geom/math.hpp>
#include <algorithm> #include <algorithm>
@ -39,7 +40,6 @@ namespace psemek::geom
static matrix scalar(T const & s); static matrix scalar(T const & s);
}; };
template <typename T, std::size_t R, std::size_t C> template <typename T, std::size_t R, std::size_t C>
matrix<T, R, C> matrix<T, R, C>::zero() matrix<T, R, C> matrix<T, R, C>::zero()
{ {
@ -269,6 +269,22 @@ namespace psemek::geom
return result; return result;
} }
template <typename T, std::size_t R, std::size_t C>
T frobenius_norm_sqr(matrix<T, R, C> const & m)
{
T r{0};
for (std::size_t i = 0; i < R; ++i)
for (std::size_t j = 0; j < C; ++j)
r += sqr(m[i][j]);
return r;
}
template <typename T, std::size_t R, std::size_t C>
T frobenius_norm(matrix<T, R, C> const & m)
{
return std::sqrt(frobenius_norm_sqr(m));
}
template <typename Stream, typename T, std::size_t R, std::size_t C> template <typename Stream, typename T, std::size_t R, std::size_t C>
Stream & operator << (Stream & os, matrix<T, R, C> const & m) Stream & operator << (Stream & os, matrix<T, R, C> const & m)
{ {