PCBDDC: move change of basis on pressure during PreSolve

Added check for b(v_I,p_0) = 0 for all v_I (o on Gamma and \partial\Omega)

PCBDDC: extend computation of coarse basis functions and coarse matrix for the benign trick

debug output temporary suppressed,
the code is still not working.
- GMRES with BDDC setup and dummy applic

PCBDDC: extend computation of coarse basis functions and coarse matrix for the benign trick

debug output temporary suppressed,
the code is still not working.
- GMRES with BDDC setup and dummy application (PCApply as VecCopy) does not converge
- still need to extract p0 from residual

show more ...

PCBDDC: added support for the benign trick in PCBDDCSubSchurs

The setup code of SubSchurs need to be expanded to symmetric indefinite factorizations support from MUMPS

PCBDDC: need to change iteration matrix in case the benign trick has been requested

PCBDDC: added function to Pop or Push B0 from local mat

PCBDDC: handle change of basis on local pressures

B0 is extracted from the local matrices before setup. Still needs to be added back when setup is done.

Replaced PetscViewerASCIISynchronizedAllow() with PetscViewerASCIIPushSynchronized() PetscViewerASCIIPopSynchronized()

Merge branch 'master' into knepley/feature-snes-deflation

* master: (206 commits)
Bib: Added ref
fix the output file for the ts example ex20adj
Bugfix: setting post-step function in ztsf.c
f

Merge branch 'master' into knepley/feature-snes-deflation

* master: (206 commits)
Bib: Added ref
fix the output file for the ts example ex20adj
Bugfix: setting post-step function in ztsf.c
fix makefile format error with 'all-legacy'
doc: handle missing dev.html in a release tarball
doc: add in petscistypes.h to doc sourcefile list
Increase patchlevel to 3.6.1
update some example output files
Doc+Vec: Fixed VecSetValuesBlocked() Reported by: Shao-Ching Huang <huangsc@gmail.com>
clarify that PETSc's "thread safe" code only works with OpenMP
allow adjoint integration over split time intervals
DT: Goddman it. Why are people stepping on my commits
Bib: Added ref
fix the bug in the endpoint version of discrete adjoint theta and correct the affected example outputs
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
fix: buggy change in 33da457b864643cfcac2768e47197c286471a5a4
updated petscthreads.html with information on thread per object model
fix for ex56
removed the dead links to threads.html
send petsc repository messages to Karl no longer from his mcs address
...

Conflicts:
config/builder.py
src/dm/impls/plex/examples/tests/ex3.c
src/dm/impls/plex/examples/tests/makefile
src/dm/impls/plex/makefile

show more ...

Merge branch 'master' into knepley/feature-plex-point-location

* master: (202 commits)
Bib: Added ref
fix makefile format error with 'all-legacy'
Increase patchlevel to 3.6.1
update some exa

Merge branch 'master' into knepley/feature-plex-point-location

* master: (202 commits)
Bib: Added ref
fix makefile format error with 'all-legacy'
Increase patchlevel to 3.6.1
update some example output files
Doc+Vec: Fixed VecSetValuesBlocked() Reported by: Shao-Ching Huang <huangsc@gmail.com>
clarify that PETSc's "thread safe" code only works with OpenMP
allow adjoint integration over split time intervals
DT: Goddman it. Why are people stepping on my commits
Bib: Added ref
fix the bug in the endpoint version of discrete adjoint theta and correct the affected example outputs
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
fix: buggy change in 33da457b864643cfcac2768e47197c286471a5a4
updated petscthreads.html with information on thread per object model
fix for ex56
removed the dead links to threads.html
send petsc repository messages to Karl no longer from his mcs address
Plex: Fix for all the hidebound morons living with c89. Good luck doing computational science in the 80s!
DT: Fix tanh-sinh - Removed wrong var
SNES: Added SNESQN types to Fortran Reported-by: Chris Bradley <c.bradley@auckland.ac.nz>
Bib: More refs
...

Conflicts:
config/builder.py

show more ...

PCBDDC: compute column sum of zero block

PCBDDC: add command line switch for detecting saddle point

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (111 commits)
add filters for website tests that removes warnings about MPI type tag mismatches due to const
remove

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (111 commits)
add filters for website tests that removes warnings about MPI type tag mismatches due to const
remove Peter from blame-email list.
don't warn about MPI const long * versus long long *
Configure: Added GMP an MPRF - These allow arbitrary precision arithmetic
SNES ex69: SolCx finally works - Fixed order of parameters - Made different versions of pointwise functions
SNES ex69: Fixed bug where SolKx and SolCx have different numbers of parameters
SNES ex69: Small bug fixes - Initialize coordinates for SolCx - Fixed declaration length
I believe this call on line 357: ierr = DMSetFromOptions(dau); should actually set `dak`. The DM `dau` is already set on line 349.
SNESSolve_KSPOnly() did not compute function at new solution unless monitoring is turned on but SNESSolve_FAS() internals assume that each subsnes solve it calls does compute the new function value so broke when -fas_levels_snes_monitor was not turned on but -fas_levels_snes_type kspmonitor was used
Have SNESVIGetInactiveSet return the *actual* current inactive set.
One now passes the nullspace in the matrix for a linear system as the first Amat argument to KSPSetOperators, SNESSetOperators, TSSetIJacobian() not the second pmat argument
Remove change documentation of Removed MatMFFDAddNullSpace() just use MatSetNullSpace() as it is now backported in maint (36.html)
MatMFFD doesn't need to have its own null space management cherry-picked from 39601f4 (barry/remove-matmffdaddnullspace/master)
Removed trailing '\\' in the authors line.
add note to PetscOptionsInsertFile() that it should not be used for PETSc wide options
Add Fortran stub for SNESSetUpdate()
improve manual page slightly for clarity
PCGASM.
Basic PCGASM description added to manual.
add DMDAVecGet/RestoreArrayReadF90()
...

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc

Merge branch 'stefano_zampini/matis-feature-localtoglobalmap-handling'

PCBDDC: raise an error if the preconditioning matrix is not square

MATIS: enable rectangular shape for MATIS

Merge branch 'master' into knepley/feature-snes-deflation

* master: (239 commits)
Bib: Moved entry
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in c

Merge branch 'master' into knepley/feature-snes-deflation

* master: (239 commits)
Bib: Moved entry
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure
release: set v3.6 strings
doc: add PCGAMGSetReuseInterpolation() to changes file
fix: workarround static initialization issue (perhaps compiler bug) on opensolaris c/complex build
fix: ex40,ex41 errors in src/ts/examples/tutorials/ on bsd with c++/complex Perhaps this code nees a separate section when 'enorm == 0) - as pow(0,-0.33) = inf.
manual: fix format
more changes..
PCGASM/MatGetSubMatricesMPI() changes in dev.html.
changes: more updates
Add size checks for static arrays of assumed max size [-Warray-bounds]
Doc: Update to Changes
configure: default to useThreads= 0as the useThreads=1 code currently prints output in jumbled order [producing hard to read configure.log]
Delete PETSc viewer attribute if none are left on the linked list instead of putting an empty attribute
MatGetSubMatricesMPI_MPIAIJ(): guard for an empty B-matrix when delogging memory in MatSetSeqMats().
PCBDDC: fix typos in the specification of MPI operation
MPIU_SUM was used with MPIU_INT when it should only be used with MPIU_REAL or MPIU_SCALAR
fixed bad PetscPrintf() formatting in two examples
Plex: In pointwise callbacks, use embedding dimension
updates for generating Apple installations due to changes in locations of include files
...

show more ...

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (219 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure

Merge branch 'master' into mesgarnejad/feature-plex-globaltonatural

* master: (219 commits)
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html - enable cuda in configure
release: set v3.6 strings
doc: add PCGAMGSetReuseInterpolation() to changes file
fix: workarround static initialization issue (perhaps compiler bug) on opensolaris c/complex build
fix: ex40,ex41 errors in src/ts/examples/tutorials/ on bsd with c++/complex Perhaps this code nees a separate section when 'enorm == 0) - as pow(0,-0.33) = inf.
manual: fix format
more changes..
PCGASM/MatGetSubMatricesMPI() changes in dev.html.
changes: more updates
Add size checks for static arrays of assumed max size [-Warray-bounds]
Doc: Update to Changes
configure: default to useThreads= 0as the useThreads=1 code currently prints output in jumbled order [producing hard to read configure.log]
Delete PETSc viewer attribute if none are left on the linked list instead of putting an empty attribute
MatGetSubMatricesMPI_MPIAIJ(): guard for an empty B-matrix when delogging memory in MatSetSeqMats().
PCBDDC: fix typos in the specification of MPI operation
MPIU_SUM was used with MPIU_INT when it should only be used with MPIU_REAL or MPIU_SCALAR
fixed bad PetscPrintf() formatting in two examples
Plex: In pointwise callbacks, use embedding dimension
updates for generating Apple installations due to changes in locations of include files
configure: do not fail when checking MPICH version [similar to 2cadc8fc42fbad3c229d43b955581b8d29e02ad6 ]
...

show more ...

Merge branch 'master' of bitbucket.org:petsc/petsc

Merge remote-tracking branch 'petsc/master' into matpartitioning-hierach

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

Merge branch 'stefano_zampini/pcbddc-adaptive-selection'

Conflicts:
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcscalingbasic.c

PCBDDC: updated man pages

PCBDDC: added option to aggregate the adjacency between subdomains on a subset of the active processes

non contiguos option removed from MatISGetSubassemblingPattern

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran inclu

Merge branch 'master' into knepley/solkx

* master: (6933 commits)
Bib: Added reference
Includegraph: Updated for compatibility with new folder layout.
Add SNESLineSearchReason to fortran includes Also truncate names to fit 32 char fortran limit
fix compile warnings /usr/home/balay/petsc.clone-2/src/vec/vec/interface/vector.c:1944: warning: division by zero in '1.0e+0 / 0.' /usr/home/balay/petsc.clone-2/src/ksp/ksp/interface/itfunc.c:508: warning: 'vec_rhs' may be used uninitialized in this function
updated output for new trust region initial size
initial size of the trust region is set as a percent of the norm of the initial guess, it should not be related to the initial norm of the function (which is kind of nuts).
The default bounds for SNESVISetVariableBounds() in the manual page were reversed from correct values
Bib: Added refs
Bib: Update MPICH webpage
Plex: Forgot to close file
removed nonexistent ex46f from makefile test rule
Revert "fixed bad merge into master"
fixed bad merge into master
fix memory leak in mkl_pardiso fix some formatting in mkl_pardiso code
mv runex111 from TESTEXAMPLES_C to TESTEXAMPLES_DATAFILESPATH
Mat+Doc: More Fortran docs
Plex ex3: Run tests with the correct number of field components
Bib: Added TetGen cite
parmetis: make sure 'ldd libparmetis.so' points to the correct libmetis.so that its linked with.
small fix as reported by the nightly test
...

Conflicts:
config/builder.py

show more ...