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_dtgevc_work.c
Go to the documentation of this file.
1
/*****************************************************************************
2
Copyright (c) 2011, 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 dtgevc
30
* Author: Intel Corporation
31
* Generated November, 2011
32
*****************************************************************************/
33
34
#include "
lapacke_utils.h
"
35
36
lapack_int
LAPACKE_dtgevc_work
(
int
matrix_order,
char
side,
char
howmny,
37
const
lapack_logical
* select,
lapack_int
n,
38
const
double
* s,
lapack_int
lds,
39
const
double
* p,
lapack_int
ldp,
double
* vl,
40
lapack_int
ldvl,
double
* vr,
lapack_int
ldvr,
41
lapack_int
mm,
lapack_int
* m,
double
* work )
42
{
43
lapack_int
info = 0;
44
if
( matrix_order ==
LAPACK_COL_MAJOR
) {
45
/* Call LAPACK function and adjust info */
46
LAPACK_dtgevc
( &side, &howmny, select, &n, s, &lds, p, &ldp, vl, &ldvl,
47
vr, &ldvr, &mm, m, work, &info );
48
if
( info < 0 ) {
49
info = info - 1;
50
}
51
}
else
if
( matrix_order ==
LAPACK_ROW_MAJOR
) {
52
lapack_int
ldp_t =
MAX
(1,n);
53
lapack_int
lds_t =
MAX
(1,n);
54
lapack_int
ldvl_t =
MAX
(1,n);
55
lapack_int
ldvr_t =
MAX
(1,n);
56
double
* s_t = NULL;
57
double
* p_t = NULL;
58
double
* vl_t = NULL;
59
double
* vr_t = NULL;
60
/* Check leading dimension(s) */
61
if
( ldp < n ) {
62
info = -9;
63
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
64
return
info;
65
}
66
if
( lds < n ) {
67
info = -7;
68
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
69
return
info;
70
}
71
if
( ldvl < mm ) {
72
info = -11;
73
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
74
return
info;
75
}
76
if
( ldvr < mm ) {
77
info = -13;
78
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
79
return
info;
80
}
81
/* Allocate memory for temporary array(s) */
82
s_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * lds_t *
MAX
(1,n) );
83
if
( s_t == NULL ) {
84
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
85
goto
exit_level_0;
86
}
87
p_t = (
double
*)
LAPACKE_malloc
(
sizeof
(
double
) * ldp_t *
MAX
(1,n) );
88
if
( p_t == NULL ) {
89
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
90
goto
exit_level_1;
91
}
92
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'l'
) ) {
93
vl_t = (
double
*)
94
LAPACKE_malloc
(
sizeof
(
double
) * ldvl_t *
MAX
(1,mm) );
95
if
( vl_t == NULL ) {
96
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
97
goto
exit_level_2;
98
}
99
}
100
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'r'
) ) {
101
vr_t = (
double
*)
102
LAPACKE_malloc
(
sizeof
(
double
) * ldvr_t *
MAX
(1,mm) );
103
if
( vr_t == NULL ) {
104
info =
LAPACK_TRANSPOSE_MEMORY_ERROR
;
105
goto
exit_level_3;
106
}
107
}
108
/* Transpose input matrices */
109
LAPACKE_dge_trans
( matrix_order, n, n, s, lds, s_t, lds_t );
110
LAPACKE_dge_trans
( matrix_order, n, n, p, ldp, p_t, ldp_t );
111
if
( (
LAPACKE_lsame
( side,
'l'
) ||
LAPACKE_lsame
( side,
'b'
) ) &&
112
LAPACKE_lsame
( howmny,
'b'
) ) {
113
LAPACKE_dge_trans
( matrix_order, n, mm, vl, ldvl, vl_t, ldvl_t );
114
}
115
if
( (
LAPACKE_lsame
( side,
'r'
) ||
LAPACKE_lsame
( side,
'b'
) ) &&
116
LAPACKE_lsame
( howmny,
'b'
) ) {
117
LAPACKE_dge_trans
( matrix_order, n, mm, vr, ldvr, vr_t, ldvr_t );
118
}
119
/* Call LAPACK function and adjust info */
120
LAPACK_dtgevc
( &side, &howmny, select, &n, s_t, &lds_t, p_t, &ldp_t,
121
vl_t, &ldvl_t, vr_t, &ldvr_t, &mm, m, work, &info );
122
if
( info < 0 ) {
123
info = info - 1;
124
}
125
/* Transpose output matrices */
126
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'l'
) ) {
127
LAPACKE_dge_trans
(
LAPACK_COL_MAJOR
, n, mm, vl_t, ldvl_t, vl,
128
ldvl );
129
}
130
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'r'
) ) {
131
LAPACKE_dge_trans
(
LAPACK_COL_MAJOR
, n, mm, vr_t, ldvr_t, vr,
132
ldvr );
133
}
134
/* Release memory and exit */
135
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'r'
) ) {
136
LAPACKE_free
( vr_t );
137
}
138
exit_level_3:
139
if
(
LAPACKE_lsame
( side,
'b'
) ||
LAPACKE_lsame
( side,
'l'
) ) {
140
LAPACKE_free
( vl_t );
141
}
142
exit_level_2:
143
LAPACKE_free
( p_t );
144
exit_level_1:
145
LAPACKE_free
( s_t );
146
exit_level_0:
147
if
( info ==
LAPACK_TRANSPOSE_MEMORY_ERROR
) {
148
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
149
}
150
}
else
{
151
info = -1;
152
LAPACKE_xerbla
(
"LAPACKE_dtgevc_work"
, info );
153
}
154
return
info;
155
}
lapacke
src
lapacke_dtgevc_work.c
Generated on Mon Dec 30 2013 16:09:25 for LAPACK by
1.8.1.2