aboutsummaryrefslogtreecommitdiff
path: root/src
Commit message (Collapse)AuthorAge
* Peter Diener's RK65 and RK87 adaptive timestep integrators.hawke2006-01-23
| | | | | | | | | As yet not altered to do grid arrays. As with RK45, adaptive timestepping does not work with mesh refinement. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@106 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
* Handle CCTK_SCALAR in the same way as CCTK_ARRAY. Previously,schnetter2005-12-11
| | | | | | | | | | CCTK_SCALAR variables were rejected. This is most likely just an oversight, since no one so far wanted to integrate CCTK_SCALAR variables. The change is trivial. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@103 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
* Correct critical error in the RK45 time integrator.schnetter2005-12-11
| | | | | | | | The current time and the time step sizes were not set correctly. That led to first order convergence instead of fifth order convergence. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@101 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Replace "MoL" in CCTK_VWarn by CCTK_THORNSTRING.schnetter2005-12-11
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@100 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Replace "MoL" in CCTK_VWarn by CCTK_THORNSTRING.schnetter2005-12-11
| | | | | | | Correct typo in warning message. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@99 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct error in calculating the number of grid points for the Countschnetter2005-11-22
| | | | | | | variable. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@97 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Small optimisation: substitute division by multiplication with the inverse.tradke2005-08-29
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@91 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Check that a constrained GF may have more than one timelevel beforehawke2005-05-04
| | | | | | | | | | complaining that it doesn't when copying data (if a constrained GF only has one timelevel then such a copy is unnecessary). Fixes PR/1921. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@89 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct a bug in the registration routines; if a variable washawke2005-03-30
| | | | | | | | | registered as SandR and then as evolved it was not correctly removed from the SandR list, leading to incorrect evolution. Thanks to J Thornburg for finding this one. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@86 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Make the warning message about calling the registration too early morehawke2005-02-01
| | | | | | | verbose, so that you have some idea which thorn did it. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@85 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Catch the case where people that didn't read the documentationhawke2005-01-31
| | | | | | | | properly try to register, e.g., evolved function with MoL before the index arrays have been set up. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@83 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Adaptive step size control using RK45. Due to Erik Schnetter.hawke2005-01-27
| | | | | | | | | | | Note that if you want to use this with Carpet you currently have to use the development (darcs) version of Carpet together with the parameter carpet::adaptive_stepsize = "yes". git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@82 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Initialize some of the grid scalars as early as possible. These arehawke2005-01-07
| | | | | | | | not reset until evolution so NaNChecker was finding NaNs after INITIAL in certain cases (according to Bela Szilagyi). git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@81 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add missing #include <stdio.h> and #include <stddef.h> statementsschnetter2004-12-29
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@80 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Tighten error checking when changing variable type. Certain no-op operations ↵hawke2004-10-15
| | | | | | | | | | | are now safe and potential segfaults caught. Correct warning messages. Fixes PRs 1865,1866. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@77 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct some warning statementsyye002004-10-11
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@76 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* committing the patch sent by Ian Hawke to ChangeType.cyye002004-10-05
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@75 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Runge-Kutta (Fehlberg) 45 with error estimation.hawke2004-09-06
| | | | | | | | | | | | | Fourth order accurate evolution with an additional fifth order step for error estimation. How much sense the error makes is unclear to me, but hey. For the moment the error is stored in an internal MoL array ErrorEstimate; there is one per evolved variable. At a later point this may be moved out to user thorns who can register their own etc. As the implementation uses 6 evaluations of the RHS (necessary) and 6 levels of scratch space (one more than necessary - laziness kicked in) then this is very expensive. This is a partial fix for PR/1840. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@74 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
* If the new keyword parameter 'verbose' is set to 'register' thenhawke2004-06-21
| | | | | | | report the names of all registered variables. Closes PR 1771. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@69 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Fix memory leak. Hint to self: NEVER use group names when trying tohawke2004-04-21
| | | | | | | extract group information. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@64 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Switch the order of warning messages and arrays accesses so that you get the ↵hawke2004-04-02
| | | | | | WARN(0,...) before the segmentation fault... git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@61 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Reduce warning for not provided aliased function to level 4.hawke2004-04-02
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@60 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Forgot to remove inclusion of Carpet header files.hawke2004-03-30
| | | | | | | Ensure that storage is not switched on for empty ArrayScratchSpace group. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@58 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Use the aliased function to en/dis able prolongation instead of linking ↵hawke2004-03-29
| | | | | | | | | direct through the header file. Add parameter to say if we want to use this function (by default we do, so nothing changes). git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@56 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Include Carpet header fileschnetter2004-03-23
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@55 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Initialise the RHS with zero before scheduling the physics routines.schnetter2004-03-23
| | | | | | | | This is necessary because MoL switches off boundary prolongation when there is mesh refinement, leaving certain grid points undefined. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@54 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Convert to int before printing with %dschnetter2004-03-23
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@52 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
* Nasty bugs in the save and restore variables.hawke2004-01-16
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@47 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Erik's changes to do general dissipation adding.hawke2004-01-12
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@45 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct the way that the FillAllTimelevels routine works, now thathawke2003-12-09
| | | | | | | | | | QueryGroupStorage changed. Should have done this months ago Of course, it would be better if the initial data thorns set up all timelevels themselves, but hey. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@44 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Correct all the warning messages in the registration functions.hawke2003-12-05
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@43 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Check that we don't try to allocate zero size arrays.hawke2003-11-19
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@42 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Stupid bug in change variables to constrained type.hawke2003-11-11
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@41 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Remove "const" that shouldn't have been there.schnetter2003-10-30
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@38 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add missing #include statement.schnetter2003-10-30
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@37 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Add the ability to NaNcheck the RHS GFs.hawke2003-09-04
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@34 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
* Protect against buffer overruns...hawke2003-07-22
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@30 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
* Switch off prolongation in Carpet during evolution - patch from Erik Schnetter.hawke2003-07-21
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@26 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Formatting changes and extra debugging statements only.hawke2003-07-18
| | | | git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@25 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Erik Schnetter's implementation of ICN with averaging, so the intermediate ↵hawke2003-07-18
| | | | | | steps are always at t+dt. git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@24 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
* Optionally don't fill all timelevels after initial data is computed (changes ↵hawke2003-07-17
| | | | | | the default but should have no effect). git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@23 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-29 22:46:07 +0000