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