25 #if ! defined (octave_sparse_lu_h) 26 #define octave_sparse_lu_h 1 28 #include "octave-config.h" 46 template <
typename lu_type>
55 : Lfact (), Ufact (), Rfact (), cond (0), P (),
Q () { }
62 bool FixedQ =
false,
double droptol = -1.0,
63 bool milu =
false,
bool udiag =
false);
66 : Lfact (
a.Lfact), Ufact (
a.Ufact), Rfact (), cond (
a.cond),
86 lu_type
L (
void)
const {
return Lfact; }
88 lu_type
U (
void)
const {
return Ufact; }
92 lu_type Y (
void)
const;
110 double rcond (
void)
const {
return cond; }
lu_type::element_type lu_elt_type
MArray< octave_idx_type > Q
SparseMatrix R(void) 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
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
sparse_lu(const sparse_lu &a)
void scale(Matrix &m, double x, double y, double z)
const octave_idx_type * row_perm(void) const
const octave_idx_type * col_perm(void) const
MArray< octave_idx_type > P