DM: DMCreateSubDM() field argument should be const

Section: Added PetscSectionCreateSupersection()
- This merges several sections together

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

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

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 'knepley/fix-plex-example-logging' into knepley/feature-plex-torus

* knepley/fix-plex-example-logging: (151 commits)
Plex ex12: Added logging stages
Plex ex1: Added logging stages

Merge branch 'knepley/fix-plex-example-logging' into knepley/feature-plex-torus

* knepley/fix-plex-example-logging: (151 commits)
Plex ex12: Added logging stages
Plex ex1: Added logging stages
update to latest hypre tag on github
fix 'inf' vs 'Inf' diff wrt vec_vec_utils_tagger_tutorials-ex1_8,9,18,19
MatSetOption(mat,MAT_IGNORE_ZERO_ENTRIES,PETSC_TRUE) no longer applies to zero values set or added to the matrix diagonal, they are now not ignored and are included in the matrix nonzero structure. This is useful for Jacobians passed to TS.
replace use/mention of depreciated -log_summary in many places
add extchem.py which was written years ago but never put in the distribution (hence it is broken).
petsc4py: update to snapshot with MatGetSubMatrices -> MatCreateSubMatrices etc API changes
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes
add help menu for MatCreateSubMatricesMPI()
PetscCMPLX: assume type punning is safe
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
...

Conflicts:
src/dm/impls/plex/examples/tests/ex1.c
src/dm/impls/plex/examples/tests/ex12.c

show more ...

Merge branch 'knepley/fix-plex-example-logging' into knepley/feature-plex-tests-new

* knepley/fix-plex-example-logging: (151 commits)
Plex ex12: Added logging stages
Plex ex1: Added logging stag

Merge branch 'knepley/fix-plex-example-logging' into knepley/feature-plex-tests-new

* knepley/fix-plex-example-logging: (151 commits)
Plex ex12: Added logging stages
Plex ex1: Added logging stages
update to latest hypre tag on github
fix 'inf' vs 'Inf' diff wrt vec_vec_utils_tagger_tutorials-ex1_8,9,18,19
MatSetOption(mat,MAT_IGNORE_ZERO_ENTRIES,PETSC_TRUE) no longer applies to zero values set or added to the matrix diagonal, they are now not ignored and are included in the matrix nonzero structure. This is useful for Jacobians passed to TS.
replace use/mention of depreciated -log_summary in many places
add extchem.py which was written years ago but never put in the distribution (hence it is broken).
petsc4py: update to snapshot with MatGetSubMatrices -> MatCreateSubMatrices etc API changes
Minor fixes for tchem examples due to changes in PETSc infrastructure
update webpage about API changes
add help menu for MatCreateSubMatricesMPI()
PetscCMPLX: assume type punning is safe
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
...

Conflicts:
src/dm/impls/plex/examples/tests/ex12.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 'tisaac/feature-vec-tagger'

* tisaac/feature-vec-tagger: (31 commits)
VecTagger: documentation
TS tutorials ex11: use updated VecTagger interface
VecTagger: rename Intervals and B

Merge branch 'tisaac/feature-vec-tagger'

* tisaac/feature-vec-tagger: (31 commits)
VecTagger: documentation
TS tutorials ex11: use updated VecTagger interface
VecTagger: rename Intervals and Boxes, declare VecTaggerBox
PetscCMPLX: fix inline declaration
VecTagger: silence some compiler complaints
VecTagger: 64-bit indices fixes
TS tutorials ex11: use VecTagger
VecTagger: cxx fixes
VecTagger: complex fixes and tests
VecTagger: fix view of andor types
VecTagger: fixes for complex
VecTagger: Implementation of union and intersection taggers
VecTagger: cdf interval implementation
VecTagger: create regression tests based on ex1
VecTagger: complex/c++ fixes
VecTagger: added tutorial ex1
VecTagger: simple backend should check that interval is set at setup
VecTagger: relax strict block size restrictions
VecTagger: comment out unimplemented registrations
VecTagger: fix header function misnames
...

show more ...

Merge remote-tracking branch 'origin/pr628/mapdes/scalable-isgtol/master'

IS: Add ISIntersect()

Add initial proper support for multiple implementations of ISGlobalToLocalMappingApply()

Provides two subtypes basic and hash. Selects hash automatically if the memory needed for basic
(end - start)

Add initial proper support for multiple implementations of ISGlobalToLocalMappingApply()

Provides two subtypes basic and hash. Selects hash automatically if the memory needed for basic
(end - start) is greater than 4 times the size of the local to global mapping (n) and 1,000,000
Utilizes the hash support provided by Lawrence Mitchell

Thanks-to: Lawrence Mitchell

show more ...

ISLtoG: Allow controlling GtoL behaviour from options

Use -ltogm_use_hash_table to use a memory scalable approach for
applying the global to local mapping.

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/allow-late-matsetblocksizes'

Merge branch 'master' into tisaac/dmplex-feature-real-to-ref

I need updated PetscFE interfaces in what I'm going to do

* master: (603 commits)
parmetis: fix typo in testsuite
Added PETSC_HAVE_A

Merge branch 'master' into tisaac/dmplex-feature-real-to-ref

I need updated PetscFE interfaces in what I'm going to do

* master: (603 commits)
parmetis: fix typo in testsuite
Added PETSC_HAVE_ATTRIBUTEALIGNED flag to makefile
PCSetUp_SVD() used incorrect matrix for matrix size in parallel, also need to skip LAPACK call if matrix size is zero
strumpack: update to latest version 1.0.4
Remove valgrind warning about MPI_Allreduce() with uninitialized data
sed for test output that has atol convergence with quad precision
updated ex1_nest tests
fixed incorrect key value
fix: bad merge a69d34a
test: move runex78 to TESTEXAMPLES_C_NOTSINGLE
build: fix warnings
sys: update test to produce same result - irrespective of debug/opt build - and presence of ~/.petscrc
sys: src/sys/examples/tests/ex30.c -> src/sys/examples/tests/ex32.c to avoid conflict in master
Sed output, since quad precision reduces below atol on these tests
remove automatically generated matrix names from test output files since can be different
Fixes for MatMPIXXXSetPreallocation(), they were not reentrant, since they kept the information used for matrix-vector products
GAMG: fix PCView_GAMG() format - and update output/ex56_0.out
superlu_dist: 64-bit-indices fix https://github.com/xiaoyeli/superlu_dist/issues/4
Add SNESSetDivergenceTolerance / -snes_divergence_tolerance.
fixed unit read in Plex
...

show more ...

Added ISLocalToGlobalMappingDuplicate

if refcnt > 0 we cannot change the block size of the l2gmap in-place

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

First attempt to allow changing the block size of a matrix (and any associated l2gmap)

Merge branch 'barry/add-matfdcoloring-ghosted-example'

Changed IS_COLORING_GHOSTED to IS_COLORING_LOCAL, added test case, fixed code to work properly

Time: 5 hours
Reported-by: Hong Zhang <hzhang@mcs.anl.gov>

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

fixed merge with master

Merge branch 'tisaac/dualspace-feature-symmetry'

In finite element assembly, to compute a residual or assemble a Jacobian
on an element, the dual space functionals of a reference finite element are

Merge branch 'tisaac/dualspace-feature-symmetry'

In finite element assembly, to compute a residual or assemble a Jacobian
on an element, the dual space functionals of a reference finite element are related
(i.e., pushed forward) to the functionals that make of the degrees of
freedom of the discretized system.

In doing so, we have to draw a correspondence between the boundary
points of the reference element (its edges and vertices) and the
boundary points of each cell in the mesh. Sometimes, this
correspondence is not perfect, because the orientations of the boundary
points of a cell may be different than the orientations of the boundary
points of the reference element. For example, an edge that runs
clockwise around the reference element may map onto an edge that runs
counter-clockwise around the real cell.

For some finite elements (such as the most common P1 Lagrange finite
elements) this doesn't matter. But for others---high-order elements,
H-div and H-curl conforming elements---a difference in orientation
affects the mapping between reference functionals and degrees of
freedom.

The simplest example is a high-order nodal Lagrange finite element in
2D: the degrees of freedom on an edge of the reference element are
numbered left-to-right, but if the orientation of the real edge is
reverse, then we must read/write degrees of freedom right-to-left.

Things get more complicated in 3D: a quadrilateral face between
hexahedra may have any of 8 different orientations (that make of the
dihedral symmetries), and each one corresponds to a different reordering
of the 2D grid of nodes supported on that quadrilateral.

Things also get more complicated for H-div elements. Functionals on the
interfaces between cells represent the flux from one cell to another,
and thus have direction. Reversing the orientation of an edge changes
the sign of the mapping between the reference functional and the degree
of freedom.

To assemble generic finite elements, we thus need to accommodate
transforms between reference functionals and degrees of freedom that
encompass (1) permutations, and (2) scalar multipliers.

This branch accomplishes this in two steps:

- We add PetscDualSpaceGetSymmetries() to the interface, to allow a dual
space to describe how symmetries of the mesh points affect
referent-to-real maps of functionals.

- We add PetscSectionSym, an object that encapsulates symmetries of the
degrees of freedom whose layout is described by a PetscSection.

When a section is created from a DM that has a finite element
discretization in its PetscDS, it will automatically construct the
appropriate PetscSectionSym when it creates a PetscSection to describe a
discretized function space.

If a PetsFE finite element discretization is not used, the user can
create a PetscSectionSym for their degrees of freedom for themselves.
We give an example of this for spectal elements in
src/dm/impls/plex/examples/tutoerials/ex6.c.

* tisaac/dualspace-feature-symmetry: (22 commits)
cast memzeros to void to avoid compiler complaints about const qualifiers
DMPlex: fix stray PetscScalar to PetscReal
DMPlex: fix stray PetscScalar to PetscReal
Plex tutorials ex6: use symmetries in SEM example
Plex tests ex3: added high order tests
Plex tests ex3: add matrix-free near null space test
DMPlex: fix PetscInt %D PetscPrintf character
DMPlexTree: use symmetries when computing constraint matrices
DMPlex: use symmetries in DMPlexGetIndicesPoint{Fields}_Internal() and DMPlexAnchorsModifyMat()
DMPlex: use PetscSectionGet/RestorePointSyms() in DMPlexVecGet/SetClosure()
DMPlex: added DMPlexGet/RestoreCompressedClosure()
DMPlex: added DMPlexGetPointDualSpaceFEM()
DMPlex: set PetscSectionSym in DMPlexCreateSectionInitial()
DMLabel: added PETSCSECTIONSYMLABEL
PetscSection: add PetscSectionSym
PetscDualSpace: added PetscDualSpaceGetSymmetries()
PetscDualSpace_Lagrange: setup recursively on dimension
PetscDualSpace_Lagrange: order functions now go after tuple used
PetscDualSpace_Lagrange: use lexicographic order for nodes
DMPlex: added DMPlexComputePointGeometry_Internal()
...

show more ...