Mercurial > jwe > octave
view liboctave/array/dDiagMatrix.h @ 31170:06c76eca5f95 default tip @
maint: merge stable to default
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Fri, 29 Jul 2022 12:02:42 -0400 |
parents | 796f54d4ddbf |
children |
line wrap: on
line source
//////////////////////////////////////////////////////////////////////// // // Copyright (C) 1994-2022 The Octave Project Developers // // See the file COPYRIGHT.md in the top-level directory of this // distribution or <https://octave.org/copyright/>. // // This file is part of Octave. // // Octave is free software: you can redistribute it and/or modify it // under the terms of the GNU General Public License as published by // the Free Software Foundation, either version 3 of the License, or // (at your option) any later version. // // Octave is distributed in the hope that it will be useful, but // WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // GNU General Public License for more details. // // You should have received a copy of the GNU General Public License // along with Octave; see the file COPYING. If not, see // <https://www.gnu.org/licenses/>. // //////////////////////////////////////////////////////////////////////// #if ! defined (octave_dDiagMatrix_h) #define octave_dDiagMatrix_h 1 #include "octave-config.h" #include "DET.h" #include "MDiagArray2.h" #include "dColVector.h" #include "dRowVector.h" #include "mx-defs.h" class OCTAVE_API DiagMatrix : public MDiagArray2<double> { public: typedef Matrix full_matrix_type; DiagMatrix (void) = default; DiagMatrix (const DiagMatrix& a) = default; DiagMatrix& operator = (const DiagMatrix& a) = default; ~DiagMatrix (void) = default; DiagMatrix (octave_idx_type r, octave_idx_type c) : MDiagArray2<double> (r, c) { } DiagMatrix (octave_idx_type r, octave_idx_type c, double val) : MDiagArray2<double> (r, c, val) { } DiagMatrix (const MDiagArray2<double>& a) : MDiagArray2<double> (a) { } template <typename U> DiagMatrix (const DiagArray2<U>& a) : MDiagArray2<double> (a) { } explicit DiagMatrix (const Array<double>& a) : MDiagArray2<double> (a) { } DiagMatrix (const Array<double>& a, octave_idx_type r, octave_idx_type c) : MDiagArray2<double> (a, r, c) { } OCTAVE_API bool operator == (const DiagMatrix& a) const; OCTAVE_API bool operator != (const DiagMatrix& a) const; OCTAVE_API DiagMatrix& fill (double val); OCTAVE_API DiagMatrix& fill (double val, octave_idx_type beg, octave_idx_type end); OCTAVE_API DiagMatrix& fill (const ColumnVector& a); OCTAVE_API DiagMatrix& fill (const RowVector& a); OCTAVE_API DiagMatrix& fill (const ColumnVector& a, octave_idx_type beg); OCTAVE_API DiagMatrix& fill (const RowVector& a, octave_idx_type beg); DiagMatrix transpose (void) const { return MDiagArray2<double>::transpose (); } DiagMatrix abs (void) const; friend OCTAVE_API DiagMatrix real (const ComplexDiagMatrix& a); friend OCTAVE_API DiagMatrix imag (const ComplexDiagMatrix& a); // resize is the destructive analog for this one OCTAVE_API Matrix extract (octave_idx_type r1, octave_idx_type c1, octave_idx_type r2, octave_idx_type c2) const; // extract row or column i. OCTAVE_API RowVector row (octave_idx_type i) const; OCTAVE_API RowVector row (char *s) const; OCTAVE_API ColumnVector column (octave_idx_type i) const; OCTAVE_API ColumnVector column (char *s) const; OCTAVE_API DiagMatrix inverse (void) const; OCTAVE_API DiagMatrix inverse (octave_idx_type& info) const; OCTAVE_API DiagMatrix pseudo_inverse (double tol = 0.0) const; // other operations ColumnVector extract_diag (octave_idx_type k = 0) const { return MDiagArray2<double>::extract_diag (k); } OCTAVE_API DET determinant (void) const; OCTAVE_API double rcond (void) const; // i/o friend OCTAVE_API std::ostream& operator << (std::ostream& os, const DiagMatrix& a); }; OCTAVE_API DiagMatrix real (const ComplexDiagMatrix& a); OCTAVE_API DiagMatrix imag (const ComplexDiagMatrix& a); // diagonal matrix by diagonal matrix -> diagonal matrix operations OCTAVE_API DiagMatrix operator * (const DiagMatrix& a, const DiagMatrix& b); MDIAGARRAY2_FORWARD_DEFS (MDiagArray2, DiagMatrix, double) #endif