Mercurial > hg > octave-nkf

comparison liboctave/dbleSCHUR.cc @ 1253:bb67a902760b

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1995-04-11 16:35:23 by jwe]
author jwe
date Tue, 11 Apr 1995 16:35:23 +0000
parents 97eac19837dc
children f93b7fa5e113
comparison
equal deleted inserted replaced
1252:ccb22498f289 1253:bb67a902760b
30 #include "lo-error.h" 30 #include "lo-error.h"
31 #include "f77-uscore.h" 31 #include "f77-uscore.h"
32 32
33 extern "C" 33 extern "C"
34 { 34 {
35 int F77_FCN (dgeesx) (const char*, const char*, 35 int F77_FCN (dgeesx, DGEESX) (const char*, const char*,
36 int (*)(const double&, const double&), 36 int (*)(const double&, const double&),
37 const char*, const int&, double*, const int&, 37 const char*, const int&, double*,
38 int&, double*, double*, double*, const int&, 38 const int&, int&, double*, double*,
39 double&, double&, double*, const int&, int*, 39 double*, const int&, double&, double&,
40 const int&, int*, int&, long, long); 40 double*, const int&, int*, const int&,
41 int*, int&, long, long);
41 } 42 }
42 43
43 static int 44 static int
44 select_ana (const double& a, const double& b) 45 select_ana (const double& a, const double& b)
45 { 46 {
98 bwork = new int [n]; 99 bwork = new int [n];
99 } 100 }
100 101
101 if (*ord == 'A' || *ord == 'a') 102 if (*ord == 'A' || *ord == 'a')
102 { 103 {
103 F77_FCN (dgeesx) (jobvs, sort, select_ana, sense, n, s, n, 104 F77_FCN (dgeesx, DGEESX) (jobvs, sort, select_ana, sense, n, s,
104 sdim, wr, wi, q, n, rconde, rcondv, work, 105 n, sdim, wr, wi, q, n, rconde, rcondv,
105 lwork, iwork, liwork, bwork, info, 1L, 1L); 106 work, lwork, iwork, liwork, bwork,
107 info, 1L, 1L);
106 } 108 }
107 else if (*ord == 'D' || *ord == 'd') 109 else if (*ord == 'D' || *ord == 'd')
108 { 110 {
109 F77_FCN (dgeesx) (jobvs, sort, select_dig, sense, n, s, n, 111 F77_FCN (dgeesx, DGEESX) (jobvs, sort, select_dig, sense, n, s,
110 sdim, wr, wi, q, n, rconde, rcondv, work, 112 n, sdim, wr, wi, q, n, rconde, rcondv,
111 lwork, iwork, liwork, bwork, info, 1L, 1L); 113 work, lwork, iwork, liwork, bwork,
114 info, 1L, 1L);
112 115
113 } 116 }
114 else 117 else
115 { 118 {
116 F77_FCN (dgeesx) (jobvs, sort, (void *) 0, sense, n, s, 119 F77_FCN (dgeesx, DGEESX) (jobvs, sort, (void *) 0, sense, n, s,
117 n, sdim, wr, wi, q, n, rconde, rcondv, 120 n, sdim, wr, wi, q, n, rconde, rcondv,
118 work, lwork, iwork, liwork, bwork, info, 121 work, lwork, iwork, liwork, bwork,
119 1L, 1L); 122 info, 1L, 1L);
120 } 123 }
121 124
122 schur_mat = Matrix (s, n, n); 125 schur_mat = Matrix (s, n, n);
123 unitary_mat = Matrix (q, n, n); 126 unitary_mat = Matrix (q, n, n);
124 127