16 #ifndef dealii_lapack_full_matrix_h 17 #define dealii_lapack_full_matrix_h 20 #include <deal.II/base/config.h> 22 #include <deal.II/base/smartpointer.h> 23 #include <deal.II/base/table.h> 24 #include <deal.II/base/thread_management.h> 26 #include <deal.II/lac/lapack_support.h> 27 #include <deal.II/lac/vector_memory.h> 33 DEAL_II_NAMESPACE_OPEN
36 template <
typename number>
38 template <
typename number>
40 template <
typename number>
42 template <
typename number>
58 template <
typename number>
65 using size_type = std::make_unsigned<types::blas_int>::type;
109 template <
typename number2>
119 template <
typename number2>
184 template <
typename MatrixType>
286 template <
typename MatrixType>
288 fill(
const MatrixType &src,
293 const number factor = 1.,
324 template <
typename number2>
326 vmult(Vector<number2> & w,
327 const Vector<number2> &v,
328 const bool adding =
false)
const;
336 const bool adding =
false)
const;
344 template <
typename number2>
346 vmult_add(Vector<number2> &w,
const Vector<number2> &v)
const;
365 template <
typename number2>
367 Tvmult(Vector<number2> & w,
368 const Vector<number2> &v,
369 const bool adding =
false)
const;
377 const bool adding =
false)
const;
385 template <
typename number2>
387 Tvmult_add(Vector<number2> &w,
const Vector<number2> &v)
const;
413 const bool adding =
false)
const;
422 const bool adding =
false)
const;
441 const bool adding =
false)
const;
450 const bool adding =
false)
const;
471 const bool adding =
false)
const;
490 const bool adding =
false)
const;
499 const bool adding =
false)
const;
519 const bool adding =
false)
const;
528 const bool adding =
false)
const;
672 const bool transposed)
const;
694 const bool left_eigenvectors =
false);
717 const number upper_bound,
718 const number abs_accuracy,
751 const number lower_bound,
752 const number upper_bound,
753 const number abs_accuracy,
860 const unsigned int precision = 3,
861 const bool scientific =
true,
862 const unsigned int width = 0,
863 const char * zero_string =
" ",
864 const double denominator = 1.,
865 const double threshold = 0.)
const;
872 norm(
const char type)
const;
889 mutable std::vector<number>
work;
894 mutable std::vector<types::blas_int>
iwork;
902 std::vector<types::blas_int>
ipiv;
913 std::vector<typename numbers::NumberTraits<number>::real_type>
wr;
919 std::vector<number>
wi;
924 std::vector<number>
vl;
929 std::vector<number>
vr;
935 std::unique_ptr<LAPACKFullMatrix<number>>
svd_u;
941 std::unique_ptr<LAPACKFullMatrix<number>>
svd_vt;
957 template <
typename number>
981 template <
typename number>
987 (*this)(i, j) = value;
991 template <
typename number>
998 template <
typename number>
1005 template <
typename number>
1006 template <
typename MatrixType>
1010 this->
reinit(M.m(), M.n());
1015 for (
size_type row = 0; row < M.m(); ++row)
1017 const typename MatrixType::const_iterator end_row = M.end(row);
1018 for (
typename MatrixType::const_iterator entry = M.begin(row);
1021 this->
el(row, entry->column()) = entry->value();
1029 template <
typename number>
1030 template <
typename MatrixType>
1037 const number factor,
1042 for (
size_type row = src_offset_i; row < M.m(); ++row)
1044 const typename MatrixType::const_iterator end_row = M.end(row);
1045 for (
typename MatrixType::const_iterator entry = M.begin(row);
1049 const size_type i = transpose ? entry->column() : row;
1050 const size_type j = transpose ? row : entry->column();
1052 const size_type dst_i = dst_offset_i + i - src_offset_i;
1053 const size_type dst_j = dst_offset_j + j - src_offset_j;
1055 (*this)(dst_i, dst_j) = factor * entry->value();
1063 template <
typename number>
1064 template <
typename number2>
1067 const Vector<number2> &,
1071 ExcMessage(
"LAPACKFullMatrix<number>::vmult must be called with a " 1072 "matching Vector<double> vector type."));
1076 template <
typename number>
1077 template <
typename number2>
1080 const Vector<number2> &)
const 1083 ExcMessage(
"LAPACKFullMatrix<number>::vmult_add must be called with a " 1084 "matching Vector<double> vector type."));
1088 template <
typename number>
1089 template <
typename number2>
1092 const Vector<number2> &,
1096 ExcMessage(
"LAPACKFullMatrix<number>::Tvmult must be called with a " 1097 "matching Vector<double> vector type."));
1101 template <
typename number>
1102 template <
typename number2>
1105 const Vector<number2> &)
const 1109 "LAPACKFullMatrix<number>::Tvmult_add must be called with a " 1110 "matching Vector<double> vector type."));
1114 template <
typename number>
1115 inline std::complex<number>
1124 return std::complex<number>(
wi[i]);
1126 return std::complex<number>(
wr[i],
wi[i]);
1130 template <
typename number>
1143 DEAL_II_NAMESPACE_CLOSE
LAPACKFullMatrix(const size_type size=0)
void Tmmult(LAPACKFullMatrix< number > &C, const LAPACKFullMatrix< number > &B, const bool adding=false) const
std::vector< number > work
void rank1_update(const number a, const Vector< number > &v)
number linfty_norm() const
void mTmult(LAPACKFullMatrix< number > &C, const LAPACKFullMatrix< number > &B, const bool adding=false) const
std::unique_ptr< LAPACKFullMatrix< number > > svd_vt
LAPACKFullMatrix< number > & operator/=(const number factor)
number reciprocal_condition_number() const
Contents is actually a matrix.
std::array< std::pair< Number, Tensor< 1, dim, Number > >, std::integral_constant< int, dim >::value > eigenvectors(const SymmetricTensor< 2, dim, Number > &T, const SymmetricTensorEigenvectorMethod method=SymmetricTensorEigenvectorMethod::ql_implicit_shifts)
number singular_value(const size_type i) const
std::array< Number, 1 > eigenvalues(const SymmetricTensor< 2, 1, Number > &T)
LAPACKSupport::State state
typename TableBase< 2, T >::size_type size_type
void remove_row_and_column(const size_type row, const size_type col)
std::vector< types::blas_int > ipiv
#define AssertIndexRange(index, range)
number norm(const char type) const
void vmult(Vector< number2 > &w, const Vector< number2 > &v, const bool adding=false) const
void solve(Vector< number > &v, const bool transposed=false) const
void vmult_add(Vector< number2 > &w, const Vector< number2 > &v) const
void compute_eigenvalues_symmetric(const number lower_bound, const number upper_bound, const number abs_accuracy, Vector< number > &eigenvalues, FullMatrix< number > &eigenvectors)
void reinit(const size_type size)
void Tvmult_add(Vector< number2 > &w, const Vector< number2 > &v) const
void print_formatted(std::ostream &out, const unsigned int precision=3, const bool scientific=true, const unsigned int width=0, const char *zero_string=" ", const double denominator=1., const double threshold=0.) const
void apply_lu_factorization(Vector< number > &v, const bool transposed) const
const TableIndices< N > & size() const
std::vector< types::blas_int > iwork
static::ExceptionBase & ExcState(State arg1)
void set(const size_type i, const size_type j, const number value)
static::ExceptionBase & ExcInvalidState()
static::ExceptionBase & ExcMessage(std::string arg1)
#define Assert(cond, exc)
void compute_inverse_svd_with_kernel(const unsigned int kernel_size)
void copy_from(const MatrixType &)
std::make_unsigned< types::blas_int >::type size_type
void TmTmult(LAPACKFullMatrix< number > &C, const LAPACKFullMatrix< number > &B, const bool adding=false) const
number frobenius_norm() const
void add(const number a, const LAPACKFullMatrix< number > &B)
std::vector< number > inv_work
Matrix is the inverse of a singular value decomposition.
void compute_lu_factorization()
SymmetricTensor< rank_, dim, Number > transpose(const SymmetricTensor< rank_, dim, Number > &t)
std::unique_ptr< LAPACKFullMatrix< number > > svd_u
void grow_or_shrink(const size_type size)
void scale_rows(const Vector< number > &V)
LAPACKFullMatrix< number > & operator*=(const number factor)
void compute_inverse_svd(const double threshold=0.)
void compute_cholesky_factorization()
void Tvmult(Vector< number2 > &w, const Vector< number2 > &v, const bool adding=false) const
LAPACKFullMatrix< number > & operator=(const LAPACKFullMatrix< number > &)
void fill(const MatrixType &src, const size_type dst_offset_i=0, const size_type dst_offset_j=0, const size_type src_offset_i=0, const size_type src_offset_j=0, const number factor=1., const bool transpose=false)
std::vector< typename numbers::NumberTraits< number >::real_type > wr
Matrix contains singular value decomposition,.
void set_property(const LAPACKSupport::Property property)
reference el(const size_type i, const size_type j)
void compute_generalized_eigenvalues_symmetric(LAPACKFullMatrix< number > &B, const number lower_bound, const number upper_bound, const number abs_accuracy, Vector< number > &eigenvalues, std::vector< Vector< number >> &eigenvectors, const types::blas_int itype=1)
void compute_eigenvalues(const bool right_eigenvectors=false, const bool left_eigenvectors=false)
number determinant() const
std::complex< number > eigenvalue(const size_type i) const
void mmult(LAPACKFullMatrix< number > &C, const LAPACKFullMatrix< number > &B, const bool adding=false) const
Eigenvalue vector is filled.
LAPACKSupport::Property property
static::ExceptionBase & ExcInternalError()