TS: add TSSetTransientVariable (alias for DMTSSetTransientVariable)Requested-by: Hong Zhang <hongzhang@anl.gov>
TS: add support for transient variables in BDFThis feature enables conservative time integration even when usingnonconservative state variables. The formulation, which requires onlythe definitio
TS: add support for transient variables in BDFThis feature enables conservative time integration even when usingnonconservative state variables. The formulation, which requires onlythe definition of transient (conservative) variables C(U) as afunction of the state variable U, is appropriate for linear multistepmethods for ODE and DAE.This technique could be extended to Runge-Kutta methods by creating anindex-1 DAE by augmentation with the equation C = C(U). The stageequations in diagonally implicit methods generally ask to solve for Y: F(Y, alpha Y + Z) = 0where Y = [C, U] and Z is known from prior stages. With the augmentedequation, this decouples to C - C(U) = 0 f(U, alpha C + C_z) = 0This commit includes an example demonstrating loss of conservation inchain rule primitive-variable IFunction, while preserving it forpure-conservative variables and the transient variable formulation.Resolves #547Requested-by: Gautam Bisht <gautam.bisht@pnnl.gov>Thanks-to: Hong Zhang <hongzhang@anl.gov>
show more ...
TS+ConvEst: Made suggested improvements- Renamed initialGuess to initialConditions- Made refinement factor usable settable, but still only support 2 for space
ConvEst: Generalize to use TS as well- Now the solver is a PetscObject- ts_convergence_estimate calculates the order for a TS method- Added tsconvest.c- Added initial guess and error methods that
ConvEst: Generalize to use TS as well- Now the solver is a PetscObject- ts_convergence_estimate calculates the order for a TS method- Added tsconvest.c- Added initial guess and error methods that can be overridden
Merge branch 'jczhang/fix-dm-originaldm/maint' [PR #1785]* jczhang/fix-dm-originaldm/maint:Set DMKSP originaldm when doing copy or getting write access.
Initialize DMTS originaldm when creating DMTS
SNES+Plex: Refactored DMSNESCheck() into 3 checks- Check for discretization error in solution- Check residual of solution- Check Taylor test for Jacobian (removed more simplistic Jacobian check)-
SNES+Plex: Refactored DMSNESCheck() into 3 checks- Check for discretization error in solution- Check residual of solution- Check Taylor test for Jacobian (removed more simplistic Jacobian check)- Changed DMSNESCheckFromOptions_Internal() to DMSNESCheck_Internal()
DM: Reorganized DM fields- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)- This will feed into the creation of PetscDS for every different cell system-
DM: Reorganized DM fields- A field is defined by a discretization (PetscFE or PetscFV) and a region (DMLabel or NULL)- This will feed into the creation of PetscDS for every different cell system- Changed interface for Get/SetField to add the region
Merge branch 'knepley/feature-plex-elasticity-test' [PR #1071]* knepley/feature-plex-elasticity-test:Makes a test for regular elasticity.Matt setups everything so that we can do GMG (ex56 does AM
Merge branch 'knepley/feature-plex-elasticity-test' [PR #1071]* knepley/feature-plex-elasticity-test:Makes a test for regular elasticity.Matt setups everything so that we can do GMG (ex56 does AMG).There are also crucial improvements in the check routines.
Doc: More function documentation
Plex+TS: We cannot run the check with the same vector that is the TS solution
Merge branch 'master' into knepley/feature-plex-elasticity-test* master: (1553 commits) -options_table was changed to options_view in 2014 but documentation was not updated to reflect it config
Merge branch 'master' into knepley/feature-plex-elasticity-test* master: (1553 commits) -options_table was changed to options_view in 2014 but documentation was not updated to reflect it configure: fix F2003 test Plex test ex5: Removed "future" test, as it resulted in issues in the nightly tests. Some minor typos in TS Add ADIOS2.py which I forgot to git add in the past Plex ex5: Disable test Doc: Added note for interface change KSPLSQRMonitorDefault fixed for quad precision Tao: Added alternative output for chwirut2 test. TS ex45,ex46: Removed overly refined tests for faster nightly tests. KSP ex72: Fixed needlessly long-running PCTFS tests. SNES ex56: Reduced test size to speed up nightly test suite. KSPSolve_LSQR: don't compute unused norm of rhs SNES ex13: Reduced problem sizes for nightly tests. PetscOptionsDeprecated(): fix for NULL prefix PetscOptionsDeprecated(): add option to print additional info SNES ex12: Reduced sizes for BDDC tests. PetscFileRetrieve(): Suppress output to STDERR when downloading a file. MatMatMult: Fixed unused variable pN. DMDestroy: Guard call to dm->ops->destroy * Similar guards are in place in many other methods (IS,Vec,Mat) and permit the DM to be destroyed even if DMSetType() was never called. ...
Merge branch 'knepley/fix-plex-residual-loop'* knepley/fix-plex-residual-loop: Plex+SNES: Protect Destroy Plex+SNES/TS: Moved zeroing of Jacobian out of DMPlexComputeJacobian_Internal() - Shoul
Merge branch 'knepley/fix-plex-residual-loop'* knepley/fix-plex-residual-loop: Plex+SNES: Protect Destroy Plex+SNES/TS: Moved zeroing of Jacobian out of DMPlexComputeJacobian_Internal() - Should probably go in DMLocal versions, but now its in the top-level Plex FEM versions SNES+Plex: Needed cell offset, not cell number Sys: Added a check for gnome-terminal, which never uses -display - Added some docs Plex: Made a public interface for ComputeBdResidualSingle() Plex: The first argument to a pointfunc should be the coordinate dimension
ghost points for X_t is added to TSComputeIFunction_DMDA
Merge branch 'master' into knepley/fix-plex-residual-loopWarning: This merge includes real changes* master: (188 commits) Add JuMP citation Copyedits and updates to developers manual Add mor
Merge branch 'master' into knepley/fix-plex-residual-loopWarning: This merge includes real changes* master: (188 commits) Add JuMP citation Copyedits and updates to developers manual Add more protection to intrinsics MPI: protect arguments in MPI macros Add the missing conditionals for 64-bit build Plex+VTK: Fix field type determination Fix bug in name cleanup Bib: Added ref Plex+VTK: Fix field type determination Plex: Updated test output - I am now hating this shit forgot to save file after moving prototypes to petscsys.h update Apple example project settings per Xcode Two functions should actually be public, used, for example from Apple MacOS and iOS programs includes should always be included relative to include directory Apple removed old locations of docsets SNES ex12: Updated test output SNES ex62: Set name of solution vector - Update test output revert sources to dev - reset PETSC_VERSION_RELEASE - reset changes/dev.html release: set v3.9 strings move output file to the correct location. i.e fix 305862f9cd0dc399542ae11dd744508e4b6ad2e3 ...Conflicts: src/snes/utils/dmplexsnes.c src/ts/utils/dmplexts.c
Plex+SNES/TS: Moved zeroing of Jacobian out of DMPlexComputeJacobian_Internal()- Should probably go in DMLocal versions, but now its in the top-level Plex FEM versions
DM+SNES: Changed signature for internal check function
DMPlex: rename newly created label
DMPlexSNES: caching of new geometry constructsThe old paradigm cached one geometry object for all cells or all faces,an approach that doesn't take into account that different quadraturesmay requi
DMPlexSNES: caching of new geometry constructsThe old paradigm cached one geometry object for all cells or all faces,an approach that doesn't take into account that different quadraturesmay require different geometry objects. Now the geometry objectis composed with the index set over which the geometry ranges, and the compositionstring include the object id of the quadrature. To take advantage ofthis caching, we have to be careful not to reconstruct index sets on thefly when they are known a priori. In particular, the hybrid boundshave been encapsulated as a label, which makes them accessible as ISs.
Clean up simple spectral adjoint example, add advection termAdd crosslinks for some TS example and additional monitor functionsCommit-type: testing-fix,style-fix, feature
Remove the use and definition of __FUNCT__ throughout the codeSince all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETScTime: 1.5 hoursThanks-t
Remove the use and definition of __FUNCT__ throughout the codeSince all modern C/C++ compilers provide this functionality we no longer need to provide it manually in PETScTime: 1.5 hoursThanks-to: Andreas Mang <andreas@ices.utexas.edu>
DMPlexTSComputeRHSFunctionFVM: ADD_VALUES from locF to F instead of INSERT_VALUESWe've actually designed it so that each cells contribution is only calculated once, but there are empty ghost cells
DMPlexTSComputeRHSFunctionFVM: ADD_VALUES from locF to F instead of INSERT_VALUESWe've actually designed it so that each cells contribution is only calculated once, but there are empty ghost cells on other processes with zero contributions. There's no version of INSERT_VALUES that just inserts the local values: it's essentially non-deterministic which value gets inserted. So we have to add all the zero contributions.
DMPlexTS: removed unused variables
DMPlex/DMP4est: small fixes for next warnings and errors
DMPlexTS, DMPlexSNES: cache in the dm instead of the dmts, dmsnes objects
1...<<11121314151617181920>>...22