16 #ifndef dealii__precondition_block_base_h 17 #define dealii__precondition_block_base_h 20 #include <deal.II/base/config.h> 21 #include <deal.II/base/exceptions.h> 22 #include <deal.II/base/subscriptor.h> 23 #include <deal.II/base/smartpointer.h> 24 #include <deal.II/base/memory_consumption.h> 25 #include <deal.II/lac/householder.h> 26 #include <deal.II/lac/lapack_full_matrix.h> 30 DEAL_II_NAMESPACE_OPEN
33 template <
typename number>
class Vector;
53 template <
typename number>
104 void reinit(
unsigned int nblocks, size_type blocksize,
bool compress,
145 unsigned int size()
const;
151 number
el(size_type i, size_type j)
const;
156 template <
typename number2>
162 template <
typename number2>
294 template <
typename number>
307 template <
typename number>
325 template <
typename number>
373 std::vector<FullMatrix<number> >
382 std::vector<FullMatrix<number> >
392 std::vector<LAPACKFullMatrix<number> >
404 template <
typename number>
412 template <
typename number>
413 template <
typename number2>
441 template <
typename number>
442 template <
typename number2>
470 template <
typename number>
483 template <
typename number>
496 template <
typename number>
509 template <
typename number>
524 template <
typename number>
539 template <
typename number>
554 template <
typename number>
569 template <
typename number>
584 template <
typename number>
592 template <
typename number>
600 template <
typename number>
608 template <
typename number>
616 template <
typename number>
620 deallog <<
"PreconditionBlockBase: " <<
size() <<
" blocks; ";
624 unsigned int kermin = 100000000, kermax = 0;
625 double sigmin = 1.e300, sigmax= -1.e300;
626 double kappamin = 1.e300, kappamax= -1.e300;
636 const double co = sm/s0;
638 if (kermin > k) kermin = k-1;
639 if (kermax < k) kermax = k-1;
640 if (s0 < sigmin) sigmin = s0;
641 if (sm > sigmax) sigmax = sm;
642 if (co < kappamin) kappamin = co;
643 if (co > kappamax) kappamax = co;
645 deallog <<
"dim ker [" << kermin <<
':' << kermax
646 <<
"] sigma [" << sigmin <<
':' << sigmax
647 <<
"] kappa [" << kappamin <<
':' << kappamax <<
']' << std::endl;
663 template <
typename number>
668 std::size_t mem =
sizeof(*this);
677 DEAL_II_NAMESPACE_CLOSE
void inverse_Tvmult(size_type i, Vector< number2 > &dst, const Vector< number2 > &src) const
number singular_value(const size_type i) const
bool store_diagonals() const
types::global_dof_index size_type
#define AssertIndexRange(index, range)
unsigned int size() const
static::ExceptionBase & ExcIndexRange(int arg1, int arg2, int arg3)
void inverses_computed(bool are_they)
std::vector< FullMatrix< number > > var_diagonal
FullMatrix< number > & inverse(size_type i)
PreconditionBlockBase(bool store_diagonals=false, Inversion method=gauss_jordan)
unsigned int global_dof_index
#define Assert(cond, exc)
number el(size_type i, size_type j) const
#define DeclException0(Exception0)
bool same_diagonal() const
std::vector< FullMatrix< number > > var_inverse_full
void reinit(unsigned int nblocks, size_type blocksize, bool compress, Inversion method=gauss_jordan)
void log_statistics() const
unsigned int n_diagonal_blocks
std_cxx11::enable_if< std_cxx11::is_fundamental< T >::value, std::size_t >::type memory_consumption(const T &t)
std::size_t memory_consumption() const
static::ExceptionBase & ExcDiagonalsNotStored()
std::vector< LAPACKFullMatrix< number > > var_inverse_svd
LAPACKFullMatrix< number > & inverse_svd(size_type i)
Householder< number > & inverse_householder(size_type i)
static::ExceptionBase & ExcNotImplemented()
std::vector< Householder< number > > var_inverse_householder
bool inverses_ready() const
FullMatrix< number > & diagonal(size_type i)
unsigned int n_cols() const
static::ExceptionBase & ExcInverseNotAvailable()
void inverse_vmult(size_type i, Vector< number2 > &dst, const Vector< number2 > &src) const