Mercurial > octave
view liboctave/array/MDiagArray2.h @ 22197:e43d83253e28
refill multi-line macro definitions
Use the Emacs C++ mode style for line continuation markers in
multi-line macro definitions.
* make_int.cc, __dsearchn__.cc, __magick_read__.cc, besselj.cc,
bitfcns.cc, bsxfun.cc, cellfun.cc, data.cc, defun-dld.h, defun-int.h,
defun.h, det.cc, error.h, find.cc, gcd.cc, graphics.cc, interpreter.h,
jit-ir.h, jit-typeinfo.h, lookup.cc, ls-mat5.cc, max.cc, mexproto.h,
mxarray.in.h, oct-stream.cc, ordschur.cc, pr-output.cc, profiler.h,
psi.cc, regexp.cc, sparse-xdiv.cc, sparse-xpow.cc, tril.cc, txt-eng.h,
utils.cc, variables.cc, variables.h, xdiv.cc, xpow.cc, __glpk__.cc,
ov-base.cc, ov-base.h, ov-cell.cc, ov-ch-mat.cc, ov-classdef.cc,
ov-complex.cc, ov-cx-mat.cc, ov-cx-sparse.cc, ov-float.cc, ov-float.h,
ov-flt-complex.cc, ov-flt-cx-mat.cc, ov-flt-re-mat.cc,
ov-int-traits.h, ov-lazy-idx.h, ov-perm.cc, ov-re-mat.cc,
ov-re-sparse.cc, ov-scalar.cc, ov-scalar.h, ov-str-mat.cc,
ov-type-conv.h, ov.cc, ov.h, op-class.cc, op-int-conv.cc, op-int.h,
op-str-str.cc, ops.h, lex.ll, Array.cc, CMatrix.cc, CSparse.cc,
MArray.cc, MArray.h, MDiagArray2.cc, MDiagArray2.h, MSparse.h,
Sparse.cc, dMatrix.cc, dSparse.cc, fCMatrix.cc, fMatrix.cc,
idx-vector.cc, f77-fcn.h, quit.h, bsxfun-decl.h, bsxfun-defs.cc,
lo-specfun.cc, oct-convn.cc, oct-convn.h, oct-norm.cc, oct-norm.h,
oct-rand.cc, Sparse-op-decls.h, Sparse-op-defs.h, mx-inlines.cc,
mx-op-decl.h, mx-op-defs.h, mach-info.cc, oct-group.cc, oct-passwd.cc,
oct-syscalls.cc, oct-time.cc, data-conv.cc, kpse.cc, lo-ieee.h,
lo-macros.h, oct-cmplx.h, oct-glob.cc, oct-inttypes.cc,
oct-inttypes.h, oct-locbuf.h, oct-sparse.h, url-transfer.cc,
oct-conf-post.in.h, shared-fcns.h: Refill macro definitions.
author | John W. Eaton <jwe@octave.org> |
---|---|
date | Mon, 01 Aug 2016 12:40:18 -0400 |
parents | 1473547f50f5 |
children | bac0d6f07a3e |
line wrap: on
line source
// Template array classes with like-type math ops /* Copyright (C) 1996-2015 John W. Eaton Copyright (C) 2010 VZLU Prague 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 <http://www.gnu.org/licenses/>. */ #if ! defined (octave_MDiagArray2_h) #define octave_MDiagArray2_h 1 #include "octave-config.h" #include "DiagArray2.h" #include "MArray.h" template <typename T> class MDiagArray2; template <typename T> MDiagArray2<T> operator + (const MDiagArray2<T>&); template <typename T> MDiagArray2<T> operator - (const MDiagArray2<T>&); template <typename T> MDiagArray2<T> operator * (const MDiagArray2<T>&, const T&); template <typename T> MDiagArray2<T> operator / (const MDiagArray2<T>&, const T&); template <typename T> MDiagArray2<T> operator * (const T&, const MDiagArray2<T>&); template <typename T> MDiagArray2<T> operator + (const MDiagArray2<T>&, const MDiagArray2<T>&); template <typename T> MDiagArray2<T> operator - (const MDiagArray2<T>&, const MDiagArray2<T>&); template <typename T> MDiagArray2<T> product (const MDiagArray2<T>&, const MDiagArray2<T>&); //! Template for two dimensional diagonal array with math operators. template <typename T> class MDiagArray2 : public DiagArray2<T> { public: MDiagArray2 (void) : DiagArray2<T> () { } MDiagArray2 (octave_idx_type r, octave_idx_type c) : DiagArray2<T> (r, c) { } MDiagArray2 (octave_idx_type r, octave_idx_type c, const T& val) : DiagArray2<T> (r, c, val) { } MDiagArray2 (const MDiagArray2<T>& a) : DiagArray2<T> (a) { } MDiagArray2 (const DiagArray2<T>& a) : DiagArray2<T> (a) { } template <typename U> MDiagArray2 (const DiagArray2<U>& a) : DiagArray2<T> (a) { } explicit MDiagArray2 (const Array<T>& a) : DiagArray2<T> (a) { } MDiagArray2 (const Array<T>& a, octave_idx_type r, octave_idx_type c) : DiagArray2<T> (a, r, c) { } ~MDiagArray2 (void) { } MDiagArray2<T>& operator = (const MDiagArray2<T>& a) { DiagArray2<T>::operator = (a); return *this; } MArray<T> array_value () const { return DiagArray2<T>::array_value (); } octave_idx_type nnz (void) const { octave_idx_type retval = 0; const T *d = this->data (); octave_idx_type nel = this->length (); for (octave_idx_type i = 0; i < nel; i++) { if (d[i] != T ()) retval++; } return retval; } MArray<T> diag (octave_idx_type k = 0) const { return DiagArray2<T>::extract_diag (k); } MDiagArray2<T> transpose (void) const { return DiagArray2<T>::transpose (); } MDiagArray2<T> hermitian (T (*fcn) (const T&) = 0) const { return DiagArray2<T>::hermitian (fcn); } bool is_multiple_of_identity (T val) const; // Currently, the OPS functions don't need to be friends, but that // may change. friend MDiagArray2<T> operator + <> (const MDiagArray2<T>&); friend MDiagArray2<T> operator - <> (const MDiagArray2<T>&); friend MDiagArray2<T> operator * <> (const MDiagArray2<T>&, const T&); friend MDiagArray2<T> operator / <> (const MDiagArray2<T>&, const T&); friend MDiagArray2<T> operator * <> (const T&, const MDiagArray2<T>&); friend MDiagArray2<T> operator + <> (const MDiagArray2<T>&, const MDiagArray2<T>&); friend MDiagArray2<T> operator - <> (const MDiagArray2<T>&, const MDiagArray2<T>&); friend MDiagArray2<T> product <> (const MDiagArray2<T>&, const MDiagArray2<T>&); }; #define MDIAGARRAY2_FORWARD_DEFS(B, R, T) \ inline R \ operator + (const R& x) \ { \ return R (operator + (dynamic_cast<const B<T>&> (x))); \ } \ inline R \ operator - (const R& x) \ { \ return R (operator - (dynamic_cast<const B<T>&> (x))); \ } \ inline R \ operator * (const R& x, const T& y) \ { \ return R (operator * (dynamic_cast<const B<T>&> (x), (y))); \ } \ inline R \ operator / (const R& x, const T& y) \ { \ return R (operator / (dynamic_cast<const B<T>&> (x), (y))); \ } \ inline R \ operator * (const T& x, const R& y) \ { \ return R (operator * ( (x), dynamic_cast<const B<T>&> (y))); \ } \ inline R \ operator + (const R& x, const R& y) \ { \ return R (operator + (dynamic_cast<const B<T>&> (x), \ dynamic_cast<const B<T>&> (y))); \ } \ inline R \ operator - (const R& x, const R& y) \ { \ return R (operator - (dynamic_cast<const B<T>&> (x), \ dynamic_cast<const B<T>&> (y))); \ } \ inline R \ product (const R& x, const R& y) \ { \ return R (product (dynamic_cast<const B<T>&> (x), \ dynamic_cast<const B<T>&> (y))); \ } #endif