PCBDDC: rearrange symmetry, positive-definitess stuffs

Added pc_bddc_symmetric command line options to shift to unsymmetric basis computation

PCBDDC: no need for eigensolver is S == St

PCBDDC: when n_subs == 1 and size_schur == size_active_schur reuse already inverted block

MatSOR takes real shift not complex

Add interface for changing Euclid smoother parameters in BoomerAMG

The number of levels k, drop tolerance and whether a Block-Jacobi smoother is used in the Euclid ILU(k) smoother can now be set wit

Add interface for changing Euclid smoother parameters in BoomerAMG

The number of levels k, drop tolerance and whether a Block-Jacobi smoother is used in the Euclid ILU(k) smoother can now be set with the following PETSc options

-pc_hypre_boomeramg_eu_level (default = 0)
-pc_hypre_boomeramg_eu_droptolerance (default=0)
-pc_hypre_boomeramg_eu_bj (default = not set)

show more ...

Default number of complex smoother levels set to 25

Set the default number of levels on which the more complex smoothers are used to 25 (otherwise it will be 0, and the more complex smoothers are no

Default number of complex smoother levels set to 25

Set the default number of levels on which the more complex smoothers are used to 25 (otherwise it will be 0, and the more complex smoothers are not used).

show more ...

PCASM: Added multiplicative local composition
- Added PCASMSet/GetLocalType()

PCBDDC: fixes for compiler warnings

PCBDDC: fix for complexes

hermitian and positive definitess checks should be replaced by something which is more robust

PCBDDC: complete support for MUMPS Schur solver reusage

is_R and is_B should be created during SubSchursSetUp
Added preprocessor guards

PCBDDC: added support for transpose application of the BDDC when reusing MUMPS solver

PCBDDC: reuse Schur solver provided by MUMPS when possible

MATIS: n_sends should stay there!!!!

if n_sends is not part of the logic, we cannot overlap levels!

Allow setting the number of levels on which the more complex smoothers are used

By using the option -pc_hypre_boomeramg_smooth_num_levels the number of levels on which the more complex smoothers are

Allow setting the number of levels on which the more complex smoothers are used

By using the option -pc_hypre_boomeramg_smooth_num_levels the number of levels on which the more complex smoothers are used is set via calling HYPRE_BoomerAMGSetSmoothNumLevels. The default is 25, which is the default number of multigrid levels.

show more ...

Allow selection of more complex smoothers for hype's BoomerAMG

By using the option -pc_hypre_boomeramg_smooth_type more complex smoothers can be set for hypre via the method HYPRE_BoomerAMGSetSmooth

Allow selection of more complex smoothers for hype's BoomerAMG

By using the option -pc_hypre_boomeramg_smooth_type more complex smoothers can be set for hypre via the method HYPRE_BoomerAMGSetSmoothType()

show more ...

Zero out destroyed subdomain lists to prevent weird bugs(hence, pass by pointer).

PCBDDC: optimize for reusage of Schur solver (if any)

R_to_B and (optionally) R_to_D scatters are different
reusage is not possible when change of basis has been requested

Merge branch 'mark/gamg-square'

Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED inste

Infrastructure that allows failures in PCSetUp(), PCApply(), MatMult() etc due to, for example a zero pivot or function evaluation outside its domain, to propagate up and become a KSP_DIVERGED instead of generating an error that stops the program. In response to Issue 96. This includes failures in MatCreateSNESMF() applications due to domain errors. The mechanism to propagate some errors is by setting Info or Nan into the output vector and using the norm or inner product reductions in SNES or KSP to propagate the error condition to all processes and then handling them immediately after the norm or inner product.

This allows, for example, ODE integrators to try again with a smaller time-step if the PCSetUp failed instead of requiring a complete restart of the run with other options.

Currently some error conditions, such as function domain error in a line search may not get progated up using the correct SNESConvergedReason

See src/snes/examples/tests/ex69.c for the handling of several conditions

show more ...

PCBDDC: reuse as much as possible the MUMPS Schur solver

optimizations still need to be applied

Merge remote-tracking branch 'origin/stefano_zampini/mattransposedense-ldafix' into stefano_zampini/pcbddc-adaptive-selection

Merge branch 'jose/upgrade-cusp-0.5.1'

jose/upgrade-cusp-0.5.1:
Fixes build problems with CUSP 0.5.0 and 0.5.1 as well as CUDA 7.

Fixed function declaration/definition mismatch: PetscMPIInt —> PetscInt.

PCBDDC: bugfixes

- change of basis falg should be set for each cc (not primal constraint)
- work should be set to zero when importing boundary values of A_RRmA_RV to avoid undefined values
- A_VV sh

PCBDDC: bugfixes

- change of basis falg should be set for each cc (not primal constraint)
- work should be set to zero when importing boundary values of A_RRmA_RV to avoid undefined values
- A_VV should be converted to BAIJ when working with SeqBAIJ

show more ...

changed GAMG square graph parameter to int from bool