backend/mplapack/binary80.hpp

Go to the documentation of this file.

#pragma once
 
#include <span>
#include <uni20/core/runtime_precision.hpp>
#include <uni20/core/types.hpp>
#include <vector>
 
namespace uni20::mplapack
{
void gemm(std::size_t m, std::size_t n, std::size_t k, float80 const& alpha, std::span<float80 const> a,
          std::span<float80 const> b, float80 const& beta, std::span<float80> c);
void gemm(std::size_t m, std::size_t n, std::size_t k, complex160 const& alpha, std::span<complex160 const> a,
          std::span<complex160 const> b, complex160 const& beta, std::span<complex160> c);
std::size_t getrf(std::size_t n, std::span<float80> a, std::vector<std::size_t>& pivots);
#if UNI20_HAS_MPLAPACK_FLOAT80_COMPLEX
std::size_t getrf(std::size_t n, std::span<complex160> a, std::vector<std::size_t>& pivots);
#endif
void getrs(std::size_t n, std::size_t nrhs, std::span<float80 const> a, std::span<std::size_t const> pivots,
           std::span<float80> b);
#if UNI20_HAS_MPLAPACK_FLOAT80_COMPLEX
void getrs(std::size_t n, std::size_t nrhs, std::span<complex160 const> a, std::span<std::size_t const> pivots,
           std::span<complex160> b);
#endif
void qr(std::size_t m, std::size_t n, std::span<float80> a, std::span<float80> q, std::span<float80> triangular,
        Precision p);
 
void lq(std::size_t m, std::size_t n, std::span<float80> a, std::span<float80> q, std::span<float80> triangular,
        Precision p);
 
void eigh(std::size_t n, std::span<float80> a, std::span<float80> values, bool vectors, char triangle, Precision p);
 
void svd(std::size_t m, std::size_t n, std::span<float80> a, std::span<float80> values, std::span<float80> u,
         std::span<float80> vh, char jobu, char jobvt, Precision p);
 
void eigen(std::size_t n, std::span<float80> a, std::span<complex160> values, std::span<complex160> vectors,
           bool compute_vectors, Precision p);
 
#if UNI20_HAS_MPLAPACK_FLOAT80_COMPLEX
void qr(std::size_t m, std::size_t n, std::span<complex160> a, std::span<complex160> q,
        std::span<complex160> triangular, Precision p);
 
void lq(std::size_t m, std::size_t n, std::span<complex160> a, std::span<complex160> q,
        std::span<complex160> triangular, Precision p);
 
void eigh(std::size_t n, std::span<complex160> a, std::span<float80> values, bool vectors, char triangle, Precision p);
 
void svd(std::size_t m, std::size_t n, std::span<complex160> a, std::span<float80> values, std::span<complex160> u,
         std::span<complex160> vh, char jobu, char jobvt, Precision p);
 
void eigen(std::size_t n, std::span<complex160> a, std::span<complex160> values, std::span<complex160> vectors,
           bool compute_vectors, Precision p);
#endif
 
std::size_t gelss(std::size_t m, std::size_t n, std::size_t nrhs, std::span<float80> a, std::span<float80> b,
                  std::span<float80> values, float80 const& cutoff, std::size_t& rank, Precision p);
float80 gecon(std::size_t n, std::span<float80 const> a, float80 const& anorm, char norm, Precision p);
#if UNI20_HAS_MPLAPACK_FLOAT80_COMPLEX
std::size_t gelss(std::size_t m, std::size_t n, std::size_t nrhs, std::span<complex160> a, std::span<complex160> b,
                  std::span<float80> values, float80 const& cutoff, std::size_t& rank, Precision p);
float80 gecon(std::size_t n, std::span<complex160 const> a, float80 const& anorm, char norm, Precision p);
#endif
} // namespace uni20::mplapack