aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetReduce
Commit message (Collapse)AuthorAge
...
* CarpetReduce: Make helper functions min/max/etc. privateErik Schnetter2011-12-14
|
* CarpetReduce: Replace vector.at by vector.ATErik Schnetter2011-12-14
|
* CarpetReduce: Handle current time correctly when interpolating in timeErik Schnetter2011-12-14
|
* CarpetReduce: Set up mask from new dh class entriesErik Schnetter2011-12-14
| | | | | Simplify mask setup code and improve its performance by using the new dh classes and their new entries.
* CarpetReduce: Use nan handling from CarpetLibErik Schnetter2011-12-14
|
* CarpetReduce: Update commentErik Schnetter2011-12-14
|
* Re-organise time level handlingErik Schnetter2011-12-14
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Store the current Cactus time (and not a fake Carpet time) in the th "time hiearchy". This removes the now redundant "leveltimes" data structure in Carpet. Add past time levels to th, so that it can store the time for past time levels instead of assuming the time step size is constant. This allows changing the time step size during evolution. Share the time hierarchy between all maps, instead of having one time hierarchy per map. Simplify the time level cycling and time stepping code used during evolution. Improve structure of the code that loops over time levels for certain schedule bins. Introduce a new Carpet variable "timelevel", similar to "reflevel". This also makes it possible to avoid time interpolation for the past time levels during regridding. The past time levels of the fine grid then remain aligned (in time) with the past time levels of the coarse grid. This is controlled by a new parameter "time_interpolation_during_regridding", which defaults to "yes" for backwards compatibility. Simplify the three time level initialisation. Instead of initialising all three time levels by taking altogether three time steps (forwards and backwards), initialise only one past time level by taking one time step backwards. The remaining time level is initialised during the first time step of the evolution, which begins by cycling time levels, which drops the non-initialised last time level anyway. Update Carpet and the mode handling correspondingly. Update the CarpetIOHDF5 checkpoint format correspondingly. Update CarpetInterp, CarpetReduce, and CarpetRegrid2 correspondingly. Update CarpetJacobi and CarpetMG correspondingly.
* CarpetReduce: Remove debugging codeErik Schnetter2011-12-14
|
* CarpetReduce: Correct multi-processor error in initialising the reduction maskErik Schnetter2011-12-14
|
* CarpetReduce: Correct handling buffer zonesErik Schnetter2011-12-14
| | | | | | Omit buffer zones when reducing. Clean up code.
* Combine CarpetLib's INSTANTIATE and Carpet's TYPECASE mechanism into aErik Schnetter2011-12-14
| | | | | | single mechanism provided by CarpetLib. Use this mechanism everywhere.
* Import CarpetErik Schnetter2011-12-14
| | | | Ignore-this: 309b4dd613f4af2b84aa5d6743fdb6b3
* Carpet/*/test/: also activate thorn LoopControl (required by CarpetReduce) ↵Thomas Radke2008-10-30
| | | | in testsuite parfiles
* CarpetReduce: Use LoopControlErik Schnetter2008-10-03
|
* CarpetReduce: Add test cases for weight on boundariesErik Schnetter2008-09-28
|
* CarpetReduce: Correct weight on staggered boundariesErik Schnetter2008-09-26
| | | | | | | Correct the weight on staggered boundaries. Correct empty domains consisting of a single internal grid point, where this grid point should have a weight of 1.
* CarpetReduce: Use C++ cast operatorsErik Schnetter2008-07-15
|
* CarpetReduce: Handle multi-patch boundaries correctlyErik Schnetter2008-07-15
| | | | | Use MultiPatch_GetBoundarySpecification to determine the boundary specifications if it is available; otherwise, use GetBoundarySpecification.
* CarpetReduce: Don't double-count grid pointsErik Schnetter2008-02-26
| | | | darcs-hash:20080226141117-dae7b-9ce593fa0c737c838211bb689937593cb32480ca.gz
* CarpetReduce: Modernise codeErik Schnetter2008-02-19
| | | | | | Declare variables where they are used for the first time. darcs-hash:20080219052703-dae7b-bed8b44832012e1abab2e92dd4b50d79506e8638.gz
* CarpetReduce: Correct error in reduction operationsErik Schnetter2007-08-28
| | | | | | | | | | | | | Make the function initialise and reduce in the reduction operators handle the count argument as well. Introduce "combine" functions into the reduction operators. These combine results from several independent reduction operations. Use the above new functions to combine reductions, correcting an error. darcs-hash:20070828044751-dae7b-02070b6aa39d25fb91481ab6a9b45cf37b10fe2f.gz
* CarpetReduce: Take weights into account when reducing in level modeErik Schnetter2007-08-27
| | | | darcs-hash:20070827201932-dae7b-2129f2816c6dcf448bf2706cf5d49cfc83e8d540.gz
* CarpetReduce: Add support for OpenMPErik Schnetter2007-08-21
| | | | darcs-hash:20070821190958-dae7b-1ed180793dd5d6482f48aa0bbf13bd5523a4b40b.gz
* CarpetReduce: Correct schedulingErik Schnetter2007-08-21
| | | | darcs-hash:20070821190829-dae7b-9efe9944356b64e8c171e6f3917f03dee4ba1deb.gz
* CarpetReduce: Implement pointwise reduction operatorsErik Schnetter2007-05-26
| | | | | | | | | | | | | | | | | | | | | | | | Implement a new set of reduction operators for Carpet. These have the same definitions as the old operators, except that they give the same weight to all grid points. The existing operators give a larger weight to coarser grid points, so that physical reduction operators work correctly, e.g. for calculating the total mass. The new operators have the same name as the old operators plus a prefix "i" (the lower case letter "eye"). For example, the new L2 norm is called "inorm2", while the old norm continues to be called "norm2". The physical (old) operators converge to a result which is independent of the grid structure. The pointwise (new) operators still converge, but to a value which depends on the grid structure. In order to make it possible to compare situations with and without symmetries, boundary points still carry a weight of 1/2 or whatever is appropriate. That means that there exists a simple relationship between the pointwise norms with and without a symmetry boundary. darcs-hash:20070526195016-dae7b-6a7e2a39bebd9c0a6377f908a4d7f446f8288191.gz
* CarpetReduce: Set up the weight function after recoveringErik Schnetter2007-05-16
| | | | | | | Set up the weight function after recovering, since the weight function is not checkpointed. darcs-hash:20070516222541-dae7b-a62fe7ade7eb4ef4a381c5710feee91321605658.gz
* CarpetReduce: Make parameter verbose steerableErik Schnetter2007-04-21
| | | | darcs-hash:20070421040012-dae7b-17a08c5094234b1b07bf2ad275d92ab2c11a4fc1.gz
* CarpetReduce: Allow negative weightsErik Schnetter2007-04-19
| | | | darcs-hash:20070419021335-dae7b-fc08da29dd6e24cb828bac615392dd9635241a8c.gz
* CarpetReduce: Update to recent changesErik Schnetter2007-04-19
| | | | darcs-hash:20070419021320-dae7b-9bdc102a1d6dd1a7f998dc6b8cb446572abb2b7c.gz
* CarpetReduce: Add second argument to static_assertErik Schnetter2007-02-28
| | | | darcs-hash:20070228023136-dae7b-86633994d6887c2dd91d7c57d56bc1db620ed2da.gz
* CarpetReduce: Disable an incorrect special case for the reduction maskErik Schnetter2007-01-12
| | | | | | | | Disable the special case that is supposed to handle empty domains for the reduction mask. The implementation of the test is flawed, not the idea; the test should be re-implemented correctly. darcs-hash:20070112224810-dae7b-b3000e90e618f81b9ec188dea2a3223c01889b10.gz
* CarpetReduce: Correct definition of mymin and mymax; introduce myisnanErik Schnetter2007-01-12
| | | | darcs-hash:20070112224438-dae7b-38b62080071060c00a7bf2be053a584af739991b.gz
* CarpetReduce: Adapt to region_t changesErik Schnetter2007-01-12
| | | | darcs-hash:20070112224413-dae7b-a6a68d471d22ef1f0f1b2a1561ef4b302e017bc1.gz
* CarpetReduce: Schedule also in postregridinitial binErik Schnetter2006-12-06
| | | | darcs-hash:20061206163437-dae7b-c2647071f4257c4a2cee820a470ff6b80c20e106.gz
* CarpetReduce: Formatting changesErik Schnetter2007-01-12
| | | | darcs-hash:20070112224318-dae7b-68b8b4b34cba9b14fff76fc6c7fdbb4032297430.gz
* CarpetReduce: Add grouptype argument when calling mode handling functionsErik Schnetter2006-10-03
| | | | darcs-hash:20061003234056-dae7b-c72aaf74712cba7006df16c3382a6ba54e0dad5a.gz
* CarpetReduce: Clean up source codeErik Schnetter2006-09-25
| | | | darcs-hash:20060925220745-dae7b-eb145a85ef330dbff56b7d1c8844449dff85dd0b.gz
* CarpetReduce: Use HAVE_CCTK_REAL* instead of CCTK_REAL*Erik Schnetter2006-09-25
| | | | | | | Use the macros HAVE_CCTK_REAL* instead of checking whether the macros CCTK_REAL* are defined. darcs-hash:20060925220641-dae7b-008b16676bf9b39fd14d718ddeadb009709b6096.gz
* CarpetReduce: Rename local variables imin/imax to bmin/bmaxErik Schnetter2006-09-25
| | | | | | | When storing the extent of the boundary, use the local variable names bmin/bmax instead imin/imax. darcs-hash:20060925220452-dae7b-7a12b1227486a66e0e4689554627c3f713b9d36b.gz
* CarpetReduce: Explicitly convert from/to CCTK_INT and from/to CCTK_REALErik Schnetter2006-06-13
| | | | | | | Explicitly convert between int and CCTK_INT, and between double and CCTK_REAL, where necessary. darcs-hash:20060613171604-dae7b-6cdecab56f7600ef471621e808847d32e45baa75.gz
* CarpetReduce: Correct signature of scheduled functionErik Schnetter2006-06-07
| | | | darcs-hash:20060607225507-dae7b-31f5f3ab59b713a68644e8d574dc7ae9f9db55aa.gz
* CarpetReduce: Handle groups without storage gracefullyErik Schnetter2006-04-13
| | | | | | | Do not abort when some refinement levels of some groups have no storage. Instead, return an error code to the called. darcs-hash:20060413193236-dae7b-63053a270f4790c9f5608b4fe284a6ad4ab33604.gz
* CarpetReduce: Check whether input and output types agreeErik Schnetter2006-04-13
| | | | darcs-hash:20060413201723-dae7b-59318803a1ad41a05efde70bc212cd9a37c260f8.gz
* CarpetReduce: Cache the grid function index of the weight functionErik Schnetter2005-11-19
| | | | | | Obtaining the weight function's index showed up high in a profile. darcs-hash:20051119213646-dae7b-c0f463975c782b1a6854aa66a84169609efe2443.gz
* CarpetReduce: Add new norm sum_abs_squaredErik Schnetter2005-11-19
| | | | | | | This norm calculates sum_i |x_i|^2, which is an ingredient for the L2 norm. darcs-hash:20051119213505-dae7b-e59c1b2fdb5e24455331c20053cb248cca4e898b.gz
* CarpetReduce: Handle sqrt() of integer values in a more C++ like fashionErik Schnetter2005-11-19
| | | | | | | | | Use static_cast<> instead of C style casts. Do not use the STL's sqrt() function for integer variables because e.g. Intel 7 does not support this. darcs-hash:20051119213126-dae7b-5d33bbb6cfc786b980dd93332289528af6e5977f.gz
* CarpetReduce: Use static_assert to ensure that dim==3Erik Schnetter2005-11-19
| | | | darcs-hash:20051119213612-dae7b-2f590847df48bdaec3e05b248adcdab8f564d647.gz
* CarpetReduce: tag CarpetReduce::weight with 'checkpoint="no"'Thomas Radke2005-09-13
| | | | | | The weight is set up at basegrid so it doesn't need to be checkpointed/recovered. darcs-hash:20050913163020-776a0-a171ae3372b9e38ad521a7db85470a8741088b46.gz
* CarpetReduce: Change to new style C++ header filesErik Schnetter2005-07-26
| | | | darcs-hash:20050726110949-891bb-fe7d442ee41c585b7078fd6746840cd68de1d770.gz
* CarpetReduce: Correct error in calculating the past levels' current timesErik Schnetter2005-07-04
| | | | darcs-hash:20050704162500-891bb-589be0634ace55be5a3a5e0b633b86435c3f847b.gz