Merge remote-tracking branch 'origin/hongzh/adjoint-parallel'

Merge remote-tracking branch 'origin/barry/update-petscfileretrieve'

Merge remote-tracking branch 'origin/barry/fix-use-objectid-not-object-for-comparison-in-ts'

Merge remote-tracking branch 'origin/barry/fix/even-huger-flaw-in-ts'

Merge remote-tracking branch 'origin/barry/fix-huge-flaw-in-ts'

Merge remote-tracking branch 'origin/barry/feature-ts-monitor-solution'

Merge remote-tracking branch 'origin/barry/consistent-formating-ascii-solver-viewers'

Make ASCII PC/KSP/SNES/TSView() code and output in a standard style

The code that prints the ASCII view for solvers was occasionally inconsistent:
1) TSView printed the type specific information aft

Make ASCII PC/KSP/SNES/TSView() code and output in a standard style

The code that prints the ASCII view for solvers was occasionally inconsistent:
1) TSView printed the type specific information after the general info while
all other solvers printed it before (right after the name of the type is printed)
2) KSPView consistently printed the name of the subtype on EACH line of the subtypes
output as did a few of the PC and SNES viewers. Since they are all printed indented
directly below the subtype name there is no reason to print this information on each line
3) TSView printed output about number of linear solvers and SNES even for explicit
methods, which is goofy.
4) a few other minor formatting and consistency issues where fix.

Commit-type: style-fix, feature
Time: 16 hours

show more ...

fix problems caused by interface change and TSSetExactFinalTime

stability_9bus examples

Merge branch 'knepley/test-process-arg' into knepley/feature-adapt-interface

* knepley/test-process-arg: (49 commits)
Test: Can now use make test NP=<num procs>
Testing: Fix error reporting
I

Merge branch 'knepley/test-process-arg' into knepley/feature-adapt-interface

* knepley/test-process-arg: (49 commits)
Test: Can now use make test NP=<num procs>
Testing: Fix error reporting
I didn't realize there was a hacky special case for dense B matrix
Fix example to support -snes_mf_operator
cleanup src/ksp/ksp/examples/tutorials/network/makefile
add missing input in the help menu of PetscOptionsView() Reported-by: Richard Katz <richard.katz@earth.ox.ac.uk>
Fix docs for MatMatMult() that C matrix can only be created by user for sequential case.
protect '-' in macros
fix compile error with --download-petsc4py
Plex: Removed useless CHKERRQ
Configure: Replace use of PETSc modules
Plex ex7: Damn, quad precision is named __float128
TSSetFromOptions: Give subtypes a chance to set the default adapt type
TS ex11.c: fix corrupt filename initialization
small fix
Plex ex1: Update output
DT: Fixed Fortran module - Add DT module to DM and Petsc modules - Fix makefile to generate DT F90 module - Fix PetscQuadratureGet/Restore/SetData() F90 defs
Plex ex19: Updated output
Plex ex7: Update output
p4est: Update to new API
...

show more ...

Merged in dalcinl/fix-tssetfromopts-adapt (pull request #691)

TSSetFromOptions: Give subtypes a chance to set the default adapt type

Approved-by: Jed Brown <jed@59a2.org>
Approved-by: BarryFSmith <

Merged in dalcinl/fix-tssetfromopts-adapt (pull request #691)

TSSetFromOptions: Give subtypes a chance to set the default adapt type

Approved-by: Jed Brown <jed@59a2.org>
Approved-by: BarryFSmith <bsmith@mcs.anl.gov>

show more ...

TSSetFromOptions: Give subtypes a chance to set the default adapt type

Merge branch 'master' into knepley/feature-plex-wedge

* master: (61 commits)
Added ignore
Plex ex7: Fixed mesh name
Plex ex1/ex2: Fixed requirements
PCGAMG: fix spelling in option docs
Pet

Merge branch 'master' into knepley/feature-plex-wedge

* master: (61 commits)
Added ignore
Plex ex7: Fixed mesh name
Plex ex1/ex2: Fixed requirements
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
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: Handle unused variable
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
Mat: Added example for MatIS
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
...

Conflicts:
src/dm/impls/plex/plexcreate.c

show more ...

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

* master: (56 commits)
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
TSMonitorDrawSo

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

* master: (56 commits)
PCGAMG: fix spelling in option docs
PetscDS: Fix print
Plex ex19: Output fixes for proper destruction
TSMonitorDrawSolutionPhase: Add missing draw pause
fix warning with gcc-7.1
Plex ex7: Added 600 cell test
Plex: Added 600-cell for S^3
Math: Add the golden ratio
Plex: Handle unused variable
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
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
Plex: Fixing insertion of boundary values - If we have a dmAux that is a boundary, we need minHeight = 1, since faces can be mapped from submesh to mesh, but not cells (since they have incomplete cones) - However, in this case the effectiveHeight should be h, whereas when using a boundary as the mesh as well, the effectiveHeight should be h - minHeight - Obviously, we need more thinking here to make the iteration make sense in general
PetscFE: Derivatives for Point class are all zero
Small fixes
TS ex46: Fix sin and cos
TS ex18: Fix sin and cos
...

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

TSMonitorDrawSolutionPhase: Add missing draw pause

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

Force users to provide Jacp matrix

Users need to make sure that the JacobianP matrix can be multiplied by the mu vectors from right.

Change TSSetCostIntegrand interface

Originally PETSc creates a sequential vec to hold the integral values for users. This is ok for common senarios where there is just one integral to calculate. Whe

Change TSSetCostIntegrand interface

Originally PETSc creates a sequential vec to hold the integral values for users. This is ok for common senarios where there is just one integral to calculate. When there are multiple integrals (e.g. coming from multiple cost functions), one might want to do it in parallel. The new interface allows users to provide their own vec to support parallelization. If the users do not provide one, a sequential vec will be created.

show more ...

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

Fix TSAdjointMonitorGradient -> TSAdjointMonitorSensi

Commit-type: bug-fix

Change TSAdjointMonitorGradient() to TSAdjointMonitorSensi() to match other usage in TS

Commit-type: style-fix
Funded-by: IMEX
Project: adjoints
Reported-by: Hong Zhang <hongzang@anl.gov>

Do not use pointer values to compare two objects

One of the objects may have been destroyed and then recreated at the same memory
location. PetscObjectIds are unique so compare with them instead.

C

Do not use pointer values to compare two objects

One of the objects may have been destroyed and then recreated at the same memory
location. PetscObjectIds are unique so compare with them instead.

Commit-type: bug-fix
Funded-by: ECP
Time: .2 hours

show more ...

Add simple TSAdjointMonitorGradient() and improve manual pages for possible RK schemes

Commit-type: feature,documentation
Funded-by: IMEX
Project: ajoints
Time: .5 hours

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