LAPACK
3.5.0
LAPACK: Linear Algebra PACKage
Main Page
Data Types List
Files
File List
File Members
All
Classes
Files
Functions
Variables
Typedefs
Macros
lapacke_dtpmqrt_work.c
Go to the documentation of this file.
1
/*****************************************************************************
2
Copyright (c) 2010, Intel Corp.
3
All rights reserved.
4
5
Redistribution and use in source and binary forms, with or without
6
modification, are permitted provided that the following conditions are met:
7
8
* Redistributions of source code must retain the above copyright notice,
9
this list of conditions and the following disclaimer.
10
* Redistributions in binary form must reproduce the above copyright
11
notice, this list of conditions and the following disclaimer in the
12
documentation and/or other materials provided with the distribution.
13
* Neither the name of Intel Corporation nor the names of its contributors
14
may be used to endorse or promote products derived from this software
15
without specific prior written permission.
16
17
THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
18
AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
19
IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
20
ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
21
LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
22
CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
23
SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
24
INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
25
CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
26
ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF
27
THE POSSIBILITY OF SUCH DAMAGE.
28
******************************************************************************
29
* Contents: Native middle-level C interface to LAPACK function dtpmqrt
30
* Author: Intel Corporation
31
* Generated November, 2011
32
*****************************************************************************/
33
34
#include "
lapacke_utils.h
"
35
36
lapack_int
LAPACKE_dtpmqrt_work
(
int
matrix_order,
char
side,
char
trans,
37
lapack_int
m,
lapack_int
n,
lapack_int
k,
38
lapack_int
l,
lapack_int
nb,
const
double
* v,
39
lapack_int
ldv,
const
double
* t,
40
lapack_int
ldt,
double
* a,
lapack_int
lda,
41
double
*
b
,
lapack_int
ldb,
double
* work )
42
{
43
lapack_int
info = 0;
44
if
( matrix_order ==
LAPACK_COL_MAJOR
) {
45
/* Call LAPACK function and adjust info */
46
LAPACK_dtpmqrt
( &side, &trans, &m, &n, &k, &l, &nb, v, &ldv, t, &ldt, a,
47
&lda, b, &ldb, work, &info );
48
if
( info < 0 ) {
49
info = info - 1;
50
}
51
}
else
if
( matrix_order ==
LAPACK_ROW_MAJOR
) {
52
lapack_int
r =
LAPACKE_lsame
( side,
'r'
) ? k : k;
53
lapack_int
lda_t =
MAX
(1,k);
54
lapack_int
ldb_t =
MAX
(1,m);
55
lapack_int
ldt_t =
MAX
(1,ldt);
56
lapack_int
ldv_t =
MAX
(1,ldv);
57
double
* v_t = NULL;
58
double
* t_t = NULL;
59
double
* a_t = NULL;
60
double
* b_t = NULL;
61
/* Check leading dimension(s) */
62
if
( lda < m ) {
63
info = -14;
64
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
65
return
info;
66
}
67
if
( ldb < n ) {
68
info = -16;
69
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
70
return
info;
71
}
72
if
( ldt < nb ) {
73
info = -12;
74
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
75
return
info;
76
}
77
if
( ldv < k ) {
78
info = -10;
79
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
80
return
info;
81
}
82
/* Allocate memory for temporary array(s) */
83
v_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * ldv_t *
MAX
(1,k) );
84
if
( v_t == NULL ) {
85
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
86
goto
exit_level_0;
87
}
88
t_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * ldt_t *
MAX
(1,nb) );
89
if
( t_t == NULL ) {
90
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
91
goto
exit_level_1;
92
}
93
a_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * lda_t *
MAX
(1,m) );
94
if
( a_t == NULL ) {
95
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
96
goto
exit_level_2;
97
}
98
b_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * ldb_t *
MAX
(1,n) );
99
if
( b_t == NULL ) {
100
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
101
goto
exit_level_3;
102
}
103
/* Transpose input matrices */
104
LAPACKE_dge_trans
( matrix_order, ldv, k, v, ldv, v_t, ldv_t );
105
LAPACKE_dge_trans
( matrix_order, ldt, nb, t, ldt, t_t, ldt_t );
106
LAPACKE_dge_trans
( matrix_order, k, m, a, lda, a_t, lda_t );
107
LAPACKE_dge_trans
( matrix_order, m, n, b, ldb, b_t, ldb_t );
108
/* Call LAPACK function and adjust info */
109
LAPACK_dtpmqrt
( &side, &trans, &m, &n, &k, &l, &nb, v_t, &ldv_t, t_t,
110
&ldt_t, a_t, &lda_t, b_t, &ldb_t, work, &info );
111
if
( info < 0 ) {
112
info = info - 1;
113
}
114
/* Transpose output matrices */
115
LAPACKE_dge_trans
(
LAPACK_COL_MAJOR
, k, m, a_t, lda_t, a, lda );
116
LAPACKE_dge_trans
(
LAPACK_COL_MAJOR
, m, n, b_t, ldb_t, b, ldb );
117
/* Release memory and exit */
118
LAPACKE_free
( b_t );
119
exit_level_3:
120
LAPACKE_free
( a_t );
121
exit_level_2:
122
LAPACKE_free
( t_t );
123
exit_level_1:
124
LAPACKE_free
( v_t );
125
exit_level_0:
126
if
( info ==
LAPACK_TRANSPOSE_MEMORY_ERROR
) {
127
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
128
}
129
}
else
{
130
info = -1;
131
LAPACKE_xerbla
(
"LAPACKE_dtpmqrt_work"
, info );
132
}
133
return
info;
134
}
lapacke
src
lapacke_dtpmqrt_work.c
Generated on Mon Dec 30 2013 16:09:25 for LAPACK by
1.8.1.2