24 #if defined (HAVE_CONFIG_H) 44 if (len !=
a.numel ())
60 if (r < 0 || r + a_len >
numel ())
61 (*current_liboctave_error_handler) (
"range error for insert");
95 if (r1 < 0 || r2 < 0 || r1 >= len || r2 >= len)
96 (*current_liboctave_error_handler) (
"range error for fill");
131 return do_mx_unary_map<double, double, std::abs> (*this);
184 F77_INT a_len = octave::to_f77_int (
a.numel ());
197 double *
y =
retval.fortran_vec ();
199 F77_XFCN (dgemv, DGEMV, (F77_CONST_CHAR_ARG2 (
"N", 1),
200 nr, nc, 1.0, m.
data (), nr,
201 a.data (), 1, 0.0,
y, 1
202 F77_CHAR_ARG_LEN (1)));
219 F77_INT a_len = octave::to_f77_int (
a.numel ());
224 if (nr == 0 || nc == 0)
249 double res =
elem (0);
265 double res =
elem (0);
279 os <<
a.elem (
i) <<
"\n";
octave_idx_type rows(void) const
ColumnVector & fill(double val)
ColumnVector operator*(const Matrix &m, const ColumnVector &a)
const double * data(void) const
identity matrix If supplied two scalar respectively For allows like xample val
std::istream & operator>>(std::istream &is, ColumnVector &a)
void mx_inline_real(size_t n, T *r, const std::complex< T > *x)
double & elem(octave_idx_type n)
octave_idx_type rows(void) const
MArray< T > transpose(void) const
#define F77_XFCN(f, F, args)
T elem(octave_idx_type r, octave_idx_type c) const
void mx_inline_imag(size_t n, T *r, const std::complex< T > *x)
octave_idx_type cols(void) const
octave_value resize(const dim_vector &dv, bool fill=false) const
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
ColumnVector abs(void) const
octave_idx_type cols(void) const
void err_nonconformant(const char *op, octave_idx_type op1_len, octave_idx_type op2_len)
ColumnVector extract(octave_idx_type r1, octave_idx_type r2) const
bool operator!=(const ColumnVector &a) const
With real return the complex result
double & xelem(octave_idx_type n)
ColumnVector extract_n(octave_idx_type r1, octave_idx_type n) const
octave_f77_int_type F77_INT
ColumnVector stack(const ColumnVector &a) const
ColumnVector & insert(const ColumnVector &a, octave_idx_type r)
bool operator==(const ColumnVector &a) const
RowVector transpose(void) const
the element is set to zero In other the statement xample y
ColumnVector imag(const ComplexColumnVector &a)
std::ostream & operator<<(std::ostream &os, const ColumnVector &a)
octave_idx_type numel(void) const
Number of elements in the array.
ColumnVector real(const ComplexColumnVector &a)
bool mx_inline_equal(size_t n, const T1 *x, const T2 *y)
write the output to stdout if nargout is