PetscDualSpace: Small cleanup of parameter name

PetscSpace: Added height subspaces for PetscSpace
- Added PetscSapceGetHeightSubspace() and member function
- Added subspaces array

Merge remote-tracking branch 'origin/knepley/feature-convergence-estimator'

* origin/knepley/feature-convergence-estimator:
ConvEst: Restore the SNES to the original condition after estimation
C

Merge remote-tracking branch 'origin/knepley/feature-convergence-estimator'

* origin/knepley/feature-convergence-estimator:
ConvEst: Restore the SNES to the original condition after estimation
ConvEst: Discretization can change, so move name extraction into loop
SNES ex13: Fixed test dependencies
SNES ex13: Update to new API
ConvEst: Initialize variables
Fortran: Fixed Fortran bindings - Need custom binding for all View() to handle default viewers
ConvEst: I hate Fortran
ConvEst: Damn damn damn
ConvEst: Added Fortran type
ConvEst: Fixed spelling, fixed type
ConvEst: Fixing LAPACK type
SNES ex13: Now estimate convergence using built-in test
SNES: Allow convergence estimation using -snes_convergence_estimate
PetscConvEst: Added an object which estimates the convergence rate of a discretization
PetscDT: Turn on Fortran binding
SNES ex13: Added test output
SNES ex13: Looking at automatic convergence estimation
PetscDS: Added exact solution functions - We will use these to evaluate convergence using MMS

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

show more ...

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

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

PetscDS: Fix print

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

* master:
Docs+PetscFE: Evaluate() operates in reference space
Mat: Added example for MatIS
Modified DMNetworkDistribute functi

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

* master:
Docs+PetscFE: Evaluate() operates in reference space
Mat: Added example for MatIS
Modified DMNetworkDistribute function call

show more ...

fix warning with gcc-7.1

/home/balay/petsc/src/dm/dt/interface/dtfe.c: In function ‘PetscFEGetTabulation’:
/home/balay/petsc/src/dm/dt/interface/dtfe.c:3789:5: warning: this ‘if’ clause does not gua

fix warning with gcc-7.1

/home/balay/petsc/src/dm/dt/interface/dtfe.c: In function ‘PetscFEGetTabulation’:
/home/balay/petsc/src/dm/dt/interface/dtfe.c:3789:5: warning: this ‘if’ clause does not guard... [-Wmisleading-indentation]
if (B) *B = NULL; if (D) *D = NULL; if (H) *H = NULL;
^~

show more ...

Merge branch 'knepley/feature-plex-boundary-values'

* knepley/feature-plex-boundary-values:
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix upd

Merge branch 'knepley/feature-plex-boundary-values'

* knepley/feature-plex-boundary-values:
PetscSpace: Handle 0 point case for evaluation
PetscSpacePoint: Fix memory handling
PetscDS: Fix update fund for new API
PetscFE: Derivatives for Point class are all zero
Small fixes
TS ex46: Fix sin and cos
TS ex18: Fix sin and cos
PetscFE: Fixing Point subclass
Plex: Fixed bug with refinement, fixed composed functions for Plex - When we refine, we were calling DMDestroy_Plex(), which wiped out the composed funds, but not calling DMInitialize_Plex() to reinsert them in to the refined mesh
PetscSpacePoint: Fix handling of points
PetscDS: Added update() function for fields - This is intended to directly update the values at points
Plex: DMPlexInsertBoundaryValues() can now be overridden by the user

Conflicts:
include/petsc/private/petscdsimpl.h
include/petscds.h

show more ...

PetscFE: Added Fortran bindings

Docs+PetscFE: Evaluate() operates in reference space

Fortran: Fixed Fortran bindings
- Need custom binding for all View() to handle default viewers

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 'knepley/feature-fe-kernel-scalars'

* knepley/feature-fe-kernel-scalars:
SNES ex56: Update test output
Plex ex3: Update to new point function API
PetscDS: Turn off Fortran binding

Merge branch 'knepley/feature-fe-kernel-scalars'

* knepley/feature-fe-kernel-scalars:
SNES ex56: Update test output
Plex ex3: Update to new point function API
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

Conflicts:
src/dm/dt/interface/dtfe.c

show more ...

PetscSpace: Handle 0 point case for evaluation

PetscSpacePoint: Fix memory handling

PetscDT: Turn on Fortran binding

PetscDS: Added exact solution functions
- We will use these to evaluate convergence using MMS

PetscDS: Fix update fund for new API

PetscDS: Turn off Fortran bindings for Constants

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

PetscFE: Initialize temp array before calling f1

PetscFE: Derivatives for Point class are all zero

Small fixes

PetscFE: Fixing Point subclass