|
Reference documentation for deal.II version 9.2.0
|
\(\newcommand{\dealvcentcolon}{\mathrel{\mathop{:}}}\)
\(\newcommand{\dealcoloneq}{\dealvcentcolon\mathrel{\mkern-1.2mu}=}\)
\(\newcommand{\jump}[1]{\left[\!\left[ #1 \right]\!\right]}\)
\(\newcommand{\average}[1]{\left\{\!\left\{ #1 \right\}\!\right\}}\)
Go to the documentation of this file.
16 #ifndef dealii_scalapack_h
17 #define dealii_scalapack_h
21 #ifdef DEAL_II_WITH_SCALAPACK
76 template <
typename NumberType>
96 const std::shared_ptr<const Utilities::MPI::ProcessGrid> &process_grid,
112 const std::shared_ptr<const Utilities::MPI::ProcessGrid> &process_grid,
130 const std::string & filename,
131 const std::shared_ptr<const Utilities::MPI::ProcessGrid> &process_grid,
152 const std::shared_ptr<const Utilities::MPI::ProcessGrid> &process_grid,
166 const std::shared_ptr<const Utilities::MPI::ProcessGrid> &process_grid,
209 const unsigned int rank);
266 const std::pair<unsigned int, unsigned int> &offset_A,
267 const std::pair<unsigned int, unsigned int> &offset_B,
268 const std::pair<unsigned int, unsigned int> &submatrix_size)
const;
294 const NumberType a = 0.,
295 const NumberType
b = 1.,
296 const bool transpose_B =
false);
343 mult(
const NumberType
b,
347 const bool transpose_A =
false,
348 const bool transpose_B =
false)
const;
370 const bool adding =
false)
const;
392 const bool adding =
false)
const;
414 const bool adding =
false)
const;
437 const bool adding =
false)
const;
460 save(
const std::string & filename,
461 const std::pair<unsigned int, unsigned int> &
chunk_size =
477 load(
const std::string &filename);
525 std::vector<NumberType>
527 const std::pair<unsigned int, unsigned int> &index_limits,
528 const bool compute_eigenvectors);
538 std::vector<NumberType>
540 const std::pair<NumberType, NumberType> &value_limits,
541 const bool compute_eigenvectors);
559 std::vector<NumberType>
561 const std::pair<unsigned int, unsigned int> &index_limits,
562 const bool compute_eigenvectors);
574 std::vector<NumberType>
576 const std::pair<NumberType, NumberType> &value_limits,
577 const bool compute_eigenvectors);
607 std::vector<NumberType>
755 local_el(
const unsigned int loc_row,
const unsigned int loc_column)
const;
761 local_el(
const unsigned int loc_row,
const unsigned int loc_column);
773 template <
class InputVector>
787 template <
class InputVector>
816 std::vector<NumberType>
818 const bool compute_eigenvectors,
819 const std::pair<unsigned int, unsigned int> &index_limits =
822 const std::pair<NumberType, NumberType> &value_limits =
823 std::make_pair(std::numeric_limits<NumberType>::quiet_NaN(),
824 std::numeric_limits<NumberType>::quiet_NaN()));
845 std::vector<NumberType>
847 const bool compute_eigenvectors,
848 const std::pair<unsigned int, unsigned int> &index_limits =
851 const std::pair<NumberType, NumberType> &value_limits =
852 std::make_pair(std::numeric_limits<NumberType>::quiet_NaN(),
853 std::numeric_limits<NumberType>::quiet_NaN()));
861 const std::pair<unsigned int, unsigned int> &
chunk_size)
const;
876 const std::pair<unsigned int, unsigned int> &
chunk_size)
const;
902 std::shared_ptr<const Utilities::MPI::ProcessGrid>
grid;
942 mutable std::vector<NumberType>
work;
995 template <
typename NumberType>
998 const unsigned int loc_column)
const
1000 return (*
this)(loc_row, loc_column);
1005 template <
typename NumberType>
1008 const unsigned int loc_column)
1010 return (*
this)(loc_row, loc_column);
1014 template <
typename NumberType>
1023 template <
typename NumberType>
1032 template <
typename NumberType>
1036 return n_local_rows;
1041 template <
typename NumberType>
1045 return n_local_columns;
1053 #endif // DEAL_II_WITH_SCALAPACK
NumberType reciprocal_condition_number(const NumberType a_norm) const
std::vector< NumberType > eigenpairs_symmetric_by_index(const std::pair< unsigned int, unsigned int > &index_limits, const bool compute_eigenvectors)
void compute_lu_factorization()
~ScaLAPACKMatrix() override=default
void scale_rows(const InputVector &factors)
NumberType l1_norm() const
void compute_cholesky_factorization()
ScaLAPACKMatrix< NumberType > & operator=(const FullMatrix< NumberType > &)
std::shared_ptr< const Utilities::MPI::ProcessGrid > grid
LAPACKSupport::State state
const int first_process_row
void copy_from(const LAPACKFullMatrix< NumberType > &matrix, const unsigned int rank)
std::vector< NumberType > eigenpairs_symmetric_by_value(const std::pair< NumberType, NumberType > &value_limits, const bool compute_eigenvectors)
NumberType norm_general(const char type) const
ScaLAPACKMatrix(const size_type n_rows, const size_type n_columns, const std::shared_ptr< const Utilities::MPI::ProcessGrid > &process_grid, const size_type row_block_size=32, const size_type column_block_size=32, const LAPACKSupport::Property property=LAPACKSupport::Property::general)
void save_parallel(const std::string &filename, const std::pair< unsigned int, unsigned int > &chunk_size) const
void mult(const NumberType b, const ScaLAPACKMatrix< NumberType > &B, const NumberType c, ScaLAPACKMatrix< NumberType > &C, const bool transpose_A=false, const bool transpose_B=false) const
void save_serial(const std::string &filename, const std::pair< unsigned int, unsigned int > &chunk_size) const
void copy_to(FullMatrix< NumberType > &matrix) const
SymmetricTensor< 2, dim, Number > C(const Tensor< 2, dim, Number > &F)
void reinit(const size_type n_rows, const size_type n_columns, const std::shared_ptr< const Utilities::MPI::ProcessGrid > &process_grid, const size_type row_block_size=32, const size_type column_block_size=32, const LAPACKSupport::Property property=LAPACKSupport::Property::general)
NumberType norm_symmetric(const char type) const
std::vector< NumberType > eigenpairs_symmetric_by_value_MRRR(const std::pair< NumberType, NumberType > &value_limits, const bool compute_eigenvectors)
void TmTmult(ScaLAPACKMatrix< NumberType > &C, const ScaLAPACKMatrix< NumberType > &B, const bool adding=false) const
void copy_transposed(const ScaLAPACKMatrix< NumberType > &B)
void least_squares(ScaLAPACKMatrix< NumberType > &B, const bool transpose=false)
@ symmetric
Matrix is symmetric.
void save(const std::string &filename, const std::pair< unsigned int, unsigned int > &chunk_size=std::make_pair(numbers::invalid_unsigned_int, numbers::invalid_unsigned_int)) const
void set_property(const LAPACKSupport::Property property)
void mmult(ScaLAPACKMatrix< NumberType > &C, const ScaLAPACKMatrix< NumberType > &B, const bool adding=false) const
LAPACKSupport::Property property
NumberType local_el(const unsigned int loc_row, const unsigned int loc_column) const
unsigned int global_row(const unsigned int loc_row) const
#define DEAL_II_NAMESPACE_OPEN
SymmetricTensor< 2, dim, Number > b(const Tensor< 2, dim, Number > &F)
void load_serial(const std::string &filename)
@ matrix
Contents is actually a matrix.
void mTmult(ScaLAPACKMatrix< NumberType > &C, const ScaLAPACKMatrix< NumberType > &B, const bool adding=false) const
unsigned int global_column(const unsigned int loc_column) const
@ general
No special properties.
std::vector< NumberType > eigenpairs_symmetric_MRRR(const bool compute_eigenvectors, const std::pair< unsigned int, unsigned int > &index_limits=std::make_pair(numbers::invalid_unsigned_int, numbers::invalid_unsigned_int), const std::pair< NumberType, NumberType > &value_limits=std::make_pair(std::numeric_limits< NumberType >::quiet_NaN(), std::numeric_limits< NumberType >::quiet_NaN()))
LAPACKSupport::State get_state() const
void scale_columns(const InputVector &factors)
void load(const std::string &filename)
unsigned int local_m() const
void Tmmult(ScaLAPACKMatrix< NumberType > &C, const ScaLAPACKMatrix< NumberType > &B, const bool adding=false) const
void add(const ScaLAPACKMatrix< NumberType > &B, const NumberType a=0., const NumberType b=1., const bool transpose_B=false)
unsigned int pseudoinverse(const NumberType ratio)
static const unsigned int invalid_unsigned_int
unsigned int local_n() const
std::vector< NumberType > eigenpairs_symmetric(const bool compute_eigenvectors, const std::pair< unsigned int, unsigned int > &index_limits=std::make_pair(numbers::invalid_unsigned_int, numbers::invalid_unsigned_int), const std::pair< NumberType, NumberType > &value_limits=std::make_pair(std::numeric_limits< NumberType >::quiet_NaN(), std::numeric_limits< NumberType >::quiet_NaN()))
const int first_process_column
void Tadd(const NumberType b, const ScaLAPACKMatrix< NumberType > &B)
NumberType frobenius_norm() const
std::vector< NumberType > compute_SVD(ScaLAPACKMatrix< NumberType > *U=nullptr, ScaLAPACKMatrix< NumberType > *VT=nullptr)
std::vector< NumberType > work
NumberType linfty_norm() const
std::vector< NumberType > eigenpairs_symmetric_by_index_MRRR(const std::pair< unsigned int, unsigned int > &index_limits, const bool compute_eigenvectors)
const TableIndices< N > & size() const
#define DEAL_II_NAMESPACE_CLOSE
const int submatrix_column
void load_parallel(const std::string &filename)
LAPACKSupport::Property get_property() const