aboutsummaryrefslogtreecommitdiff
path: root/src/GenericRK.c
Commit message (Collapse)AuthorAge
* MoL Updateeschnett2013-01-22
| | | | | | | | | | | | | | | | | | New integrator Euler. This is an explicit, first-order method, mostly useful for debugging. New integrators AB (Adams-Bashforth) with various orders. These are explicit integrators using several past timelevels to provide higher-order integration with a single RHS evaluation each. Introduce LinearCombination, a generic routine to calculate linear combinations of grid functions. This simplifies existing code, and can be overloaded if MoL should run on a device (e.g. with OpenCL). Replace cctk_lsh with cctk_ash where necessary. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@190 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Remove unused variable declarationshinder2012-11-26
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@187 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Parallelise loops with OpenMP.schnetter2008-09-23
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@129 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Modification for evolving array variables; change to localhawke2006-08-01
| | | | | | | | | arrays. This gets around the problems with the driver allocating the scratch space differently on multiprocessors from the individual arrays. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@116 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Fix for evolving arrays on multi-processors with the generic RKhawke2006-07-28
| | | | | | | | method. This fix should use minimal memory, and works even when there is no storage for the array on the local processor. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@115 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Remember the variable index of the scratch variables.schnetter2005-12-11
| | | | | | | | | | This saves calling CCTK_FirstVarIndex during each MoL iteration. This function is expensive, as it performs case insensitive string comparisons, and it showed up high on a profile of some of the Mexico tests, which integrate small domains over long times. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@105 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add const and restrict qualifiers to the pointers.schnetter2005-12-11
| | | | | | | | | | Simplify some complex arithmetic. Initialise the error variable in the RK45 integrator only after the last iteration. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@104 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Replace "MoL" by CCTK_THORNSTRING in calls to CCTK_VWarn.schnetter2005-12-11
| | | | | | | | | Correct some typos in warnings. Remove some commented out variable declarations. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@102 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add the Classic RK3 method (as a generic method, so use Generic_Type).hawke2004-07-01
| | | | | | | | | | | | Agrees with other RK3's to floating point round off (except at boundaries) for linear case. Uses more storage and is slower than standard RK3 so I don't recommend it. This showed up (so I fixed) a bug with the generic methods when used with Carpet. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@72 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Update CCTK_FILEVERSION macros to CactusBase.hawke2004-07-01
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@71 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Make MoL independent of dimensions for GFs.hawke2004-06-23
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@70 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct a bug in storing the scratch space levels. This will fix RK4.hawke2004-02-28
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@49 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Formatting changes only.hawke2003-07-24
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@32 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add two methods:hawke2003-07-22
| | | | | | | | | | | RK3. The optimized version of the TVD RK3 solver. Requires no scratch space so is about as efficient as ICN, but third order. Generic method from a parameter table. By specifying the number of intermediate steps and the alpha and beta arrays, create your own method at parameter time. Not well (or at all) documented because it doesn't seem to work correctly at the moment. Some tidying of extraneous code as well. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@29 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Take the ScratchSpace pointers from VarIndex instead of making assumptions ↵hawke2003-07-17
| | | | | | about memory layout. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@22 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correctly (?!) evolve the array variables when using the generic RK methods.hawke2003-06-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@17 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Comment or ifdef out complex variable support as Carpet doesn't like it.hawke2003-05-24
| | | | | | | To use the complex variable stuff right now you need to uncomment the appropriate groups in the interface.ccl and #def MOLDOESCOMPLEX in the appropriate files. It still probably won't work. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@13 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add support for evolving complex GFs and (real and complex) GAs.hawke2003-05-21
| | | | | | | | | Only works with ICN or RK2 for now - in fact this commit may break the generic RK methods temporarily. Note the documentation isn't quite right - there's no longer a seperate function for each different type... git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@12 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Change the commented out debugging statements to ifdef'd out.hawke2003-04-25
| | | | | | | Remove old, buggy commented out stuff. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@6 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Change the fileversion strings and comments from MoL2->MoL (cosmetic).hawke2003-04-24
| | | | | | | Change the grid scalar used in the schedule.ccl from MoL2->MoL (necessary). git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@4 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* The Method of Lines thorn (version 2 - see below).hawke2003-04-23
MoL provides generic integration methods for multiple thorns simultaneously. By providing a layer between the driver and evolution thorns, this should mean that some technical issues to do with mesh refinement can be ignored. It also allows you to choose different evolution methods (in time, at least). But the primary purpose is to unambiguously evolve models in different thorns at the same time. This is version 2 - the one that will work with mesh refinement. It's a straight copy of HawkeCVS/Public/CactusMoL2/MoL2 and I haven't checked that it will work "as is" with the new name (one of the reasons it goes into Alpha). At the moment the only evolution method guaranteed to work with mesh refinement is ICN, although RK2 should. The generic RK methods will do something that will be subtly wrong... Note that the "old" way of registering variables (through the functions declared in header files) is still there. As soon as function aliasing settles down, this will be removed. Also to be done: Better documentation Tidy up the code (especially the debugging statements) Optimize Add various useful time evolution methods. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@2 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b