MatIS: added implementation for PtAP operation

supported P are {Seq|MPI}(B)AIJ and HYPRE
Interface for {MPI|Seq}BAIJ and HYPRE code works, however:
- local PtAP ops are missing for BAIJ
- MatCreateS

MatIS: added implementation for PtAP operation

supported P are {Seq|MPI}(B)AIJ and HYPRE
Interface for {MPI|Seq}BAIJ and HYPRE code works, however:
- local PtAP ops are missing for BAIJ
- MatCreateSubMatr* ops are missing for HYPRE (will be added in a separate branch)

If requested, carry over the l2l representation of P by attaching it to the resulting mat
This is needed by the MFEM constructor for BDDC

added few more tests to ex23.c
L2L storage is disabled because of a bug in MatDestroySubMatrices that will be fixed in a separate branch

show more ...

Merge branch 'master' into knepley/feature-error-logging

* master: (221 commits)
dashboard: count and list example compile failures like:
fix fix-pr879-jczhang-optimize-vecscatter-multi-copy
a

Merge branch 'master' into knepley/feature-error-logging

* master: (221 commits)
dashboard: count and list example compile failures like:
fix fix-pr879-jczhang-optimize-vecscatter-multi-copy
add manual page cross links for SNESLineSearchSetPostCheck() and related routines
fix help manual of MatCreateSubmatrices() Reported-by: Tim Steinhoff <kandanovian@gmail.com>
Add manual pages for PetscRealPart() and PetscImaginaryPart()
fix compile error with missing MPI_C_DOBULE_COMPLEX
fix typo Reported-by: Febrian Setianto <febrian.setianto@gmail.com>
plex tutorial ex5: updated test requirements
DMPlexLoad_HDF5_Xdmf_Internal: fix for complex PetscScalar.
DMPlexLoad_HDF5_Xdmf_Internal: fix for 64-bit PetscInt.
DMPlexBuildFromCellList_Internal: use DMPlexInvertCell_Internal (with PetscInt* arg)
petscviewer.h: removed notice which is no longer relevant
DMPlexCreateFromFile: document option
check only if PETSC_USE_DEBUG
DMPlexBuildFromCellList_Internal optionally inverts cells.
ex5: new test 4 with -compare without repartitioning
ex5: new test 3 using -compare
ex5: HDF5 write+read moved to separate function
ex5: option -format processed in ProcessOptions()
ex5: added missing hdf5 requirement
...

Conflicts:
include/petsclog.h

show more ...

Merge branch 'maint'

Merge remote-tracking branch 'origin/hongzh/guard-fma-intrinsics/maint' into maint

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead

Merge branch 'master' into knepley/feature-snes-ex13-spectral

* master: (2350 commits)
migrate warning fix from next-tmp
Revert "Remove dead code."
clarifying comment in the code
Remove dead code.
Add JuMP citation
Copyedits and updates to developers manual
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build
Plex+VTK: Fix field type determination
Fix bug in name cleanup
Bib: Added ref
Plex+VTK: Fix field type determination
Plex: Updated test output - I am now hating this shit
forgot to save file after moving prototypes to petscsys.h
update Apple example project settings per Xcode
Two functions should actually be public, used, for example from Apple MacOS and iOS programs
includes should always be included relative to include directory
Apple removed old locations of docsets
SNES ex12: Updated test output
...

Conflicts:
src/snes/examples/tutorials/ex13.c

show more ...

Merge branch 'master' into knepley/fix-plex-residual-loop
Warning: This merge includes real changes

* master: (188 commits)
Add JuMP citation
Copyedits and updates to developers manual
Add mor

Merge branch 'master' into knepley/fix-plex-residual-loop
Warning: This merge includes real changes

* master: (188 commits)
Add JuMP citation
Copyedits and updates to developers manual
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build
Plex+VTK: Fix field type determination
Fix bug in name cleanup
Bib: Added ref
Plex+VTK: Fix field type determination
Plex: Updated test output - I am now hating this shit
forgot to save file after moving prototypes to petscsys.h
update Apple example project settings per Xcode
Two functions should actually be public, used, for example from Apple MacOS and iOS programs
includes should always be included relative to include directory
Apple removed old locations of docsets
SNES ex12: Updated test output
SNES ex62: Set name of solution vector - Update test output
revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html
release: set v3.9 strings
move output file to the correct location. i.e fix 305862f9cd0dc399542ae11dd744508e4b6ad2e3
...

Conflicts:
src/snes/utils/dmplexsnes.c
src/ts/utils/dmplexts.c

show more ...

Add guard macros for FMA intrinsics

mm256 FMA instructions do not belong to AVX2 (but AVX512 contains mm512 FMA). There are processors that supports AVX2 but no FMA, e.g. i7-4800MQ.

Funded-by:
Proj

Add guard macros for FMA intrinsics

mm256 FMA instructions do not belong to AVX2 (but AVX512 contains mm512 FMA). There are processors that supports AVX2 but no FMA, e.g. i7-4800MQ.

Funded-by:
Project: PETSc for Theta
Time: 0.5 hours
Reported-by: Randall Mackie <rlmackie862@gmail.com>
Thanks-to:

show more ...

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

# By Hong Zhang (2) and Jed Brown (1)
# Via Satish Balay
* 'master' of https://bitbucket.org/petsc/petsc:
Add more protection to intrinsi

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

# By Hong Zhang (2) and Jed Brown (1)
# Via Satish Balay
* 'master' of https://bitbucket.org/petsc/petsc:
Add more protection to intrinsics
MPI: protect arguments in MPI macros
Add the missing conditionals for 64-bit build

show more ...

Merge branch 'maint'

Merge remote-tracking branch 'origin/hongzh/fix-64bit-knl-build/maint' into maint

Add more protection to intrinsics

Funded-by:
Project: PETSc for Theta
Time: 0.5 hours
Reported-by: "Nourgaliev, Robert Nr" <nourgaliev1@llnl.gov>
Thanks-to:

Merge remote-tracking branch 'origin/knepley/fix-snes-ex69' into tisaac/feature-dmfield

* origin/knepley/fix-snes-ex69: (311 commits)
PetscFECreateDefault: fix custom Fortran interface to use MPI_

Merge remote-tracking branch 'origin/knepley/fix-snes-ex69' into tisaac/feature-dmfield

* origin/knepley/fix-snes-ex69: (311 commits)
PetscFECreateDefault: fix custom Fortran interface to use MPI_Comm
SNES ex69: Updated test output
SNES ex69: Removed all mantle stuff and change fields to constants - Also use dsymutil for executables
SNES ex69: Removed special partition
SNES ex69: Fixed null space creation - Null space locks its vectors
PetscFE: PetscFECreateDefault() now takes MPI_Comm instead of DM
cleanup very ugly PCGAMG manual page
minor clarifications for VecScatterCreate() manual page
gcov: PETSC_ARCH is required onlyfor -run_gcov - and not for -merge_gcov [where it doesn't exist]
make clearer in help statements what --with-external-packages-dir is for.
Website: Added news entry for PETSc User Meeting 2018 on main page.
CUDA: Further test output fixes for mat_tests-ex2
Turn off elemental test for parmetis running since this does not work on batch systems.
3.8 Changes: KSPCHEBYSHEV
CUDA: Fixed most test errors due to MatMultTranspose*
website: update urls
website: fix Jason's name to match pic. Also add Todd's pic
Dev manual: copyedit
DMPlexDistributeLabels: Propagate the output flag
cuda/cusp: update minimum supported versions of cuda to 7.5 and cusp to 0.5.0
...

show more ...

Merge branch 'pr862/dkokron/IntelIntrinsicsForBAIJ9by9MDotMaxpy/maint'

Merge branch 'barry/remove-mdiag-baij'

update #define conditions used to turn on Intel intrinsics

Commit-type: portability-fix, style-fix

Remove unused mdiag[] array from BAIJ that was filled but never used

Commit-type: style-fix
Reported-by: Xiangdong <epscodes@gmail.com>

Added case specific intrinsic kernels for VecMdot_Seq, VecMAXPY_Seq, and BAIJ bocksize 9 by 9

Merge commit 'd47bf9aaf1e5266cc3f0ff499f934c85788965a9' into jed/fix-matcreatempibaij/maint

Obtain changed (v3.5) handling of MatMPIBAIJSetPreallocationCSR to match
documentation.

Merge branch 'pr746/irinasok/feature-baijmkl/master'

* pr746/irinasok/feature-baijmkl/master:
Implementation of BAIJMKL matrix type similar to AIJMKL.
Supports MKL calls for MatMult (add and transpo

Merge branch 'pr746/irinasok/feature-baijmkl/master'

* pr746/irinasok/feature-baijmkl/master:
Implementation of BAIJMKL matrix type similar to AIJMKL.
Supports MKL calls for MatMult (add and transposed).
In case if used MKL version does not support mkl_sparse_optimize() (MKL 11.2 and earlier) PETSc default kernels are used.
Test added on configuration step to check if used version of MKL supports column major layout with zero based indexing.
All released versions of MKL do not support such cross case and MatSeqBAIJMKL_create_mkl_handle() converts matrix arrays to one based indexing before passing them to mkl routines.

show more ...

Added new matrix classes to include files
Added define MKL_SUPPORTS_BAIJ_ZERO_BASED to suuport different versions
of MKL (older verions support only one based indexing for column major
blocks in baij

Added new matrix classes to include files
Added define MKL_SUPPORTS_BAIJ_ZERO_BASED to suuport different versions
of MKL (older verions support only one based indexing for column major
blocks in baij format)

show more ...

Merge remote-tracking branch 'origin/hzhang/reorg-matdestroysubmats'

for baij and dummy matrices

Merge branch 'master' into knepley/fix-fe-default

* master: (297 commits)
Bib: Added refs and fixed one
Doc: Fixed formatting problems for MatCreateAIJ() manual page.
Add to PCSOR manual page

Merge branch 'master' into knepley/fix-fe-default

* master: (297 commits)
Bib: Added refs and fixed one
Doc: Fixed formatting problems for MatCreateAIJ() manual page.
Add to PCSOR manual page information about convergence testing when using KSPRICHARDSON
SNES ex5: Restore old output since this breaks an example I have
Noted that some Cray XC systems will use 'aprun' instead of 'srun'. Thanks to Herr Doktor Hong Zhang for this suggestion.
Added example configure script for Cray XC40 with KNL processors (e.g., NERSC Cori).
Fix possible sprintf -Wformat-overflow
PetscMalloc2-7: fix -Wint-in-bool-context for GCC 7
SNES ex5: make boundary condition handling consistent for symmetry
SNES ex5: use L2 norm and *_monitor_short for output
SNES ex5: the L^2 norm scales with the l^2 norm over sqrt(N)
SNES ex5: refactor MMS capability to avoid extreme code duplication
Documentation fix: typo in PetscLogStageRegister manual page.
nodes -> vertices
rename DMNetworkGetComponentTypeOffset -> DMNetworkGetComponentKeyOffset
rename DMNetworkGetConnectedNodes() -> DMNetworkGetConnectedVertices()
Jacobians for adjoints for explicit methods and constant Jacobians never got initialized
fix complex build error
fix build failure with fortran interface
fix compile error with cdc753b64c658d532d7813e942002e67025667a6
...

Conflicts:
src/snes/examples/tutorials/ex13.c
src/snes/examples/tutorials/ex62.c
src/snes/examples/tutorials/ex69.c

show more ...

MatCopy implementations need to update the state since values have changed

Commit-type: bug-fix
Development Tools: Vim, Emacs, Eclipse

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warni

Merge branch 'master' into stefano_zampini/feature-gmsh-periodic

* master: (1016 commits)
metis: setself.downloaddirnames for change in 6d90219d8fa4f8b2fd242cc6ff9008358d87c26f
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ may be used uninitialized in this function [-Wmaybe-uninitialized] ierr = DMInterpolationSetUp(interpolator, dm, pointsAllProcs);CHKERRQ(ierr);
migrate 'build requires: !mpiuni' to new format
migrate 'build requires: !complex' to new format
workarround 'division by zero' compiler warning on windows C:\cygwin64\home\petsc\PETSC~3.CLO\src\ts\examples\TUTORI~1\ex11.c(416): warning #39: division by zero if (uL->h < 0 || uR->h < 0) {for (i=0; i<1+dim; i++) flux[i] = 0./0.; return;} /* SETERRQ(PETSC_COMM_SELF,PETSC_ERR_ARG_OUTOFRANGE,"Reconstructed thickness is negative"); */
dos2unix to prevent warnings on opensolaris "/export/home/petsc/petsc.clone/src/ksp/ksp/examples/tests/ex53.c", line 6: warning: invalid white space character in directive
fix output file - corresponding to code changes in efd4aadf157bf1ba2d80c2be092fcf4247860003
disable test on windows
update log from opensolaris so that the diff is clean on linux aswell
update grep in sync with change in efd4aadf157bf1ba2d80c2be092fcf4247860003
Fixed code that did not correctly set current color for block insertion
Update output files that did not get fixed after previous merges
Change MPIUni to not immediately abort on operations it cannot perform, instead return error code
Add filter for example that can produce 1 different iteration count
Example add system calloc, bad formating of %d in print, commented out code
Always use (void) for functions that take no arguments, never ()
always need to caste PetscReal variables to (double) in print statements
updates to manual page for clarity
update help menu for MatRARt() -- MatRARt_mpiaij_mpiaij is supported now, but not efficient.
fix nightlybuild error with c1cbb1def3f8ddca04c3998b9c06618fd1ffef35 on arch-linux-IPL64 src/ksp/pc/examples/tests ex5_1
...

Conflicts:
include/petscdm.h
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/dm/impls/plex/plexhdf5.c

show more ...