32 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
40 lu_type Yout (nr, nc, Lfact.nnz () + Ufact.nnz () - (nr<nz?nr:nz));
48 Yout.xridx (ii) = Ufact.ridx (i);
49 Yout.xdata (ii++) = Ufact.data (i);
55 i < Lfact.cidx (j +1); i++)
57 Yout.xridx (ii) = Lfact.ridx (i);
58 Yout.xdata (ii++) = Lfact.data (i);
61 Yout.xcidx (j + 1) = ii;
67 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
74 p_type Pout (nr, nr, nr);
79 Pout.ridx (P (i)) = i;
87 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
97 Pout.
xelem (i) =
static_cast<double> (P(i) + 1);
102 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
109 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
115 p_type Pout (nc, nc, nc);
120 Pout.ridx (i) =
Q (i);
128 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
138 Pout.
xelem (i) =
static_cast<double> (
Q(i) + 1);
143 template <
class lu_type,
class lu_elt_type,
class p_type,
class p_elt_type>
T & xelem(octave_idx_type n)
F77_RET_T const octave_idx_type const octave_idx_type const octave_idx_type double const octave_idx_type double const octave_idx_type double * Q