spblas.hpp

Go to the documentation of this file.

/*
 * $Id$
 * Copyright (C) 2010, 2018, 2019, 2020 John D Lamb
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or (at
 * your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 */
 
#ifndef CCGSL_SPBLAS
#define CCGSL_SPBLAS
 
#include<cmath>
#include<gsl/gsl_spblas.h>
#include"vector.hpp"
#include"spmatrix.hpp"
#include"blas.hpp"
 
namespace gsl {
  namespace spblas {
    inline int dgemv( CBLAS_TRANSPOSE_t const TransA, double const alpha, spmatrix const& A,
              vector const& x, double const beta, vector& y ){
      return gsl_spblas_dgemv( TransA, alpha, A.get(), x.get(), beta, y.get() ); } 
    inline int dgemm( double const alpha, spmatrix const& A, spmatrix const& B, spmatrix& C ){
      return gsl_spblas_dgemm( alpha, A.get(), B.get(), C.get() ); } 
#ifndef DOXYGEN_SKIP
    inline size_t scatter( spmatrix const& A, size_t const j, double const alpha, int* w,
               double* x, int const mark, spmatrix& C, size_t nz ){
      return gsl_spblas_scatter( A.get(), j, alpha, w, x, mark, C.get(), nz ); }
#endif // DOXYGEN_SKIP
  }
}
 
#endif