Scroll to navigation

RsbMatrix< NT >(3) Library Functions Manual RsbMatrix< NT >(3)

NAME

RsbMatrix - Represent a sparse matrix in RSB format by means of librsb.

SYNOPSIS

#include <rsb.hpp>

Public Types


enum RsbSym { IsGen = RSB_FLAG_NOFLAGS, IsHer = RSB_FLAG_HERMITIAN, IsSym = RSB_FLAG_SYMMETRIC, IsTri = RSB_FLAG_TRIANGULAR }

Public Member Functions


RsbMatrix (rsb_coo_idx_t nrA, rsb_coo_idx_t ncA, const RsbSym sym=IsGen)
RsbMatrix (rsb_coo_idx_t nrA, const rsb_coo_idx_t *RP, const rsb_coo_idx_t *JA, const NT *VA, const RsbSym sym=IsGen)
RsbMatrix (const rsb_coo_idx_t *IA, const rsb_coo_idx_t *JA, const NT *VA, rsb_nnz_idx_t nnzA, const rsb_flags_t flagsA=RSB_FLAG_NOFLAGS)
RsbMatrix (const rsb_char_t *filename, const RsbSym sym=IsGen)
RsbMatrix (const RsbMatrix &A_Rsb, bool do_trans=false, rsb_flags_t flagsA=RSB_FLAG_NOFLAGS)
RsbMatrix (RsbMatrix &&other)
~RsbMatrix (void)
RSBP_RVT RSBP_DEPRECATED Err_t _add (rsb_coo_idx_t i, rsb_coo_idx_t j, NT val)
RSBP_RVT Err_t close (void)
RSBP_RVT RSBP_DEPRECATED Err_t _close (void)
RSBP_RVT RSBP_DEPRECATED Err_t spmv (rsb_trans_t transA, const NT *alphap, const NT *Xp, rsb_coo_idx_t incX, const NT *betap, NT *Yp, rsb_coo_idx_t incY) const
RSBP_RVT Err_t spmv (rsb_trans_t transA, const NT alpha, const NT *Xp, rsb_coo_idx_t incX, const NT beta, NT *Yp, rsb_coo_idx_t incY) const
RSBP_RVT Err_t spmv (rsb_trans_t transA, const NT alpha, const NT *Xp, const NT beta, NT *Yp) const
RSBP_RVT Err_t spmv (NT *y, const NT *x, bool do_trans=false) const
RSBP_RVT RSBP_DEPRECATED Err_t spmm (rsb_trans_t transA, const NT *alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, rsb_nnz_idx_t ldB, const NT *betap, NT *Cp, rsb_nnz_idx_t ldC) const
RSBP_RVT Err_t spmm (rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, rsb_nnz_idx_t ldB, const NT beta, NT *Cp, rsb_nnz_idx_t ldC) const
RSBP_RVT Err_t spmm (rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, const NT beta, NT *Cp) const
RSBP_RVT RSBP_DEPRECATED Err_t spsm (rsb_trans_t transT, const NT *alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *betap, const NT *Bp, rsb_nnz_idx_t ldB, NT *Cp, rsb_nnz_idx_t ldC) const
RSBP_RVT Err_t spsm (rsb_trans_t transT, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT beta, const NT *Bp, rsb_nnz_idx_t ldB, NT *Cp, rsb_nnz_idx_t ldC) const
RSBP_RVT Err_t spsm (rsb_trans_t transT, const NT alpha, rsb_coo_idx_t nrhs, const NT *Bp, NT *Cp) const
RSBP_RVT Err_t spsm (NT *y, const NT *x, rsb_coo_idx_t nrhs, bool do_trans=false) const
RSBP_RVT Err_t spsm (NT *y, rsb_coo_idx_t nrhs, bool do_trans=false) const
RSBP_RVT RSBP_DEPRECATED Err_t spsv (rsb_trans_t transT, const NT *alphap, const NT *Xp, rsb_coo_idx_t incX, NT *Yp, rsb_coo_idx_t incY) const
RSBP_RVT Err_t spsv (rsb_trans_t transT, const NT alpha, const NT *Xp, NT *Yp) const
RSBP_RVT Err_t spsv (NT *y, const NT *x, bool do_trans=false) const
RSBP_RVT Err_t spsv (NT *y, bool do_trans=false) const
size_t get_info_size_t (enum rsb_mif_t mif) const
rsb_flags_t get_info_rsb_flags_t (enum rsb_mif_t mif) const
rsb_blk_idx_t get_info_blk_t (enum rsb_mif_t mif) const
rsb_nnz_idx_t get_info_nnz_t (enum rsb_mif_t mif) const
rsb_flags_t get_flags_t (enum rsb_mif_t mif) const
rsb_type_t get_type_t (enum rsb_mif_t mif) const
rsb_coo_idx_t get_info_coo_t (enum rsb_mif_t mif) const
size_t _get_index_storage_bytes (void) const
size_t _get_storage_bytes (void) const
rsb_nnz_idx_t nnz (void) const
rsb_blk_idx_t blocks (void) const
rsb_coo_idx_t rows (void) const
rsb_coo_idx_t cols (void) const
RSBP_RVT Err_t get_vals (NT *VA, const rsb_coo_idx_t *IA, const rsb_coo_idx_t *JA, rsb_nnz_idx_t nnz, rsb_flags_t flags) const
NT get_val (const rsb_coo_idx_t i, const rsb_coo_idx_t j, rsb_flags_t flags=RSB_FLAG_NOFLAGS) const
RSBP_RVT Err_t set_val (const NT val, const rsb_coo_idx_t i, const rsb_coo_idx_t j, rsb_flags_t flags=RSB_FLAG_NOFLAGS)
RSBP_RVT Err_t set_vals (const NT *VA, const rsb_coo_idx_t *IA, const rsb_coo_idx_t *JA, rsb_nnz_idx_t nnz, rsb_flags_t flags)
RSBP_RVT Err_t get_vec (NT *Dp, enum rsb_extff_t flags) const
RSBP_RVT RSBP_DEPRECATED Err_t get_coo (rsb_trans_t transA, NT *VA, rsb_coo_idx_t *IA, rsb_coo_idx_t *JA, rsb_flags_t flags) const
RSBP_RVT RSBP_DEPRECATED Err_t get_csr (rsb_trans_t transA, NT *VA, rsb_coo_idx_t *RP, rsb_coo_idx_t *JA, rsb_flags_t flags) const
RSBP_RVT RSBP_DEPRECATED Err_t get_rows_sparse (rsb_trans_t transA, const NT *alphap, NT *VA, rsb_coo_idx_t *IA, rsb_coo_idx_t *JA, rsb_coo_idx_t frA, rsb_coo_idx_t lrA, rsb_nnz_idx_t *rnzp, rsb_flags_t flags) const
RSBP_RVT Err_t upd_vals (enum rsb_elopf_t elop_flags, const NT &omega)
RSBP_RVT Err_t upd_vals (enum rsb_elopf_t elop_flags, const NT *omegap)
RSBP_RVT Err_t get_nrm (NT *Np, enum rsb_extff_t flags) const
rsb_type_t rsbtype (void) const
rsb_flags_t rsbflags (void) const
rsb_string_t get_info_str (const char *key) const
RSBP_RVT Err_t get_info (enum rsb_mif_t miflags, void *minfop) const
rsb_string_t _info (void) const
RSBP_RVT RSBP_DEPRECATED Err_t tune_spsm_threads (rsb_real_t *sfp=RSBP_NULL, rsb_int_t *tnp=RSBP_NULL, rsb_int_t maxr=0, rsb_time_t maxt=0, rsb_trans_t transA=RSB_TRANSPOSITION_N, const NT *alphap=RSBP_NULL, rsb_coo_idx_t nrhs=1, rsb_flags_t order=RSB_FLAG_WANT_COLUMN_MAJOR_ORDER, const NT *Bp=RSBP_NULL, rsb_nnz_idx_t ldB=0, const NT *betap=RSBP_NULL, NT *Cp=RSBP_NULL, rsb_nnz_idx_t ldC=0) const
RSBP_RVT RSBP_DEPRECATED Err_t tune_spmm_threads (rsb_real_t *sfp=RSBP_NULL, rsb_int_t *tnp=RSBP_NULL, rsb_int_t maxr=0, rsb_time_t maxt=0, rsb_trans_t transA=RSB_TRANSPOSITION_N, const NT *alphap=RSBP_NULL, rsb_coo_idx_t nrhs=1, rsb_flags_t order=RSB_FLAG_WANT_COLUMN_MAJOR_ORDER, const NT *Bp=RSBP_NULL, rsb_nnz_idx_t ldB=0, const NT *betap=RSBP_NULL, NT *Cp=RSBP_NULL, rsb_nnz_idx_t ldC=0) const
RSBP_RVT RSBP_DEPRECATED Err_t tune_spmm (rsb_real_t *sfp, rsb_int_t *tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT *alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, rsb_nnz_idx_t ldB, const NT *betap, NT *Cp, rsb_nnz_idx_t ldC)
RSBP_RVT RSBP_DEPRECATED Err_t tune_spmm (rsb_real_t *sfp, rsb_int_t *tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, rsb_nnz_idx_t ldB, const NT beta, NT *Cp, rsb_nnz_idx_t ldC)
RSBP_RVT Err_t tune_spmm (rsb_real_t &sf, rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, const NT beta, NT *Cp)
RSBP_RVT RSBP_DEPRECATED Err_t tune_spsm (rsb_real_t *sfp, rsb_int_t *tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT *alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT *Bp, rsb_nnz_idx_t ldB, const NT *betap, NT *Cp, rsb_nnz_idx_t ldC)
RSBP_RVT Err_t file_save (const rsb_char_t *filename=RSBP_NULL) const
RsbMatrix & operator= (const RsbMatrix &A_Rsb)
bool _is_complex (void) const
bool operator== (const RsbMatrix &B_Rsb) const
bool operator!= (const RsbMatrix &B_Rsb) const
NT normOne (void) const
NT normInf (void) const
RSBP_RVT Err_t rndr (const rsb_char_t *filename=RSBP_NULL, rsb_coo_idx_t pmWidth=512, rsb_coo_idx_t pmHeight=512, rsb_marf_t rflags=RSB_MARF_EPS) const

Detailed Description

template<RSBP_Scalar_t NT>

class RsbMatrix< NT >"Represent a sparse matrix in RSB format by means of librsb.

Manage construction, destruction, and std::move of numerical matrices.
Most of the member functions here translate directly to a single function call to librsb (rsb.h), and pass the parameters as they are, so the error checking is done by librsb.
While most of librsb C functions use void* pointers instead of numerical data, RsbMatrix is templated by a type parameter. This introduces type safety at compile time.

Users of member functions can choose among several overloads. So in additional to the more direct overloads passing e.g. $ alpha $ and $ beta $ by reference, here a user can pass them by value.

Parameters

NT the numerical type, at least for the four canonical ones (float, double, std::complex<float>, std::complex<double>); see matrix_supported_numerical_types_section and rsb_type_t for more.

Note

Default error propagation is by exception throw for all constructors and most member functions.
Functions declared to return Err_t can be specialized in rsb_err_t so not to throw exceptions, but to return an error code instead.
Exceptions thrown by member functions (not constructors) can be deactivated at build time by defining RSBP_NOTHROW before including <rsb.hpp>.

One may turn on return error value as default at build time by defining RSBP_WANT_REV=1 .

Warning

The error model is work in progress and subject to change.

Todo

While the rsb.h interface is stable, the rsb.hpp interface is neither stable, nor complete: early users' feedback is very welcome.
Shall all $ alpha $ and $ beta $ be passed by values only? This is natural and fits C++. But rsb_tune_spmm() / rsb_tune_spsm() have more parameters with a nullptr-like `default' value: how to deal with them in RsbMatrix::tune_spmm() / RsbMatrix::tune_spsm() without introducing too many special cases ?
Similarly for the consistency of RsbMatrix::get_flags_t(), RsbMatrix::get_type_t(), RsbMatrix::get_info_coo_t(), and similar: shall they throw an exception if given a flag not matching the return type ?
Or maybe keeping only RsbMatrix::get_info(), with per-type reference overloads ?
Furthermore: if working under C++20, shall RsbLib avoid pointer-based interfaces completely (using std::span only)?

Member Enumeration Documentation

template<RSBP_Scalar_t NT> enum RsbMatrix::RsbSym

Matrix structure: either general, symmetric, hermitian, or triangular (also see rsb_flags_t).

Enumerator

General matrix, no triangle structure or symmetry assumed.
Hermitian ( $ A == A^H $). Please pass only lower/upper triangle.
Symmetric ( $ A == A^T $). Please pass only lower/upper triangle.
Triangular (required for spsv/ spsm).

Constructor & Destructor Documentation

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (rsb_coo_idx_t nrA, rsb_coo_idx_t ncA, const RsbSym sym = IsGen) [inline]

Begin assembling a sparse matrix of given dimensions and type.
Then you can use set_val() or set_vals() repeatedly to populate the matrix.
After populating the matrix, use close() to terminate its assembly.
Example snip from examples/assemble.cpp:



  const rsb_coo_idx_t nrA { 4 }, ncA { 4 };


  RsbMatrix<double> mtx(nrA,ncA); // begin matrix assembly


  // insert elements of a tridiagonal matrix, one by one


  for (auto i = 0; i < nrA; ++i )


    for (auto j = i-1; j <= i+1; ++j )


      if ( i >= 0 && i < nrA )


        if ( j >= 0 && j < ncA )


          mtx.set_val((i+1)*100+(j+1),i,j); // add entry


  mtx.close(); // finish matrix assembly


  assert(mtx.nnz() == 3 * nrA - 2);

Note

Directly based on C function rsb_mtx_alloc_from_coo_begin() from <rsb.h>.

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (rsb_coo_idx_t nrA, const rsb_coo_idx_t * RP, const rsb_coo_idx_t * JA, const NT * VA, const RsbSym sym = IsGen) [inline]

Assemble a sparse matrix given CSR input.
Note

Directly based on C function rsb_mtx_alloc_from_csr_const() from <rsb.h>.

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (const rsb_coo_idx_t * IA, const rsb_coo_idx_t * JA, const NT * VA, rsb_nnz_idx_t nnzA, const rsb_flags_t flagsA = RSB_FLAG_NOFLAGS) [inline]

Assemble a sparse matrix given COO input.
Example snip from examples/misc.cpp:



    const rsb_nnz_idx_t nnzA { 7 };


    const rsb_coo_idx_t nrA { 6 }, ncA { 6 }, nrhs { 1 };


    const std::vector<rsb_coo_idx_t> IA {0,1,2,3,4,5,1}, JA {0,1,2,3,4,5,0};


    const std::vector<double> VA {1,1,1,1,1,1,2}, X(ncA,1);


    RsbMatrix<double> mtx(IA.data(),JA.data(),VA.data(),nnzA);

Note

Directly based on C function rsb_mtx_alloc_from_coo_const() from <rsb.h>.

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (const rsb_char_t * filename, const RsbSym sym = IsGen) [inline]

Assemble a sparse matrix given filename input.
Note

Directly based on C function rsb_file_mtx_load() from <rsb.h>.

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (const RsbMatrix< NT > & A_Rsb, bool do_trans = false, rsb_flags_t flagsA = RSB_FLAG_NOFLAGS) [inline]

Copy a sparse matrix given example input.
Can either clone it, or transpose it or change flags (structure) in the process.

template<RSBP_Scalar_t NT> RsbMatrix< NT >::RsbMatrix (RsbMatrix< NT > && other) [inline]

Move constructor.
The moved matrix object will be invalid afterwards.
Example snip from examples/misc.cpp:



  assert( mtx1.nnz() == nnzA );


  RsbMatrix<double> mtx3 { std::move(mtx1) };


  assert( mtx3.nnz() == nnzA );

template<RSBP_Scalar_t NT> RsbMatrix< NT >::~RsbMatrix (void) [inline]

Destructor.
Frees matrix object memory.
Note

Directly based on C function rsb_mtx_free() from <rsb.h>.

Member Function Documentation

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::_add (rsb_coo_idx_t i, rsb_coo_idx_t j, NT val) [inline]

Deprecated

Use set_val() and set_vals() instead.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::_close (void) [inline]

Deprecated

Use close() instead.

template<RSBP_Scalar_t NT> size_t RsbMatrix< NT >::_get_index_storage_bytes (void) const [inline]

Warning

The name of this member function is expected to change.

template<RSBP_Scalar_t NT> size_t RsbMatrix< NT >::_get_storage_bytes (void) const [inline]

Warning

The name of this member function is expected to change.

template<RSBP_Scalar_t NT> rsb_string_t RsbMatrix< NT >::_info (void) const [inline]

Warning

The name of this member function is expected to change.

template<RSBP_Scalar_t NT> bool RsbMatrix< NT >::_is_complex (void) const [inline]

Warning

The name of this member function is expected to change.

template<RSBP_Scalar_t NT> rsb_blk_idx_t RsbMatrix< NT >::blocks (void) const [inline]

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::close (void) [inline]

Terminate assembly of a previously started and populated matrix.
Shall be called once.

Example snip from examples/assemble.cpp:



  const rsb_coo_idx_t nrA { 4 }, ncA { 4 };


  RsbMatrix<double> mtx(nrA,ncA); // begin matrix assembly


  // insert elements of a tridiagonal matrix, one by one


  for (auto i = 0; i < nrA; ++i )


    for (auto j = i-1; j <= i+1; ++j )


      if ( i >= 0 && i < nrA )


        if ( j >= 0 && j < ncA )


          mtx.set_val((i+1)*100+(j+1),i,j); // add entry


  mtx.close(); // finish matrix assembly


  assert(mtx.nnz() == 3 * nrA - 2);

See also

RsbMatrix::RsbMatrix(rsb_coo_idx_t nrA, rsb_coo_idx_t ncA, const RsbSym sym = IsGen );

Note

Directly based on C function rsb_mtx_alloc_from_coo_end() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_coo_idx_t RsbMatrix< NT >::cols (void) const [inline]

Example snip from examples/twonnz.cpp:



        std::cout << "# Matrix sized " << mtx.rows() << "x" << mtx.cols() << ", " << nnzA << " nnz  built in " << dt << " s and occupies " << mtxocc << " bytes " << std::endl;

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::file_save (const rsb_char_t * filename = RSBP_NULL) const [inline]

Note

Directly based on C function rsb_file_mtx_save() from <rsb.h>.

Example snip from examples/misc.cpp:



    mtx.file_save(); // print to stdout

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::get_coo (rsb_trans_t transA, NT * VA, rsb_coo_idx_t * IA, rsb_coo_idx_t * JA, rsb_flags_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_coo() from <rsb.h>.

Warning

Only transA=RSB_TRANSPOSITION_N currently supported.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::get_csr (rsb_trans_t transA, NT * VA, rsb_coo_idx_t * RP, rsb_coo_idx_t * JA, rsb_flags_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_csr() from <rsb.h>.

Warning

Only transA=RSB_TRANSPOSITION_N currently supported.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> rsb_flags_t RsbMatrix< NT >::get_flags_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::get_info (enum rsb_mif_t miflags, void * minfop) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_blk_idx_t RsbMatrix< NT >::get_info_blk_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_coo_idx_t RsbMatrix< NT >::get_info_coo_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_nnz_idx_t RsbMatrix< NT >::get_info_nnz_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_flags_t RsbMatrix< NT >::get_info_rsb_flags_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> size_t RsbMatrix< NT >::get_info_size_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_string_t RsbMatrix< NT >::get_info_str (const char * key) const [inline]

Note

Directly based on C function rsb_mtx_get_info_str() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::get_nrm (NT * Np, enum rsb_extff_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_nrm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::get_rows_sparse (rsb_trans_t transA, const NT * alphap, NT * VA, rsb_coo_idx_t * IA, rsb_coo_idx_t * JA, rsb_coo_idx_t frA, rsb_coo_idx_t lrA, rsb_nnz_idx_t * rnzp, rsb_flags_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_rows_sparse() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> rsb_type_t RsbMatrix< NT >::get_type_t (enum rsb_mif_t mif) const [inline]

Note

Directly based on C function rsb_mtx_get_info() from <rsb.h>.

template<RSBP_Scalar_t NT> NT RsbMatrix< NT >::get_val (const rsb_coo_idx_t i, const rsb_coo_idx_t j, rsb_flags_t flags = RSB_FLAG_NOFLAGS) const [inline]

Note

Directly based on C function rsb_mtx_get_vals() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::get_vals (NT * VA, const rsb_coo_idx_t * IA, const rsb_coo_idx_t * JA, rsb_nnz_idx_t nnz, rsb_flags_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_vals() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::get_vec (NT * Dp, enum rsb_extff_t flags) const [inline]

Note

Directly based on C function rsb_mtx_get_vec() from <rsb.h>.

template<RSBP_Scalar_t NT> rsb_nnz_idx_t RsbMatrix< NT >::nnz (void) const [inline]

template<RSBP_Scalar_t NT> NT RsbMatrix< NT >::normInf (void) const [inline]

template<RSBP_Scalar_t NT> NT RsbMatrix< NT >::normOne (void) const [inline]

template<RSBP_Scalar_t NT> bool RsbMatrix< NT >::operator!= (const RsbMatrix< NT > & B_Rsb) const [inline]

Example snip from examples/misc.cpp:



  assert(   mtx1 == mtx2  );


  assert( !(mtx1 != mtx2) );

See also

RsbMatrix::operator==(const RsbMatrix & B_Rsb) const;

template<RSBP_Scalar_t NT> RsbMatrix & RsbMatrix< NT >::operator= (const RsbMatrix< NT > & A_Rsb) [inline]

A copy constructor. Will clone the input matrix contents.

template<RSBP_Scalar_t NT> bool RsbMatrix< NT >::operator== (const RsbMatrix< NT > & B_Rsb) const [inline]

Deep comparison: compare if the two matrices have same dimensions, nonzeroes count, nonzeroes pattern and value. Meant for very sporadic use. Inefficient: it can involve matrices copying.

Example snip from examples/misc.cpp:



  assert(   mtx1 == mtx2  );


  assert( !(mtx1 != mtx2) );

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::rndr (const rsb_char_t * filename = RSBP_NULL, rsb_coo_idx_t pmWidth = 512, rsb_coo_idx_t pmHeight = 512, rsb_marf_t rflags = RSB_MARF_EPS) const [inline]

Note

Directly based on C function rsb_mtx_rndr() from <rsb.h>.

Example snip from examples/render.cpp:



                mtx.rndr(psfilename.c_str());

template<RSBP_Scalar_t NT> rsb_coo_idx_t RsbMatrix< NT >::rows (void) const [inline]

Example snip from examples/twonnz.cpp:



        std::cout << "# Matrix sized " << mtx.rows() << "x" << mtx.cols() << ", " << nnzA << " nnz  built in " << dt << " s and occupies " << mtxocc << " bytes " << std::endl;

template<RSBP_Scalar_t NT> rsb_flags_t RsbMatrix< NT >::rsbflags (void) const [inline]

template<RSBP_Scalar_t NT> rsb_type_t RsbMatrix< NT >::rsbtype (void) const [inline]

Example snip from examples/twonnz.cpp:



                std::cout << "# type=" << mtx.rsbtype() << " nt=1," << rnt << " n=" << n << " nrhs=" << nrhs << " order=" << oc << " alpha=" << alpha << " beta=" << beta << " dt=" << dta[0] << ".." << dta[1] << " spmm-scalability=" << dta[0]/dta[1] << " nnz/s=" << nnzA/dta[0] << ".." << nnzA/dta[1] << " flops=" << flops_c/dta[0] << ".." << flops_c/dta[1] << " occ.=" << opocc << " " << std::endl;

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::set_val (const NT val, const rsb_coo_idx_t i, const rsb_coo_idx_t j, rsb_flags_t flags = RSB_FLAG_NOFLAGS) [inline]

Note

Directly based on C function rsb_mtx_set_vals() from <rsb.h>.

Example snip from examples/assemble.cpp:



  const rsb_coo_idx_t nrA { 4 }, ncA { 4 };


  RsbMatrix<double> mtx(nrA,ncA); // begin matrix assembly


  // insert elements of a tridiagonal matrix, one by one


  for (auto i = 0; i < nrA; ++i )


    for (auto j = i-1; j <= i+1; ++j )


      if ( i >= 0 && i < nrA )


        if ( j >= 0 && j < ncA )


          mtx.set_val((i+1)*100+(j+1),i,j); // add entry


  mtx.close(); // finish matrix assembly


  assert(mtx.nnz() == 3 * nrA - 2);

See also

set_vals().

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::set_vals (const NT * VA, const rsb_coo_idx_t * IA, const rsb_coo_idx_t * JA, rsb_nnz_idx_t nnz, rsb_flags_t flags) [inline]

Add a single entry during the assembly of a matrix created empty.
Use close() to terminate matrix assembly.

Example snip from examples/assemble.cpp:



  const rsb_coo_idx_t nrA { 4 }, ncA { 4 };


  RsbMatrix<double> mtx(nrA,ncA); // begin matrix assembly


  // insert elements of a tridiagonal matrix, one by one


  for (auto i = 0; i < nrA; ++i )


    for (auto j = i-1; j <= i+1; ++j )


      if ( i >= 0 && i < nrA )


        if ( j >= 0 && j < ncA )


          mtx.set_val((i+1)*100+(j+1),i,j); // add entry


  mtx.close(); // finish matrix assembly


  assert(mtx.nnz() == 3 * nrA - 2);

Note

Directly based on C function rsb_mtx_set_vals() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::spmm (rsb_trans_t transA, const NT * alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, rsb_nnz_idx_t ldB, const NT * betap, NT * Cp, rsb_nnz_idx_t ldC) const [inline]

Note

Directly based on C function rsb_spmm() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spmm (rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, const NT beta, NT * Cp) const [inline]

Note

Indirectly based on C function rsb_spmm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spmm (rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, rsb_nnz_idx_t ldB, const NT beta, NT * Cp, rsb_nnz_idx_t ldC) const [inline]

Example snip from examples/bench.cpp:



                        mtx.spmm(transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC);

Note

Indirectly based on C function rsb_spmm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spmv (NT * y, const NT * x, bool do_trans = false) const [inline]

Note

Indirectly based on C function rsb_spmv() from <rsb.h>.

Example snip from examples/span.cpp:



  mtx.tune_spmm(nullptr,&tn,0,0.0,RSB_TRANSPOSITION_N,alpha,nrhs,RSB_FLAG_WANT_COLUMN_MAJOR_ORDER,X,ncA,beta,Y,nrA);


  mtx.spmv(RSB_TRANSPOSITION_N, alpha, X, beta, Y);

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::spmv (rsb_trans_t transA, const NT * alphap, const NT * Xp, rsb_coo_idx_t incX, const NT * betap, NT * Yp, rsb_coo_idx_t incY) const [inline]

Note

Directly based on C function rsb_spmv() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spmv (rsb_trans_t transA, const NT alpha, const NT * Xp, const NT beta, NT * Yp) const [inline]

Note

Indirectly based on C function rsb_spmv() from <rsb.h>.

Example snip from examples/span.cpp:



  mtx.tune_spmm(nullptr,&tn,0,0.0,RSB_TRANSPOSITION_N,alpha,nrhs,RSB_FLAG_WANT_COLUMN_MAJOR_ORDER,X,ncA,beta,Y,nrA);


  mtx.spmv(RSB_TRANSPOSITION_N, alpha, X, beta, Y);

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spmv (rsb_trans_t transA, const NT alpha, const NT * Xp, rsb_coo_idx_t incX, const NT beta, NT * Yp, rsb_coo_idx_t incY) const [inline]

Note

Indirectly based on C function rsb_spmv() from <rsb.h>.

Example snip from examples/span.cpp:



  mtx.tune_spmm(nullptr,&tn,0,0.0,RSB_TRANSPOSITION_N,alpha,nrhs,RSB_FLAG_WANT_COLUMN_MAJOR_ORDER,X,ncA,beta,Y,nrA);


  mtx.spmv(RSB_TRANSPOSITION_N, alpha, X, beta, Y);

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsm (NT * y, const NT * x, rsb_coo_idx_t nrhs, bool do_trans = false) const [inline]

Note

Indirectly based on C function rsb_spsm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsm (NT * y, rsb_coo_idx_t nrhs, bool do_trans = false) const [inline]

Note

Indirectly based on C function rsb_spsm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::spsm (rsb_trans_t transT, const NT * alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * betap, const NT * Bp, rsb_nnz_idx_t ldB, NT * Cp, rsb_nnz_idx_t ldC) const [inline]

Note

Directly based on C function rsb_spsm() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsm (rsb_trans_t transT, const NT alpha, rsb_coo_idx_t nrhs, const NT * Bp, NT * Cp) const [inline]

Note

Indirectly based on C function rsb_spsm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsm (rsb_trans_t transT, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT beta, const NT * Bp, rsb_nnz_idx_t ldB, NT * Cp, rsb_nnz_idx_t ldC) const [inline]

Note

Indirectly based on C function rsb_spsm() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsv (NT * y, bool do_trans = false) const [inline]

Note

Indirectly based on C function rsb_spsv() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsv (NT * y, const NT * x, bool do_trans = false) const [inline]

Note

Indirectly based on C function rsb_spsv() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::spsv (rsb_trans_t transT, const NT * alphap, const NT * Xp, rsb_coo_idx_t incX, NT * Yp, rsb_coo_idx_t incY) const [inline]

Note

Directly based on C function rsb_spsv() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::spsv (rsb_trans_t transT, const NT alpha, const NT * Xp, NT * Yp) const [inline]

Note

Indirectly based on C function rsb_spsv() from <rsb.h>.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::tune_spmm (rsb_real_t & sf, rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, const NT beta, NT * Cp) [inline]

Note

Indirectly based on C function rsb_tune_spmm() from <rsb.h>.

Example snip from examples/autotune.cpp:



                        tt = -rsb_time();


                        mtx.tune_spmm(&sf,&tn,maxr,tmax,transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC);


                        tt += rsb_time();


                        auto nnsmA {mtx.blocks()};


                        std::cout << "Tuning took " << tt << " s ( " << tt / dt << " ops ) and changed " << nsmA << " to " << nnsmA << " blocks" << std::endl;


                        mtx.spmm(transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC); // caches warmup

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::tune_spmm (rsb_real_t * sfp, rsb_int_t * tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT * alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, rsb_nnz_idx_t ldB, const NT * betap, NT * Cp, rsb_nnz_idx_t ldC) [inline]

Note

Directly based on C function rsb_tune_spmm() from <rsb.h>.

Example snip from examples/autotune.cpp:



                        tt = -rsb_time();


                        mtx.tune_spmm(&sf,&tn,maxr,tmax,transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC);


                        tt += rsb_time();


                        auto nnsmA {mtx.blocks()};


                        std::cout << "Tuning took " << tt << " s ( " << tt / dt << " ops ) and changed " << nsmA << " to " << nnsmA << " blocks" << std::endl;


                        mtx.spmm(transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC); // caches warmup

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::tune_spmm (rsb_real_t * sfp, rsb_int_t * tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT alpha, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, rsb_nnz_idx_t ldB, const NT beta, NT * Cp, rsb_nnz_idx_t ldC) [inline]

Note

Indirectly based on C function rsb_tune_spmm() from <rsb.h>.

Example snip from examples/autotune.cpp:



                        tt = -rsb_time();


                        mtx.tune_spmm(&sf,&tn,maxr,tmax,transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC);


                        tt += rsb_time();


                        auto nnsmA {mtx.blocks()};


                        std::cout << "Tuning took " << tt << " s ( " << tt / dt << " ops ) and changed " << nsmA << " to " << nnsmA << " blocks" << std::endl;


                        mtx.spmm(transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC); // caches warmup

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::tune_spmm_threads (rsb_real_t * sfp = RSBP_NULL, rsb_int_t * tnp = RSBP_NULL, rsb_int_t maxr = 0, rsb_time_t maxt = 0, rsb_trans_t transA = RSB_TRANSPOSITION_N, const NT * alphap = RSBP_NULL, rsb_coo_idx_t nrhs = 1, rsb_flags_t order = RSB_FLAG_WANT_COLUMN_MAJOR_ORDER, const NT * Bp = RSBP_NULL, rsb_nnz_idx_t ldB = 0, const NT * betap = RSBP_NULL, NT * Cp = RSBP_NULL, rsb_nnz_idx_t ldC = 0) const [inline]

Note

Directly based on C function rsb_tune_spmm() from <rsb.h>.

Example snip from examples/autotune.cpp:



                        tt = -rsb_time();


                        mtx.tune_spmm(&sf,&tn,maxr,tmax,transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC);


                        tt += rsb_time();


                        auto nnsmA {mtx.blocks()};


                        std::cout << "Tuning took " << tt << " s ( " << tt / dt << " ops ) and changed " << nsmA << " to " << nnsmA << " blocks" << std::endl;


                        mtx.spmm(transA,&alpha,nrhs,order,B.data(),ldB,&beta,C.data(),ldC); // caches warmup

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::tune_spsm (rsb_real_t * sfp, rsb_int_t * tnp, rsb_int_t maxr, rsb_time_t maxt, rsb_trans_t transA, const NT * alphap, rsb_coo_idx_t nrhs, rsb_flags_t order, const NT * Bp, rsb_nnz_idx_t ldB, const NT * betap, NT * Cp, rsb_nnz_idx_t ldC) [inline]

Note

Directly based on C function rsb_tune_spsm() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT RSBP_DEPRECATED Err_t RsbMatrix< NT >::tune_spsm_threads (rsb_real_t * sfp = RSBP_NULL, rsb_int_t * tnp = RSBP_NULL, rsb_int_t maxr = 0, rsb_time_t maxt = 0, rsb_trans_t transA = RSB_TRANSPOSITION_N, const NT * alphap = RSBP_NULL, rsb_coo_idx_t nrhs = 1, rsb_flags_t order = RSB_FLAG_WANT_COLUMN_MAJOR_ORDER, const NT * Bp = RSBP_NULL, rsb_nnz_idx_t ldB = 0, const NT * betap = RSBP_NULL, NT * Cp = RSBP_NULL, rsb_nnz_idx_t ldC = 0) const [inline]

Note

Directly based on C function rsb_tune_spsm() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::upd_vals (enum rsb_elopf_t elop_flags, const NT & omega) [inline]

Note

Directly based on C function rsb_mtx_upd_vals() from <rsb.h>.

Example snip from examples/bench.cpp:



        mtx.upd_vals(RSB_ELOPF_POW,nt_t{0.0}); // set matrix values to ones

template<RSBP_Scalar_t NT> RSBP_RVT Err_t RsbMatrix< NT >::upd_vals (enum rsb_elopf_t elop_flags, const NT * omegap) [inline]

Note

Directly based on C function rsb_mtx_upd_vals() from <rsb.h>.

Deprecated

This function overload is deprecated: use the corresponding overload without the C-style pointers-to-scalar, or request one if non existing yet.

Author

Generated automatically by Doxygen for librsb from the source code.

Version 1.3.0.2 librsb