Fix MatMult/MatMultTranspose SeqAIJCUSPARSE

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
ini

Merge branch 'master' into knepley/feature-plex-sphere-mesh

* master: (2633 commits)
Bib: Added ref
SNES ex56: Update test output
add a citation to petsc.bib
Plex: Initialize variables
initialize ap as complained by the nightly tests
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
SNES ex12: Mix misnaming Nf --> Nc
update tests
add tests
add MatView_SeqXAIJ_ASCII_structonly(); enable A->structure_only in MatSetValuesBlocked_SeqBAIJ
typo
fill-in some references
Bib: Fixed ref
Plex: Added PLY to file support
Plex: Adding PLY support
PetscDS: Turn off Fortran bindings for Constants
Plex tests: Updated to new point function API
SNES ex56: Moved tests to new format
DM: Allow field projection functions to take an array of scalars - This will allow us to pass in user data to kernels without void *, more efficiency than aux fields
...

Conflicts:
config/builder.py
include/petscdmplex.h
src/dm/impls/plex/examples/tests/ex1.c
src/dm/impls/plex/plexcreate.c
src/dm/impls/plex/plexgmsh.c
src/snes/examples/tutorials/ex63.c

show more ...

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes

Merge branch 'master' into knepley/feature-plex-pragmatic-parallel

* master: (1379 commits)
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes
add help menu for MatCreateSubMatricesMPI()
rename struct Mat_SubMat -> Mat_SubSppt
add a little information to the manual about using tchem
rename MatCreateSubMatrixComposite -> MatCreateSubMatrixVirtual; Mat_SubMatrix -> Mat_SubVirtual
need to protect the \ passed on the makefile on my Mac
Bib: Updated ref
update src/mat/examples/tests/output/ex23_#.out
rename MatGetSubMat -> MatCreateSubMat
Mark MatCreateMPIAIJWithSplitArrays() and MatMPIAIJSetPreallocationCSR() with MAT_NO_OFF_PROC_ENTRIES
superlu: libray is installed in PREFIX/lib64 - fix this to use PREFIX/lib
backport fixes cc105cc3 2b35c90a - required for 3bf428d - that was backported from master - 80c96bb
'echo -n' fails on opensolaris/osx - so avoid using it.
fix diffs in nightlybuilds due to options in ~/.petscrc
add Kelley03 to src/docs/tex/petsc.bib
build: determine moudles, mpi.h used during buildtime and print in 'info'
fix analyzer warnings
fix build failure on opensolaris - where -lm from flibs is removed due to -lm in cxxlibs. But then - only flib gets used for some configure tests [blaslapack, parms]
config: prefer __restrict over C99 restrict keyword
...

show more ...

Merge branch 'barry/remove-funct'

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-t

Remove the use and definition of __FUNCT__ throughout the code

Since all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETSc
Time: 1.5 hours
Thanks-to: Andreas Mang <andreas@ices.utexas.edu>

show more ...

Merge branch 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/featu

Merge branch 'stefano_zampini/feature-convert-matnest-matis' into stefano_zampini/feature-pcbddc-saddlepoint

I have to merge it here, since there's an API addition for MatIS in stefano_zampini/feature-pcbddc-saddlepoint
that prevents me to smoothly merge stefano_zampini/feature-convert-matnest-matis into next.

This will also facilitate merging stefano_zampini/feature-pcbddc-saddlepoint into master

Conflicts:
src/docs/website/documentation/changes/dev.html

show more ...

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in m

Merge branch 'master' into barry/fix-successful-numerical-factorization-after-failure/master

The maint patch for the given problem cannot merge cleanly into next due to changes
in the same code in master (the introduction of MatFactorGetError()) therefor need
to merge in master and use the new construct before merging to next

show more ...

Merge branch 'master' into psanan/pc-telescope-nearnullspace

Merge branch 'knepley/fix-plex-output' into mark/snes-ex56

* knepley/fix-plex-output:
Quadrature: Fix creation init
Plex: Fix output optimization - We obviously have no test for this logic
All

Merge branch 'knepley/fix-plex-output' into mark/snes-ex56

* knepley/fix-plex-output:
Quadrature: Fix creation init
Plex: Fix output optimization - We obviously have no test for this logic
Allow the Richardson scale factor to be reset at each iteration
prevent a possible segfault when m is 0 (or less)
when configured --with-errorchecking=0, CHKERRV is missing
essl: fix PetscOptionsGetReal() usage
Fix code comment. Routine actually does not add arguments to database it just reads them in
Fix performance issue in MatAIJCUSPARSE after 34d6c7a
Plex ex6: fix PetscInt print codes
DMPlex: fix PetscPartitionerPartition() to deal with overlap

show more ...

Merge branch 'master' into mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increa

Merge branch 'master' into mark/snes-ex56

* master: (93 commits)
Web: Added a tutorial
typo
rename X -> U
sowing: update to version that creates manpage URLs for 'enum' types/values
Increase patchlevel to 3.7.3
move runex3_3 to NOTSINGLE [ partial backport 9ed0bd76d5c329eb91401ebabfbd851446214e20 ]
DMPlex: local coordinate vector should be a sequential vector
DM: Name coordinate vectors in subdomain and restrict hooks
tests: disable a couple of TS examples from single precision tests
reduce TS ex12_2 output
fix PetscFormatRealArray() to cast PetscReal in print statement to double so works for __float128 precision
Vec: enable Fortran interfaces for VecGetLocalVectorRead and VecRestoreLocalVectorRead
Bib: Added refs
PCShell: man page fixes
update handling of latex in example code to use https
fix typo in isltog comment
test: Avoid conversion of object <-> function pointer type [-Wpedantic]
fortran: enable interface for VecGetLocalVector(), VecRestoreLocalVector()
update citation for cmz2016
PCBDDC: silent compilation warnings from nightlybuilds
...

# Please enter a commit message to explain why this merge is necessary,
# especially if it merges an updated upstream into a topic branch.
#
# Lines starting with '#' will be ignored, and an empty message aborts
# the commit.

show more ...

Merge branch 'jose/fix-cusparse-flag'

* jose/fix-cusparse-flag:
In MatSeqAIJCUSPARSECopyToGPU the case of simplified copy must also set the GPU consistency flag
PR #514

Fix performance issue in MatAIJCUSPARSE after 34d6c7a

In MatSeqAIJCUSPARSECopyToGPU the case of simplified copy must also
set the GPU consistency flag

Merge branch 'jose/veccuda-tests'

Conflicts:
config/PETSc/Regression.py

* jose/veccuda-tests:
This is an attempt to modify configure to run regression tests with VECCUDA.
The tests for CUSP ar

Merge branch 'jose/veccuda-tests'

Conflicts:
config/PETSc/Regression.py

* jose/veccuda-tests:
This is an attempt to modify configure to run regression tests with VECCUDA.
The tests for CUSP are left unchanged.
The new tests include CUSPARSE and work with real and complex scalars.
See PR #490 for discussion.

show more ...

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

Conflicts:
src/docs/website/documentation/changes/dev.html
src/ksp/ksp/examples/tutorials/ex59.c
src/ksp/pc/impls/bddc/bddc.c
src/ksp/pc/impls/bddc/bddcfetidp.c
src/ksp/pc/impls/bddc/bddcprivate.c
src/ksp/pc/impls/bddc/bddcschurs.c
src/mat/impls/is/matis.c
src/mat/impls/is/matis.h

show more ...

Miscellaneous changes in AIJCUSPARSE code

- Some operations need synchronization of data from CPU to GPU, e.g. MatScale
- Add assembly operation, so that lvec is set to type CUDA
- In MatSeqAIJCUSPA

Miscellaneous changes in AIJCUSPARSE code

- Some operations need synchronization of data from CPU to GPU, e.g. MatScale
- Add assembly operation, so that lvec is set to type CUDA
- In MatSeqAIJCUSPARSECopyToGPU do not recompute indices every time
- Operations MatMultAdd* need to get the array with ReadWrite
- In MatScale, do not change coherence flag if alpha=1.0
- VecDuplicate_MPICUDA must make sure that the GPU array exists
- Fix null pointer in VecScatterCUDAIndicesCreate_PtoP

show more ...

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary c

Merge branch 'master' into tisaac/ts-ex11-amr

* master: (141 commits)
chombo: fix dirname in tarball and update to use snapshot with gcc-6.1 compile fix
petsc4py: fix tarball url [for arbitary commitid] - and dirname in tarball
Plex: Moved ReconstructGradients_Internal() and ApplyLimiter_Internal() to Plex from SNES
CUDA: Added external declaration for cuBLAS handle.
docs: fix reference to KSP_Create in parent commit
Add information on PETSC_EXTERN and PETSC_INTERN to developers manual
Plex: Symbol is shared between DM and SNES, so must be EXTERN
Plex: Remove unused variables
Configure: Fix check for spaces in PETSC_DIR
Plex: Updated examples to new DMPlexCreateBoxMesh()
sick of failure due to incorrect master branch
SNES ex12: Fixed unstructured GMG test
SNES ex12: Added another unstructured GMG test
petsc4py: use the snapshot with the fix for dmplex change
Examples that require Fortran to compile but were listed in tests that did not require Fortran
Plex ex19: Added DMSetFromOptions()
Plex: Store a flag indicating whether the boundary should change during remeshing - Set using -dm_plex_remesh_bd - Changed Remesh_Internal() signature
Always coarsen boundaries.
Plex: Fixed option name
Plex: Protect unused variables
...

Conflicts:
src/dm/impls/plex/makefile
src/snes/utils/dmplexsnes.c

show more ...

Merge branch 'jose/veccuda-witharray'

jose/veccuda-witharray:
The aim of this PR is to be able to create a VecCUDA vector
where the user provides the array on the GPU side.

Also, this PR (#466) inc

Merge branch 'jose/veccuda-witharray'

jose/veccuda-witharray:
The aim of this PR is to be able to create a VecCUDA vector
where the user provides the array on the GPU side.

Also, this PR (#466) includes a few fixes:
- Fix use of not initialized variable
- Fix CUDA_VERSION check in configure
- Increase the required CUDA_VERSION >= 7.5 for complex builds
(CUDA version 7.0 contains thrust/complex.h but it produces compilation errors)
- Stop using CUDA_SUCCESS and CUBLAS_STATUS_SUCCESS in CHKERRCUDA macro (otherwise the user must include cuda.h)

show more ...

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()

Merge branch 'master' into knepley/feature-parallel-mesh-load

* master: (180 commits)
Bib: Added ref
-random_view caused error because PetscRandomView() was not provided in PetscRandomCreate()
Update the xSDKTrilinos test message to use interactive session and ctest
gcc6: fix comple warnings with -Wmisleading-indentation
fix compile warning
bugfix for MatTransposeMatMultSymbolic_MPIAIJ_MPIAIJ_nonscalable()
minor formatting cleanup of gamg.c
add ml example to tutorials page for xsdk
bugfix of MatDuplicate for mpiaij matrix product Reported-by: Miorelli, Federico <Federico.Miorelli@cgg.com>
The deprecated -log_summary_exclude was still being used when it should have been -log_exclude
cpardiso: requires mkl_blacs so check for that. Also common code to detect MKL for both pardiso and cpardiso
cpardiso: fix segv and memory leaks
resume MatGetSubMatrix_SeqSBAIJ() because it cannot share routine with BAIJ due to cast
fix bad relative links
replace MatSeqBAIJSetPreallocation_SeqBAIJ with MatSeqBAIJSetPreallocation to skip this call when input mat is sbaij type
fix a variety of set but not used variables in examples detected by clang static analyzer
removed declared but not used variable
ANL Report formatting requires listing authors by affiliation
DMP4est: uncrustify
DMP4est: remove memset()
...

show more ...

Add const to VecCUDAGetArrayRead() to match with VecGetArrayRead()

Fix use of not initialized variable

Merge branch 'master' into stefano_zampini/feature-pcbddc-saddlepoint

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

Merge branch 'jose/veccuda-complex'

jose/veccuda-complex:
Adds/reenables complex arithmetic for VECCUDA.

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make

Merge branch 'tisaac/sf-fix-multi-sf-leaves' into tisaac/dmp4est-feature-injection

* tisaac/sf-fix-multi-sf-leaves: (174 commits)
PetscSF: fix PetscSFGetMultiSF() for sparse leaves
PetscSF: make regression test that fails
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release version
Updates for xSDKTrilinos builds
Added --download-xsdktrilinos
Communicator passed to VecCreateSeq should be PETSC_COMM_SELF
PCMG man page: fix typo in options name
Bib: Added refs
DMLabel: Removed unused variable
Metis tries to use backtrace which requires -lexecinfo on some systems such as freebsd Remove -lexecinfo from freebsd configure files since it is only needed by metis
add alternative output that occurs on some systems due to use of random number generator in partitioner
Plex: When creating cohesive submeshes, fix memory leak
Plex ex11: Updated output for new distribution algorithm
Plex: When explicitly distributing the depth label, we must recreate empty strata
Plex: When partitioning, the cell numbering should include hybrid points
Plex: Allow cell and vertex numberings to include hybrid points
need to show diffs of ex22 if it fails
Support two different output for test example, since due to roundoff this example can produce slightly different convergence history
...

show more ...

Merge branch 'master' into sanderarens/fix-plex-neumann-bc

* master: (28 commits)
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release

Merge branch 'master' into sanderarens/fix-plex-neumann-bc

* master: (28 commits)
few more clang static analyzer fixes
many issues detected by clang static analyzer
Update hypre to its release version
Updates for xSDKTrilinos builds
Added --download-xsdktrilinos
DMLabel: Removed unused variable
Plex: When creating cohesive submeshes, fix memory leak
Plex ex11: Updated output for new distribution algorithm
Plex: When explicitly distributing the depth label, we must recreate empty strata
Plex: When partitioning, the cell numbering should include hybrid points
Plex: Allow cell and vertex numberings to include hybrid points
Plex ex1: Fixed output for empty strata
DMLabel: Added DMLabelHasStratum()
Minor changes and fix to bugs introduced in 82f73ecaa
Modified makefile to fix aijcusparse tests
Add missing code for the case PETSC_HAVE_VECCUDA
Fix bugs introduced in 82f73ecaa
Use VECCUDA with MATAIJCUSPARSE
CUDA: Fixed visibility and disabled Fortran bindings for VecScatters on GPU.
Replace kernels with thrust
...

show more ...

Fix VECCUDA compilation with complex and single precision