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 FLENS_BLAS_LEVEL2_MV_H
34 #define FLENS_BLAS_LEVEL2_MV_H 1
35
36 #include <cxxblas/cxxblas.h>
37 #include <flens/matrixtypes/matrixtypes.h>
38 #include <flens/typedefs.h>
39 #include <flens/vectortypes/vectortypes.h>
40
41 namespace flens { namespace blas {
42
43 //== GeneralMatrix - Vector products ===========================================
44
45 //-- gemv
46 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
47 void
48 mv(Transpose trans,
49 const ALPHA &alpha, const GeMatrix<MA> &A, const DenseVector<VX> &x,
50 const BETA &beta, DenseVector<VY> &y);
51
52
53 //== TriangularMatrix - Vector products ========================================
54
55 //-- trmv
56 template <typename MA, typename VX>
57 void
58 mv(Transpose trans, const TrMatrix<MA> &A, DenseVector<VX> &x);
59
60
61 //== SymmetricMatrix - Vector products =========================================
62
63 //-- symv
64 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
65 void
66 mv(const ALPHA &alpha, const SyMatrix<MA> &A, const DenseVector<VX> &x,
67 const BETA &beta, DenseVector<VY> &y);
68
69
70 //== HermitianMatrix - Vector products =========================================
71
72 //-- hemv
73 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
74 void
75 mv(const ALPHA &alpha, const HeMatrix<MA> &A, const DenseVector<VX> &x,
76 const BETA &beta, DenseVector<VY> &y);
77
78
79 //== forwarding ================================================================
80
81 //-- GeneralMatrix - Vector products
82 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
83 typename RestrictTo<IsGeneralMatrix<MA>::value &&
84 !IsClosure<MA>::value && !IsClosure<VX>::value &&
85 IsSame<VY, typename VY::Impl>::value,
86 void>::Type
87 mv(Transpose trans,
88 const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
89
90 //-- TriangularMatrix - Vector products
91 template <typename MA, typename VX>
92 typename RestrictTo<IsTriangularMatrix<MA>::value &&
93 !IsClosure<MA>::value &&
94 IsSame<VX, typename VX::Impl>::value,
95 void>::Type
96 mv(Transpose trans, const MA &A, VX &&x);
97
98 //-- Symmetric Matrix - Vector products
99 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
100 typename RestrictTo<IsSymmetricMatrix<MA>::value &&
101 !IsClosure<MA>::value && !IsClosure<VX>::value &&
102 IsSame<VY, typename VY::Impl>::value,
103 void>::Type
104 mv(const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
105
106 //-- Hermitian Matrix - Vector products
107 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
108 typename RestrictTo<IsHermitianMatrix<MA>::value &&
109 !IsClosure<MA>::value && !IsClosure<VX>::value &&
110 IsSame<VY, typename VY::Impl>::value,
111 void>::Type
112 mv(const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
113
114
115
116 } } // namespace blas, flens
117
118 #endif // FLENS_BLAS_LEVEL2_MV_H
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 FLENS_BLAS_LEVEL2_MV_H
34 #define FLENS_BLAS_LEVEL2_MV_H 1
35
36 #include <cxxblas/cxxblas.h>
37 #include <flens/matrixtypes/matrixtypes.h>
38 #include <flens/typedefs.h>
39 #include <flens/vectortypes/vectortypes.h>
40
41 namespace flens { namespace blas {
42
43 //== GeneralMatrix - Vector products ===========================================
44
45 //-- gemv
46 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
47 void
48 mv(Transpose trans,
49 const ALPHA &alpha, const GeMatrix<MA> &A, const DenseVector<VX> &x,
50 const BETA &beta, DenseVector<VY> &y);
51
52
53 //== TriangularMatrix - Vector products ========================================
54
55 //-- trmv
56 template <typename MA, typename VX>
57 void
58 mv(Transpose trans, const TrMatrix<MA> &A, DenseVector<VX> &x);
59
60
61 //== SymmetricMatrix - Vector products =========================================
62
63 //-- symv
64 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
65 void
66 mv(const ALPHA &alpha, const SyMatrix<MA> &A, const DenseVector<VX> &x,
67 const BETA &beta, DenseVector<VY> &y);
68
69
70 //== HermitianMatrix - Vector products =========================================
71
72 //-- hemv
73 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
74 void
75 mv(const ALPHA &alpha, const HeMatrix<MA> &A, const DenseVector<VX> &x,
76 const BETA &beta, DenseVector<VY> &y);
77
78
79 //== forwarding ================================================================
80
81 //-- GeneralMatrix - Vector products
82 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
83 typename RestrictTo<IsGeneralMatrix<MA>::value &&
84 !IsClosure<MA>::value && !IsClosure<VX>::value &&
85 IsSame<VY, typename VY::Impl>::value,
86 void>::Type
87 mv(Transpose trans,
88 const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
89
90 //-- TriangularMatrix - Vector products
91 template <typename MA, typename VX>
92 typename RestrictTo<IsTriangularMatrix<MA>::value &&
93 !IsClosure<MA>::value &&
94 IsSame<VX, typename VX::Impl>::value,
95 void>::Type
96 mv(Transpose trans, const MA &A, VX &&x);
97
98 //-- Symmetric Matrix - Vector products
99 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
100 typename RestrictTo<IsSymmetricMatrix<MA>::value &&
101 !IsClosure<MA>::value && !IsClosure<VX>::value &&
102 IsSame<VY, typename VY::Impl>::value,
103 void>::Type
104 mv(const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
105
106 //-- Hermitian Matrix - Vector products
107 template <typename ALPHA, typename MA, typename VX, typename BETA, typename VY>
108 typename RestrictTo<IsHermitianMatrix<MA>::value &&
109 !IsClosure<MA>::value && !IsClosure<VX>::value &&
110 IsSame<VY, typename VY::Impl>::value,
111 void>::Type
112 mv(const ALPHA &alpha, const MA &A, const VX &x, const BETA &beta, VY &&y);
113
114
115
116 } } // namespace blas, flens
117
118 #endif // FLENS_BLAS_LEVEL2_MV_H