Home
Browse Files
  1 /*
  2  *   Copyright (c) 2009, Michael Lehn
  3  *
  4  *   All rights reserved.
  5  *
  6  *   Redistribution and use in source and binary forms, with or without
  7  *   modification, are permitted provided that the following conditions
  8  *   are met:
  9  *
 10  *   1) Redistributions of source code must retain the above copyright
 11  *      notice, this list of conditions and the following disclaimer.
 12  *   2) Redistributions in binary form must reproduce the above copyright
 13  *      notice, this list of conditions and the following disclaimer in
 14  *      the documentation and/or other materials provided with the
 15  *      distribution.
 16  *   3) Neither the name of the FLENS development group nor the names of
 17  *      its contributors may be used to endorse or promote products derived
 18  *      from this software without specific prior written permission.
 19  *
 20  *   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 21  *   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 22  *   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
 23  *   A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
 24  *   OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 25  *   SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
 26  *   LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
 27  *   DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
 28  *   THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 29  *   (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
 30  *   OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 31  */
 32 
 33 #ifndef CXXBLAS_LEVEL1_SCAL_TCC
 34 #define CXXBLAS_LEVEL1_SCAL_TCC 1
 35 
 36 namespace cxxblas {
 37 
 38 template <typename IndexType, typename ALPHA, typename Y>
 39 void
 40 scal_generic(IndexType n, const ALPHA &alpha, Y *y, IndexType incY)
 41 {
 42     CXXBLAS_DEBUG_OUT("scal_generic");
 43 
 44     for (IndexType i=0, iY=0; i<n; ++i, iY+=incY) {
 45         y[iY] *= alpha;
 46     }
 47 }
 48 
 49 template <typename IndexType, typename ALPHA, typename Y>
 50 void
 51 scal(IndexType n, const ALPHA &alpha, Y *y, IndexType incY)
 52 {
 53     if (incY<0) {
 54         y -= incY*(n-1);
 55     }
 56     scal_generic(n, alpha, y, incY);
 57 }
 58 
 59 #ifdef HAVE_CBLAS
 60 
 61 // sscal
 62 template <typename IndexType>
 63 typename If<IndexType>::isBlasCompatibleInteger
 64 scal(IndexType n, float alpha, float *x, IndexType incX)
 65 {
 66     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_sscal");
 67 
 68     cblas_sscal(n, alpha, x, incX);
 69 }
 70 
 71 // dscal
 72 template <typename IndexType>
 73 typename If<IndexType>::isBlasCompatibleInteger
 74 scal(IndexType n, double alpha, double *x, IndexType incX)
 75 {
 76     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_dscal");
 77 
 78     cblas_dscal(n, alpha, x, incX);
 79 }
 80 
 81 // cscal
 82 template <typename IndexType>
 83 typename If<IndexType>::isBlasCompatibleInteger
 84 scal(IndexType n, const ComplexFloat *alpha, ComplexFloat *x, IndexType incX)
 85 {
 86     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_cscal");
 87 
 88     cblas_cscal(n, reinterpret_cast<const float *>(alpha),
 89                    reinterpret_cast<const float *>(x), incX);
 90 }
 91 
 92 // zscal
 93 template <typename IndexType>
 94 typename If<IndexType>::isBlasCompatibleInteger
 95 scal(IndexType n, const ComplexDouble *alpha, ComplexDouble *x, IndexType incX)
 96 {
 97     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_zscal");
 98 
 99     cblas_zscal(n, reinterpret_cast<const double *>(alpha),
100                    reinterpret_cast<const double *>(x), incX);
101 }
102 
103 // csscal
104 template <typename IndexType>
105 typename If<IndexType>::isBlasCompatibleInteger
106 scal(IndexType n, float alpha, ComplexFloat *x, IndexType incX)
107 {
108     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_csscal");
109 
110     cblas_csscal(n, alpha, reinterpret_cast<const float *>(x), incX);
111 }
112 
113 // zdscal
114 template <typename IndexType>
115 typename If<IndexType>::isBlasCompatibleInteger
116 scal(IndexType n, double alpha, ComplexDouble *x, IndexType incX)
117 {
118     CXXBLAS_DEBUG_OUT("[" BLAS_IMPL "] cblas_zdscal");
119 
120     cblas_zdscal(n, alpha, reinterpret_cast<const double *>(x), incX);
121 }
122 
123 #endif // HAVE_CBLAS
124 
125 } // namespace cxxblas
126 
127 #endif // CXXBLAS_LEVEL1_SCAL_TCC
Home