25 #if ! defined (octave_MDiagArray2_h)
26 #define octave_MDiagArray2_h 1
28 #include "octave-config.h"
93 const T *
d = this->
data ();
128 const MDiagArray2<T>&);
132 #define MDIAGARRAY2_FORWARD_DEFS(B, R, T) \
134 operator + (const R& x) \
136 return R (operator + (dynamic_cast<const B<T>&> (x))); \
139 operator - (const R& x) \
141 return R (operator - (dynamic_cast<const B<T>&> (x))); \
144 operator * (const R& x, const T& y) \
146 return R (operator * (dynamic_cast<const B<T>&> (x), (y))); \
149 operator / (const R& x, const T& y) \
151 return R (operator / (dynamic_cast<const B<T>&> (x), (y))); \
154 operator * (const T& x, const R& y) \
156 return R (operator * ( (x), dynamic_cast<const B<T>&> (y))); \
159 operator + (const R& x, const R& y) \
161 return R (operator + (dynamic_cast<const B<T>&> (x), \
162 dynamic_cast<const B<T>&> (y))); \
165 operator - (const R& x, const R& y) \
167 return R (operator - (dynamic_cast<const B<T>&> (x), \
168 dynamic_cast<const B<T>&> (y))); \
171 product (const R& x, const R& y) \
173 return R (product (dynamic_cast<const B<T>&> (x), \
174 dynamic_cast<const B<T>&> (y))); \
MDiagArray2< T > hermitian(T(*fcn)(const T &)=0) const
MDiagArray2< T > product(const MDiagArray2< T > &, const MDiagArray2< T > &)
octave_idx_type nnz(void) const
MDiagArray2(const DiagArray2< U > &a)
identity matrix If supplied two scalar respectively For allows like xample val
MDiagArray2< T > operator*(const MDiagArray2< T > &, const T &)
MDiagArray2< T > & operator=(const MDiagArray2< T > &a)
MDiagArray2(const Array< T > &a)
Array< T > extract_diag(octave_idx_type k=0) const
const T * data(void) const
MDiagArray2< T > operator-(const MDiagArray2< T > &)
Template for N-dimensional array classes with like-type math operators.
DiagArray2< T > hermitian(T(*fcn)(const T &)=0) const
MDiagArray2(octave_idx_type r, octave_idx_type c, const T &val)
F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_REAL const F77_REAL F77_REAL &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T F77_DBLE &F77_RET_T F77_REAL &F77_RET_T F77_REAL &F77_RET_T F77_DBLE &F77_RET_T const F77_DBLE F77_DBLE &F77_RET_T const F77_REAL F77_REAL &F77_RET_T F77_REAL F77_REAL &F77_RET_T F77_DBLE F77_DBLE &F77_RET_T const F77_DBLE const F77_DBLE F77_DBLE * d
calling an anonymous function involves an overhead quite comparable to the overhead of an m file function Passing a handle to a built in function is because the interpreter is not involved in the internal loop For a
Template for two dimensional diagonal array with math operators.
DiagArray2< T > & operator=(const DiagArray2< T > &a)
MDiagArray2< T > operator/(const MDiagArray2< T > &, const T &)
the sparsity preserving column transformation such that that defines the pivoting threshold can be given in which case it defines the c
MDiagArray2(octave_idx_type r, octave_idx_type c)
MArray< T > diag(octave_idx_type k=0) const
MDiagArray2(const DiagArray2< T > &a)
Array< T > array_value(void) const
MArray< T > array_value() const
N Dimensional Array with copy-on-write semantics.
bool is_multiple_of_identity(T val) const
MDiagArray2(const Array< T > &a, octave_idx_type r, octave_idx_type c)
=val(i)}if ode{val(i)}occurs in table i
MDiagArray2< T > operator+(const MDiagArray2< T > &)
DiagArray2< T > transpose(void) const
MDiagArray2< T > transpose(void) const
MDiagArray2(const MDiagArray2< T > &a)