PCMG Galerkin: use A_coarse = R*A*P when R != P^T

Formerly, Galerkin coarse operators were always defined as P^T A P which
is not consistent when the user has chosen R != P^T.

Reported-by: Torquil

PCMG Galerkin: use A_coarse = R*A*P when R != P^T

Formerly, Galerkin coarse operators were always defined as P^T A P which
is not consistent when the user has chosen R != P^T.

Reported-by: Torquil Macdonald Sørensen <torquil@gmail.com>
Reported-by: Boris Kaus <kaus@uni-mainz.de>

show more ...

Merge branch 'master' into barry/remove-big-g

Doc: Fixed MG Galerkin docs to note that different restriction and interpolation are possible
Reported-by: Torquil Macdonald Sørensen <torquil@gmail.com>

completed some percentage of changes from %G format in prints

Sys: drop explicit type arguments from PetscNew() and PetscNewLog()

git grep -l 'PetscNew\(Log\)\?(' -- $1 | xargs perl -pi -e '
s@PetscNew\([^,;()]+ *, *@PetscNew(@;
s@PetscNewLog\(([^,;()]+) *,[^,

Sys: drop explicit type arguments from PetscNew() and PetscNewLog()

git grep -l 'PetscNew\(Log\)\?(' -- $1 | xargs perl -pi -e '
s@PetscNew\([^,;()]+ *, *@PetscNew(@;
s@PetscNewLog\(([^,;()]+) *,[^,;()]+, *@PetscNewLog($1,@'

show more ...

Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(

Sys: add PetscMalloc1 macro, array allocation without redundant types

The type is inferred from the pointer return type. This patch is
automated via the following script:

git grep -l 'PetscMalloc(.*sizeof' src | xargs perl -pi -e 's@PetscMalloc\(([^,;]*[^,; ]) *\* *sizeof\([^,;()]+\),@PetscMalloc1($1,@'

This commit contains an additional bug-fix in csrperm.c, fixing pointer
arity. The code was introduced in 2006, but the allocation could not
have been correct at any time. This probably means that
MatDuplicate_SeqAIJPERM has never been tested.

a54129beb540034ba105796c682d589e7e1111f2
Richard Tran Mills <rmills@ornl.gov>

Added MATSEQCSRPERM support for MatDuplicate() and conversion to/from
MATSEQAIJ. Note that these changes are not quite debugged.

show more ...

Merge branch 'madams/mat-residual'

* madams/mat-residual:
refactored residuals by adding a 'residual' method to Mat

Merge branch 'jed/converged-default-and-skip'

This process started with SNESDefaultConverged -> SNESConvergedDefault
in v3.4, but the other parts have been added piecemeal. Deprecation
support in C

Merge branch 'jed/converged-default-and-skip'

This process started with SNESDefaultConverged -> SNESConvergedDefault
in v3.4, but the other parts have been added piecemeal. Deprecation
support in C interface.

* jed/converged-default-and-skip:
SNES and KSP: add deprecation for *DefaultConverged and *SkipConverged
KSP: rename KSPDefaultConverged* to KSPConvergedDefault*
KSP: rename KSPSkipConverged -> KSPConvergedSkip
SNES: document SNESSkipConverged -> SNESConvergedSkip
snes->ttol to be used only in convergence tests.
Renamed SNESKSPEW_Pre/PostSolve --> KSPPre/PostSolve_SNESEW.
Renamed SNESSkipConverged --> SNESConvergedSkip.

Conflicts:
include/petscsnes.h
src/docs/website/documentation/changes/dev.html

show more ...

refactored residuals by adding a 'residual' method to Mat

KSP: rename KSPSkipConverged -> KSPConvergedSkip