Mercurial > octave

comparison liboctave/CMatrix.cc @ 1195:0da1cc2c1947

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
[project @ 1995-03-31 16:11:43 by jwe]
author jwe
date Fri, 31 Mar 1995 16:16:34 +0000
parents b6360f2d4fa6
children 8302fab9fe24
comparison
equal deleted inserted replaced
1194:621fef7bcca1 1195:0da1cc2c1947
932 Complex *z = new Complex [nr]; 932 Complex *z = new Complex [nr];
933 Complex *tmp_data = dup (data (), len); 933 Complex *tmp_data = dup (data (), len);
934 934
935 F77_FCN (zgeco) (tmp_data, &nr, &nc, ipvt, &rcond, z); 935 F77_FCN (zgeco) (tmp_data, &nr, &nc, ipvt, &rcond, z);
936 936
937 volatile double tmp_rcond = rcond; 937 volatile double rcond_plus_one = rcond + 1.0;
938 if (tmp_rcond + 1.0 == 1.0) 938 if (rcond_plus_one == 1.0)
939 { 939 {
940 info = -1; 940 info = -1;
941 copy (tmp_data, data (), len); // Restore contents. 941 copy (tmp_data, data (), len); // Restore contents.
942 } 942 }
943 else 943 else
1203 Complex *z = new Complex [nr]; 1203 Complex *z = new Complex [nr];
1204 Complex *tmp_data = dup (data (), length ()); 1204 Complex *tmp_data = dup (data (), length ());
1205 1205
1206 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); 1206 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z);
1207 1207
1208 volatile double tmp_rcond = rcond; 1208 volatile double rcond_plus_one = rcond + 1.0;
1209 if (tmp_rcond + 1.0 == 1.0) 1209 if (rcond_plus_one == 1.0)
1210 { 1210 {
1211 info = -1; 1211 info = -1;
1212 retval = ComplexDET (); 1212 retval = ComplexDET ();
1213 } 1213 }
1214 else 1214 else
1285 Complex *z = new Complex [nr]; 1285 Complex *z = new Complex [nr];
1286 Complex *tmp_data = dup (data (), length ()); 1286 Complex *tmp_data = dup (data (), length ());
1287 1287
1288 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); 1288 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z);
1289 1289
1290 volatile double tmp_rcond = rcond; 1290 volatile double rcond_plus_one = rcond + 1.0;
1291 if (tmp_rcond + 1.0 == 1.0) 1291 if (rcond_plus_one == 1.0)
1292 { 1292 {
1293 info = -2; 1293 info = -2;
1294 } 1294 }
1295 else 1295 else
1296 { 1296 {
1348 Complex *z = new Complex [nr]; 1348 Complex *z = new Complex [nr];
1349 Complex *tmp_data = dup (data (), length ()); 1349 Complex *tmp_data = dup (data (), length ());
1350 1350
1351 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z); 1351 F77_FCN (zgeco) (tmp_data, &nr, &nr, ipvt, &rcond, z);
1352 1352
1353 volatile double tmp_rcond = rcond; 1353 volatile double rcond_plus_one = rcond + 1.0;
1354 if (tmp_rcond + 1.0 == 1.0) 1354 if (rcond_plus_one == 1.0)
1355 { 1355 {
1356 info = -2; 1356 info = -2;
1357 } 1357 }
1358 else 1358 else
1359 { 1359 {