Mercurial > octave

comparison liboctave/CmplxLU.cc @ 1251:97eac19837dc

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1995-04-11 15:58:32 by jwe]
author jwe
date Tue, 11 Apr 1995 15:58:32 +0000
parents b6360f2d4fa6
children bb67a902760b
comparison
equal deleted inserted replaced
1250:5cca5ae20299 1251:97eac19837dc
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 (zgesv) (const int*, const int*, Complex*, const int*, 35 int F77_FCN (zgesv) (const int&, const int&, Complex*, const int&,
36 int*, Complex*, const int*, int*); 36 int*, Complex*, const int&, int&);
37 } 37 }
38 38
39 ComplexLU::ComplexLU (const ComplexMatrix& a) 39 ComplexLU::ComplexLU (const ComplexMatrix& a)
40 { 40 {
41 int a_nr = a.rows (); 41 int a_nr = a.rows ();
50 50
51 int *ipvt = new int [n]; 51 int *ipvt = new int [n];
52 int *pvt = new int [n]; 52 int *pvt = new int [n];
53 Complex *tmp_data = dup (a.data (), a.length ()); 53 Complex *tmp_data = dup (a.data (), a.length ());
54 int info = 0; 54 int info = 0;
55 int zero = 0; 55 Complex *b;
56 Complex b;
57 56
58 F77_FCN (zgesv) (&n, &zero, tmp_data, &n, ipvt, &b, &n, &info); 57 F77_FCN (zgesv) (n, 0, tmp_data, n, ipvt, b, n, info);
59 58
60 ComplexMatrix A_fact (tmp_data, n, n); 59 ComplexMatrix A_fact (tmp_data, n, n);
61 60
62 int i; 61 int i;
63 62