octave-nkf: liboctave/array/CRowVector.h annotate

annotate liboctave/array/CRowVector.h @ 20633:ffc6cdcd02c5 stable

Fix segfault when complex double matrix calls ZGETRF (bug #45577). * CMatrix.cc (finverse, determinant, rcond, fsolve): Calculate norm of matrix and if it is NaN, skip calling ZGETRF in LAPACK and set info to non-zero value to signal an error.

author	Rik <rik@octave.org>
date	Sat, 10 Oct 2015 16:46:00 -0700
parents	4197fc428c7d
children

rev	line source
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	1 /*
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	2
19731 4197fc428c7d maint: Update copyright notices for 2015. John W. Eaton <jwe@octave.org> parents: 17822 diff changeset	3 Copyright (C) 1994-2015 John W. Eaton
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	4
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	5 This file is part of Octave.
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	6
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	7 Octave is free software; you can redistribute it and/or modify it
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	8 under the terms of the GNU General Public License as published by the
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	9 Free Software Foundation; either version 3 of the License, or (at your
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	10 option) any later version.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	11
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	12 Octave is distributed in the hope that it will be useful, but WITHOUT
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	13 ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	14 FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	15 for more details.
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	16
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	17 You should have received a copy of the GNU General Public License
7016 93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	18 along with Octave; see the file COPYING. If not, see
93c65f2a5668 [project @ 2007-10-12 06:40:56 by jwe] jwe parents: 6108 diff changeset	19 <http://www.gnu.org/licenses/>.
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	20
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	21 */
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	22
17822 ebb3ef964372 maint: Use common #define syntax "octave_filename_h" in h_files. Rik <rik@octave.org> parents: 17769 diff changeset	23 #if !defined (octave_CRowVector_h)
ebb3ef964372 maint: Use common #define syntax "octave_filename_h" in h_files. Rik <rik@octave.org> parents: 17769 diff changeset	24 #define octave_CRowVector_h 1
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	25
1214 0bf4d2b7def4 [project @ 1995-04-06 02:33:59 by jwe] jwe parents: 1205 diff changeset	26 #include "MArray.h"
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	27 #include "dRowVector.h"
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	28
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	29 #include "mx-defs.h"
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	30
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	31 class
6108 143b556ce725 [project @ 2006-10-27 01:45:54 by jwe] jwe parents: 5307 diff changeset	32 OCTAVE_API
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	33 ComplexRowVector : public MArray<Complex>
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	34 {
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	35 friend class ComplexColumnVector;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	36
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	37 public:
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	38
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	39 ComplexRowVector (void) : MArray<Complex> (dim_vector (1, 0)) { }
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	40
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	41 explicit ComplexRowVector (octave_idx_type n)
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	42 : MArray<Complex> (dim_vector (1, n)) { }
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	43
9612 66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	44 explicit ComplexRowVector (const dim_vector& dv) : MArray<Complex> (dv) { }
66970dd627f6 further liboctave design improvements Jaroslav Hajek <highegg@gmail.com> parents: 8920 diff changeset	45
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	46 ComplexRowVector (octave_idx_type n, const Complex& val)
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	47 : MArray<Complex> (dim_vector (1, n), val) { }
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	48
d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	49 ComplexRowVector (const ComplexRowVector& a) : MArray<Complex> (a) { }
d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	50
11570 57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	51 ComplexRowVector (const MArray<Complex>& a)
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	52 : MArray<Complex> (a.as_row ()) { }
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	53
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	54 ComplexRowVector (const Array<Complex>& a)
57632dea2446 attempt better backward compatibility for Array constructors John W. Eaton <jwe@octave.org> parents: 11523 diff changeset	55 : MArray<Complex> (a.as_row ()) { }
3585 d9803711e047 [project @ 2000-02-08 04:35:39 by jwe] jwe parents: 3582 diff changeset	56
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	57 explicit ComplexRowVector (const RowVector& a) : MArray<Complex> (a) { }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	58
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	59 ComplexRowVector& operator = (const ComplexRowVector& a)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	60 {
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	61 MArray<Complex>::operator = (a);
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	62 return *this;
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	63 }
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	64
2386 4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	65 bool operator == (const ComplexRowVector& a) const;
4fc9fd1424a9 [project @ 1996-10-12 18:31:34 by jwe] jwe parents: 1993 diff changeset	66 bool operator != (const ComplexRowVector& a) const;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	67
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	68 // destructive insert/delete/reorder operations
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	69
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	70 ComplexRowVector& insert (const RowVector& a, octave_idx_type c);
23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	71 ComplexRowVector& insert (const ComplexRowVector& a, octave_idx_type c);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	72
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	73 ComplexRowVector& fill (double val);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	74 ComplexRowVector& fill (const Complex& val);
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	75 ComplexRowVector& fill (double val, octave_idx_type c1, octave_idx_type c2);
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	76 ComplexRowVector& fill (const Complex& val,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	77 octave_idx_type c1, octave_idx_type c2);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	78
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	79 ComplexRowVector append (const RowVector& a) const;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	80 ComplexRowVector append (const ComplexRowVector& a) const;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	81
7789 82be108cc558 First attempt at single precision tyeps David Bateman <dbateman@free.fr> parents: 7503 diff changeset	82 ComplexColumnVector hermitian (void) const;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	83 ComplexColumnVector transpose (void) const;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	84
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	85 friend ComplexRowVector conj (const ComplexRowVector& a);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	86
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	87 // resize is the destructive equivalent for this one
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	88
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	89 ComplexRowVector extract (octave_idx_type c1, octave_idx_type c2) const;
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	90
5275 23b37da9fd5b [project @ 2005-04-08 16:07:35 by jwe] jwe parents: 5105 diff changeset	91 ComplexRowVector extract_n (octave_idx_type c1, octave_idx_type n) const;
4316 236c10efcde2 [project @ 2003-02-10 21:57:15 by jwe] jwe parents: 4192 diff changeset	92
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	93 // row vector by row vector -> row vector operations
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	94
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	95 ComplexRowVector& operator += (const RowVector& a);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	96 ComplexRowVector& operator -= (const RowVector& a);
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	97
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	98 // row vector by matrix -> row vector
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	99
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	100 friend ComplexRowVector operator * (const ComplexRowVector& a,
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	101 const ComplexMatrix& b);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	102
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	103 friend ComplexRowVector operator * (const RowVector& a,
10312 cbc402e64d83 untabify liboctave header files John W. Eaton <jwe@octave.org> parents: 10158 diff changeset	104 const ComplexMatrix& b);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	105
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	106 // other operations
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	107
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	108 Complex min (void) const;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	109 Complex max (void) const;
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	110
1359 a6994c934a50 [project @ 1995-09-05 21:30:26 by jwe] jwe parents: 1315 diff changeset	111 // i/o
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	112
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	113 friend std::ostream& operator << (std::ostream& os,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	114 const ComplexRowVector& a);
3504 5eef8a2294bd [project @ 2000-02-01 10:06:51 by jwe] jwe parents: 2847 diff changeset	115 friend std::istream& operator >> (std::istream& is, ComplexRowVector& a);
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	116
14616 13cc11418393 improve handling of default resize fill value for arrays John W. Eaton <jwe@octave.org> parents: 14138 diff changeset	117 void resize (octave_idx_type n, const Complex& rfv = Complex (0))
11574 a83bad07f7e3 attempt better backward compatibility for Array resize functions John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	118 {
a83bad07f7e3 attempt better backward compatibility for Array resize functions John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	119 Array<Complex>::resize (dim_vector (1, n), rfv);
a83bad07f7e3 attempt better backward compatibility for Array resize functions John W. Eaton <jwe@octave.org> parents: 11570 diff changeset	120 }
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	121
12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	122 void clear (octave_idx_type n)
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	123 { Array<Complex>::clear (1, n); }
10350 12884915a8e4 merge MArray classes & improve Array interface Jaroslav Hajek <highegg@gmail.com> parents: 10312 diff changeset	124
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	125 };
38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	126
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	127 // row vector by column vector -> scalar
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	128
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	129 Complex OCTAVE_API operator * (const ComplexRowVector& a,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	130 const ColumnVector& b);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	131
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	132 Complex OCTAVE_API operator * (const ComplexRowVector& a,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	133 const ComplexColumnVector& b);
1205 8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	134
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	135 // other operations
8302fab9fe24 [project @ 1995-04-04 02:05:01 by jwe] jwe parents: 1196 diff changeset	136
17769 49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	137 OCTAVE_API ComplexRowVector linspace (const Complex& x1, const Complex& x2,
49a5a4be04a1 maint: Use GNU style coding conventions for code in liboctave/ Rik <rik@octave.org> parents: 17744 diff changeset	138 octave_idx_type n);
1196 c01e881ccbb1 [project @ 1995-03-31 17:12:19 by jwe] jwe parents: 1195 diff changeset	139
3573 6ae6f1180e62 [project @ 2000-02-04 09:01:59 by jwe] jwe parents: 3504 diff changeset	140 MARRAY_FORWARD_DEFS (MArray, ComplexRowVector, Complex)
6ae6f1180e62 [project @ 2000-02-04 09:01:59 by jwe] jwe parents: 3504 diff changeset	141
458 38cb88095913 [project @ 1994-06-06 00:41:10 by jwe] jwe parents: diff changeset	142 #endif

Mercurial > octave-nkf

annotate liboctave/array/CRowVector.h @ 20633:ffc6cdcd02c5 stable