#include <schur.hh>
Inheritance diagram for SchurComplement:
Collaboration diagram for SchurComplement:
Public Member Functions | |
| SchurComplement (IS ia, Distribution *ds) | |
| SchurComplement (SchurComplement &other) | |
| LinSys * | get_system () const |
| Distribution * | get_distribution () const |
| ~SchurComplement () | |
| void | form_rhs () |
| void | set_complement (LinSys_PETSC *ls) |
| Set complement LinSys object. More... | |
| Distribution * | make_complement_distribution () |
| get distribution of complement object if complement is defined More... | |
| double | get_solution_precision () override |
| get precision of solving More... | |
| int | solve () override |
 Public Member Functions inherited from LinSys_PETSC | |
| LinSys_PETSC (const Distribution *rows_ds) | |
| LinSys_PETSC (LinSys_PETSC &other) | |
| const Distribution * | get_ds () |
| const Mat & | get_matrix () |
| const Vec & | get_rhs () |
| PetscErrorCode | set_matrix (Mat &matrix, MatStructure str) |
| PetscErrorCode | set_rhs (Vec &rhs) |
| PetscErrorCode | mat_zero_entries () |
| PetscErrorCode | rhs_zero_entries () |
| void | start_allocation () |
| void | start_add_assembly () |
| void | start_insert_assembly () |
| void | mat_set_values (int nrow, int *rows, int ncol, int *cols, double *vals) |
| void | rhs_set_values (int nrow, int *rows, double *vals) |
| void | preallocate_values (int nrow, int *rows, int ncol, int *cols) |
| void | preallocate_matrix () |
| void | finish_assembly () |
| void | finish_assembly (MatAssemblyType assembly_type) |
| void | apply_constrains (double scalar=1.) |
| void | set_initial_guess_nonzero (bool set_nonzero=true) |
| double | get_absolute_accuracy () |
| void | get_whole_solution (std::vector< double > &globalSolution) |
| void | view () |
| void | set_from_input (const Input::Record in_rec) |
| ~LinSys_PETSC () | |
| Public Member Functions inherited from LinSys | |
| LinSys (const Distribution *rows_ds) | |
| LinSys (LinSys &other) | |
| unsigned int | size () |
| unsigned int | vec_lsize () |
| void | set_matrix_changed () |
| void | set_rhs_changed () |
| bool | is_matrix_changed () |
| bool | is_rhs_changed () |
| const Vec & | get_solution () |
| void | set_solution (double *sol_array) |
| double * | get_solution_array () |
| virtual void | set_whole_solution (std::vector< double > &globalSolution) |
| void | mat_set_value (int row, int col, double val) |
| void | rhs_set_value (int row, double val) |
| void | set_values (int nrow, int *rows, int ncol, int *cols, PetscScalar *mat_vals, PetscScalar *rhs_vals) |
| Set values in the system matrix and values in the right-hand side vector on corresponding rows. More... | |
| void | set_values (std::vector< int > &row_dofs, std::vector< int > &col_dofs, const arma::mat &matrix, const arma::vec &rhs, const arma::vec &row_solution, const arma::vec &col_solution) |
| void | add_constraint (int row, double value) |
| double | get_residual_norm () |
| double | get_relative_accuracy () |
| void | set_symmetric (bool flag=true) |
| bool | is_symmetric () |
| void | set_positive_definite (bool flag=true) |
| void | set_negative_definite (bool flag=true) |
| bool | is_positive_definite () |
| bool | is_negative_definite () |
| bool | is_new () |
| bool | is_preallocated () |
| void | set_spd_via_symmetric_general (bool flag=true) |
| bool | is_spd_via_symmetric_general () |
| virtual | ~LinSys () |
| LinSys (unsigned int lsize, double *sol_array=NULL) | |
| unsigned int | size () |
| unsigned int | vec_lsize () |
| const Distribution & | ds () |
| const Mat & | get_matrix () |
| const Vec & | get_rhs () |
| const Vec & | get_solution () |
| double * | get_solution_array () |
| void | start_add_assembly () |
| void | start_insert_assembly () |
| void | partial_assembly () |
| void | finish_assembly () |
| virtual void | mat_set_values (int nrow, int *rows, int ncol, int *cols, PetscScalar *vals)=0 |
| void | mat_set_value (int row, int col, PetscScalar val) |
| void | rhs_set_values (int nrow, int *rows, PetscScalar *vals) |
| void | rhs_set_value (int row, PetscScalar val) |
| void | set_values (int nrow, int *rows, int ncol, int *cols, PetscScalar *mat_vals, PetscScalar *rhs_vals, std::vector< bool > &constrains_row_mask=std::vector(0), double *constrain_values=NULL) |
| void | add_constrain (int row, double value) |
| void | apply_constrains (std::vector< Constrain > &constraints) |
| void | set_symmetric (bool flag=true) |
| bool | is_symmetric () |
| void | set_positive_definite (bool flag=true) |
| bool | is_positive_definite () |
| void | view (std::ostream output_stream, int *output_mapping=NULL) |
| virtual | ~LinSys () |
Protected Member Functions | |
| void | create_inversion_matrix () |
| create IA matrix More... | |
| void | form_schur () |
| void | resolve () |
Protected Attributes | |
| Mat | IA |
| Mat | B |
| Mat | Bt |
| Mat | xA |
| Mat | IAB |
| int | loc_size_A |
| int | loc_size_B |
| IS | IsA |
| IS | IsB |
| Vec | RHS1 |
| Vec | RHS2 |
| Vec | Sol1 |
| Vec | Sol2 |
| SchurState | state |
| int | orig_lsize |
| Size of local vector part of original system. More... | |
| LinSys_PETSC * | Compl |
| Distribution * | ds_ |
| Protected Attributes inherited from LinSys_PETSC | |
| std::string | params_ |
| command-line-like options for the PETSc solver More... | |
| bool | init_guess_nonzero |
| flag for starting from nonzero guess More... | |
| Mat | matrix_ |
| Petsc matrix of the problem. More... | |
| Vec | rhs_ |
| PETSc vector constructed with vx array. More... | |
| double * | v_rhs_ |
| local RHS array pointing to Vec rhs_ More... | |
| Vec | on_vec_ |
| Vectors for counting non-zero entries in diagonal block. More... | |
| Vec | off_vec_ |
| Vectors for counting non-zero entries in off-diagonal block. More... | |
| double | solution_precision_ |
| Protected Attributes inherited from LinSys | |
| double | r_tol_ |
| double | a_tol_ |
| int | max_it_ |
| MPI_Comm | comm_ |
| SetValuesMode | status_ |
| Set value status of the linear system. More... | |
| const unsigned | lsize_ |
| local number of matrix rows (non-overlapping division of rows) More... | |
| unsigned | size_ |
| global number of matrix rows, i.e. problem size More... | |
| const Distribution * | rows_ds_ |
| final distribution of rows of MH matrix More... | |
| bool | symmetric_ |
| bool | positive_definite_ |
| bool | negative_definite_ |
| bool | spd_via_symmetric_general_ |
| bool | matrix_changed_ |
| true if the matrix was changed since the last solve More... | |
| bool | rhs_changed_ |
| true if the right hand side was changed since the last solve More... | |
| Vec | solution_ |
| PETSc vector constructed with vb array. More... | |
| double * | v_solution_ |
| local solution array pointing into Vec solution_ More... | |
| bool | own_solution_ |
| Indicates if the solution array has been allocated by this class. More... | |
| double | residual_norm_ |
| local solution array pointing into Vec solution_ More... | |
| ConstraintVec_ | constraints_ |
| std::vector< double > | globalSolution_ |
| global solution in numbering for linear system More... | |
| Input::Record | in_rec_ |
| Distribution | vec_ds |
| Distribution of continuous blocks of system rows among the processors. More... | |
| bool | symmetric |
| Flag for the symmetric system. More... | |
| bool | positive_definite |
| Flag for positive definite system. More... | |
| bool | own_solution |
| Indicates if the solution array has been allocated by this class. More... | |
| SetValuesMode | status |
| Set value status of the linear system. More... | |
| Mat | matrix |
| Petsc matrix of the problem. More... | |
| Vec | rhs |
| PETSc vector constructed with vx array. More... | |
| Vec | solution |
| PETSc vector constructed with vb array. More... | |
| double * | v_rhs |
| RHS vector. More... | |
| double * | v_solution |
| Vector of solution. More... | |
| int * | subdomain_indices |
| Remember indices which created mapping. More... | |
| Mat | local_matrix |
| local matrix of subdomain (used in LinSys_MATIS) More... | |
Additional Inherited Members | |
| Public Types inherited from LinSys | |
| enum | SetValuesMode { INSERT =INSERT_VALUES, ADD =ADD_VALUES, ALLOCATE, DONE, NONE, INSERT =INSERT_VALUES, ADD =ADD_VALUES, ALLOCATE, DONE, NONE } |
| enum | SetValuesMode { INSERT =INSERT_VALUES, ADD =ADD_VALUES, ALLOCATE, DONE, NONE, INSERT =INSERT_VALUES, ADD =ADD_VALUES, ALLOCATE, DONE, NONE } |
| enum | LinSysType { PETSC_MPIAIJ_preallocate_by_assembly, PETSC_MPIAIJ_assembly_by_triples, BDDC, PETSC_schur_complement } |
| Public Attributes inherited from LinSys | |
| LinSysType | type |
| Particular type of the linear system. More... | |
| Static Public Attributes inherited from LinSys_PETSC | |
| static Input::Type::Record | input_type |
| Protected Types inherited from LinSys | |
| typedef std::pair< unsigned, double > | Constraint_ |
| typedef std::vector< Constraint_ > | ConstraintVec_ |
Detailed Description
Constructor & Destructor Documentation
| SchurComplement::SchurComplement | ( | IS | ia, |
| Distribution * | ds | ||
| ) |
Constructor
Gets linear system with original matrix A and creates its inversion (IA matrix)
In current implementation the index set IsA has to be continuous sequence at the beginning of the local block of indices.
Create Schur complement system.
- Parameters
-
[in] orig : original system [in] inv_a : inversion of the A block [in] ia : index set of the A block, default continuous given by inv_a: proc 1 2 3
Orig: ****** ****** **** IA : *** ** ***
| SchurComplement::SchurComplement | ( | SchurComplement & | other) |
| SchurComplement::~SchurComplement | ( | ) |
Destructor. In particular it also delete complement linear system if it was passed in through the set_complement() method.
SCHUR COMPLEMENT destructor
Member Function Documentation
|
protected |
| void SchurComplement::form_rhs | ( | ) |
Compute only right hand side. This is useful when you change only rhs of the original system. TODO: We should ask original system if the matrix has changed (using LazyDependency) and possibly call only form_rhs, then this can be protected
Definition at line 246 of file schur.cc.
Here is the caller graph for this function:
|
protected |
COMPUTE A SCHUR COMPLEMENT OF A PETSC MATRIX
given symmetric original matrix Orig has form A B x_1 RHS_1 B' C * x_2 = RHS_2 where the first block is given by index set IsA, and the second block by IsB user has to provide inverse IA of the A-block we suppose that original matrix have non-zero pattern for the schur complement
we return: Shur - schur complement, ShurRHS - RHS of the complemented system: (B' * IA * B - C) * x_2 = (B' * IA * RHS_1 - RHS_2) IAB - a matrix to compute eliminated part of the solution: x_1 = IA * RHS_1 - IAB * x_2
Actually as B' is stored separetly, the routine can be used also for nonsymetric original system
Definition at line 141 of file schur.cc.
Here is the caller graph for this function:
|
inline |
Returns distribution of the original system (solved by class SchurComplement).
Definition at line 97 of file schur.hh.
Here is the caller graph for this function:
|
overridevirtual |
|
inline |
| Distribution * SchurComplement::make_complement_distribution | ( | ) |
|
protected |
| void SchurComplement::set_complement | ( | LinSys_PETSC * | ls) |
|
overridevirtual |
Member Data Documentation
|
protected |
|
protected |
|
protected |
|
protected |
The documentation for this class was generated from the following files:
