#ifndef HPC_ULMBLAS_SCAL_H #define HPC_ULMBLAS_SCAL_H 1 namespace hpc { namespace ulmblas { template <typename Index, typename Alpha, typename TX> void scal(Index n, const Alpha &alpha, TX *x, Index incX) { if (alpha!=Alpha(1)) { for (Index i=0; i<n; ++i) { x[i*incX] *= alpha; } } } } } // namespace ulmblas, hpc #endif // HPC_ULMBLAS_SCAL_H |