File: lapmr.tcc

package info (click to toggle)
xtensor-blas 0.23.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 13,860 kB
  • sloc: cpp: 98,000; makefile: 201; perl: 178; python: 153
file content (124 lines) | stat: -rw-r--r-- 3,918 bytes parent folder | download | duplicates (3)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
/*
 *   Copyright (c) 2012, Michael Lehn, Klaus Pototzky
 *
 *   All rights reserved.
 *
 *   Redistribution and use in source and binary forms, with or without
 *   modification, are permitted provided that the following conditions
 *   are met:
 *
 *   1) Redistributions of source code must retain the above copyright
 *      notice, this list of conditions and the following disclaimer.
 *   2) Redistributions in binary form must reproduce the above copyright
 *      notice, this list of conditions and the following disclaimer in
 *      the documentation and/or other materials provided with the
 *      distribution.
 *   3) Neither the name of the FLENS development group nor the names of
 *      its contributors may be used to endorse or promote products derived
 *      from this software without specific prior written permission.
 *
 *   THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
 *   "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
 *   LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
 *   A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT
 *   OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
 *   SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
 *   LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
 *   DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
 *   THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
 *   (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
 *   OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 */

#ifndef CXXLAPACK_INTERFACE_LAPMR_TCC
#define CXXLAPACK_INTERFACE_LAPMR_TCC 1

#include <iostream>
#include "xflens/cxxlapack/interface/interface.h"
#include "xflens/cxxlapack/netlib/netlib.h"

namespace cxxlapack {

template <typename IndexType>
void
lapmr(bool                  forwrd,
      IndexType             m,
      IndexType             n,
      float                 *X,
      IndexType             ldX,
      IndexType             &k)
{
    CXXLAPACK_DEBUG_OUT("slapmr");

    IndexType forwrd_ = forwrd;
    LAPACK_IMPL(slapmr)(&forwrd_,
                        &m,
                        &n,
                        X,
                        &ldX,
                        &k);
}

template <typename IndexType>
void
lapmr(bool                  forwrd,
      IndexType             m,
      IndexType             n,
      double                *X,
      IndexType             ldX,
      IndexType             &k)
{
    CXXLAPACK_DEBUG_OUT("dlapmr");

    IndexType forwrd_ = forwrd;
    LAPACK_IMPL(dlapmr)(&forwrd_,
                        &m,
                        &n,
                        X,
                        &ldX,
                        &k);
}

template <typename IndexType>
void
lapmr(bool                  forwrd,
      IndexType             m,
      IndexType             n,
      std::complex<float>   *X,
      IndexType             ldX,
      IndexType             &k)
{
    CXXLAPACK_DEBUG_OUT("clapmr");

    IndexType forwrd_ = forwrd;
    LAPACK_IMPL(clapmr)(&forwrd_,
                        &m,
                        &n,
                        reinterpret_cast<float  *>(X),
                        &ldX,
                        &k);
}

template <typename IndexType>
void
lapmr(bool                  forwrd,
      IndexType             m,
      IndexType             n,
      std::complex<double>  *X,
      IndexType             ldX,
      IndexType             &k)
{
    CXXLAPACK_DEBUG_OUT("zlapmr");

    IndexType forwrd_ = forwrd;
    LAPACK_IMPL(zlapmr)(&forwrd_,
                        &m,
                        &n,
                        reinterpret_cast<double *>(X),
                        &ldX,
                        &k);
}

} // namespace cxxlapack

#endif // CXXLAPACK_INTERFACE_LAPMR_TCC