aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetSlab
Commit message (Collapse)AuthorAge
* global: Remove $Header$ keywords from Carpet arrangementErik Schnetter2005-01-01
| | | | darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
* CarpetSlab: iterate over length[hdim], not over length[dim]tradke2004-11-17
| | | | darcs-hash:20041117143412-3fd61-afe895a742cf18f818a020ca856750497550a608.gz
* CarpetSlab: Switch to singlemap mode before accessing cctk_gshErik Schnetter2004-11-17
| | | | darcs-hash:20041117024250-891bb-aa469ef874902146d300c229edf6f55ae06096f2.gz
* CarpetSlab: Allow hyperslabbing in level mode when there is only one patchErik Schnetter2004-11-15
| | | | darcs-hash:20041115210251-891bb-6b732110e24423efd30cee728478d883ee0c1c33.gz
* Split hyperslabbing routines into multiple filesErik Schnetter2004-09-18
| | | | darcs-hash:20040918133657-891bb-69e4ea212109578a1e65b2e8b8632bb3dee276ce.gz
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* Change indentation.schnetter2004-08-19
| | | | | | | Change indentation. Introduce a type for the conversion function pointers. darcs-hash:20040819043536-07bb3-9e71a221023019f059d6d8cf6133d2db15937374.gz
* Use "requires thorns" to ensure that the thorns are activated in theschnetter2004-06-08
| | | | | | | Use "requires thorns" to ensure that the thorns are activated in the correct order. darcs-hash:20040608190925-07bb3-70b0ac8ee8e570c0b7fff50b02395a288482257d.gz
* Try to catch insanely large hyperslab specifications earlyschnetter2004-05-27
| | | | darcs-hash:20040527102633-07bb3-0a998692efd2d9d8f926653e4d6988a88753e95b.gz
* no changeschnetter2004-05-21
| | | | darcs-hash:20040521161353-07bb3-88f01ca8c5fb75607ec40b9229f8ea170c1cd144.gz
* Undo interface/capability update, because now Carpet conflicts withschnetter2004-05-04
| | | | | | | | Undo interface/capability update, because now Carpet conflicts with PUGH. Need to sort things out, maybe switch temporarily back to the "requires thorns" mechanism. darcs-hash:20040504205239-07bb3-771c6ec96a2b454d63c87a6dc7a1f44f62cd8fd7.gz
* Begin to clean up inheritance and requires relationshipsschnetter2004-05-04
| | | | darcs-hash:20040504200954-07bb3-cd9f17c981ce37247fa2122858b10b46425f678a.gz
* Use "requires" instead of "requires thorns".schnetter2004-04-18
| | | | darcs-hash:20040418120959-07bb3-05050a71e475171cac6b293de35118e4af8df5ba.gz
* Some more cleanup and bug fixes. Hyperslabbing probably never workedschnetter2004-04-16
| | | | | | | Some more cleanup and bug fixes. Hyperslabbing probably never worked with Carpet. darcs-hash:20040416120226-07bb3-816cca419723e96b25dd1d17218e4068738d889d.gz
* Allow the downsampling argument to be NULL.schnetter2004-04-16
| | | | darcs-hash:20040416095218-07bb3-19ec29d50cb450c06a9d8ca276f3e3fd706cd6c2.gz
* Undo last change. Was a stupid idea.schnetter2004-04-01
| | | | darcs-hash:20040401070546-07bb3-86d0c146a8a5007f00edc360d5d416ff40a8d7d5.gz
* Provide the capability Hyperslab.schnetter2004-04-01
| | | | darcs-hash:20040401070253-07bb3-61fecc73140ff36fb597d3dd9694c2cdb410e582.gz
* Use vector::at() instead of vector::operator[] to access vector elements.schnetter2004-03-23
| | | | darcs-hash:20040323183014-07bb3-4562b1a7b1dda25fbcccc7689c19bbd432b8a1ff.gz
* Import the recently announced changes:schnetter2004-01-25
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Import the recently announced changes: 1. Carpet has now an infrastructure for multiple maps (aka "grid patches"). Instead of a single grid hierarchy there can now be several. This is largely untested, because the remainder of Cactus cannot handle multiple coordinate systems. 2. The order in which the schedule bins are called has changed. As Ian Hawke pointed out, the previous order during time evolution was inconsistent. The initial data ordering did not allow for recovering and was not usable for progressively solving elliptic equations for initial data. 3. Carpet now supports convergence levels. The convergence level specifies by how many factors of two the resolution in the parameter file should be coarsened (or refined, if negative). This should make convergence tests and test runs much easier. It is, in principle, also possible to run several convergence levels at once. This has not been tested because the remainder of Cactus cannot handle multiple resolutions. This will be necessary for a multigrid solver, and also for having a shadow hierarchy to determine where to refine adaptively. 4. Carpet works together with the new CoordBase domain specification parameters. Without these, using convergence levels will lead to very strange results. 5. The "modes" have changed. There are now: meta mode: the whole simulation global mode: one convergence level level mode: one refinement level singlemap mode: one map on one refinement level local mode: as previously The whole mode handling has been cleaned up. 6. The regridding thorn has been cleaned up. 7. The kind of prolongation stencil is now determined in Carpet, i.e. at a fairly hight level, instead of in CarpetLib. 8. The low-order prolongation operators have been made much more efficient (as have previously the higher-order ones). 9. Assorted smaller changes. For Carpet users, there should be no major incompatibilities. The major improvements are 3 and 4 combined. Here is an example: CoordBase::domainsize = extent CoordBase::spacing = gridspacing CoordBase::zero_origin_x = yes CoordBase::zero_origin_y = yes CoordBase::zero_origin_z = yes CoordBase::xextent = 20.0 CoordBase::yextent = 20.0 CoordBase::zextent = 20.0 CoordBase::dx = 1.0 CoordBase::dy = 1.0 CoordBase::dz = 1.0 CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_shiftout_y_lower = 1 CoordBase::boundary_shiftout_z_lower = 1 Carpet::domain_from_coordbase = yes Carpet::convergence_level = 0 grid::type = coordbase grid::domain = octant grid::avoid_origin = no This gives you a grid that extends from the origin ("zero_origin") up to 20.0 with a grid spacing of 1.0. Symmetry zones and boundary zones are added automatically. The "shiftout" says that there is no boundary point on the origin. The staggering parameters (not shown) default to "no". In order to change the resolution, only the convergence level has to be adjusted. Note that the old way of specifying the domain extent still works. For Carpet developers, one major change is the new mode handling. As described in 5, the looping macros (that loop over all refinement levels, or all components) have changed. darcs-hash:20040125135727-07bb3-51c9647c1b5080e7e180b52a1b81fa155cfd19e9.gz
* Import the recently announced changes:schnetter2004-01-25
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Import the recently announced changes: 1. Carpet has now an infrastructure for multiple maps (aka "grid patches"). Instead of a single grid hierarchy there can now be several. This is largely untested, because the remainder of Cactus cannot handle multiple coordinate systems. 2. The order in which the schedule bins are called has changed. As Ian Hawke pointed out, the previous order during time evolution was inconsistent. The initial data ordering did not allow for recovering and was not usable for progressively solving elliptic equations for initial data. 3. Carpet now supports convergence levels. The convergence level specifies by how many factors of two the resolution in the parameter file should be coarsened (or refined, if negative). This should make convergence tests and test runs much easier. It is, in principle, also possible to run several convergence levels at once. This has not been tested because the remainder of Cactus cannot handle multiple resolutions. This will be necessary for a multigrid solver, and also for having a shadow hierarchy to determine where to refine adaptively. 4. Carpet works together with the new CoordBase domain specification parameters. Without these, using convergence levels will lead to very strange results. 5. The "modes" have changed. There are now: meta mode: the whole simulation global mode: one convergence level level mode: one refinement level singlemap mode: one map on one refinement level local mode: as previously The whole mode handling has been cleaned up. 6. The regridding thorn has been cleaned up. 7. The kind of prolongation stencil is now determined in Carpet, i.e. at a fairly hight level, instead of in CarpetLib. 8. The low-order prolongation operators have been made much more efficient (as have previously the higher-order ones). 9. Assorted smaller changes. For Carpet users, there should be no major incompatibilities. The major improvements are 3 and 4 combined. Here is an example: CoordBase::domainsize = extent CoordBase::spacing = gridspacing CoordBase::zero_origin_x = yes CoordBase::zero_origin_y = yes CoordBase::zero_origin_z = yes CoordBase::xextent = 20.0 CoordBase::yextent = 20.0 CoordBase::zextent = 20.0 CoordBase::dx = 1.0 CoordBase::dy = 1.0 CoordBase::dz = 1.0 CoordBase::boundary_shiftout_x_lower = 1 CoordBase::boundary_shiftout_y_lower = 1 CoordBase::boundary_shiftout_z_lower = 1 Carpet::domain_from_coordbase = yes Carpet::convergence_level = 0 grid::type = coordbase grid::domain = octant grid::avoid_origin = no This gives you a grid that extends from the origin ("zero_origin") up to 20.0 with a grid spacing of 1.0. Symmetry zones and boundary zones are added automatically. The "shiftout" says that there is no boundary point on the origin. The staggering parameters (not shown) default to "no". In order to change the resolution, only the convergence level has to be adjusted. Note that the old way of specifying the domain extent still works. For Carpet developers, one major change is the new mode handling. As described in 5, the looping macros (that loop over all refinement levels, or all components) have changed. darcs-hash:20040125125727-07bb3-3368611314b2dcb8c8ae58ab3f501b683d7edb8f.gz
* Renamed some aliased hyperslab functions to comply with the maximum lengthtradke2004-01-22
| | | | | | | | | | | | | for names of Fortran routines (31 characters): Hyperslab_DefineGlobalMappingByPhys -> Hyperslab_GlobalMappingByPhys Hyperslab_DefineGlobalMappingByIndex -> Hyperslab_GlobalMappingByIndex Hyperslab_DefineLocalMappingByPhys -> Hyperslab_LocalMappingByPhys Hyperslab_DefineLocalMappingByIndex -> Hyperslab_LocalMappingByIndex Also renamed the providing functions accordingly. darcs-hash:20040122123107-1d9bf-95166b54dcc264e30c13bacbf58c094a25b34dc8.gz
* Many changes that accumulated while Cactus and Carpet diverged.schnetter2003-11-05
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Many changes that accumulated while Cactus and Carpet diverged. Add processor splitting mechanism "along-dir" that splits along a specified direction. Rename group PostRestrict to bin POSTRESTRICT. Prolongate initial data only when desired. This saves much time. Sorry, Ian. Fix bug in time level cycling of grid arrays. (Note: grid arrays should not have time levels.) Fix time_t bug on IRIX. Make sure that there is no integer overflow when there are many refinement levels. Always put parentheses around (maxreflevelfact/reflevelfact). Fix typo in Carpet verbose output. Add debug output in processor splitting. Communicate in three stages: Irecv, (work), Isend, Wait. This might be more efficient. Much more, potentially. Fix bug in processor layout of grid arrays. Sorry, Ian. Make the interpolator interpolate between time levels. Untested. Fix bug in processor communication in interpolator. Sorry, Ian. Rewrite prolongation operators to make them twice as fast. There you are, Ian. Move prolongation operator kind handling from data to gdata. Add official hyperslabbing interfaces to CarpetSlab. Adapt to new cGH * handling. darcs-hash:20031105151837-07bb3-758a87ff0355dba053269df4b7d7d79bea018669.gz
* Many changes that accumulated while Cactus and Carpet diverged.schnetter2003-11-05
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Many changes that accumulated while Cactus and Carpet diverged. Add processor splitting mechanism "along-dir" that splits along a specified direction. Rename group PostRestrict to bin POSTRESTRICT. Prolongate initial data only when desired. This saves much time. Sorry, Ian. Fix bug in time level cycling of grid arrays. (Note: grid arrays should not have time levels.) Fix time_t bug on IRIX. Make sure that there is no integer overflow when there are many refinement levels. Always put parentheses around (maxreflevelfact/reflevelfact). Fix typo in Carpet verbose output. Add debug output in processor splitting. Communicate in three stages: Irecv, (work), Isend, Wait. This might be more efficient. Much more, potentially. Fix bug in processor layout of grid arrays. Sorry, Ian. Make the interpolator interpolate between time levels. Untested. Fix bug in processor communication in interpolator. Sorry, Ian. Rewrite prolongation operators to make them twice as fast. There you are, Ian. Move prolongation operator kind handling from data to gdata. Add official hyperslabbing interfaces to CarpetSlab. Adapt to new cGH * handling. darcs-hash:20031105141838-07bb3-d348d033307d82905893e241be485c4db2ce8ce3.gz
* Make gf and data objects store a Cactus variable index that thisschnetter2003-10-14
| | | | | | | | Make gf and data objects store a Cactus variable index that this object is associated with. This can be used to access various kinds of information, e.g. variable names. darcs-hash:20031014143916-07bb3-678d114393008db7790b5ed72d6462673414c06a.gz
* Replace num_points() by size().schnetter2003-09-20
| | | | darcs-hash:20030920115318-07bb3-50458524a69bce820359042b32f6ee9696bb903e.gz
* Carpet infrastructure thorns now require the thorns they depend upontradke2003-09-04
| | | | | | | | (CarpetLib, Carpet, FlexIO) to be compiled in and activated via the "Requires thorns: ..." attribute in the configuration.ccl files. No need to inherit from those thorns/implementations anymore. darcs-hash:20030904142321-1d9bf-309be01f0950d97b8cb45d041d44edce01aff2a5.gz
* Do not include header files from other Carpet thorns directly.schnetter2003-06-18
| | | | | | | | Do not include header files from other Carpet thorns directly. Instead, use the "INCLUDES HEADER" and "USES INCLUDE HEADER" mechanism. darcs-hash:20030618162807-07bb3-a81444cde6c76e6a24516d108861fc1b5541c643.gz
* Major update after a quiet time.schnetter2003-06-18
| | | | | | | | | | | | | | | | | | | | Major update after a quiet time. Carpet: The flesh now has new cGH fields cctk_levoff[], cctk_levoffdenom[], and cctk_timefac that describe the spatial offset and temporal refinement factor between the base and the current refinement level. These fields are now set and used; they change how coordinates are handled. CarpetIOASCII: Fix bugs regarding choosing the output hyperslab and the output coordinates. ID*, *Toy*: New WaveToy examples with various formulations and different integrations methods. Currently, none of them converge to second order except the standard WaveToy formulation. These updates require the recent flesh, base thorn (and MoL) updates. darcs-hash:20030618162427-07bb3-70761f74bce6ae246b5a2943a385647657d46d34.gz
* Allow hyperslabbing grid arrays with dim<3.schnetter2003-05-21
| | | | darcs-hash:20030521123129-07bb3-6b61f6f1d8bfba3ae3929c4feec24a13185c7a04.gz
* Make code const-correct. Well, sort of.schnetter2003-05-13
| | | | darcs-hash:20030513101942-07bb3-2a9b4aa40aafc2566a13522ff7704e480e9817b6.gz
* Handle empty grids.schnetter2003-05-13
| | | | | | | | | | Handle empty grids. Fix global and level mode handling of the interpolator. Allow local mode calls for the interpolator, reduction operator, and hyperslabber, assuming that people know what they are doing. darcs-hash:20030513101400-07bb3-42af7544c2ea98dbef5a7ad6d5112bea2dcc404f.gz
* Abort when in global mode.schnetter2003-05-12
| | | | darcs-hash:20030512142511-07bb3-2931d8eaaf3f69dd17b70b1f5f6cc25d53541e7d.gz
* Replaced new[] and delete[] with a more elegant use of vector<>.schnetter2003-02-28
| | | | darcs-hash:20030228090927-07bb3-d089883a3f0b5f7704418bc4995a89b3c7ddf696.gz
* Removed dependence on alloca, which isn't standard and isn't supportedschnetter2003-02-28
| | | | | | | Removed dependence on alloca, which isn't standard and isn't supported on all systems. darcs-hash:20030228085604-07bb3-f55515e91912d6e5f696cb669d08ccbd4d7a1bef.gz
* Removed the "dimgeneric_*" data types from CarpetLib. They were neverschnetter2003-01-03
| | | | | | | | | | | | | Removed the "dimgeneric_*" data types from CarpetLib. They were never really used anyway. Changed "th" into "th<D>". Renamed "generic_*" to "g*". Propagated these changes to all other places. darcs-hash:20030103144936-07bb3-c4d74d2d94a56a898f9bb735666eefc01e2610e2.gz
* Updated parameter files.schnetter2002-10-24
| | | | darcs-hash:20021024191153-07bb3-90dd8efad24bdb5625a2365fba0a57ea96adc14d.gz
* Renamed source file from carpetslab.* to slab.*.schnetter2002-10-24
| | | | | | | Renamed source file from carpetslab.* to slab.*. Added arrangement name and final semicolon to CCTK_FILEVERSION macros. darcs-hash:20021024085347-07bb3-f25e59ad814c9143dbe65face1f5f273e1a2a3cc.gz
* Fixed quite a few include guard mishaps.schnetter2002-09-01
| | | | darcs-hash:20020901125223-07bb3-f494618480c1ffc46118365c4c957d7a313d7be9.gz
* Changed "abort()" statements to "assert(0)" statements, because theschnetter2002-08-02
| | | | | | | Changed "abort()" statements to "assert(0)" statements, because the latter do output a message to the screen. darcs-hash:20020802131112-07bb3-a89ca2582fce00c0d5b4cc82bf61473b7c8bdcb5.gz
* Added CCTK_FILEVERSION macros.schnetter2002-03-26
| | | | darcs-hash:20020326122225-07bb3-0c7746df97fb64f00918be1bee33191fb9d44b0b.gz
* Fixed a bug in CarpetSlab; was using the local instead of the globalschnetter2002-03-23
| | | | | | | | | | | | | | | | Fixed a bug in CarpetSlab; was using the local instead of the global extent. Added a more comfortable way to specify the refined region. The new method has no restrictions. Added a way to specify the location of the outer boundary explicitely. Added dummy parameter "periodic*" so that Carpet can be used as drop-in replacement of PUGH (at least when running without refinement). darcs-hash:20020323192054-07bb3-2b9208c8a6b54a55c98d4035095cd4ff6ebcd447.gz
* Corrected interface.ccl files which were missing inter-thorn dependencies.schnetter2001-07-12
| | | | | | | Corrected interface.ccl files which were missing inter-thorn dependencies. Removed legacy .deps files. darcs-hash:20010712041104-07bb3-1fe3f038149a3b61cfdf1c20f633cc741a5cc2d9.gz
* Changed handling of scalars. Scalars are now zero-dimensional arrays.schnetter2001-07-09
| | | | | | | | | Changed handling of scalars. Scalars are now zero-dimensional arrays. Now handling CCTK_GroupDynamicData correctly. Now using "include header". Added results of test case. darcs-hash:20010709070002-07bb3-ba0f4339acb8652e35a884fbdf7e8482b7236f8e.gz
* Broke the file "carpet.cc" into several files.schnetter2001-07-04
| | | | | | | | Broke the file "carpet.cc" into several files. Simplified support for arrays with a dimension different from 3. Added a testing thorn. darcs-hash:20010704102946-07bb3-82132f0c643e91f0de33bbaef93f9c693ce7227f.gz
* Second stab at having arrays with a dimension different than gridschnetter2001-07-02
| | | | | | | Second stab at having arrays with a dimension different than grid functions. Not yet tested. darcs-hash:20010702112209-07bb3-816053f629c333189a63014698d4b34c124168aa.gz
* Adapted Carpet to the new timelevel scheme in Cactus.schnetter2001-05-16
| | | | darcs-hash:20010516122945-07bb3-0ff1f074df5f2e304820f3a662e3260760b5cd4b.gz
* Modified the Carpet driver so that it uses the same new timelevelschnetter2001-04-23
| | | | | | | | | | | | | | | | | | Modified the Carpet driver so that it uses the same new timelevel cycling as does PUGH. Modified the CarpetIOASCII and CarpetIOFlexIO output routines so that they always output the current time level. Found and fixed a bug in the time interpolation in CarpetLib. Fixed a bug where the necessary number of ghost zones for space interpolation was calculated incorrectly. Fixed bugs in the example parameter files that used an insufficient number of ghost zones. darcs-hash:20010423061011-07bb3-c69e59f7bc5f283d7632a663c0770578cdb73881.gz
* Brought in latest differences from the SGI version. This is workeschnett2001-03-22
| | | | | | | Brought in latest differences from the SGI version. This is work towards a code that compiles on both architectures. darcs-hash:20010322174200-f6438-23ab5f26cf84d2666312791c6bdb5a0fc1d0390a.gz
* There is a difference between Cactus timelevels and Carpet timelevels.eschnett2001-03-17
| | | | darcs-hash:20010317151452-f6438-567639b4b41992db94ec23419f077aa4f3e8d9ac.gz
* Carpet now passes the Cart3d test suite. This required a rewrite ofeschnett2001-03-10
| | | | | | | | | | | | | | | | | Carpet now passes the Cart3d test suite. This required a rewrite of CarpetSlabe, which now uses the data<> class instead of doing things its own way. This in turn required some updates to data<>. CarpetSlab now ignores (i. e. does not use) boundaries, except outer boundaries. For that to work, the dh<> class has been given the notion of an outer boundary. In order to increase performance, the cctk_bbox[] information in Carpet is not set correctly. As an approximation, it just reflects whether the current grid component lies at the outer boundary, as calculated using cctk_lbnd, cctk_ubnd, and cctk_gsh. darcs-hash:20010310195459-f6438-8ee874d722de4df717e5eb3002722ceb8d4a9bb7.gz