trace.h

Go to the documentation of this file.

#ifndef FUSED_ALGORITHMS_TRACE_H
#define FUSED_ALGORITHMS_TRACE_H
 
#include "config.h"
#include "fused/fused_matrix.h"
 
namespace matrix_algorithms
{
 
    template <typename T, my_size_t N>
    T trace(const FusedMatrix<T, N, N> &A)
    {
        T sum = T(0);
 
        for (my_size_t i = 0; i < N; ++i)
        {
            sum += A(i, i);
        }
 
        return sum;
    }
 
} // namespace matrix_algorithms
 
#endif // FUSED_ALGORITHMS_TRACE_H