23 #if defined (HAVE_CONFIG_H)
37 template <
typename MT>
48 template <
typename MT>
53 if (U.is_square () && fact.
regular ())
141 bool issparse = (nargin > 0 &&
args(0).is_sparse_type ());
143 if (nargin < 1 || (issparse && nargin > 3) || (! issparse && nargin > 2))
152 if (
args(n).is_string ())
159 error (
"lu: unrecognized string argument");
164 error (
"lu: can not define pivoting threshold THRES for full matrices");
167 if (tmp.
numel () == 1)
173 else if (tmp.
numel () == 2)
176 error (
"lu: THRES must be a 1 or 2-element vector");
235 retval.
resize (scale ? 5 : 4);
301 retval.
resize (scale ? 5 : 4);
345 retval =
ovl (fact.
Y ());
379 retval =
ovl (fact.
Y ());
416 retval =
ovl (fact.
Y ());
450 retval =
ovl (fact.
Y ());
610 if (nargin != 4 && nargin != 5)
613 bool pivoted = (nargin == 5);
624 error (
"luupdate: L, U, X, and Y must be numeric");
626 if (! check_lu_dims (argl, argu, argp))
627 error (
"luupdate: dimension mismatch");
ColumnVector Pc_vec(void) const
bool is_real_type(void) const
static octave_value get_lu_l(const octave::math::lu< MT > &fact)
octave_idx_type rows(void) const
void resize(octave_idx_type nr, octave_idx_type nc, double rfv=0)
FloatComplexMatrix float_complex_matrix_value(bool frc_str_conv=false) const
OCTAVE_EXPORT octave_value_list isa nd deftypefn *return ovl(args(0).is_integer_type())
OCTINTERP_API void print_usage(void)
octave_idx_type numel(void) const
Number of elements in the array.
PermMatrix Pr_mat(void) const
bool is_numeric_type(void) const
#define DEFUN(name, args_name, nargout_name, doc)
void error(const char *fmt,...)
PermMatrix transpose(void) const
bool is_perm_matrix(void) const
ColumnVector Pr_vec(void) const
const octave_idx_type * col_perm(void) const
octave_idx_type columns(void) const
static octave_value get_lu_u(const octave::math::lu< MT > &fact)
OCTAVE_EXPORT octave_value_list return the number of command line arguments passed to Octave If called with the optional argument the function xample nargout(@histc)
const octave_idx_type * row_perm(void) const
SparseMatrix R(void) const
nd deftypefn *octave_map m
bool is_complex_type(void) const
void update_piv(const VT &u, const VT &v)
SparseComplexMatrix sparse_complex_matrix_value(bool frc_str_conv=false) const
Matrix matrix_value(bool frc_str_conv=false) const
void err_wrong_type_arg(const char *name, const char *s)
bool is_empty(void) const
=val(i)}if ode{val(i)}occurs in table i
ComplexMatrix complex_matrix_value(bool frc_str_conv=false) const
FloatMatrix float_matrix_value(bool frc_str_conv=false) const
PermMatrix Pc_mat(void) const
the element is set to zero In other the statement xample y
void scale(Matrix &m, double x, double y, double z)
bool is_undefined(void) const
void update(const VT &u, const VT &v)
void resize(octave_idx_type n, const octave_value &rfv=octave_value())
ColumnVector P_vec(void) const
PermMatrix perm_matrix_value(void) const
SparseMatrix sparse_matrix_value(bool frc_str_conv=false) const
bool is_single_type(void) const
static PermMatrix eye(octave_idx_type n)
If this string is the system will ring the terminal sometimes it is useful to be able to print the original representation of the string
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 * x
return octave_value(v1.char_array_value().concat(v2.char_array_value(), ra_idx),((a1.is_sq_string()||a2.is_sq_string())? '\'': '"'))
charNDArray min(char d, const charNDArray &m)