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
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