17 #ifndef __deal2__sparse_direct_h
18 #define __deal2__sparse_direct_h
22 #include <deal.II/base/config.h>
24 #include <deal.II/base/subscriptor.h>
25 #include <deal.II/base/thread_management.h>
26 #include <deal.II/lac/vector.h>
27 #include <deal.II/lac/sparse_matrix.h>
28 #include <deal.II/lac/sparse_matrix_ez.h>
29 #include <deal.II/lac/block_sparse_matrix.h>
31 #ifdef DEAL_II_WITH_MUMPS
32 # include <deal.II/base/utilities.h>
33 # include <dmumps_c.h>
141 template <
class Matrix>
147 template <
class Matrix>
249 template <
class Matrix>
250 void solve (
const Matrix &matrix,
252 bool transpose =
false);
257 template <
class Matrix>
258 void solve (
const Matrix &matrix,
260 bool transpose =
false);
272 <<
"UMFPACK routine " << arg1
273 <<
" returned error status " << arg2
274 <<
". See the file <bundled/umfpack/UMFPACK/Include/umfpack.h>"
275 <<
" for a description of 'status codes'.");
285 void *numeric_decomposition;
298 template <
typename number>
301 template <
typename number>
304 template <
typename number>
310 std::vector<long int>
Ap;
311 std::vector<long int> Ai;
312 std::vector<double> Ax;
344 #ifdef DEAL_II_WITH_MUMPS
346 #endif // DEAL_II_WITH_MUMPS
349 std::vector<double> rhs;
359 template<
class Matrix>
404 template <
class Matrix>
412 template <
class Matrix>
430 DEAL_II_NAMESPACE_CLOSE
432 #endif // __deal2__sparse_direct_h
void * symbolic_decomposition
DeclException0(ExcInitializeAlreadyCalled)
types::global_dof_index size_type
void Tvmult_add(Vector< double > &dst, const Vector< double > &src) const
void Tvmult(Vector< double > &dst, const Vector< double > &src) const
void factorize(const Matrix &matrix)
void copy_solution(Vector< double > &vector)
void initialize_matrix(const Matrix &matrix)
std::vector< long int > Ap
unsigned int global_dof_index
void solve(Vector< double > &vector)
void solve(Vector< double > &rhs_and_solution, bool transpose=false) const
void vmult_add(Vector< double > &dst, const Vector< double > &src) const
std::vector< double > control
DeclException2(ExcUMFPACKError, char *, int,<< "UMFPACK routine "<< arg1<< " returned error status "<< arg2<< ". See the file <bundled/umfpack/UMFPACK/Include/umfpack.h>"<< " for a description of 'status codes'.")
void initialize(const SparsityPattern &sparsity_pattern)
void sort_arrays(const SparseMatrixEZ< number > &)
void vmult(Vector< double > &dst, const Vector< double > &src)
void initialize(const Matrix &matrix, const Vector< double > &vector)
void vmult(Vector< double > &dst, const Vector< double > &src) const
types::global_dof_index size_type