50 template <
class T1,
class T2>
69 #define INSTANTIATE_MX_LEFTDIV_CONFORM(T1, T2) \
70 template bool mx_leftdiv_conform (const T1&, const T2&)
85 template <
class T1,
class T2>
104 #define INSTANTIATE_MX_DIV_CONFORM(T1, T2) \
105 template bool mx_div_conform (const T1&, const T2&)
146 Matrix result = btmp.
solve (btyp, atmp, info, rcond,
294 template <
typename RT,
typename SM,
typename DM>
311 typedef typename DM::element_type DM_elt_type;
312 const DM_elt_type zero = DM_elt_type ();
318 const DM_elt_type s = d.dgelem (j);
320 r.xcidx (j) = k_result;
324 r.xdata (k_result) = a.data (k) / s;
325 r.xridx (k_result) = a.ridx (k);
329 r.xcidx (nc) = k_result;
331 r.maybe_compress (
true);
339 return do_rightdiv_sm_dm<SparseMatrix> (a, b);
346 return do_rightdiv_sm_dm<SparseComplexMatrix> (a, b);
353 return do_rightdiv_sm_dm<SparseComplexMatrix> (a, b);
360 return do_rightdiv_sm_dm<SparseComplexMatrix> (a, b);
561 template <
typename RT,
typename DM,
typename SM>
578 typedef typename DM::element_type DM_elt_type;
579 const DM_elt_type zero = DM_elt_type ();
586 r.xcidx (j) = k_result;
592 const DM_elt_type s = d.dgelem (i);
595 r.xdata (k_result) = a.data (k) / s;
596 r.xridx (k_result) = i;
602 r.xcidx (a_nc) = k_result;
604 r.maybe_compress (
true);
612 return do_leftdiv_dm_sm<SparseMatrix> (
d, a);
619 return do_leftdiv_dm_sm<SparseComplexMatrix> (
d, a);
626 return do_leftdiv_dm_sm<SparseComplexMatrix> (
d, a);
634 return do_leftdiv_dm_sm<SparseComplexMatrix> (
d, a);
octave_idx_type cols(void) const
void gripe_nonconformant(const char *op, octave_idx_type op1_len, octave_idx_type op2_len)
octave_idx_type rows(void) const
static void solve_singularity_warning(double rcond)
ComplexMatrix solve(MatrixType &typ, const Matrix &b) const
SparseMatrix transpose(void) const
octave_idx_type * cidx(void)
bool mx_leftdiv_conform(const T1 &a, const T2 &b)
T & elem(octave_idx_type n)
Matrix solve(MatrixType &typ, const Matrix &b) const
F77_RET_T const double const double double * d
ComplexMatrix hermitian(void) const
Matrix x_el_div(double a, const SparseMatrix &b)
MatrixType transpose(void) const
#define INSTANTIATE_MX_LEFTDIV_CONFORM(T1, T2)
Matrix transpose(void) const
Matrix xleftdiv(const SparseMatrix &a, const Matrix &b, MatrixType &typ)
void gripe_singular_matrix(double rcond)
SparseComplexMatrix hermitian(void) const
bool mx_div_conform(const T1 &a, const T2 &b)
octave_idx_type * ridx(void)
Matrix xdiv(const Matrix &a, const SparseMatrix &b, MatrixType &typ)
#define INSTANTIATE_MX_DIV_CONFORM(T1, T2)
std::complex< double > Complex
RT do_leftdiv_dm_sm(const DM &d, const SM &a)
RT do_rightdiv_sm_dm(const SM &a, const DM &d)
charNDArray min(char d, const charNDArray &m)