src/dense_cpu.c
void sgemv_wrapper ( int N, data_struct* DATA, float* x, float* y )
{
    checkpoint t0 = tic();
    // read the relevant data from the struct
    float* A = DATA->A;
    int LDA = DATA->LDA;
    // calculate y <-- Ax
    cblas_sgemv ( CblasColMajor, CblasNoTrans, N, N, 1.0, A, LDA, x, 1, 0.0, y, 1);
    DATA->sgemv_time += toc ( t0 );
    DATA->sgemv_calls += 1;
}
//!
//! Calculate the matrix-vector product  y <- Ax for the eigensystem defined by DATA. This
//! routine is called from the Fortran backend whenever ARPACK requires a matrix-vector
//! product to be calculated.
//! \param[in] N The dimension of the eigensystem.
//! \param[in,out] DATA The structure representing the eigenproblem.
//! \param[in] x The vector that must be multiplied by A.
//! \param[out] y The vector that must store the result.
Calculate the matrix-vector product  y <- Ax for the eigensystem defined by DATA. This
routine is called from the Fortran backend whenever ARPACK requires a matrix-vector
product to be calculated.
\param[in] N The dimension of the eigensystem.
\param[in,out] DATA The structure representing the eigenproblem.
\param[in] x The vector that must be multiplied by A.
\param[out] y The vector that must store the result.
Calculate the matrix - vector product y < - Ax for the eigensystem defined by DATA . This routine is called from the Fortran backend whenever ARPACK requires a matrix - vector product to be calculated . \ param [ in ] N The dimension of the eigensystem . \ param [ in out ] DATA The structure representing the eigenproblem . \ param [ in ] x The vector that must be multiplied by A . \ param [ out ] y The vector that must store the result .
void sgemv_wrapper ( int N, data_struct* DATA, float* x, float* y )
{
    checkpoint t0 = tic();
    float* A = DATA->A;
    int LDA = DATA->LDA;
    cblas_sgemv ( CblasColMajor, CblasNoTrans, N, N, 1.0, A, LDA, x, 1, 0.0, y, 1);
    DATA->sgemv_time += toc ( t0 );
    DATA->sgemv_calls += 1;
}
void sgemv_wrapper ( int N , data_struct * DATA , float * x , float * y ) { checkpoint t0 = tic ( ) ; float * A = DATA -> A ; int LDA = DATA -> LDA ; cblas_sgemv ( CblasColMajor , CblasNoTrans , N , N , 1.0 , A , LDA , x , 1 , 0.0 , y , 1 ) ; DATA -> sgemv_time += toc ( t0 ) ; DATA -> sgemv_calls += 1 ; }
Calculate the matrix-vector product  y <- Ax for the eigensystem defined by DATA.
Calculate the matrix - vector product y < - Ax for the eigensystem defined by DATA .
// read the relevant data from the struct // calculate y <-- Ax