MRK: allow more than 2 paritions in the nonsplit version
MRK: add multirate suport for more than 2 partitionsThe key idea is to allow for slow-fast paritioning recursively. For example, the first level could consist of a slow part and a fast part. Then t
MRK: add multirate suport for more than 2 partitionsThe key idea is to allow for slow-fast paritioning recursively. For example, the first level could consist of a slow part and a fast part. Then the fast part can be futher split into a slow part and a fast part. Each subts has its own rk context, but share as many working vectors as possible.
show more ...
Move multirate implementation out of rk.c
Simplify the nonsplit version of MRK- VecISCopy is more efficient- EvaluateStep_Rk can be reused
Overhaul the MRK implementation using RHSSplit
add a dense output for tsrk5dp
add multirate method to rkAlso- add partitioned RHS to mrk- add command line option -ts_rk_multirate_type
Fixes for begining of users manual and a variety of webpages including manual pagesCommit-type: bug-fix, documentation
TSGetStages: make each return argument optional
Add if (*ierr) return checks for all Fortran stubs that utilize FREECHAR()Otherwise the ierr is overwritten by the FREECHAR() thus returning a zero error code when it should notCommit-type: bug-f
Add if (*ierr) return checks for all Fortran stubs that utilize FREECHAR()Otherwise the ierr is overwritten by the FREECHAR() thus returning a zero error code when it should notCommit-type: bug-fixReported-by: Tim Steinhoff <kandanovian@gmail.com>
Remove memory leaks of XXXRegister() is called but XXX object is never usedThis means adding YYYInitialPackage() at the beginning of XXXRegister() routinesMinor fixes of missing prototypesCommi
Remove memory leaks of XXXRegister() is called but XXX object is never usedThis means adding YYYInitialPackage() at the beginning of XXXRegister() routinesMinor fixes of missing prototypesCommit-type: style-fixReported-by: Jed Brown <jed@jedbrown.org>Thanks-to: Fande Kong <fande.kong@inl.gov>
Merge branch 'maint'
Fix the wrong time used by TSXXXCostIntegral_RK()
Fix integral evaluation for RKRK resets the integral value to the previous value when timestep adaptivity is turned on and a time step is rejected. This commit fixed the bug.Reported-by: Mirko Ha
Fix integral evaluation for RKRK resets the integral value to the previous value when timestep adaptivity is turned on and a time step is rejected. This commit fixed the bug.Reported-by: Mirko Hahn <hahnm@anl.gov>
Fix for PetscFormatRealArray() due to updates to PetscVSNPrintf()Also removed duplicate copies of PetscFormatRealArray()Commit-type: bug-fixReported-by: nightly tests
TSRK: some optimizations to TSAdjointStep_RKreorder operations so that just 1 temporary vector (VecSensiTemp) is neededTemporary vector obtained from DM (remove storage in RK private data)Use Vec
TSRK: some optimizations to TSAdjointStep_RKreorder operations so that just 1 temporary vector (VecSensiTemp) is neededTemporary vector obtained from DM (remove storage in RK private data)Use VecMAXPY with VecDeltaLamSkip first stage if last b is zero (i.e. RK3BS, RK5BS and RK5DP)No need to assemble JpFix ex9adj of power_grid example
Merge remote-tracking branch 'origin/jed/variadic-malloc'
Fixes for use of PetscMallocN() with const type inputMicrosoft compilers warn about incompatiblity between pointers passed to new PetscMallocA() functionCommit-type: portability-fix, testing-fix
Fixes for use of PetscMallocN() with const type inputMicrosoft compilers warn about incompatiblity between pointers passed to new PetscMallocA() functionCommit-type: portability-fix, testing-fixReported-by: nightly tests
DMSubDomainHookRemove: add interface and use to clean up TS hooks
Refactor forward SARenamed a TS public function (TSAdjointComputeCostIntegrand->TSComputeCostIntegrand); implementations have been simplifed.
Merge remote-tracking branch 'origin/barry/feature-ts-monitor-solution'
Make ASCII PC/KSP/SNES/TSView() code and output in a standard styleThe 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 styleThe 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 line3) 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, featureTime: 16 hours
TSRK: Fix tableaux to use PetscRealConstant() for FP literalsEmacs tip: M-x query-replace-regexp \(-?[0-9]*\.[0-9]+\) -> RC(\,\1)
Fix types introduced in docs for TSRK schemes and improve docs for TSARKIMEX schemesCommit-type: documentationFunded-by: IMEXTime: .1 hour
12345678910>>...16