Some comments from Barry.
Initial support for SeqSBAIJ matrices.
Merge branch 'stefano_zampini/GPU-matdensecuda' [PR #1911]* stefano_zampini/GPU-matdensecuda:GPU: Initial implementation for SeqDense class on GPUs.
MATMUMPS: Schur complement is updated on the CPU
Merge branch 'jczhang/sf-more-opts' [PR #1567]* jczhang/sf-more-opts:Add more optimizations in SF and use it as the default for VecScatter.
Merge branch 'barry/2019-07-21/scalapack-mumps-parallel-make' [PR #1896]* barry/2019-07-21/scalapack-mumps-parallel-make:Parallel make for --download-scalapack --download-mumps
Mumps no Cholesky with Hermitian, no PetscMallocValidate() unless -malloc_*, example used PetscAbsScalar() on realMumps has no support for Cholesky for Hermitian matrix so document and error check
Mumps no Cholesky with Hermitian, no PetscMallocValidate() unless -malloc_*, example used PetscAbsScalar() on realMumps has no support for Cholesky for Hermitian matrix so document and error checkDocumentation says CHKMEMQ() and PetscMallocValidate() with -malloc_debug or -malloc_test but in fact it still checked memory even without those being set (slowing code)Example used PetscAbsScalar() on real values, also had incorrect indentation formatingCommit-type: testing-fix, style-fix, documentation, example
show more ...
Parallel make for --download-scalapack --download-mumpsRequired added a new Makefile.parallel and Makefile.objs to the scalapack git repository to handle it tag v2.0.2-pFixed spelling of controll
Parallel make for --download-scalapack --download-mumpsRequired added a new Makefile.parallel and Makefile.objs to the scalapack git repository to handle it tag v2.0.2-pFixed spelling of controller in mpishm.cRefactoring of part of MATMUMPS manual page to clafiy and simplify the discussion of hybrid MPI/OpenMP, note that itsupports both modesl for MPI/OpenMPCommit-type: configuration
MUMPS: Add comments for a TODO optimization
MUMPS: Use VecScatterCreateToZero to simplify RHS scatter code
Add PetscArraycmp(), PetscArraycpy(), PetscArrayzero(), PetscArraymove(), update code to use new formsReplace a couple multi-mallocs with a single mallocReplace backwards loops in MatSetValues_*()
Add PetscArraycmp(), PetscArraycpy(), PetscArrayzero(), PetscArraymove(), update code to use new formsReplace a couple multi-mallocs with a single mallocReplace backwards loops in MatSetValues_*() with PetscArraymove()Commit-type: style-fix, feature
MUMPS: activate Block Low-Rank parametersVersion 5.2.0 of MUMPS features BLR (Block Low-Rank) factorization.Extended ICNTL, INFO and INFOG to sizes 60, 80 and 80.Added description and control par
MUMPS: activate Block Low-Rank parametersVersion 5.2.0 of MUMPS features BLR (Block Low-Rank) factorization.Extended ICNTL, INFO and INFOG to sizes 60, 80 and 80.Added description and control parameters (ICNTL) 35, 36 and 38.
Use MatDenseGetArrayRead() rather than MatDenseGetArray() when possibleCommit-type: style-fix, feature
MUMPS: Fixed uninitialized variable warning.Follow-up to PR #1592.
rm malloc of iidx array of size nrhs*M for sparse rhs matrix.Requested by Jaysaval, Piyoosh <Piyoosh.Jaysaval@pnnl.gov>
add support of MatMatTransposeSolve() for mumps CholeskyRequested-by: Jaysaval, Piyoosh <Piyoosh.Jaysaval@pnnl.gov>
Revert the function name where it is referenced
Add the missed VecRestoreArrayVecGetArray() must have a matching VecRestoreArray()
Improve the error propagation between and within KSPs and PCs- Added KSPCheckSolve() to be called consistently by PCs to set failures if the inner solvers fail.- Added use of KSPCheckNorm() and KS
Improve the error propagation between and within KSPs and PCs- Added KSPCheckSolve() to be called consistently by PCs to set failures if the inner solvers fail.- Added use of KSPCheckNorm() and KSPCheckDot() to all the Krylov methods to propagate errors from PCs to the KSPs.Commit-type: feature
Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do.Changed old VecScatterCreate() to VecScatterCreateWithData()Completed basic refactorizati
Changed VecScatterCreate() to take a communicator and return scatter as other first class PETSc objects do.Changed old VecScatterCreate() to VecScatterCreateWithData()Completed basic refactorization of VecScatter code, ready for new implementations to be added.(There is still some code duplication between the VecScatter MPI1 and MPI3 code but it can be fixed at anytime without requiring changes to the VecScatter API.Part 5 of a refactorization of VecScatter to be treated as a first class PETSc objectCommit-type: style-fix
Replace old ifdefs with new macro PETSC_HAVE_OPENMP_SUPPORT
Allow -mat_mumps_use_omp_threads to have a default value
Update MUMPS document, change -mumps_omp_num_threads to -mat_mumps_use_omp_threads
fix nightlybuild issues - disable hypre [for openmp] and fix uninitialized variables
Use mmap to allocate the shared pthread barrier if available
12345678910>>...25