sparse-dmsolve.cc File Reference

#include <algorithm>
#include "CMatrix.h"
#include "CSparse.h"
#include "MArray.h"
#include "MSparse.h"
#include "MatrixType.h"
#include "dSparse.h"
#include "lo-error.h"
#include "oct-inttypes-fwd.h"
#include "oct-locbuf.h"
#include "oct-sort.h"
#include "oct-sparse.h"
#include "quit.h"
#include "sparse-dmsolve.h"
#include "sparse-qr.h"

Include dependency graph for sparse-dmsolve.cc:

Go to the source code of this file.

Functions
template<typename RT , typename ST , typename T >
RT	dmsolve (const ST &a, const T &b, octave_idx_type &info)
template OCTAVE_API ComplexMatrix	dmsolve< ComplexMatrix, SparseComplexMatrix, ComplexMatrix > (const SparseComplexMatrix &, const ComplexMatrix &, octave_idx_type &)
template OCTAVE_API ComplexMatrix	dmsolve< ComplexMatrix, SparseComplexMatrix, Matrix > (const SparseComplexMatrix &, const Matrix &, octave_idx_type &)
template OCTAVE_API ComplexMatrix	dmsolve< ComplexMatrix, SparseMatrix, ComplexMatrix > (const SparseMatrix &, const ComplexMatrix &, octave_idx_type &)
template OCTAVE_API Matrix	dmsolve< Matrix, SparseMatrix, Matrix > (const SparseMatrix &, const Matrix &, octave_idx_type &)
template OCTAVE_API SparseComplexMatrix	dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseComplexMatrix > (const SparseComplexMatrix &, const SparseComplexMatrix &, octave_idx_type &)
template OCTAVE_API SparseComplexMatrix	dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseMatrix > (const SparseComplexMatrix &, const SparseMatrix &, octave_idx_type &)
template OCTAVE_API SparseComplexMatrix	dmsolve< SparseComplexMatrix, SparseMatrix, SparseComplexMatrix > (const SparseMatrix &, const SparseComplexMatrix &, octave_idx_type &)
template OCTAVE_API SparseMatrix	dmsolve< SparseMatrix, SparseMatrix, SparseMatrix > (const SparseMatrix &, const SparseMatrix &, octave_idx_type &)
template<typename T >
static MArray< T >	dmsolve_extract (const MArray< T > &m, const octave_idx_type *, const octave_idx_type *, octave_idx_type r1, octave_idx_type r2, octave_idx_type c1, octave_idx_type c2)
template<typename T >
static MSparse< T >	dmsolve_extract (const MSparse< T > &A, const octave_idx_type *Pinv, const octave_idx_type *Q, octave_idx_type rst, octave_idx_type rend, octave_idx_type cst, octave_idx_type cend, octave_idx_type maxnz=-1, bool lazy=false)
template<typename T >
static void	dmsolve_insert (MArray< T > &a, const MArray< T > &b, const octave_idx_type *Q, octave_idx_type r, octave_idx_type c)
template<typename T >
static void	dmsolve_insert (MSparse< T > &a, const MSparse< T > &b, const octave_idx_type *Q, octave_idx_type r, octave_idx_type c)
template<typename T , typename RT >
static void	dmsolve_permute (MArray< RT > &a, const MArray< T > &b, const octave_idx_type *p)
template<typename T , typename RT >
static void	dmsolve_permute (MSparse< RT > &a, const MSparse< T > &b, const octave_idx_type *p)
static void	solve_singularity_warning (double)

Function Documentation

dmsolve()

template<typename RT , typename ST , typename T >

RT dmsolve	(	const ST &	a,
		const T &	b,
		octave_idx_type &	info
	)

Definition at line 368 of file sparse-dmsolve.cc.

References CXSPARSE_DNAME, d, dmsolve_extract(), dmsolve_insert(), dmsolve_permute(), MatrixType::Full, m, OCTAVE_LOCAL_BUFFER, qrsolve(), solve_singularity_warning(), to_octave_idx_type_ptr(), and to_suitesparse_intptr().

Referenced by SparseComplexMatrix::solve().

dmsolve< ComplexMatrix, SparseComplexMatrix, ComplexMatrix >()

template OCTAVE_API ComplexMatrix dmsolve< ComplexMatrix, SparseComplexMatrix, ComplexMatrix >	(	const SparseComplexMatrix &	,
		const ComplexMatrix &	,
		octave_idx_type &
	)

Referenced by SparseComplexMatrix::solve().

dmsolve< ComplexMatrix, SparseComplexMatrix, Matrix >()

template OCTAVE_API ComplexMatrix dmsolve< ComplexMatrix, SparseComplexMatrix, Matrix >	(	const SparseComplexMatrix &	,
		const Matrix &	,
		octave_idx_type &
	)

Referenced by SparseComplexMatrix::solve().

dmsolve< ComplexMatrix, SparseMatrix, ComplexMatrix >()

template OCTAVE_API ComplexMatrix dmsolve< ComplexMatrix, SparseMatrix, ComplexMatrix >	(	const SparseMatrix &	,
		const ComplexMatrix &	,
		octave_idx_type &
	)

Referenced by SparseMatrix::solve().

dmsolve< Matrix, SparseMatrix, Matrix >()

template OCTAVE_API Matrix dmsolve< Matrix, SparseMatrix, Matrix >	(	const SparseMatrix &	,
		const Matrix &	,
		octave_idx_type &
	)

Referenced by SparseMatrix::solve().

dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseComplexMatrix >()

template OCTAVE_API SparseComplexMatrix dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseComplexMatrix >	(	const SparseComplexMatrix &	,
		const SparseComplexMatrix &	,
		octave_idx_type &
	)

dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseMatrix >()

template OCTAVE_API SparseComplexMatrix dmsolve< SparseComplexMatrix, SparseComplexMatrix, SparseMatrix >	(	const SparseComplexMatrix &	,
		const SparseMatrix &	,
		octave_idx_type &
	)

Referenced by SparseComplexMatrix::solve().

dmsolve< SparseComplexMatrix, SparseMatrix, SparseComplexMatrix >()

template OCTAVE_API SparseComplexMatrix dmsolve< SparseComplexMatrix, SparseMatrix, SparseComplexMatrix >	(	const SparseMatrix &	,
		const SparseComplexMatrix &	,
		octave_idx_type &
	)

Referenced by SparseMatrix::solve().

dmsolve< SparseMatrix, SparseMatrix, SparseMatrix >()

template OCTAVE_API SparseMatrix dmsolve< SparseMatrix, SparseMatrix, SparseMatrix >	(	const SparseMatrix &	,
		const SparseMatrix &	,
		octave_idx_type &
	)

Referenced by SparseMatrix::solve().

dmsolve_extract() [1/2]

template<typename T >

static MArray<T> dmsolve_extract ( const MArray< T > &  m, const octave_idx_type *  , const octave_idx_type *  , octave_idx_type  r1, octave_idx_type  r2, octave_idx_type  c1, octave_idx_type  c2  )

static

Definition at line 143 of file sparse-dmsolve.cc.

References m, and Array< T, Alloc >::xelem().

dmsolve_extract() [2/2]

template<typename T >

static MSparse<T> dmsolve_extract ( const MSparse< T > &  A, const octave_idx_type *  Pinv, const octave_idx_type *  Q, octave_idx_type  rst, octave_idx_type  rend, octave_idx_type  cst, octave_idx_type  cend, octave_idx_type  maxnz = -1, bool  lazy = false  )

static

Definition at line 49 of file sparse-dmsolve.cc.

References A, B, OCTAVE_LOCAL_BUFFER, Q, r, and octave_sort< T >::sort().

Referenced by dmsolve().

dmsolve_insert() [1/2]

template<typename T >

static void dmsolve_insert ( MArray< T > &  a, const MArray< T > &  b, const octave_idx_type *  Q, octave_idx_type  r, octave_idx_type  c  )

static

Definition at line 173 of file sparse-dmsolve.cc.

References Array< T, Alloc >::cols(), Array< T, Alloc >::data(), Array< T, Alloc >::fortran_vec(), Q, r, and Array< T, Alloc >::rows().

Referenced by dmsolve().

dmsolve_insert() [2/2]

template<typename T >

static void dmsolve_insert ( MSparse< T > &  a, const MSparse< T > &  b, const octave_idx_type *  Q, octave_idx_type  r, octave_idx_type  c  )

static

Definition at line 200 of file sparse-dmsolve.cc.

References Sparse< T, Alloc >::cidx(), Sparse< T, Alloc >::cols(), Sparse< T, Alloc >::data(), Sparse< T, Alloc >::nnz(), OCTAVE_LOCAL_BUFFER, Q, r, Sparse< T, Alloc >::ridx(), Sparse< T, Alloc >::rows(), octave_sort< T >::sort(), Sparse< T, Alloc >::xcidx(), Sparse< T, Alloc >::xdata(), and Sparse< T, Alloc >::xridx().

dmsolve_permute() [1/2]

template<typename T , typename RT >

static void dmsolve_permute ( MArray< RT > &  a, const MArray< T > &  b, const octave_idx_type *  p  )

static

Definition at line 293 of file sparse-dmsolve.cc.

References Array< T, Alloc >::cols(), Array< T, Alloc >::data(), Array< T, Alloc >::fortran_vec(), Array< T, Alloc >::resize(), and Array< T, Alloc >::rows().

Referenced by dmsolve().

dmsolve_permute() [2/2]

template<typename T , typename RT >

static void dmsolve_permute ( MSparse< RT > &  a, const MSparse< T > &  b, const octave_idx_type *  p  )

static

Definition at line 317 of file sparse-dmsolve.cc.

References Sparse< T, Alloc >::cidx(), Sparse< T, Alloc >::cols(), Sparse< T, Alloc >::data(), Sparse< T, Alloc >::nnz(), OCTAVE_LOCAL_BUFFER, r, Sparse< T, Alloc >::ridx(), Sparse< T, Alloc >::rows(), octave_sort< T >::sort(), Sparse< T, Alloc >::xcidx(), Sparse< T, Alloc >::xdata(), and Sparse< T, Alloc >::xridx().

solve_singularity_warning()

static void solve_singularity_warning ( double  )

static

Definition at line 358 of file sparse-dmsolve.cc.

Referenced by dmsolve().