25 #if ! defined (octave_sparse_lu_h)
26 #define octave_sparse_lu_h 1
28 #include "octave-config.h"
42 template <
typename lu_type>
51 : Lfact (), Ufact (), Rfact (), cond (0),
P (),
Q () { }
58 bool FixedQ =
false,
double droptol = -1.0,
59 bool milu =
false,
bool udiag =
false);
62 : Lfact (a.Lfact), Ufact (a.Ufact), Rfact (), cond (a.cond),
82 lu_type
L (
void)
const {
return Lfact; }
84 lu_type
U (
void)
const {
return Ufact; }
88 lu_type
Y (
void)
const;
106 double rcond (
void)
const {
return cond; }
lu_type::element_type lu_elt_type
MArray< octave_idx_type > Q
Octave interface to the compression and uncompression libraries.
const octave_idx_type * col_perm(void) const
OCTAVE_EXPORT octave_value_list return the value of the option it must match the dimension of the state and the relative tolerance must also be a vector of the same length tem it must match the dimension of the state and the absolute tolerance must also be a vector of the same length The local error test applied at each integration step is xample roup calculate Y_a and Y _d item Given Y
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
const octave_idx_type * row_perm(void) const
SparseMatrix R(void) const
sparse_lu(const sparse_lu &a)
void scale(Matrix &m, double x, double y, double z)
MArray< octave_idx_type > P
F77_RET_T const F77_INT const F77_INT const F77_INT F77_DBLE const F77_INT F77_DBLE const F77_INT F77_DBLE * Q