Fixed types for PCGAMGProlongator_Classical and PCGAMGGraph_Classical, strict C fix.

Complex fixes for classical

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

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 ...

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 ...

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 ...

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)

PCGAMG: namespace PCGAMGType macros and fix Fortran type

PetscFunctionListAdd: use macro to hide (PetscVoidFunction) cast

Same rationale as the last commit:

Casting using PetscVoidFunction instead of (void(*)(void)) is needed
with C++ compilers that will

PetscFunctionListAdd: use macro to hide (PetscVoidFunction) cast

Same rationale as the last commit:

Casting using PetscVoidFunction instead of (void(*)(void)) is needed
with C++ compilers that will not implicitly convert C++ language linkage
to C language linkage. Since PetscFunctionListAdd() is meant to operate
on generic function pointers, we hide this language linkage detail
(which is inconsistently implemented) in a cast.

show more ...

PetscFunctionListFind: use macro to hide (PetscVoidFunction*) cast

Casting using PetscVoidFunction* instead of (void(**)(void)) is needed
with C++ compilers that will not implicitly convert C++ lang

PetscFunctionListFind: use macro to hide (PetscVoidFunction*) cast

Casting using PetscVoidFunction* instead of (void(**)(void)) is needed
with C++ compilers that will not implicitly convert C++ language linkage
to C language linkage. Since PetscFunctionListFind() is meant to operate
on generic function pointers, we hide this language linkage detail
(which is inconsistently implemented) in a cast.

show more ...

removed string version of function name for XXXRegister(), PetscFunctionListAdd() and PetscObjectComposeFunction()

removed path and MPI_Comm arguments from PetscFunctionListFind/Add()

small fixes

Hg-commit: 15fc3499cf02f49190f576669860cdb9b667729c

Change all PETSC_EXTERN_C to PETSC_EXTERN

Hg-commit: 8d2ebbb193fb583bccc64015e35640c4e08c3426

merge, terrible manual process with many conflicts with Jed's PETSC_INTERNAL

Hg-commit: 84df07d03c6e55bd0f27bd5ee8c1738562bd529d

visibility: use PETSC_INTERN when possible

Hg-commit: 3b4bfb5c23843293f1f08790b48437080a32eb42