1 /*
  2  *   Copyright (c) 2007, 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 FLENS_BLAS_CLOSURES_ASSIGN_TCC
 34 #define FLENS_BLAS_CLOSURES_ASSIGN_TCC 1
 35 
 36 #ifdef FLENS_DEBUG_CLOSURES
 37 #   include <flens/blas/blaslogon.h>
 38 #else
 39 #   include <flens/blas/blaslogoff.h>
 40 #endif
 41 
 42 namespace flens {
 43 
 44 //-- vector closures
 45 template <typename VX, typename VY>
 46 void
 47 assign(const Vector<VX> &x, Vector<VY> &y)
 48 {
 49     FLENS_BLASLOG_BEGIN_ASSIGNMENT(x, y);
 50 
 51     blas::copy(x.impl(), y.impl());
 52 
 53     FLENS_BLASLOG_END;
 54 }
 55 
 56 template <typename VX, typename VY>
 57 void
 58 plusAssign(const Vector<VX> &x, Vector<VY> &y)
 59 {
 60     FLENS_BLASLOG_BEGIN_PLUSASSIGNMENT(x, y);
 61 
 62     typedef typename Vector<VX>::Impl::ElementType  T;
 63     const T  One(1);
 64 
 65     blas::axpy(One, x.impl(), y.impl());
 66 
 67     FLENS_BLASLOG_END;
 68 }
 69 
 70 template <typename VX, typename VY>
 71 void
 72 minusAssign(const Vector<VX> &x, Vector<VY> &y)
 73 {
 74     FLENS_BLASLOG_BEGIN_MINUSASSIGNMENT(x, y);
 75 
 76     typedef typename Vector<VX>::Impl::ElementType  T;
 77     const T  MinusOne(-1);
 78 
 79     blas::axpy(MinusOne, x.impl(), y.impl());
 80 
 81     FLENS_BLASLOG_END;
 82 }
 83 
 84 //-- matrix closures
 85 template <typename MA, typename MB>
 86 void
 87 assign(const Matrix<MA> &A, Matrix<MB> &B)
 88 {
 89     FLENS_BLASLOG_BEGIN_ASSIGNMENT(A, B);
 90 
 91     blas::copy(NoTrans, A.impl(), B.impl());
 92 
 93     FLENS_BLASLOG_END;
 94 }
 95 
 96 template <typename MA, typename MB>
 97 void
 98 plusAssign(const Matrix<MA> &A, Matrix<MB> &B)
 99 {
100     FLENS_BLASLOG_BEGIN_PLUSASSIGNMENT(A, B);
101 
102     typedef typename Matrix<MA>::Impl::ElementType  T;
103     const T  One(1);
104 
105     blas::axpy(NoTrans, One, A.impl(), B.impl());
106 
107     FLENS_BLASLOG_END;
108 }
109 
110 template <typename MA, typename MB>
111 void
112 minusAssign(const Matrix<MA> &A, Matrix<MB> &B)
113 {
114     FLENS_BLASLOG_BEGIN_MINUSASSIGNMENT(A, B);
115 
116     typedef typename Matrix<MA>::Impl::ElementType  T;
117     const T  MinusOne(-1);
118 
119     blas::axpy(NoTrans, MinusOne, A.impl(), B.impl());
120 
121     FLENS_BLASLOG_END;
122 }
123 
124 // namespace flens
125 
126 #endif // FLENS_BLAS_CLOSURES_ASSIGN_TCC