rank_update.h

Go to the documentation of this file.

#ifndef FUSED_ALGORITHMS_RANK_UPDATE_H
#define FUSED_ALGORITHMS_RANK_UPDATE_H
 
#include "config.h"
#include "fused/fused_matrix.h"
 
namespace matrix_algorithms
{
 
    template <typename T, my_size_t N>
    FusedMatrix<T, N, N> symmetric_rank_k_update(
        const FusedMatrix<T, N, N> &F,
        const FusedMatrix<T, N, N> &P,
        const FusedMatrix<T, N, N> &Q)
    {
        FusedMatrix<T, N, N> result;
 
        auto tmp = FusedMatrix<T, N, N>::matmul(F, P);
        result = FusedMatrix<T, N, N>::matmul(tmp, F.transpose_view()) + Q;
        return result;
    }
 
    template <typename T, my_size_t N>
    FusedMatrix<T, N, N> symmetric_rank_k_update(
        const FusedMatrix<T, N, N> &F,
        const FusedMatrix<T, N, N> &P)
    {
        FusedMatrix<T, N, N> result;
 
        auto tmp = FusedMatrix<T, N, N>::matmul(F, P);
        result = FusedMatrix<T, N, N>::matmul(tmp, F.transpose_view());
        return result;
    }
 
} // namespace matrix_algorithms
 
#endif // FUSED_ALGORITHMS_RANK_UPDATE_H