Revert "replace 'const MatFactorInfo' to 'MatFactorInfo' in MatFactor routi"

This reverts commit 927faaad4e7a9177b263c19fc2f4e1f9dae8f67a.

Revert "replace 'const MatFactorInfo' -> 'MatFactorInfo'"

This reverts commit af48bdc5724c729ddc55b84e8afd48f344de3945.

add "Mat fact" as an input to MatPivotCheck() for setting errortype into matrix fact

replace 'const MatFactorInfo' -> 'MatFactorInfo'

added many checks for integer overflow in symbolic factorization and symbolic matrix matrix multiplication

replace 'const MatFactorInfo' to 'MatFactorInfo' in MatFactor routi
nes

Merge branch 'master' into barry/petscoptions-object

Merge remote-tracking branch 'origin/hzhang/matmatmult_dense'

Complete update to new PetscOptions interface

add MAT_INPLACE_MATRIX; replace MAT_REUSE_MATRIX with it for MatConvert()

Merge branch 'barry/add-mpiu_allreduce'

Merge branch 'barry/fix-ts-matshift-error'

Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()
This is a bug that can be difficult to track down.

Suggested-by

Introduce MPIU_Allreduce() to detect cases where all processes that share an MPI comm do not call the same sequence of MPI_Allreduce()
This is a bug that can be difficult to track down.

Suggested-by: Eric Chamberland <Eric.Chamberland@giref.ulaval.ca>
Suggested-by: Patrick Lacasse <placasse@giref.ulaval.ca>

show more ...

Merge branch 'maint'

If only some processes had no local entries in the matrix the MatSeqXAIJPreallocation() changed the value of aij->nonew on only those
processes thus resulting in gridlock in MatAssemblyEnd_MPIXAIJ()

If only some processes had no local entries in the matrix the MatSeqXAIJPreallocation() changed the value of aij->nonew on only those
processes thus resulting in gridlock in MatAssemblyEnd_MPIXAIJ() since some processes wanted to update the matrix nonzero state while others did not

Reported-by: Eric Chamberland <Eric.Chamberland@giref.ulaval.ca>

show more ...

fix violations of PETSc style guide: Usage of SETERRQ and NULL

Added missing MatMissingDiagonal() implementations

MatHeaderReplace() corrupted the -objects_dump array

Reported-by: Torquil Macdonald Sørensen <torquil@gmail.com>

MatHeaderReplace() and MatHeaderMerge() destroy the second matrix argument therefor

MatHeaderReplace() corrupted the -objects_dump array

Reported-by: Torquil Macdonald Sørensen <torquil@gmail.com>

MatHeaderReplace() and MatHeaderMerge() destroy the second matrix argument therefor pass it as
a pointer to it can be zeroed and not mistakenly reused.

show more ...

Merge branch 'jed/mat-assembly-perf' of bitbucket:petsc/petsc

VecAssembly and MatAssembly now use a scalable exchange pattern based on
PetscCommBuildTwoSided. This feature can be controlled with th

Merge branch 'jed/mat-assembly-perf' of bitbucket:petsc/petsc

VecAssembly and MatAssembly now use a scalable exchange pattern based on
PetscCommBuildTwoSided. This feature can be controlled with the options

-vec_assembly_bts 0 or 1 (default 0)
-matstash_bts 0 or 1 (default 0)

The rationale is that the new implementation with scalable data
structures can be slightly slower than the old version at small process
counts. The default here could be changed to depend on the process
count (leading to possibly-confusing scaling performance diagnostics) or
the implementation could learn to take a fast path.

* 'jed/mat-assembly-perf' of bitbucket:petsc/petsc: (49 commits)
MatStash: fix -Wsign-compare by using size_t for loop index when max is also size_t
Sys BuildTwoSided test: fix for non-POD std::complex
MatStash: cast to satisfy non-structural MPI type tag check
mpiuni: fix compile error /sandbox/petsc/petsc.clone-2/arch-linux-uni/lib/libpetsc.so: undefined reference to `MPI_Type_create_resized'
Vec: Silence compiler warning
Vec: fix typo in comment
Sys: fix C89 compiler warning
VecStash BTS: fix block stash InsertMode accounting
VecStash BTS: fix indexing bug counting sends to rank 0
MatStash BTS: fix memory leak on MAT_SUBSET_OFF_PROC_ENTRIES
VecAssemblyEnd_MPI_BTS: fix donotstash code path
VecAssemblyEnd_MPI_BTS: fix C++ conversion to InsertMode
Sys: fix datatypes test using MPI_Type_create_resized
MatStash BTS: work around lack of offsetof() for non-POD (std::complex)
MatAssembly: move check for InsertMode consistency into MatStashScatter impls
MatStash BTS: fix memory leak on reassembly with MAT_SUBSET_OFF_PROC_ENTRIES
MatStash BTS: small simplification to row ownership calculation
MatStash BTS: add MAT_SUBSET_OFF_PROC_ENTRIES, impl with neighbor-only comm
MatStash: initial BTS (BuildTwoSided) implementation
MatStash: add extension point for new BTS implementation
...

show more ...