matrix_decompositions.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <cstdint>
 
#include "pulsatrix/tensor.hpp"
 
namespace pulsatrix {
 
struct EigenResult {
    Tensor values;
    Tensor vectors;
};
 
struct DominantEigenResult {
    float value;
    Tensor vector;
    int64_t iterations;
    bool converged;
};
 
struct QRResult {
    Tensor q;
    Tensor r;
};
 
struct SVDResult {
    Tensor u;
    Tensor s;
    Tensor v;
};
 
[[nodiscard]] EigenResult SymmetricEigen(const Tensor& a);
 
[[nodiscard]] DominantEigenResult PowerIteration(const Tensor& a, int64_t max_iterations = 1000,
                                                 float tolerance = 1e-6f, uint64_t seed = 0);
 
[[nodiscard]] QRResult QR(const Tensor& a);
 
[[nodiscard]] SVDResult SVD(const Tensor& a);
 
}  // namespace pulsatrix