aboutsummaryrefslogtreecommitdiff
path: root/Carpet
Commit message (Collapse)AuthorAge
* 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
* Modified Carpet2ygraphCat. Works with progressive and moving meshes. Wastes ↵hawke2004-01-23
| | | | | | a lot more memory and disk space than it used to. darcs-hash:20040123074040-58737-eef0cd6ebc85e62611a2f15ab5eb0b927a4bc98a.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
* Undo my previous changes, which did not compile.schnetter2004-01-21
| | | | darcs-hash:20040121153204-07bb3-2737dc30c6c1b867713208ee5dfe49085ba1b625.gz
* Fix another (maybe) bug in calculating the prolongation regions.schnetter2004-01-21
| | | | darcs-hash:20040121132535-07bb3-8221de1d3146f586eb13b0fe3acb12836b06efd5.gz
* Check the stencil size before calling the minmod prolongationschnetter2004-01-21
| | | | | | | Check the stencil size before calling the minmod prolongation operators. darcs-hash:20040121132509-07bb3-7eb951d88ae3b66a3ad2338fecc632c8c77cb312.gz
* Use .at() instead of [] to access vectors.schnetter2004-01-20
| | | | | | | | Use .at() instead of [] to access vectors. Fix a bug in the calculation of the restriction regions. Remove unused code. darcs-hash:20040120122505-07bb3-4e4c8c59f712ebe671d261ef1bea3640ab077ada.gz
* Slight wording changesschnetter2004-01-18
| | | | darcs-hash:20040118102930-07bb3-1ed8463535f2d2e32616e985598c3d588a1f70a3.gz
* Import Carpet web pages into the Carpet repositoryschnetter2004-01-17
| | | | darcs-hash:20040117084347-07bb3-96e4d27e49955470253973d7bd8a7642fa8d5228.gz
* While restricting, don't assume that there are buffer zones atschnetter2004-01-16
| | | | | | | While restricting, don't assume that there are buffer zones at interprocessor boundaries. darcs-hash:20040116163644-07bb3-7eef930c1b8bb461228d194dcb9ecddd5534d9ad.gz
* EnableGroupStorage now returns the number of active timelevels before the ↵hawke2004-01-16
| | | | | | call. Note that this is always the maximum number of timelevels as Carpet always switches on all storage. darcs-hash:20040116095704-58737-0b01f91957109646209f596a26c1dd4a2031389d.gz
* Changed interface to aliased function "RegridLevel" - added parameter ↵cott2004-01-15
| | | | | | maxreflevels darcs-hash:20040115084558-19929-5187725fb081ea4da02a6805f2cc65ecbfdf4244.gz
* Correct location of 3 timelevel initialisationschnetter2004-01-14
| | | | darcs-hash:20040114081408-07bb3-65a7482594808f8bca77df73d41ff1763177c541.gz
* The scheduling tree for PUGH and Carpet for comparison. For PUGH,schnetter2004-01-13
| | | | | | | | The scheduling tree for PUGH and Carpet for comparison. For PUGH, this is the real (and not the documented) tree. For Carpet, however, this is an ideal tree and not the real one... darcs-hash:20040113150542-07bb3-7647db7e87e969da49608ca07211d5b70fca3af6.gz
* Added STEERABLE = ALWAYS to parameter output_all_timelevelscott2004-01-13
| | | | darcs-hash:20040113142658-19929-bdda00a3cb2bc7e14c7111108e8e471892ca99cb.gz
* Move the regridding function to a forward loop where the iteration %hawke2004-01-13
| | | | | | | Move the regridding function to a forward loop where the iteration % regrid_every is zero instead of 1... darcs-hash:20040113125119-58737-8847453265d1b17dcd8343149839a270afc8a95e.gz
* Changes to CarpetRegrid.hawke2004-01-13
| | | | | | | | | | Changes to CarpetRegrid. - Does progressive MR correctly again. - Will take a criteria for the number of levels to activate from an aliased function if requested, allowing for adaptive refinement in time. darcs-hash:20040113125005-58737-15008754eaae7b272c5a144333241be8ab6bbde2.gz
* Clarify error message.schnetter2004-01-09
| | | | darcs-hash:20040109122120-07bb3-ec22ca801d35b1de50585080dd2ab0e8d31da4d5.gz
* Fix typo in parameter descriptionschnetter2003-12-19
| | | | darcs-hash:20031219150155-07bb3-41d0f7f3df186e02d5a3a98ad823906252564440.gz
* test new mailing listschnetter2003-12-19
| | | | darcs-hash:20031219104155-07bb3-c74a7c6ed129c05fe6687c4ac627874a93cda3d5.gz
* Check new mailing list.schnetter2003-12-19
| | | | darcs-hash:20031219103309-07bb3-023267a17f4d7bcbd2399cedefc787cc228a9539.gz
* undo testing changeschnetter2003-12-19
| | | | darcs-hash:20031219094352-07bb3-20d9e41ca0bb17d441f83f83edcab34739ef1fde.gz
* Follow C++ standard.schnetter2003-11-21
| | | | darcs-hash:20031121125546-07bb3-60338023244fbd0fab9d8663adff4cdf78e31141.gz
* Round manually specified coordinates to the nearest grid point, butschnetter2003-11-21
| | | | | | | Round manually specified coordinates to the nearest grid point, but don't round strides. darcs-hash:20031121115156-07bb3-651b99788e8699892ef5e87ece189d4fe95484d0.gz
* Missing break in TVD interpolator case statement.hawke2003-11-20
| | | | darcs-hash:20031120072851-58737-e1aef776341af758afc99b27f175463b3ad508c2.gz
* Correct dependencies.schnetter2003-11-19
| | | | darcs-hash:20031119161203-07bb3-48f4411b33dbdf6dffc24fd74dd61a5fb4634f4f.gz
* Only allow MG level switches in global mode.schnetter2003-11-19
| | | | darcs-hash:20031119130607-07bb3-1a512c150eba19d5447181d614c4ae0d63b7f4f2.gz
* Allow reflevel=-1 in set_mglevel if we're going to global mode.hawke2003-11-19
| | | | darcs-hash:20031119084938-58737-5affcb34e6214f88010808e33f6cd9e0d9ffb443.gz
* More error checking when converting real coordinates to integer grid points.schnetter2003-11-14
| | | | darcs-hash:20031114115047-07bb3-e2a2aa933f412bf4b9b6d9334d3774203a80b569.gz
* Return the local interpolator's error code if something goes wrong.schnetter2003-11-13
| | | | darcs-hash:20031113162947-07bb3-267cc4954cf1498f2a1f058eb4c23904976e92e9.gz
* Implement specifying refined boxes by coordinates instead of by gridschnetter2003-11-13
| | | | | | | Implement specifying refined boxes by coordinates instead of by grid point numbers. darcs-hash:20031113150437-07bb3-3889e9e5e4f97233034f57ab4f968b1476bd515f.gz
* Allow bboxes of CCTK_REALs.schnetter2003-11-13
| | | | darcs-hash:20031113150358-07bb3-427886aab7c861fe53bf7fef96ec2a31e1d155f5.gz
* Output a real error message if a file cannot be opened for writing.schnetter2003-11-13
| | | | darcs-hash:20031113150338-07bb3-d614c53c250f86c87ead3aa8a7b97c74439b0c04.gz
* Allow arrays with fewer elements than the number of processors.schnetter2003-11-13
| | | | darcs-hash:20031113094917-07bb3-9890909505d88e3f6d0d18b109065d05098d5a4b.gz
* Fix bug in the maximum reduction operator.schnetter2003-11-12
| | | | darcs-hash:20031112162930-07bb3-eb9a8b4a6c167c98bfeca86bc0c7b10baed922b5.gz
* Output better error message.schnetter2003-11-11
| | | | darcs-hash:20031111132341-07bb3-6fdb95d03d6c9d1fb857557afd17b6418a17dc7e.gz
* Use vector<> instead of dynamic [] arrays.schnetter2003-11-07
| | | | darcs-hash:20031107122106-07bb3-a21c30bb42714054caed0490f79121cfef385697.gz
* Return local interpolation errors to the caller instead of aborting.schnetter2003-11-06
| | | | darcs-hash:20031106133319-07bb3-f59f24cdc8a5f5b170693e56b1b8f2a14b1adf0e.gz
* Fix core dump.schnetter2003-11-06
| | | | darcs-hash:20031106125023-07bb3-bc17ef32006900de77a8f185a7424322df71f9db.gz
* Update doxygen configuration.schnetter2003-11-05
| | | | darcs-hash:20031105190135-07bb3-0ee85486318563b8fe42e7a5c8de220725d2b7fd.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
* Fixed minmod interpolation loopcott2003-10-30
| | | | darcs-hash:20031030125459-19929-40e83a3467ca581eb7e493404657b56b3a78de2a.gz
* Patch the minmod prolongation routines in a fashion that really appears to ↵hawke2003-10-20
| | | | | | work for a change. darcs-hash:20031020083924-58737-0808ad70950041cf90593f8e63fc5dd8a552b17d.gz
* Fixed bug in case PROLONG_TVD (missing break statement in case 1).cvs_anon2003-10-17
| | | | | | | Fixed bug in case PROLONG_TVD (missing break statement in case 1). - Christian darcs-hash:20031017061441-62b91-caf11a681cb9c5d39858e666c1d281efb56d00b2.gz
* Fix bug in automatic domain decomposition.schnetter2003-10-16
| | | | darcs-hash:20031016150004-07bb3-9c843c9710adcce8d2b3243b0c02b798d83bab59.gz
* A much corrected version of the TVD prolongation. Looks nicer in the first ↵hawke2003-10-15
| | | | | | tests. Also simpler. darcs-hash:20031015095155-58737-cdbc9b0c881bc9970298ad5c95d5375a1eacaaab.gz
* The prolongation method can now be controlled from the tags table defined in ↵hawke2003-10-15
| | | | | | | | | | | | | | | | | | the interface.ccl. The prolongation method defaults to "Lagrange" and the other options are "None" and "TVD" (which I'm now sure doesn't work properly). The prolongation method can now be controlled from the tags table defined in the interface.ccl. The prolongation method defaults to "Lagrange" and the other options are "None" and "TVD" (which I'm now sure doesn't work properly). As an example, from CartGrid3D: REAL coordinates type=GF tags='Prolongation="None"' { x, y, z, r } switches off prolongation for the coordinates. Also, a one off opportunity to marvel at C++ code that combines the worst features of Fortran with the worst features of C... darcs-hash:20031015051401-58737-870c23e54e58164a45ffe3438e393587b1cc670b.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