Mercurial > octave

comparison liboctave/CmplxAEPBAL.cc @ 1360:7eb93d12654c

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1995-09-05 21:51:54 by jwe]
author jwe
date Tue, 05 Sep 1995 21:51:54 +0000
parents 611d403c7f3d
children a744f4d0ba59
comparison
equal deleted inserted replaced
1359:a6994c934a50 1360:7eb93d12654c
49 ComplexAEPBALANCE::init (const ComplexMatrix& a, const char *balance_job) 49 ComplexAEPBALANCE::init (const ComplexMatrix& a, const char *balance_job)
50 { 50 {
51 51
52 int n = a.cols (); 52 int n = a.cols ();
53 53
54 // Parameters for balance call. 54 // Parameters for balance call.
55 55
56 int info; 56 int info;
57 int ilo; 57 int ilo;
58 int ihi; 58 int ihi;
59 double *scale = new double [n]; 59 double *scale = new double [n];
60 60
61 // Copy matrix into local structure. 61 // Copy matrix into local structure.
62 62
63 balanced_mat = a; 63 balanced_mat = a;
64 64
65 F77_FCN (zgebal, ZGEBAL) (balance_job, n, 65 F77_FCN (zgebal, ZGEBAL) (balance_job, n,
66 balanced_mat.fortran_vec (), n, ilo, ihi, 66 balanced_mat.fortran_vec (), n, ilo, ihi,
67 scale, info, 1L, 1L); 67 scale, info, 1L, 1L);
68 68
69 // Initialize balancing matrix to identity. 69 // Initialize balancing matrix to identity.
70 70
71 balancing_mat = Matrix (n, n, 0.0); 71 balancing_mat = Matrix (n, n, 0.0);
72 for (int i = 0; i < n; i++) 72 for (int i = 0; i < n; i++)
73 balancing_mat (i, i) = 1.0; 73 balancing_mat (i, i) = 1.0;
74 74