| Commit message (Collapse) | Author | Age |
|
|
|
|
| |
Remove CarpetLib's "equals" method.
Implement equal_to, less, greater, greater_equal, less_equal, not_equal for vect and bbox classes.
|
| |
|
| |
|
| |
|
| |
|
|
|
|
| |
Change count() from counting all elements, which is already done by size(), to counting true (non-zero) elements, which is what Fortran's count() does as well.
|
|
|
|
| |
vect::map is already called vect::vmap.
|
|
|
|
| |
and a scalar
|
|
|
|
| |
Replace calls to abs() with floating arguments by calls to fabs(), to ensure the results are not accidentally truncated to integer.
|
| |
|
| |
|
|
|
|
| |
Remove Carpet's mechanisms to deal with how to call isnan. Instead, expect isnan in std::, and rely on Cactus to correct things if this is not the case.
|
| |
|
| |
|
| |
|
|
|
|
|
|
|
| |
Provide idiv/imod functions for signed integer division that always
round towards negative infinity. Also provide vector version of these.
Use these in the bbox output routines to ensure proper rounding.
|
| |
|
| |
|
|
|
|
|
| |
Capture the system's nan functions into the namespace CarpetLib::good,
then offer #defines to access these.
|
| |
|
| |
|
| |
|
|
|
|
| |
Ignore-this: 309b4dd613f4af2b84aa5d6743fdb6b3
|
| |
|
|
|
|
|
| |
Introduce a new function vect<T,D>::index(lsh,ind), which calculates the
linear grid point index given the array size and an index triplet.
|
|
|
|
|
|
|
| |
Add an overloaded function memoryof which measures the size in bytes
of CarpetLib's data structures.
darcs-hash:20080202154302-dae7b-de41c79cb04617327695fae27928c05c5d431ee8.gz
|
|
|
|
|
|
|
|
|
| |
Define vect<T,D> functions with macros.
Define many vect<T,D> functions also for vect<vect<T,D>,E> and for
combinations with scalars.
darcs-hash:20070419014657-dae7b-b3c62cc55963872a39cd48439c2057918578eda3.gz
|
|
|
|
|
|
|
|
| |
Add some new functions to make the vect class more convenient.
Remove some unused functions.
Instantiate some more templates.
darcs-hash:20070112205505-dae7b-12a0ea160b9740cdbe60d89d9e9100fc06a1faa7.gz
|
|
|
|
|
|
|
|
|
| |
Add __attribute__((const)) and __attribute__((pure)) to vect<T>
functions where appropriate. These attributes are not active, because
it is not clear that this really improves performance, and errors
could lead to wrong code.
darcs-hash:20060731152218-dae7b-20c523f60890b2ad611157764e0b993bfb185c8d.gz
|
|
|
|
|
|
|
|
|
| |
Move the arithmetic operations (+-*/ etc.) out of the vect class.
Instead of member functions, they are now standalone functions. This
makes it possible to write "int * vect<int>" with automatic
conversion.
darcs-hash:20060605220401-dae7b-eee6ed411e3d0e5850dd8e6bd85d432eed4752e8.gz
|
|
|
|
| |
darcs-hash:20050501171344-891bb-d608668bef718cc9f4c6244b721b3d43d7c33ae0.gz
|
|
|
|
| |
darcs-hash:20050209223915-891bb-4c4001d07890086e95de8d5a91deffc66c32e469.gz
|
|
|
|
| |
darcs-hash:20050209222027-891bb-e8501dfa40575303af1338c1d2d4528d08ea273c.gz
|
|
|
|
| |
darcs-hash:20050127104948-891bb-214d1033924b4db7a215a624eb337f8a619ac82d.gz
|
|
|
|
| |
darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
|
|
|
|
| |
darcs-hash:20041225183147-891bb-84276981b1deab37af4c3b541f993e62d63200e6.gz
|
|
|
|
|
|
| |
Replace all CVS header tags with the standard "$Header:$".
darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
|
|
|
|
|
|
|
| |
Explain the harmless compiler warnings about array indices out of
range.
darcs-hash:20040814054125-07bb3-232abf6b8906138d6b7425e2ec05feca8eac2368.gz
|
|
|
|
| |
darcs-hash:20040527102609-07bb3-8a80ea8a94dc7739da7e482ae925beb4a01585c3.gz
|
|
|
|
| |
darcs-hash:20040504201207-07bb3-b8a9db665cc225198475179abc2d2e6bcd7f2b2f.gz
|
|
|
|
| |
darcs-hash:20040418112556-07bb3-87cca0e1b0860665e56b2a729a06baba55fa5580.gz
|
|
|
|
| |
darcs-hash:20040311110413-07bb3-71275d869ab76640be38effb5f04522f0dcbef6c.gz
|
|
|
|
| |
darcs-hash:20040227092523-07bb3-e4cc3ad72fab0773cc5f42d1af3a27cf61d01545.gz
|
|
|
|
| |
darcs-hash:20040218141045-07bb3-0e60f56d44aa61e0e3190d5d27b9a0bfab3f2d38.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
|
|
|
|
| |
darcs-hash:20031113150358-07bb3-427886aab7c861fe53bf7fef96ec2a31e1d155f5.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
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
|
|
|
|
|
|
|
| |
Fix bug.
Remove duplicate functions.
darcs-hash:20030828190727-07bb3-ec3abefc3438a4822d11d3ff43587a0968e9a2a9.gz
|
|
|
|
| |
darcs-hash:20030815073254-07bb3-d964703599463eb0657c0d6c6a00d3450fc3c31d.gz
|
|
|
|
|
|
|
| |
Add doxygen documentation.
Doxygen is great. Forget grdoc. Forget you even ever heard its name.
darcs-hash:20030803142621-07bb3-e18d0b14080698cd6ddbe94deb678ce0a29fd339.gz
|