| Commit message (Collapse) | Author | Age |
|
|
|
| |
darcs-hash:20040207142140-07bb3-e5e38c10925f303f4ca0d67f4c52142e78403d8d.gz
|
|
|
|
|
|
|
| |
Initialise the convergence factor.
Use the cGroup tagstable entry.
darcs-hash:20040205235936-07bb3-b7806ebd48e57aa5f3e7085727461cbfe0ce7276.gz
|
|
|
|
| |
darcs-hash:20040205151120-07bb3-cfe9235ac5b627bd71559cd5c2db694b4cc895fe.gz
|
|
|
|
|
|
|
| |
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 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 grid
arrays.
darcs-hash:20040203154613-07bb3-ede6b0469cf4513d4316b6a050879191b03ca936.gz
|
|
|
|
| |
darcs-hash:20040203154501-07bb3-067b7d7d448f1f23b6bb9751aef0aaf8399f6a55.gz
|
|
|
|
|
|
|
| |
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 been
set to "none".
darcs-hash:20040203133322-07bb3-7fd7e5fc85efe04ecd39864cba0385a4570bf46a.gz
|
|
|
|
| |
darcs-hash:20040203133302-07bb3-6468aba36ff7aceb42c9d34f1d859ce44cd9e860.gz
|
|
|
|
| |
darcs-hash:20040201135456-07bb3-9c7b217df919037f256454b49062f7231a68fd00.gz
|
|
|
|
| |
darcs-hash:20040130095246-07bb3-3418eda831f610244ed8477a2e77e4ca3595cc20.gz
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
|
|
| |
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 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:
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:
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
|
|
|
|
|
|
| |
a lot more memory and disk space than it used to.
darcs-hash:20040123074040-58737-eef0cd6ebc85e62611a2f15ab5eb0b927a4bc98a.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
darcs-hash:20040121153204-07bb3-2737dc30c6c1b867713208ee5dfe49085ba1b625.gz
|
|
|
|
| |
darcs-hash:20040121132535-07bb3-8221de1d3146f586eb13b0fe3acb12836b06efd5.gz
|
|
|
|
|
|
|
| |
Check the stencil size before calling the minmod prolongation
operators.
darcs-hash:20040121132509-07bb3-7eb951d88ae3b66a3ad2338fecc632c8c77cb312.gz
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
darcs-hash:20040118102930-07bb3-1ed8463535f2d2e32616e985598c3d588a1f70a3.gz
|
|
|
|
| |
darcs-hash:20040117084347-07bb3-96e4d27e49955470253973d7bd8a7642fa8d5228.gz
|
|
|
|
|
|
|
| |
While restricting, don't assume that there are buffer zones at
interprocessor boundaries.
darcs-hash:20040116163644-07bb3-7eef930c1b8bb461228d194dcb9ecddd5534d9ad.gz
|
|
|
|
|
|
| |
call. Note that this is always the maximum number of timelevels as Carpet always switches on all storage.
darcs-hash:20040116095704-58737-0b01f91957109646209f596a26c1dd4a2031389d.gz
|
|
|
|
|
|
| |
maxreflevels
darcs-hash:20040115084558-19929-5187725fb081ea4da02a6805f2cc65ecbfdf4244.gz
|
|
|
|
| |
darcs-hash:20040114081408-07bb3-65a7482594808f8bca77df73d41ff1763177c541.gz
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
| |
darcs-hash:20040113142658-19929-bdda00a3cb2bc7e14c7111108e8e471892ca99cb.gz
|
|
|
|
| |
darcs-hash:20040113134652-19929-a9114fc6d8c5ebc0101c4991df99032093c3f35a.gz
|
|
|
|
|
|
|
| |
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.
- 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
|
|
|
|
|
|
|
|
| |
Mostly code beautification. Incorporated Erik's new WriteAttribute
functions and changed the way CarpetIOFlexIO writes dataset attributes.
All the functionality provided in CarpetIOFlexIO is now provided.
darcs-hash:20040112095031-19929-f76eb19112c1ff26b232ab05ba03b0b1beb00687.gz
|
|
|
|
|
|
|
|
| |
Mostly code beautification. Incorporated Erik's new WriteAttribute
functions and changed the way CarpetIOFlexIO writes dataset attributes.
All the functionality provided in CarpetIOFlexIO is now provided.
darcs-hash:20040112085031-19929-ce03ec0298a148bc14855d0db33cfa50ce59e681.gz
|
|
|
|
| |
darcs-hash:20040111112039-07bb3-7ffb487a26951fb48a0b090f5aaa6d160a4620f0.gz
|
|
|
|
| |
darcs-hash:20040111102039-07bb3-aa2db5eb23fb28d9175c9e2718658a4b6f868a86.gz
|
|
|
|
| |
darcs-hash:20040109144346-19929-67096355895264e0861c3792c7bd3893ee932801.gz
|
|
|
|
| |
darcs-hash:20040109122120-07bb3-ec22ca801d35b1de50585080dd2ab0e8d31da4d5.gz
|
|
|
|
|
|
|
|
|
| |
Okay. Here we go. Recovery goes through without any crashes on a very simple
test case. No test for consistency yet.
The code is still a huge construction site...
darcs-hash:20040108184333-19929-b0d1df37d71a59413ba5893926fef4e5862ea7a7.gz
|
|
|
|
| |
darcs-hash:20040107153049-19929-d817dfcb9fcaae5f794d4648530bc781a2fa0982.gz
|
|
|
|
| |
darcs-hash:20040107121404-19929-279507a275e8bcb9b2ae29c5eb0e87d86b2c34cd.gz
|
|
|
|
|
|
|
| |
Getting there... Skeleton of reovery routine for data is in place.
Need to modify scheduling for IOUtil_RecoverGH: OPTIONS: global
darcs-hash:20040107115756-19929-6e1501955bbe1a8811de898e2803d3e22875265d.gz
|
|
|
|
| |
darcs-hash:20040106075639-19929-06e2883236b3a6710cabb2229a18e62219bd5452.gz
|
|
|
|
|
|
| |
parameter distribution
darcs-hash:20040105215416-19929-b5fa95175f2996b0b9effd25433de64c4cb914af.gz
|
|
|
|
| |
darcs-hash:20031219150155-07bb3-41d0f7f3df186e02d5a3a98ad823906252564440.gz
|
|
|
|
| |
darcs-hash:20031219104155-07bb3-c74a7c6ed129c05fe6687c4ac627874a93cda3d5.gz
|
|
|
|
| |
darcs-hash:20031219103309-07bb3-023267a17f4d7bcbd2399cedefc787cc228a9539.gz
|
|
|
|
| |
darcs-hash:20031219094352-07bb3-20d9e41ca0bb17d441f83f83edcab34739ef1fde.gz
|