PCBDDC: minor changes to bddcprivate

PCBDDC: add subassembling_init and use_coarse_estimates to private bddc data struct for future use

PCBDDC: coarse_subassembling uninitialized

PCBDDC: dbg viewer is the default viewer on the communicator

PCFieldSplit: Schur can use lumping with MatSchurComplement.

PCFieldSplit implements SELFP preconditioning for Schur using MatSchurComplementGetPmat to assemble Sp.

spai: update spai to using -DSPAI_USE_MPI [instead of -DMPI]
This fixes clang compile error with -DMPI & mpi.h [from MPICH]

spai: remove un-needed 'num_ptr' and simplify the loops [and remove dead loop code]

spai: fix regression/valgrind-errors with spai test suite. The issue
must have started after MatRestoreRow() behavior was changed in 4b4408bb

Merge branch 'barry/remove-hmpi'

Conflicts:
src/sys/objects/mpinit.c

Merge branch 'barry/fix-singleton-viewer'

Merge branch 'barry/fix-mpiint-petscoptionsgetint'

Conflicts:
src/ksp/pc/impls/redistribute/redistribute.c
src/mat/impls/aij/mpi/mpiaij.c
src/mat/impls/dense/mpi/mpidense.c
src/mat/utils/matstas

Merge branch 'barry/fix-mpiint-petscoptionsgetint'

Conflicts:
src/ksp/pc/impls/redistribute/redistribute.c
src/mat/impls/aij/mpi/mpiaij.c
src/mat/impls/dense/mpi/mpidense.c
src/mat/utils/matstash.c
src/snes/examples/tutorials/ex28.c
src/sys/objects/pinit.c
src/vec/is/ao/impls/memscalable/aomemscalable.c

show more ...

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 PetscCalloc[1-7]

Replace lots of places where PetscMemzero is called immediately after
allocation.

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

PetscMalloc[2-7]: remove type arguments, infer from pointer type

This change can be performed almost entirely by this script:

git grep -l PetscMalloc[2-7] $1 | xargs perl -pi \
-e '
s@PetscMalloc

PetscMalloc[2-7]: remove type arguments, infer from pointer type

This change can be performed almost entirely by this script:

git grep -l PetscMalloc[2-7] $1 | xargs perl -pi \
-e '
s@PetscMalloc2\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc2($1,$3,$4,$6)@;
s@PetscMalloc3\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc3($1,$3,$4,$6,$7,$9)@;
s@PetscMalloc4\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc4($1,$3,$4,$6,$7,$9,$10,$12)@;
s@PetscMalloc5\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc5($1,$3,$4,$6,$7,$9,$10,$12,$13,$15)@;
s@PetscMalloc6\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc6($1,$3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18)@;
s@PetscMalloc7\(([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+),([^,;()]+)\)@PetscMalloc7($1,$3,$4,$6,$7,$9,$10,$12,$13,$15,$16,$18,$19,$21)@;
'

show more ...

fixed many examples that worked incorrectly for 64 bit indices

fixed PetscViewerASCIIPrintf() for a single viewer so that the singletons ASCII output is actually preserved

removed all HMPI code
- didn't work any more because of randomized memory messing up function pointers
- complicated code to maintain that no one used and they would prefer threads anyways

added support for PCSVD for complex numbers

Merge branch 'jed/fix-pcview-fieldsplit-early'

* jed/fix-pcview-fieldsplit-early:
PCView_FieldSplit_Schur: fix crash on view before setup

added FILE* argument to PetscSynchronizedFlush() so it knows where to flush the off process values

Previously if process 0 did not call PetscSynchronizedFPrintf() then it output the off process valu

added FILE* argument to PetscSynchronizedFlush() so it knows where to flush the off process values

Previously if process 0 did not call PetscSynchronizedFPrintf() then it output the off process values
to stdout even if the file was not stdout. Also affected PetscViewerASCIISynchronizeFlush()

Reported-by: Randall Mackie <rlmackie862@gmail.com>

show more ...

PCFieldSplit: fix documentation of full factorization and Schur

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

since external packages use a variety of different types for 64 bit indices
(though all have the same bit values) casts are needed for some compilers to
prevent warning messages when using 64 bit ind

since external packages use a variety of different types for 64 bit indices
(though all have the same bit values) casts are needed for some compilers to
prevent warning messages when using 64 bit indices

show more ...

PCBDDC: move creation of temporary MATIS matrix before coarse ksp is created.

This way, we can temporary set the coarse MATIS into coarse ksp.
This allows to call pc_bddc_coarse_ksp_fischer_guess <>

PCBDDC: move creation of temporary MATIS matrix before coarse ksp is created.

This way, we can temporary set the coarse MATIS into coarse ksp.
This allows to call pc_bddc_coarse_ksp_fischer_guess <> from option

show more ...