| Commit message (Collapse) | Author | Age |
|
|
|
|
|
|
|
|
|
|
|
|
| |
Make the function initialise and reduce in the reduction operators
handle the count argument as well.
Introduce "combine" functions into the reduction operators. These
combine results from several independent reduction operations.
Use the above new functions to combine reductions, correcting an
error.
darcs-hash:20070828044751-dae7b-02070b6aa39d25fb91481ab6a9b45cf37b10fe2f.gz
|
|
|
|
| |
darcs-hash:20070827201932-dae7b-2129f2816c6dcf448bf2706cf5d49cfc83e8d540.gz
|
|
|
|
| |
darcs-hash:20070821190958-dae7b-1ed180793dd5d6482f48aa0bbf13bd5523a4b40b.gz
|
|
|
|
| |
darcs-hash:20070821190829-dae7b-9efe9944356b64e8c171e6f3917f03dee4ba1deb.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Implement a new set of reduction operators for Carpet. These have the
same definitions as the old operators, except that they give the same
weight to all grid points. The existing operators give a larger
weight to coarser grid points, so that physical reduction operators
work correctly, e.g. for calculating the total mass.
The new operators have the same name as the old operators plus a
prefix "i" (the lower case letter "eye"). For example, the new L2
norm is called "inorm2", while the old norm continues to be called
"norm2".
The physical (old) operators converge to a result which is independent
of the grid structure. The pointwise (new) operators still converge,
but to a value which depends on the grid structure.
In order to make it possible to compare situations with and without
symmetries, boundary points still carry a weight of 1/2 or whatever is
appropriate. That means that there exists a simple relationship
between the pointwise norms with and without a symmetry boundary.
darcs-hash:20070526195016-dae7b-6a7e2a39bebd9c0a6377f908a4d7f446f8288191.gz
|
|
|
|
|
|
|
| |
Set up the weight function after recovering, since the weight function
is not checkpointed.
darcs-hash:20070516222541-dae7b-a62fe7ade7eb4ef4a381c5710feee91321605658.gz
|
|
|
|
| |
darcs-hash:20070421040012-dae7b-17a08c5094234b1b07bf2ad275d92ab2c11a4fc1.gz
|
|
|
|
| |
darcs-hash:20070419021335-dae7b-fc08da29dd6e24cb828bac615392dd9635241a8c.gz
|
|
|
|
| |
darcs-hash:20070419021320-dae7b-9bdc102a1d6dd1a7f998dc6b8cb446572abb2b7c.gz
|
|
|
|
| |
darcs-hash:20070228023136-dae7b-86633994d6887c2dd91d7c57d56bc1db620ed2da.gz
|
|
|
|
|
|
|
|
| |
Disable the special case that is supposed to handle empty domains for
the reduction mask. The implementation of the test is flawed, not the
idea; the test should be re-implemented correctly.
darcs-hash:20070112224810-dae7b-b3000e90e618f81b9ec188dea2a3223c01889b10.gz
|
|
|
|
| |
darcs-hash:20070112224438-dae7b-38b62080071060c00a7bf2be053a584af739991b.gz
|
|
|
|
| |
darcs-hash:20070112224413-dae7b-a6a68d471d22ef1f0f1b2a1561ef4b302e017bc1.gz
|
|
|
|
| |
darcs-hash:20061206163437-dae7b-c2647071f4257c4a2cee820a470ff6b80c20e106.gz
|
|
|
|
| |
darcs-hash:20070112224318-dae7b-68b8b4b34cba9b14fff76fc6c7fdbb4032297430.gz
|
|
|
|
| |
darcs-hash:20061003234056-dae7b-c72aaf74712cba7006df16c3382a6ba54e0dad5a.gz
|
|
|
|
| |
darcs-hash:20060925220745-dae7b-eb145a85ef330dbff56b7d1c8844449dff85dd0b.gz
|
|
|
|
|
|
|
| |
Use the macros HAVE_CCTK_REAL* instead of checking whether the macros
CCTK_REAL* are defined.
darcs-hash:20060925220641-dae7b-008b16676bf9b39fd14d718ddeadb009709b6096.gz
|
|
|
|
|
|
|
| |
When storing the extent of the boundary, use the local variable names
bmin/bmax instead imin/imax.
darcs-hash:20060925220452-dae7b-7a12b1227486a66e0e4689554627c3f713b9d36b.gz
|
|
|
|
|
|
|
| |
Explicitly convert between int and CCTK_INT, and between double and
CCTK_REAL, where necessary.
darcs-hash:20060613171604-dae7b-6cdecab56f7600ef471621e808847d32e45baa75.gz
|
|
|
|
| |
darcs-hash:20060607225507-dae7b-31f5f3ab59b713a68644e8d574dc7ae9f9db55aa.gz
|
|
|
|
|
|
|
| |
Do not abort when some refinement levels of some groups have no
storage. Instead, return an error code to the called.
darcs-hash:20060413193236-dae7b-63053a270f4790c9f5608b4fe284a6ad4ab33604.gz
|
|
|
|
| |
darcs-hash:20060413201723-dae7b-59318803a1ad41a05efde70bc212cd9a37c260f8.gz
|
|
|
|
|
|
| |
Obtaining the weight function's index showed up high in a profile.
darcs-hash:20051119213646-dae7b-c0f463975c782b1a6854aa66a84169609efe2443.gz
|
|
|
|
|
|
|
| |
This norm calculates sum_i |x_i|^2, which is an ingredient for the L2
norm.
darcs-hash:20051119213505-dae7b-e59c1b2fdb5e24455331c20053cb248cca4e898b.gz
|
|
|
|
|
|
|
|
|
| |
Use static_cast<> instead of C style casts.
Do not use the STL's sqrt() function for integer variables because
e.g. Intel 7 does not support this.
darcs-hash:20051119213126-dae7b-5d33bbb6cfc786b980dd93332289528af6e5977f.gz
|
|
|
|
| |
darcs-hash:20051119213612-dae7b-2f590847df48bdaec3e05b248adcdab8f564d647.gz
|
|
|
|
|
|
| |
The weight is set up at basegrid so it doesn't need to be checkpointed/recovered.
darcs-hash:20050913163020-776a0-a171ae3372b9e38ad521a7db85470a8741088b46.gz
|
|
|
|
| |
darcs-hash:20050726110949-891bb-fe7d442ee41c585b7078fd6746840cd68de1d770.gz
|
|
|
|
| |
darcs-hash:20050704162500-891bb-589be0634ace55be5a3a5e0b633b86435c3f847b.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Empty simulation domains occur e.g. with Cartoon2D, where the total
volume of the domain is zero: All evolved grid points lie on the
boundary. Before, CarpetReduce did not handle this special case
correctly.
Instead of returning 0 for all reductions, set up the grid point
weights to perform a 2D reduction on the data. That is, the domain is
treated as if it was one grid point wide in the direction in which it
has the extent zero.
darcs-hash:20050603220813-891bb-83088cbd570d967aee95170915c54d905ea5bf71.gz
|
|
|
|
|
|
|
|
|
|
| |
Add support for varying refinement factors. The spatial refinement
factors can be different in different directions, can be different
from the time refinement factor, and can be different on each level.
(However, the underlying spatial transport operators do currently not
handle any factors except two.)
darcs-hash:20050501205010-891bb-8d3a74abaad55ee6c77ef18d51fca2a2b69740de.gz
|
|
|
|
|
|
| |
When there are zero reductions, allow the output pointers to be NULL.
darcs-hash:20050430154439-891bb-95ea930518b19378040601526504ceb6be039a97.gz
|
|
|
|
|
|
|
| |
When a grid function has only one time level, emit a level 1 warning
instead of aborting with an assertion failure.
darcs-hash:20050428151159-891bb-53d3330b499a68f0304e99f86aa176af2ca00ecf.gz
|
|
|
|
| |
darcs-hash:20050303102622-891bb-8d89ed158189c376fa1235ece02a5d8fef87fcae.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Change the way in which the grid hierarchy is stored. The new hierarchy is
map
mglevel
reflevel
component
timelevel
i.e., mglevel moved from the bottom to almost the top. This is
because mglevel used to be a true multigrid level, but is now meant to
be a convergence level.
Do not allocate all storage all the time. Allow storage to be
switched on an off per refinement level (and for a single mglevel,
which prompted the change above). Handle storage management with
CCTK_{In,De}creaseGroupStorage instead of
CCTK_{En,Dis}ableGroupStorage.
darcs-hash:20050201225827-891bb-eae3b6bd092ae8d6b5e49be84c6f09f0e882933e.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Turn most of the templates in CarpetLib, which used to have the form
template<int D> class XXX
into classes, i.e., into something like
class XXX
by setting D to the new global integer constant dim, which in turn is set to 3.
The templates gf and data, which used to be of the form
template<typename T, int D> class XXX
are now of the form
template<typename T> class XXX
The templates vect, bbox, and bboxset remain templates.
This change simplifies the code somewhat.
darcs-hash:20050101182234-891bb-c3063528841f0d078b12cc506309ea27d8ce730d.gz
|
|
|
|
| |
darcs-hash:20050127105110-891bb-078f5704db1a1af7004e6cbde2dc1b7577d3f669.gz
|
|
|
|
| |
darcs-hash:20050110165353-891bb-5918670062d3b1d2ac2f3381a00bd3dc8e0867e9.gz
|
|
|
|
| |
darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
|
|
|
|
| |
darcs-hash:20041225201857-891bb-ec4076f15ac19f9aef88b0bf2ae728bd9e31f7d0.gz
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
| |
Made all variables in CarpetLib:gh to be const or non-public, to avoid
confusion and insanity.
Toward this end, created four member accessors for gh:
const rexts & extents() const
const rbnds & outer_boundaries() const
const rprocs & processors() const
const vector<vector<ibbox> > & bases() const
This involved a few changes in several files throughout the repository:
M ./Carpet/Carpet/src/Initialise.cc -2 +2
M ./Carpet/Carpet/src/Recompose.cc -8 +8
M ./Carpet/Carpet/src/SetupGH.cc -2 +2
M ./Carpet/Carpet/src/modes.cc -3 +3
M ./Carpet/CarpetIOHDF5/src/Recover.cc -4 +4
M ./Carpet/CarpetInterp/src/interp.cc -5 +5
M ./Carpet/CarpetLib/src/dh.cc -3 +3
M ./Carpet/CarpetLib/src/gh.cc -33 +33
M ./Carpet/CarpetLib/src/gh.hh -10 +28
M ./Carpet/CarpetReduce/src/mask_carpet.cc -4 +4
M ./Carpet/CarpetRegrid/src/automatic.cc -2 +2
darcs-hash:20041210131459-32473-ff0835ff0e57f6693fbe1ec23b350d1b1d18e3bc.gz
|
|
|
|
|
|
|
|
| |
Made all Carpet thorns to include the "typecase" file from the Carpet
thorn using the "include header: .. in .." mechanism, getting rid of
explicit file paths.
darcs-hash:20041201161810-32473-75f4594a4d9c19a11dbe9d07741678a84bd6f694.gz
|
|
|
|
| |
darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
|
|
|
|
| |
darcs-hash:20041008214326-891bb-92f62ca9e16ba5acb5d4b79b1865057fb345711d.gz
|
|
|
|
| |
darcs-hash:20040922204131-891bb-e4070a4c2cee6417e961ca1aa283aeb98bc695d8.gz
|
|
|
|
|
|
| |
Do not set weight to zero on ghost zones. Ghost zones are ignored anyway.
darcs-hash:20040922204017-891bb-26f1608179b6415cebf161275751ea767cf8d8f8.gz
|
|
|
|
| |
darcs-hash:20040922203956-891bb-9d66f72ca19242dc7eda1208732a4acd6c5ee923.gz
|
|
|
|
| |
darcs-hash:20040922203920-891bb-59f30143f8d685db70da4626c5129e1908484eae.gz
|
|
|
|
|
|
|
| |
Use relative instead of absolute error when checking that the times
are aligned during reduction.
darcs-hash:20040921161026-891bb-361d8f344ec91f7d1392101267fa94114d1d18e5.gz
|