aboutsummaryrefslogtreecommitdiff
path: root/Carpet
Commit message (Collapse)AuthorAge
* Use new bbox iterators in all places (missed some before).schnetter2004-02-27
| | | | darcs-hash:20040227144040-07bb3-a89704193701abcae95b27068427adc63cbc76b3.gz
* Use new bbox iterators in all places (I missed some before).schnetter2004-02-27
| | | | darcs-hash:20040227144018-07bb3-9fe22a2e0c96216fe30daae38b21e4a474d03e40.gz
* Add a dummy array element when D==0.schnetter2004-02-27
| | | | darcs-hash:20040227092523-07bb3-e4cc3ad72fab0773cc5f42d1af3a27cf61d01545.gz
* Remove #warning statement.schnetter2004-02-18
| | | | | | | Remove #warning statement. Use cGH entries instead of Carpet global variables. darcs-hash:20040218141838-07bb3-1644927d8d2948c7742c4cac9b3ab7d93193858f.gz
* Add 0-dimensional ("point") output. This is necessary to output gridschnetter2004-02-18
| | | | | | | | | | | Add 0-dimensional ("point") output. This is necessary to output grid scalars. Make the output triggering on all refinement levels be independent. Honour do_global_mode when outputting grid scalars and grid arrays. darcs-hash:20040218141229-07bb3-87d51a3ac71175ba7d190f45dd0d735e81ff0350.gz
* Add functions for element-wise floor, ceil, and pow.schnetter2004-02-18
| | | | darcs-hash:20040218141045-07bb3-0e60f56d44aa61e0e3190d5d27b9a0bfab3f2d38.gz
* Change definition for the bbox iterators, so that they work withschnetter2004-02-18
| | | | | | | | | Change definition for the bbox iterators, so that they work with 0-dimensional bboxes as well. Instantiate the templates for 0 dimensions. darcs-hash:20040218141001-07bb3-df6baa90f8b8ae69a788d863ef663d55587d1d25.gz
* Instantiate vectors for 0 dimensions as well.schnetter2004-02-18
| | | | darcs-hash:20040218131017-07bb3-f84c3cccf65509ed31d1951791368bc77c8845de.gz
* Print the grid structure in both grid points and coordinate rangesschnetter2004-02-15
| | | | | | | Print the grid structure in both grid points and coordinate ranges after regridding. darcs-hash:20040215135837-07bb3-17e0217fdff9984b9299490f96e1d4bc2472fa7d.gz
* Run global mode functions at the last possible point when they areschnetter2004-02-15
| | | | | | | | | | Run global mode functions at the last possible point when they are called after the recursion. Set the global time by a direct multiplication instead of incremental addition. This is necessary to avoid accumulation of roundoff errors. darcs-hash:20040215102751-07bb3-f2eeb8581239b1d30ed9286b8110451fa576105b.gz
* Run global mode functions at the last possible point when they areschnetter2004-02-15
| | | | | | | Run global mode functions at the last possible point when they are called after the recursion. darcs-hash:20040215102718-07bb3-edff79388446558ced60f3f3a725d29cc8eba0ab.gz
* When interpolating in time, call the local interpolator many timesschnetter2004-02-15
| | | | | | | | | | When interpolating in time, call the local interpolator many times instead of only once. This is slower, but it avoids having to rewrite all input arguments -- and there are many input arguments, most of them hidden in the options table. In fact, Carpet can not even know how many input arguments there are. darcs-hash:20040215101516-07bb3-fac5be0b7f5557733f81b66cf70cea86e79f67ec.gz
* Fix array size allocation.schnetter2004-02-15
| | | | darcs-hash:20040215093414-07bb3-1b9b2a4fd94a77be876fef35aa15a2eef7d55941.gz
* Use (rl-minrl) instead of rl to index into the count arrays (that areschnetter2004-02-14
| | | | | | | | Use (rl-minrl) instead of rl to index into the count arrays (that are used to communicate which processor works on how many grid points on what levels). darcs-hash:20040214140929-07bb3-e0305bc753a075f79d1d476c22b5773e7e6edd88.gz
* Report a warning instead of aborting when the RegridLevel functionschnetter2004-02-09
| | | | | | | Report a warning instead of aborting when the RegridLevel function returns an illegal number of levels. darcs-hash:20040209154843-07bb3-f92482392a8aa2b703a120e088a7b668cf08abae.gz
* Add support for copying (i.e. synchronising) complex variables.schnetter2004-02-09
| | | | darcs-hash:20040209135807-07bb3-0cd2c04df4f1ef2fa2e5d31f956c8b57b38a38b2.gz
* Fig bug in output: In C, quotes are not quoted by "", but by \".schnetter2004-02-09
| | | | darcs-hash:20040209135646-07bb3-2d7814ee1b34bbb2b8b3751732ffe3e8dcac8922.gz
* Print an error message when trying to synchronise an unsupported dataschnetter2004-02-09
| | | | | | | Print an error message when trying to synchronise an unsupported data type. darcs-hash:20040209135507-07bb3-eedaedfeab593e73a8c96e2d08745c5f2dc3d08f.gz
* Allow output for grid variables with just one grid point in certainschnetter2004-02-09
| | | | | | | | Allow output for grid variables with just one grid point in certain directions. This grid point will get the lower boundary of the coordinate system. darcs-hash:20040209135302-07bb3-6c022e1c5e263535fa22fe6b8f368a9e9b1889ec.gz
* Output an error message when a DISTRIB=constant group has ghost zones.schnetter2004-02-09
| | | | darcs-hash:20040209121033-07bb3-380bfaa116ba85908947c5b0f19c9d840c105db6.gz
* Always output the new grid structure after recomposing, and output itschnetter2004-02-09
| | | | | | | Always output the new grid structure after recomposing, and output it in a nicer format. darcs-hash:20040209120835-07bb3-5a740a76e4e7611461cf0f8d4b187e2f2aade9ea.gz
* Put group names in quotes when printing them to the screen.schnetter2004-02-09
| | | | darcs-hash:20040209120605-07bb3-ec1889c9298509d1807433fcb295041221c3991e.gz
* Initialise always from coarse to fine grid, instead of recursively.schnetter2004-02-09
| | | | | | | | Initialise always from coarse to fine grid, instead of recursively. Going from coarse to fine grid is necessary so that one can prolongate in the postinitial and poststep bin. darcs-hash:20040209120231-07bb3-8001c3c8a205be0ec4b52a2d0b061ec50de1b1f8.gz
* Do not output warning messages twice.schnetter2004-02-09
| | | | darcs-hash:20040209115953-07bb3-daa126dc4bb5befa183d9d22688e622ea67ccdb5.gz
* Enter and leave global and level mode only when this is necessary atschnetter2004-02-09
| | | | | | | | Enter and leave global and level mode only when this is necessary at the current iteration. This speeds up Carpet, especially when less than maxreflevels levels are used. darcs-hash:20040209115800-07bb3-9f29de569eab9b34185fe707ca25753b9c97e886.gz
* Put group names into quotes when printed to screen.schnetter2004-02-09
| | | | darcs-hash:20040209114316-07bb3-d5041d29971f88aa96a099326ca31ee0af8f4115.gz
* Change spelling.schnetter2004-02-09
| | | | darcs-hash:20040209104203-07bb3-7e86c59c4d39515b9c6ccaaa3b77bbba1f03330b.gz
* Do not accidentally forbid output for grid variables with only oneschnetter2004-02-07
| | | | | | | Do not accidentally forbid output for grid variables with only one grid point. darcs-hash:20040207152111-07bb3-60db780519b86a7d6321f1a233d9b071fe45646d.gz
* Initialise the convergence factor.schnetter2004-02-05
| | | | | | | Initialise the convergence factor. Use the cGroup tagstable entry. darcs-hash:20040205235936-07bb3-b7806ebd48e57aa5f3e7085727461cbfe0ce7276.gz
* Shorten keyword option.schnetter2004-02-05
| | | | darcs-hash:20040205151120-07bb3-cfe9235ac5b627bd71559cd5c2db694b4cc895fe.gz
* Introduce a keyword parameter that selects between fixed and variableschnetter2004-02-05
| | | | | | | Introduce a keyword parameter that selects between fixed and variable progressive mesh refinement. darcs-hash:20040205085315-07bb3-08164f32105ec814e53691afac7d88ea8e2099fc.gz
* More error checking to make sure that global mode routines are indeedschnetter2004-02-03
| | | | | | | | | | | More error checking to make sure that global mode routines are indeed called. Make the "verbose" output a bit more verbose. Move Regridding to before advancing the time. darcs-hash:20040203154807-07bb3-c946590cd7c83081fe506ed0921175d3e57147dd.gz
* Do not select a prolongation operator for grid scalars and gridschnetter2004-02-03
| | | | | | | Do not select a prolongation operator for grid scalars and grid arrays. darcs-hash:20040203154613-07bb3-ede6b0469cf4513d4316b6a050879191b03ca936.gz
* Slight wording change in info output.schnetter2004-02-03
| | | | darcs-hash:20040203154501-07bb3-067b7d7d448f1f23b6bb9751aef0aaf8399f6a55.gz
* Do call boundary prolongation independent of the prolongationschnetter2004-02-03
| | | | | | | Do call boundary prolongation independent of the prolongation operator, because CarpetLib now catches the case "operator == none". darcs-hash:20040203133434-07bb3-7f2db8305c5fd443cbc3b1cdf47a62d37a0cdefa.gz
* Do not restrict nor prolongate when the prolongation operator has beenschnetter2004-02-03
| | | | | | | Do not restrict nor prolongate when the prolongation operator has been set to "none". darcs-hash:20040203133322-07bb3-7fd7e5fc85efe04ecd39864cba0385a4570bf46a.gz
* Some sanity checking to fend off evil C++ compilersschnetter2004-02-03
| | | | darcs-hash:20040203133302-07bb3-6468aba36ff7aceb42c9d34f1d859ce44cd9e860.gz
* Fix definition of mymin and mymax for complex values.schnetter2004-02-01
| | | | darcs-hash:20040201135456-07bb3-9c7b217df919037f256454b49062f7231a68fd00.gz
* Fix stupid bugschnetter2004-01-30
| | | | darcs-hash:20040130095246-07bb3-3418eda831f610244ed8477a2e77e4ca3595cc20.gz
* Call global mode routines at each time and for the coarsest possibleschnetter2004-01-30
| | | | | | | | Call global mode routines at each time and for the coarsest possible grid there, instead of (wrongly) at each time when the level 0 grid is active there. darcs-hash:20040130094052-07bb3-552e064e7f69fb827479e404c739f62cc7570887.gz
* An additional coord_delta appeared in the local grid geometry setup. With ↵hawke2004-01-29
| | | | | | | | | this removed it finds horizons and so on again. An additional coord_delta appeared in the local grid geometry setup. With this removed it finds horizons and so on again. However, for me it seems to be an order of magnitude slower than before?! darcs-hash:20040129113915-58737-8c9c38eabdb0523434e413305dfcb6378e4f4c2b.gz
* Fixed some namespace problems by explicitly specifying namespace stdswhite2004-01-27
| | | | | | | | | | | | | | | | | | Fixed some namespace problems by explicitly specifying namespace std for complex etc (is this the best solution here) Also, got rid of some very questionable zero-argument templates mymax and mymin that were supposed to be instantiated for complex type. Beyond being unclear how the compiler should choose to instantiate such templates non-explicitly, the mathematics is very questionable. But no matter, it seems they weren't used anyway. Note: I question the code in the float versions of mymax and mymin. It might work, but I suspect the coder didn't understand something... And I question the wisdom of calling these weird complex-valued function mymin and mymax. darcs-hash:20040127154340-5d157-8bea668a378ab69adc425ed9c708e222179f82c1.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
* 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-08 04:14:32 +0000