norm.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_COMPLEX_FUNCTIONS_NORM_H
#define pRC_CORE_COMPLEX_FUNCTIONS_NORM_H
 
#include <prc/core/complex/complex.hpp>
 
namespace pRC
{
    template<Index P, Index Q = P, IsComplex T>
        requires(P == 1 && Q == 1)
    static inline constexpr auto norm(T const &a)
    {
        return abs(a.real()) + abs(a.imag());
    }
 
    template<Index P, Index Q, IsComplex T>
        requires(P == 2 && Q == 1)
    static inline constexpr auto norm(T const &a)
    {
        return a.real() * a.real() + a.imag() * a.imag();
    }
 
    template<Index P = 2, Index Q = P, IsComplex T>
        requires(P == 2 && Q == 2)
    static inline constexpr auto norm(T const &a)
    {
        return sqrt(norm<2, 1>(a));
    }
}
#endif // pRC_CORE_COMPLEX_FUNCTIONS_NORM_H