Better to use the proper error handling version of PetscViewerASCIIGetStdout() instead of PETSC_VIEWER_STDOUT_()

Merge branch 'maint'

* maint:
MatGetGhosts: add missing implementation for MPIAIJ
error in manual page from cut and paste

error in manual page from cut and paste

Reported-by: Dave May <dave.mayhem23@gmail.com>

allow clearing the NullSpace or NearNullSpace of a matrix by passing in a NULL

Merge remote-tracking branch 'petsc/jed/mat-assembly-perf' into matpartitioning-hierach

To make the interface PetscCommBuildTwoSided(MPI_Comm comm,PetscMPIInt count,MPI_Datatype dtype,PetscInt nto,c

Merge remote-tracking branch 'petsc/jed/mat-assembly-perf' into matpartitioning-hierach

To make the interface PetscCommBuildTwoSided(MPI_Comm comm,PetscMPIInt count,MPI_Datatype dtype,PetscInt nto,const PetscMPIInt *toranks,const void *todata,PetscInt *nfrom,PetscMPIInt **fromranks,void *fromdata) compatible with Jed's new changes.
DataType of 'nto' and 'nfrom' is changed to PetscMPIInt.

Conflicts:
include/petscsys.h
src/sys/utils/sorti.c
src/vec/vec/examples/tests/makefile

show more ...

Modified an interface

Mat: get rid of MatCheckPreallocated in MatCreateVecs

removed unused variable size

Fixed complier warnings and errors
and moved the implemenration of matincreaseoverlapsplit to utils

Merge remote-tracking branch 'petsc/master' into pcgasm-increaseoverlap

documentation for PetscViewerASCIIOpen() and MatView() and VecView() was WRONG for use with a
viewer that comes from PETSC_COMM_SELF. The documentation mistakenly said that the output was independent

documentation for PetscViewerASCIIOpen() and MatView() and VecView() was WRONG for use with a
viewer that comes from PETSC_COMM_SELF. The documentation mistakenly said that the output was independent for
each process when in fact the communicator of the viewer must match that of the object.

Reported-by: Longyin Cui <cuilongyin@gmail.com>

show more ...

Fixed

Implemented MatIncreaseOverlapSplit

An incomplete implementation of MatIncreaseOverlapSplit

doc: fix typo in MatMatSolve() manpage

Merge remote-tracking branch 'petsc/master' into mat-increaseoverlap

error out if invalide MatReuse option is passed in

Merged in PR312: karpeev/ksp-pcgasm-overhaul.

MatGetSubmatrix() could hang when some processes thought the entire matrix as being selected and others did not.

Reported-by: Mark Adams <mfadams@lbl.gov>

deeper cleanup of matrix when the type is reset

Merge branch 'barry/improve-nullspace-handling'

Merge branch 'barry/propagate-pcsetup-failures'

removed KSPSetNullSpace() since mat already carries that information
and having two locations for it is confusing, unnecessary and prone to error
added MatSetTransposeNullSpace() for providing the nu

removed KSPSetNullSpace() since mat already carries that information
and having two locations for it is confusing, unnecessary and prone to error
added MatSetTransposeNullSpace() for providing the null space of the transpose of the matrix and added support in KSPSolve() to remove this null space from the right handside of the linear system thus making the system consistent.

show more ...

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

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

MatGetSubMatricesParallel —> MatGetSubMatricesMPI. Implementation rewritten, tests added.