Mercurial > hg > octave-thorsten
comparison liboctave/EIG.cc @ 1365:9b2654476502
[project @ 1995-09-05 23:44:41 by jwe]
author | jwe |
---|---|
date | Tue, 05 Sep 1995 23:50:05 +0000 |
parents | 611d403c7f3d |
children | 33bb7975f866 |
comparison
equal
deleted
inserted
replaced
1364:353b06ad0194 | 1365:9b2654476502 |
---|---|
69 Matrix vr (n, n); | 69 Matrix vr (n, n); |
70 double *pvr = vr.fortran_vec (); | 70 double *pvr = vr.fortran_vec (); |
71 int lwork = 8*n; | 71 int lwork = 8*n; |
72 double *work = new double[lwork]; | 72 double *work = new double[lwork]; |
73 | 73 |
74 double *dummy; | 74 double *dummy = 0; |
75 int idummy = 1; | 75 int idummy = 1; |
76 | 76 |
77 F77_FCN (dgeev, DGEEV) ("N", "V", n, tmp_data, n, wr, wi, dummy, | 77 F77_FCN (dgeev, DGEEV) ("N", "V", n, tmp_data, n, wr, wi, dummy, |
78 idummy, pvr, n, work, lwork, info, 1L, 1L); | 78 idummy, pvr, n, work, lwork, info, 1L, 1L); |
79 | 79 |
141 | 141 |
142 int lwork = 8*n; | 142 int lwork = 8*n; |
143 Complex *work = new Complex[lwork]; | 143 Complex *work = new Complex[lwork]; |
144 double *rwork = new double[4*n]; | 144 double *rwork = new double[4*n]; |
145 | 145 |
146 Complex *dummy; | 146 Complex *dummy = 0; |
147 int idummy = 1; | 147 int idummy = 1; |
148 | 148 |
149 F77_FCN (zgeev, ZGEEV) ("N", "V", n, tmp_data, n, pw, dummy, idummy, | 149 F77_FCN (zgeev, ZGEEV) ("N", "V", n, tmp_data, n, pw, dummy, idummy, |
150 pvr, n, work, lwork, rwork, info, 1L, 1L); | 150 pvr, n, work, lwork, rwork, info, 1L, 1L); |
151 | 151 |