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