Fixed a number of bugs in Classical GAMG

Fix to the threshold calculation

Fixed PCGAMG_Classical so that thresholding decreases fill as bugs are fixed

Fixed one-processor memory corruption bug for classical GAMG.

Not the normal way those things typically go..

Namespacing in gamg/classical.c

Added classical AMG to GAMG

completed initial implementation of PetscCitationsRegister() and added proper cites for hypre, ml, mumps

-citations option causes the print out of approriate bibtex citation indices for all packages USED in the run
Most of the appropriate citations still need to be added

Merge branch 'prbrune/dmda-subdomains-per-process'

Fix PCGASM leak from not destroying subDMs from DMCreateDomainDecomposition()

gamg: change '1.e300' to 'PETSC_MAX_REAL'. Fixes single precision 'runex55_SA'

Merge branch 'jed/stack-call-blas'

Fix performance bug that calling BLAS takes an amount of time dependent
on the number of memory allocations that the program has made.

* jed/stack-call-blas:
Pe

Merge branch 'jed/stack-call-blas'

Fix performance bug that calling BLAS takes an amount of time dependent
on the number of memory allocations that the program has made.

* jed/stack-call-blas:
PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)
Stack: consolidate handling of PetscStack

show more ...

Merge branch 'madams/gamg-params'

* madams/gamg-params:
PCGAMG: change default nsmooths, reuse_interpolation, threshold

PCGAMG: change default nsmooths, reuse_interpolation, threshold

-pc_gamg_agg_nsmooths 1 (was 0)

Although unsmoothed aggregation is more robust, smoothed is faster for
the problems that AMG is

PCGAMG: change default nsmooths, reuse_interpolation, threshold

-pc_gamg_agg_nsmooths 1 (was 0)

Although unsmoothed aggregation is more robust, smoothed is faster for
the problems that AMG is most frequently used for.

-pc_gamg_reuse_interpolation FALSE (was TRUE)

Lagging interpolation by default causes GAMG to spuriously fail for
many nonlinear problems and transient problems with moving coefficient
structure, and with variable time steps.

-pc_gamg_threshold 0.0 (was .001)

ML does not use thresholding by default, so this change makes GAMG
algorithmically more similar. Increasing the threshold causes slower
coarsening, which increases grid complexity. It is critical for
anisotropic problems.

show more ...

Merge branch 'toby/pcml-nnz'

* toby/pcml-nnz:
PCML: set N_nonzeros so thatfine-grid operator complexity is reported

PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)

When CHKMEMQ occurs on every call to BLAS or Lapack, we sometimes see
unacceptable performance degradation (orders of magnitude

PetscStackCallBLAS: new macro instead of PetscStackCall (skip CHKMEMQ)

When CHKMEMQ occurs on every call to BLAS or Lapack, we sometimes see
unacceptable performance degradation (orders of magnitude).
PetscStackCallBLAS is currently identical to PetscStackCall except for
skipping CHKMEMQ, but such checking could be made optional in the
future.

show more ...

PCML: set N_nonzeros so thatfine-grid operator complexity is reported

With the version of ML in Trilinos (but not with stand-alone ML, fetched
by --download-ml), failure to set this variable results

PCML: set N_nonzeros so thatfine-grid operator complexity is reported

With the version of ML in Trilinos (but not with stand-alone ML, fetched
by --download-ml), failure to set this variable results in bogus values
(sometimes negative) being given, instead of reasonable output like:

Smoothed Aggregation : operator complexity = 1.358396e+00

show more ...

Merge branch 'madams/gamg-aggoptions'

More delegation and less method-overriding.

* madams/gamg-aggoptions:
PCGAMG: call subtype destroy from PCDestroy_GAMG
PCGAMG: move subtype pointers to _PC

Merge branch 'madams/gamg-aggoptions'

More delegation and less method-overriding.

* madams/gamg-aggoptions:
PCGAMG: call subtype destroy from PCDestroy_GAMG
PCGAMG: move subtype pointers to _PCGAMGOps.
PCGAMG: set subtype in PCSetFromOptions_GAMG instead of PCCreate_GAMG

show more ...

Merge branch 'jed/pcview-before-setup'

* jed/pcview-before-setup:
PCView_MG: avoid SEGV before PCSetUp
PCView: warn about incomplete information before PCSetUp

MatRestoreRow: zero arguments after restore to catch accidental misuse

Explicitly pass NULL when we are intentionally using the number of
columns after restoring the row.

PCView_MG: avoid SEGV before PCSetUp

PCGAMG: call subtype destroy from PCDestroy_GAMG

Having the subtype call the parent is more fragile because the
delineation of responsibilities is less clear and the call order could
be different fo

PCGAMG: call subtype destroy from PCDestroy_GAMG

Having the subtype call the parent is more fragile because the
delineation of responsibilities is less clear and the call order could
be different for different subtypes.

show more ...

PCGAMG: move subtype pointers to _PCGAMGOps.

Avoid stomping on user subtype preference in PCSetFromOptions_GAMG.

Subtype function pointers are now cleared by PCGAMGSetType to avoid
lingering from t

PCGAMG: move subtype pointers to _PCGAMGOps.

Avoid stomping on user subtype preference in PCSetFromOptions_GAMG.

Subtype function pointers are now cleared by PCGAMGSetType to avoid
lingering from the last type in cases where not all pointers get set.
Also store the type name to help with debugging and diagnostic output.

show more ...

PCGAMG: set subtype in PCSetFromOptions_GAMG instead of PCCreate_GAMG

PCGAMG: add PCGAMGInitializePackage (fixes memory leak)