Add subtype concept to SeqAIJ matrices

Introduced -mat_seqaij_type <seqaijperm, ...> MatSeqAIJSetType() to allow running code written for AIJ to work
with any subtype. See MatSeqAIJRegister() and Ma

Add subtype concept to SeqAIJ matrices

Introduced -mat_seqaij_type <seqaijperm, ...> MatSeqAIJSetType() to allow running code written for AIJ to work
with any subtype. See MatSeqAIJRegister() and MatSeqAIJRegisterAll() for registering new subtypes

Use

export PETSC_OPTIONS="-mat_seqaij_type seqaijperm -malloc_test"
make alltests

to run the entire test suite using an alternative subclass, in this case seqaijperm. This will produce output regarding the different name for
the matrix class but otherwise should produce no errors.

This was implemented to 1) allow easy testing of an MKL subtype and 2) to allow users to easily use the MKL subtype
without changing their code.

Also fixes several bugs in seqaijperm, if it isn't tested its broken.

Commit-type: testing-fix,style-fix, feature
Funded-by: ECP
Project: MKL
Time: 32 hours
Reported-by: Irina Sokolova <irina.sokolova@intel.com>

show more ...

Merge remote-tracking branch 'origin/barry/feature-simple-spectral'

Merge branch 'maint'

* maint:
MatMAIJ: generate Fortran stubs
Doc: Fixed formatting problems for MatCreateAIJ() manual page.

Get GLL code to work in complex, prevent spectral examples from running with complex

Commit-type: portability-fix
Time: .2 hours
Reported-by: nightly builds

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

Doc: Fixed formatting problems for MatCreateAIJ() manual page.

Doc: Fixed formatting problems for MatCreateAIJ() manual page.

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

Merge branch 'master' into knepley/feature-dm-bc-components

* master: (320 commits)
fix compile warning: /sandbox/petsc/petsc.clone/src/snes/examples/tests/ex2.c:250:8: warning: ‘pointsAllProcs’ m

Merge branch 'master' into knepley/feature-dm-bc-components

* master: (320 commits)
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
output has changed with e6d0a238963c2a97dd04845ea512b529992c7cdb so update the output file
...

show more ...

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

Merge branch 'hzhang/mpirart'

Add support for MatRARt_MPIAIJ()
Requested-by: Franck Houssen <franck.houssen@inria.fr>

Add support for MatRARt_MPIAIJ()
Requested-by: Franck Houssen <franck.houssen@inria.fr>

Merge branch 'hzhang/fix-blockedIS-submat'

rm unused variable

set blocksizes to submat using input isrow and iscol

fix bug for ISblocksize>1 in ISGetSeqIS_SameColDist_Private() and MatCreateSubMatrix_MPIAIJ_SameRowDist()oReported-by: Hoang Giang Bui <hgbk2008@gmail.com>
When input iscol is split/reduced to smalle

fix bug for ISblocksize>1 in ISGetSeqIS_SameColDist_Private() and MatCreateSubMatrix_MPIAIJ_SameRowDist()oReported-by: Hoang Giang Bui <hgbk2008@gmail.com>
When input iscol is split/reduced to smaller local IS, block structure might be broken, thus created local IS should not inherit blocksize from input iscol

show more ...

Merge branch 'master' into knepley/feature-ex11-simplex

* master: (63 commits)
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Docs+PetscFE: Evaluate() operates in

Merge branch 'master' into knepley/feature-ex11-simplex

* master: (63 commits)
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Docs+PetscFE: Evaluate() operates in reference space
Plex ex7: Added 600 cell test
Plex: Added 600-cell for S^3
Math: Add the golden ratio
Plex ex7: Meshing a sphere
Plex: Changed to DMPlexCreateSphereMesh() - Now can do quads or triangles
Share: Added large 3D mesh file - Its the only one I have gotten to work
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
Mat: Added example for MatIS
Plex ex3: Update to new point function API
initialize some variables as complained by nightly tests
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
update tests
...

show more ...

Merge branch 'scott/test-datafiles-fix' of https://bitbucket.org/petsc/petsc into scott/test-datafiles-fix

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-convergence-estimator

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

Merge branch 'master' into knepley/feature-convergence-estimator

* master: (27 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
update tests
add tests
add MatView_SeqXAIJ_ASCII_structonly(); enable A->structure_only in MatSetValuesBlocked_SeqBAIJ
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
PetscDS: Allow point 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
dmda matstructure_only
add 'structure_only' to seqbaij matrices -- MatView() does not work yet
PetscFE: Initialize temp array before calling f1
...

Conflicts:
include/petsc/private/petscdsimpl.h
include/petscds.h
src/dm/dt/interface/dtds.c

show more ...

Merge branch 'hzhang/mat-structureonly'

MPI_Allreduce() with PetscInt needs MPIU_INT flag, not MPI_INT

OpenMPI generates a crash but MPICH in the PETSc nightly builds does not seem to.

Commit-type: bug-fix
Reported-by: Travis CI

Merge branch 'master' into stefano_zampini/feature-fetidp-ex69