petsc-3.14.6 2021-03-30
Report Typos and Errors

PCDEFLATION

Deflation preconditioner shifts (deflates) part of the spectrum to zero or to a predefined value.

Options Database Keys

-pc_deflation_init_only <false> - if true computes only the special guess
-pc_deflation_max_lvl <0> - maximum number of levels for multilevel deflation
-pc_deflation_reduction_factor <-1> - reduction factor on bottom level coarse problem for PCTELESCOPE (default based on the size of the coarse problem)
-pc_deflation_correction <false> - if true apply coarse problem correction
-pc_deflation_correction_factor <1.0> - sets coarse problem correction factor
-pc_deflation_compute_space <haar> - compute PCDeflationSpaceType deflation space
-pc_deflation_compute_space_size <1> - size of the deflation space (corresponds to number of levels for wavelet-based deflation)

Notes

Given a (complex - transpose is always Hermitian) full rank deflation matrix W, the deflation (introduced in [1,2]) preconditioner uses projections Q = W*(W'*A*W)^{-1}*W' and P = I - Q*A, where A is pmat.

The deflation computes initial guess x0 = x_{-1} - Q*r_{-1}, which is the solution on the deflation space. If PCDeflationSetInitOnly() or -pc_deflation_init_only is set to PETSC_TRUE (InitDef scheme), the application of the preconditioner consists only of application of the additional preconditioner M^{-1}. Otherwise, the preconditioner application consists of P*M^{-1} + factor*Q. The first part of the preconditioner (PM^{-1}) shifts some eigenvalues to zero while the addition of the coarse problem correction (factor*Q) makes the preconditioner to shift some eigenvalues to the given factor. The InitDef scheme is recommended for deflation using high accuracy estimates of eigenvectors of A when it exhibits similar convergence to the full deflation but is cheaper.

The deflation matrix is by default automatically computed. The type of deflation matrix and its size to compute can be controlled by PCDeflationSetSpaceToCompute() or -pc_deflation_compute_space and -pc_deflation_compute_space_size. User can set an arbitrary deflation space matrix with PCDeflationSetSpace(). If the deflation matrix is a multiplicative MATCOMPOSITE, a multilevel deflation [3] is used. The first matrix in the composite is used as the deflation matrix, and the coarse problem (W'*A*W)^{-1} is solved by KSPFCG (if A is MAT_SPD) or KSPFGMRES preconditioned by deflation with deflation matrix being the next matrix in the MATCOMPOSITE. This scheme repeats until the maximum level is reached or there are no more matrices. If the maximum level is reached, the remaining matrices are merged (multiplied) to create the last deflation matrix. The maximum level defaults to 0 and can be set by PCDeflationSetLevels() or by -pc_deflation_levels.

The coarse problem KSP can be controlled from the command line with prefix -deflation_ for the first level and -deflation_[lvl-1] from the second level onward. You can also use PCDeflationGetCoarseKSP() to control it from code. The bottom level KSP defaults to KSPPREONLY with PCLU direct solver (MATSOLVERSUPERLU/MATSOLVERSUPERLU_DIST if available) wrapped into PCTELESCOPE. For convenience, the reduction factor can be set by PCDeflationSetReductionFactor() or -pc_deflation_recduction_factor. The default is chosen heuristically based on the coarse problem size.

The additional preconditioner can be controlled from command line with prefix -deflation_[lvl]_pc (same rules used for coarse problem KSP apply for [lvl]_ part of prefix), e.g., -deflation_1_pc_pc_type bjacobi. You can also use PCDeflationGetPC() to control the additional preconditioner from code. It defaults to PCNONE.

The coarse problem correction term (factor*Q) can be turned on by -pc_deflation_correction and the factor value can be set by pc_deflation_correction_factor or by PCDeflationSetCorrectionFactor(). The coarse problem can significantly improve convergence when the deflation coarse problem is not solved with high enough accuracy. We recommend setting factor to some eigenvalue, e.g., the largest eigenvalue so that the preconditioner does not create an isolated eigenvalue.

The options are automatically inherited from the previous deflation level.

The preconditioner supports KSPMonitorDynamicTolerance(). This is useful for the multilevel scheme for which we also recommend limiting the number of iterations for the coarse problems.

See section 3 of [4] for additional references and decription of the algorithm when used for conjugate gradients. Section 4 describes some possible choices for the deflation space.

Developer Notes

Contributed by Jakub Kruzik (PERMON), Institute of Geonics of the Czech Academy of Sciences and VSB - TU Ostrava.

Developed from PERMON code used in [4] while on a research stay with Prof. Reinhard Nabben at the Institute of Mathematics, TU Berlin.

References

[1] - A. Nicolaides. "Deflation of conjugate gradients with applications to boundary value problems", SIAM J. Numer. Anal. 24.2, 1987.
[2] - Z. Dostal. "Conjugate gradient method with preconditioning by projector", Int J. Comput. Math. 23.3-4, 1988.
[3] - Y. A. Erlangga and R. Nabben. "Multilevel Projection-Based Nested Krylov Iteration for Boundary Value Problems", SIAM J. Sci. Comput. 30.3, 2008.
[4] - J. Kruzik "Implementation of the Deflated Variants of the Conjugate Gradient Method", Master's thesis, VSB-TUO, 2018 - http://dspace5.vsb.cz/bitstream/handle/10084/130303/KRU0097_USP_N2658_2612T078_2018.pdf

See Also

PCCreate(), PCSetType(), PCType (for list of available types), PC,
PCDeflationSetInitOnly(), PCDeflationSetLevels(), PCDeflationSetReductionFactor(), PCDeflationSetCorrectionFactor(), PCDeflationSetSpaceToCompute(), PCDeflationSetSpace(), PCDeflationSpaceType, PCDeflationSetProjectionNullSpaceMat(), PCDeflationSetCoarseMat(), PCDeflationGetCoarseKSP(), PCDeflationGetPC()

Level

intermediate

Location

src/ksp/pc/impls/deflation/deflation.c
Index of all PC routines
Table of Contents for all manual pages
Index of all manual pages