Merge branch 'master' into knepley/feature-label-empty-stratum

* master: (151 commits)
moved PCTELESCOPE comment into users manual since no one will see the comment
message in makefile about hyp

Merge branch 'master' into knepley/feature-label-empty-stratum

* master: (151 commits)
moved PCTELESCOPE comment into users manual since no one will see the comment
message in makefile about hypre bugs causing test crashes
remove use of some unitialized variables, remove use of "dummy" entries in Vec and IS since Dave May thought they could not be of size zero. Now three telescope tests give no valgrind errors
code sent unitialized data via MPI and allocated too large space for creating indices
Force generateetags.py to use /usr/local/bin/ctags before using ctags because the ctags provided with Apple XCode gets confused by function usage in .tex files and does not make proper tags for TSSolve() and many many others
Optimized version of snes/examples/tutorials/ex1f.F produced incorrect results when lx_v() and lf_v() were declared as of dimension 1 with optimzed builds only. Cannot explain exactly why this happens but since we are using VecGetArray() from Fortran we are gambling with the compiler regardless
make the PetscRandom object used in GAMG to be a public object that can be controlled via the options database. This allows the same example to produce the same results with real and complex numbers by using the option -pc_gamg_random_no_imaginary_part
added parallel tester for examples. Crude but it reduces a basic alltests on my laptop from 25 minutes to 7 1/2 minutes.
remove unneccessary if
fix the problem for ex9 which involves multiple TSSolve calls
shortened the time of several test examples
decrease dramatically the time used by several test examples that take more than a minute on may machine with debug mode
Some of the vectors created in MatNullSpaceCreateRigidBody() were not normalized to 1
fix strange PCTelescope struct name of P which caused duplicate symbol error on some systems cleanup use of PetscMalloc and PetscFree for telescope for modern PETSc standards replace semiredundat to telescope in makefile that runs test examples
SNESUpdateCheckJacobian: make relative norm output syntactically correct
DMCompositeGetGlobalISs: document Fortran interface
fix warnings on bsd. Also fix usage of PetscFunctionBegin/Return
Making an even bigger mess, I cannot explain why this is not in the branch
hypre: update with "Bug fix for Parasails; src/ksp/ksp/examples/tutorials/ex10.c runex10_16"
decrease duration of ex22 by factor of 10 to decrease testing time
...

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

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of

Merge branch 'barry/fix-nonew-notcollective/maint' into jed/mat-assembly-perf

Jed this fucking 6+ month running thing in next that is not in master is a royal fucking pain in the ass and
an abuse of git. I think it is not justified to ever have anything in next for more than a few weeks at most.
Either take it out of next if it is broken or put it in master if it is not broken. Hanging around in next but
not master for ever is not a good policy.

show more ...

Merge branch 'matpartitioning-hierach' into pcgasm-increaseoverlap

Got an example ( /src/mat/examples/tests) from matpartitioning-hierach

Conflicts:
src/mat/impls/adj/mpi/mpiadj.c

Merge branch 'mat-increaseoverlap' of ssh://bitbucket.org/Fande-Kong/petsc into pr332/master/Fande-Kong/mat-increaseoverlap-scalable

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

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

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

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

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' changed to PetscMPIInt.

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

show more ...

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

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
thi

Merge branch 'master' into tisaac/dmforest

If we're going to start developing on this again, the recent changes to the
location of configuration files will make it annoying to switch to and from
this branch. Matt's and my fixes in DMPlex will also be relevant.

* master: (1557 commits)
MatLoad_MPIDense() always allocated array space even if it was previously provided
Clarified lumping options for PCFieldSplitSchurPreType SELFP.
fixed up PetscInfo in GAMG
TSErrorWeightedNorm[2|Infinity]: API changes
MatLoad_SeqDense() always allocated new space for the data even if the user already provided it
Fixed DMCompositeGetAccessArray() to support readonly vectors
Fix obviously-wrong aspects of documentation on SNESVINEWTONRSLS and SNESVINEWTONSSLS.
fixed example to compile - demonstrates a bug
var name 'len' was F90 key word
For user provided Jacobian functions they should always assembly into the second matrix, the one from which the preconditioner is built because the first may be used as matrix-free with -snes_mf_operator
Suggest --with-64-bit-indices if users mistakenly use --with-64bit-indices
Fix configure logging of integer size.
missing underscore in #define of vecgetarrayread translation of fortran stub for VECGETARRAYREAD
updated test output files to reflect differences introduced previously in PCView output
-ksp_type preonly does not support -ksp_monitor_short so remove from test line to prevent extranious warning output
missing CHKERRQ; and ugly double declaration of ierr
Plex ex3: For want of a semicolon, the kingdom was lost
MatSetValuesMPISBAIJ will now print correct global indices when error generated due to new nonzero location
completed basic use of MatSetValuesBlocked_SeqBAIJ_Inlined() by MatSetValuesBlocked_MPIBAIJ() to allow proper global indices being printed when error generated by new nonzero locations.
Added inline version of MatSetValuesBlocked_SeqBAIJ() to be used by MatSetValuesBlocked_MPIBAIJ() to print error message for nonzero block with global indices
...

show more ...

Merge branch 'master' into shri/ts-is-for-differential-variables

Conflicts:
include/petsc-private/tsimpl.h
src/ts/impls/eimex/eimex.c
src/ts/interface/ts.c

Merge branch 'master' into jed/chem
Note merging a branch with good stuff in it for a year and a half is insane!

Conflicts:
src/ts/examples/tutorials/makefile
src/ts/interface/ts.c

MatStash BTS: add MAT_SUBSET_OFF_PROC_ENTRIES, impl with neighbor-only comm

Add MatOption to imprint a communication pattern on the first assembly
and reuse on subsequent assemblies. Messages of si

MatStash BTS: add MAT_SUBSET_OFF_PROC_ENTRIES, impl with neighbor-only comm

Add MatOption to imprint a communication pattern on the first assembly
and reuse on subsequent assemblies. Messages of size 0 are sent if no
values need to be communicated to that rank. Since there are no
collectives (even Ibarrier) and no separate headers, this should be more
efficient for large process counts. The BTS implementation must
currently be activated via the -matstash_bts option.

show more ...

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c

Merge branch 'master' into barry/saws-options-merge-master

Totally confusing manual merge, could be wrong

Conflicts:
include/petscoptions.h
src/dm/impls/da/dacreate.c
src/sys/objects/aoptions.c
src/sys/threadcomm/interface/threadcomm.c

show more ...

Merge branch 'master' into pr252/dmay/PetscViewerBinarySkipHeaderVecViewLoad
Manual merge of makefile

Conflicts:
src/sys/classes/viewer/examples/tutorials/makefile

Merge branch 'master' into mlange/refactor-plex-distribute

Conflicts:
src/dm/impls/plex/examples/tests/output/ex12_1.out
src/dm/impls/plex/examples/tests/output/ex12_2.out

Merge branch 'hzhang/mat-redundant_sbaij'

rm MatGetRedundantMatrix_xxx()

mv struct 'Mat_Redundant' to '_p_Mat’ for simplify codes

compile: fixes for windows/solaris

add MatGetRedundantMatrix_MPIBAIJ() -- almost same as MatGetRedundantMatrix_MPIAIJ()

Merge branch 'maint'

Merge branch 'hzhang/mataxpy-bugfix' into maint

reduce code duplication for MatAXPYGetPreallocation_MPIXXX()

bugfix for MatAXPY_ for mpibaij and mpisbaij

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]

Merge branch 'master' into prbrune/snes-fascontinuation

* master: (855 commits)
fix for moab example had wrong blocksize of mapping
docs:remove src/docs/doxygen/manual [per Barry's instructions]
gfortran: remove 132 char limit from f90 free-form
SNES: Set up matrix-free *after* copying the user Jacobian to snes->pc
DMPlex-Gmsh: Use trueNumCells when building the coordinate section
Sys: Got guard logic reversed
Sys: The #define guard for PetscMalloc*() must match the one for PetscFree*()
PetscFV: I hate fucking complex
removed ltog scatter from DMDA since functionality can be obtained with gtol scatter and less memory
fix test of DMLocalToGlobalBegin/End and add actual test cases
improve object memory logging for DA, ISLocalToGlobal and ISBLOCK
lower high water mark of constructing parallel scatter by allocating only as much room in svalues as needed
owner in VecScatter construction is the rank, therefor can use PetscMPIInt to hold it do not PetscInt
another reuse of allocated array in 2d/3d DMDA create
improved high water mark of DMDASetups by implementing and using PETSC_USE_POINTER for ISCreateBlock()
PetscFV: Another complex fix
DMPlex+TS: Fixing types
ISIdentity() was wrong for block IS! Indices are NOT scaled by bs
PetscFV: Straighten out complex types
TS ex11: Fix casts and remove unused variable
...

show more ...