Merge branch 'maint'

Calls to MatGetDiagonal() for MATSHELL should fail with an error if the MATSHELL did not provide a MatGetDiagonal() method

Commit-type: bug-fix, testing-fix
Reported-by: Stefano Zampini <stefano.zam

Calls to MatGetDiagonal() for MATSHELL should fail with an error if the MATSHELL did not provide a MatGetDiagonal() method

Commit-type: bug-fix, testing-fix
Reported-by: Stefano Zampini <stefano.zampini@gmail.com>

show more ...

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \(

Man pages: add newlines after "Notes:"

This allows for proper formatting from sowing.

On OS X (using gsed, not the default BSD sed), from the PETSc root directory:

find src include -type f \( -name "*.c" -or -name "*.h" -or -name "*.cxx" \) | xargs gsed -i 's/Notes\s*:\s*\(\w.*\)/Notes:\n \1/'

This adds a newline and 4 spaces whenever "Notes:" is followed by any "word" character, in any .c, .h, or .cxx file in src/ or include/

show more ...

Add MatGetOperation() and fixes to MATSHELL

Missing break statement in MatShellSetOperation

missing CHKERRQ() on PetscPrintf() statements

Commit-type: style-fix
Reported-by: nightly static analyzer

Fixes doe mangled manual page

Commit-type: documentation
Reported-by: Jose E. Roman

MatCopy_Shell() needs to clean out old axpy

Commit-type: bug-fix
Reported-by: nightly tests

Missing {}

Commit-type: bug-fix
Reported-by: nightly tests

another complex number comparison to 0

Commit-type: bug-fix
Reported-by: nightly tests

Error checking and minor additional functional for MatAXPY() case for MATSHELL

Commit-type: feature

Set incorrect pointer for provided function

Commit-type: bug-fix
Reported-by: Stefano Zampini <stefano.zampini@gmail.com>

Refactorization of MATSHELL now supports matrix free version of ts ex50.c

which uses composition of matrix free matrices.

Commit-type: feature

Begin refactorization of MatShell to make the code more straightforward and extensible

Introduce MatSetOperation() for changing operations on non-MATSHELLs
Introduce two modes for MATSHELL
1) MAT

Begin refactorization of MatShell to make the code more straightforward and extensible

Introduce MatSetOperation() for changing operations on non-MATSHELLs
Introduce two modes for MATSHELL
1) MATSHELL manages scalings and shifts automatically
2) Creator of MATSHELL can provide their own code for managing scaling and shifts
This is used by MatCreateADA() and MatCreateSubMatrixFree()

The extensibility is so that a MatShell can be easily used in as many places as possible as other matrices,
for example,

/* A = diag(u) G */
ierr = MatCopy(appctx->SEMop.grad,A,SAME_NONZERO_PATTERN);CHKERRQ(ierr);
ierr = MatDiagonalScale(A,globalin,NULL);CHKERRQ(ierr);
/* A = A + diag(Gu) */
ierr = VecDuplicate(globalin,&Gglobalin);CHKERRQ(ierr);
ierr = MatMult(appctx->SEMop.grad,globalin,Gglobalin);CHKERRQ(ierr);
ierr = MatDiagonalSet(A,Gglobalin,ADD_VALUES);CHKERRQ(ierr);
ierr = VecDestroy(&Gglobalin);CHKERRQ(ierr);
/* A = K - A */
ierr = MatScale(A,-1.0);CHKERRQ(ierr);
ierr = MatAXPY(A,1.0,appctx->SEMop.keptstiff,SAME_NONZERO_PATTERN);

Commit-type: feature

show more ...

Cannot use 0 for PetscBool value

Commit-type: bug-fix, portability-fix
Reported-by: travis

Add code for automatically checking if shell matrix vector and transpose matrix vector products work with TS

Added MatShellTestMult(), MatShellTestMultTranspose(), TSRHSJacobianTest(), and TSRHSJaco

Add code for automatically checking if shell matrix vector and transpose matrix vector products work with TS

Added MatShellTestMult(), MatShellTestMultTranspose(), TSRHSJacobianTest(), and TSRHSJacobianTestTranspose()

Commit-type: feature

show more ...

Start to add 1d matrix free version of burgers spectral ODE example

Commit-type: feature, example

Merge branch 'maint'