aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetLib/src/gdata.hh
Commit message (Collapse)AuthorAge
* CarpetLib: introduce accessor in gdata to query electric_fenceRoland Haas2013-10-31
| | | | | this function effectively exposes the CarpetLib parameter to Carpet without explicit parameter sharing
* CarpetLib: add some code for electric fence like functionalityRoland Haas2013-09-27
|
* CarpetLib: Add comment in gdata classErik Schnetter2013-08-04
|
* CarpetLib: Store registered gh, dh, th, gf, data etc. via sets, not via listsErik Schnetter2013-08-04
| | | | | | Use set instead of list to store registered gh, dh, th, gf, data. Do not store an iterator as handle for de-registering. This corrects some memory corruption with the Intel compiler; the previous code may have been illegal.
* CarpetLib: Use "process" instead of "processor" in comments and outputErik Schnetter2012-11-24
|
* Allow padding in transport operatorsErik Schnetter2012-11-22
| | | | | Rewrite padding infrastructure. Add padded array extents to transport operator APIs.
* CarpetLib: Remove/correct ATTRIBUTE_CONST and ATTRIBUTE_PURE declarationsErik Schnetter2012-09-11
|
* Re-organise some refluxing internalsErik Schnetter2012-09-11
| | | | | | | | | | | Re-organise some of the internal details of refluxing. Refluxing requires restricting fluxes from fine to coarse grids. Previously, Carpet would internally store bounding boxes that were offset by 1/2 grid point, and adjust (correct) these boxes in various places in a rather ad-hoc manner. This is now cleaned up. Remove the (unused) "accumulate" operator, and the (unused) "accumulate" prolongation type.
* CarpetLib: add empty constructor to struct slab to make gcc happyRoland Haas2012-09-11
|
* CarpetLib: Extend operator API to allow hyperslabbingErik Schnetter2012-09-11
| | | | | | | Extend API of grid operators and of gdata::copy_from, gdata::transfer_from to allow hyperslabbing Implement hyperslabbing in copy_3d.cc
* CarpetLib: Introduce gdata::allocated_memory_shapeErik Schnetter2011-12-14
| | | | | | Introduce a function gdata::allocated_memory_shape that determines from a given grid size how many grid points should be allocated. This allows adding padding to grid sizes, e.g. to improve cache performance. Use this function when checking grid sizes.
* CarpetLib: Rearrange order of #include statementsErik Schnetter2011-12-14
|
* CarpetLib: Correct ggf::new_typed_dataErik Schnetter2011-12-14
|
* CarpetLib: declare local variable as constErik Schnetter2011-12-14
| | | | Ignore-this: f1efc776fad3834bdac64279ef8c38c9
* CarpetLib: Implement padding for grid variablesErik Schnetter2011-12-14
| | | | | | | | | Ignore-this: 1a389f0dd3f40a0c0edb3fdabd6e7d40 Padding grid variables means that e.g. a component of size 32x32x32 is allocated as 33x33x33 instead, but only 32x32x32 of this storage is used. This can improve cache performance considerably. This requires corresponding changes to the cGH entries.
* Import CarpetErik Schnetter2011-12-14
| | | | Ignore-this: 309b4dd613f4af2b84aa5d6743fdb6b3
* CarpetLib: Make mem<T> saferErik Schnetter2008-02-19
| | | | | | | | | | When passing in a pointer to mem<T>, also pass the size of the pointed-to memory region, so that mem<T> can check that there is enough space. Keep track of the number of allocated bytes in mempools. darcs-hash:20080219044528-dae7b-107edc6f696a35aad32ef6e58129b3281d00eb56.gz
* CarpetLib: Add function ggf::fill, which fills all time levelsErik Schnetter2008-01-14
| | | | | | | | | Add a new function ggf::fill, which fills all time levels from the current time level. Remove the function ggf::copy, which was unused. darcs-hash:20080114151255-dae7b-469b57d0393ba0ed22590c2eabe862e02702fde7.gz
* CarpetLib: Add member gdata::allocsizeErik Schnetter2007-05-26
| | | | | | | Add a member gdata::allocsize which calculates the amount of memory that a data object requires. darcs-hash:20070526170719-dae7b-622c8d34ced2a5066fbaaf56e9d9b39efeea9022.gz
* CarpetLib: Various changes to class dataErik Schnetter2007-04-19
| | | | | | | | | Remove outdated communication mechanisms from data<T> and gdata classes. Redesign communication mechanism in data<T> and gdata classes. darcs-hash:20070419014301-dae7b-ece41543dcb1e23ab81f103b04a75880ac68807f.gz
* CarpetLib: Add _storage field to gdata classErik Schnetter2007-02-04
| | | | | | | Add a field "void * _storage" to the gdata class to speed up its storage() member function. darcs-hash:20070204175843-dae7b-644ade82e1432b82384baf64f55aeb2482032fdf.gz
* CarpetLib: Reorganise prolongation and restriction operatorsErik Schnetter2007-01-12
| | | | | | | | | | | | | | | | | | | | | | | | Reorganise prolongation and restriction operators. This is a major implementation change. Most operators are now written as C++ templates instead of as Fortran 77 code. This simplifies the code, since C++ routines can be called more easily, and they also have access to CarpetLib's high-level data structures. Previously, the operators combined temporal and spatial interpolation. Now, time interpolation and space interpolation are handled separately. This may be less efficient, but simplifies the code significantly, since there are now N+M instead of N*M routines, for N time interpolation and M space interpolation methods. Remove the minmod prolongation operator, which was previously disabled. Add support for cell centering, using a method described by Simon Hern, and suggested for Carpet by Ian Hawke. darcs-hash:20070112205812-dae7b-5329795aa698e7bbc3671b1504134885dd830238.gz
* Carpet*: generalise the comm_state class for collective buffer communicationsThomas Radke2005-08-15
| | | | | | | | | CarpetLib's comm_state class (actually, it's still just a struct) has been extended to handle collective buffer communications for all possible C datatypes at the same time. This makes it unnecessary for the higher-level communication routines to loop over each individual datatype separately. darcs-hash:20050815150023-776a0-dddc1aca7ccaebae872f9f451b2c3595cd951fed.gz
* CarpetLib: clean-up of communication scheme for individual sends/recvs on ↵Thomas Radke2005-05-26
| | | | | | | | | | | | | | | | | | | | single components The default communication scheme in Carpet (which does an individual send/recv operation for each component) comes with two parameters for fine tuning: CarpetLib::use_lightweight_buffers CarpetLib::combine_recv_send the defaults of which are set to use a well-tested but also slower communication pattern (as turned out during benchmark runs). This patch cleans up the implementation of this communication scheme so that the fastest communication pattern (combined posting of send/recv; use of lightweight buffers) is now always used. The above parameters therefore became obsolete and shouldn't be used anymore in parfiles. darcs-hash:20050526114253-776a0-780933a1539a260d74da8b92522fa2f48c714964.gz
* CarpetLib: remove bool gdata::owns_storage() which isn't needed anymore in ↵Thomas Radke2005-05-26
| | | | | | this class darcs-hash:20050526113252-776a0-58480600178d1c7beae50aa34808564a443f3c92.gz
* CarpetLib: Forbid copying data<T> objectsErik Schnetter2005-04-11
| | | | darcs-hash:20050411183135-891bb-5d2ced682685fb55a00da1864560e54bd113f765.gz
* CarpetLib: Remove unused function transfer_fromErik Schnetter2005-04-11
| | | | darcs-hash:20050411182954-891bb-6f141054635439136d978f98f528af1204919199.gz
* CarpetLib, Carpet: implement and use collective communication buffersThomas Radke2005-03-30
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Collective buffers are used to gather all components' data on a processor before it gets send off to other processors in one go. This minimizes the number of outstanding MPI communications down to O(N-1) and thus improves overall efficiency as benchmarks show. Each processor allocates a pair of single send/recv buffers to communicate with all other processors. For this the class (actually, the struct) comm_state was extended by 3 more states: state_get_buffer_sizes: accumulates the sizes for the send/recv buffers state_fill_send_buffers: gathers all the data into the send buffers state_empty_recv_buffers: copies the data from the recv buffer back into the processor's components Send/recv buffers are exchanged during state_fill_send_buffers and state_empty_recv_buffers. The constructor for a comm_state struct now takes an argument <datatype> which denotes the CCTK datatype to use for the attached collective buffers. If a negative value is passed here then it falls back to using the old send/recv/wait communication scheme. The datatype argument has a default value of -1 to maintain backwards compatibility to existing code (which therefore will keep using the old scheme). The new communication scheme is chosen by setting the parameter CarpetLib::use_collective_communication_buffers to "yes". It defaults to "no" meaning that the old send/recv/wait scheme is still used. So far all the comm_state objects in the higher-level routines in thorn Carpet (restriction/prolongation, regridding, synchronization) have been enabled to use collective communication buffers. Other thorns (CarpetInterp, CarpetIO*, CarpetSlab) will follow in separate commits. darcs-hash:20050330152811-776a0-51f426887fea099d1a67b42bd79e4f786979ba91.gz
* CarpetLib: Restructure lightweight communication buffersErik Schnetter2005-01-03
| | | | | | | Restructure the lightweight communication buffers. Use lightweight communication buffers for interpolation as well. darcs-hash:20050103200712-891bb-7e42816d3b8d667916084e3f32527c8f35327d7f.gz
* CarpetLib: Add lightweight communication buffers (untested)Erik Schnetter2005-01-02
| | | | | | | | | | | Lightweight communication buffers use essentially only a vector<T> instead of a data<T> to transfer data between processors. This should reduce the computational overhead. Set the parameter "use_lightweight_buffers" to use this feature. This feature is completely untested. darcs-hash:20050102173524-891bb-6a3999cbd63e367c8520c175c8078374d294eaa8.gz
* CarpetLib: Move class commstate into its own fileErik Schnetter2005-01-01
| | | | darcs-hash:20050101193846-891bb-7bb505d29a25b04c0d23e792eea7ff404d1f4200.gz
* global: Turn CarpetLib templates into classesErik Schnetter2005-01-01
| | | | | | | | | | | | | | | | | | | 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
* CarpetLib: Move class timestat into its own fileErik Schnetter2005-01-01
| | | | darcs-hash:20050101171429-891bb-130630de8631b8f9bbe494e135662ffb089ecca0.gz
* CarpetLib: Change checking of processor number in data classErik Schnetter2004-12-30
| | | | | | | | | | | | | Remove functions gdata::lives_on_this_processor() and gdata::this_processor_is(proc). Introduce dist::rank() and dist::size() instead. Re-introduce assert statement in data::proc(). Move declaration and definition of assignment operator near the constructor and destructor. darcs-hash:20041230191026-891bb-90eeb1be4c04753c165e13e7c1e65f06847180ca.gz
* global: Remove $Header$ keywords from Carpet arrangementErik Schnetter2005-01-01
| | | | darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
* CarpetLib: Collect more timing statisticsErik Schnetter2004-12-30
| | | | | | | | | | | | | Collect more timing statistics in the data class. Print these statistics to stdout when the Cactus parameter print_timestats is set. Create a timer class "timestat". This is a timer that can be started, stopped, and it prints the total time as well as some statistics. For memory allocation statistics, count the number of objects as well as the number of bytes. darcs-hash:20041230212136-891bb-c14edfa7d539ae9b135eee76afadaad51fd0b098.gz
* CarpetLib: Use new style instead of old style C++ headersErik Schnetter2004-12-25
| | | | | | That is, say e.g. <cmath> instead of <math.h>. darcs-hash:20041225182902-891bb-2a4502d6a76ed830ab978030fc4e4765b183509f.gz
* CarpetLib: Introduce parameter to post MPI_Irecv and MPI_Isend at the same timeErik Schnetter2004-12-08
| | | | | | | | | Introduce parameter a parameter to post MPI_Irecv and MPI_Isend at the same time. Use two queues instead of one vector to store the MPI_Requests. darcs-hash:20041208222541-891bb-c7c8994a0c41b6cfb37f6dc023bc1172238f3619.gz
* CarpetLib/src/*gf.*,*data.* const public members, canonicaliseswhite2004-12-09
| | | | | | | | | | | | | | | | | | The g* classes had exposed public member variables which were never meant to be altered from outside (I imagine altering most of them would have resulted in very bad behaviour); I made most of these const. This resulted in some complaints by the compiler about the canonical constructor and assignment ops, so I provided private dummies. This patch contains the following changes: M ./Carpet/CarpetLib/src/gdata.cc +6 M ./Carpet/CarpetLib/src/gdata.hh -2 +2 M ./Carpet/CarpetLib/src/gf.hh +5 M ./Carpet/CarpetLib/src/ggf.hh -9 +15 darcs-hash:20041209172710-32473-46d2a93e11f21f305467b375fe15910d47b8e4de.gz
* CarpetLib/src/*data* CarpetIOASCII/src/ioascii.cc MPI_Comm_rank wrappersswhite2004-12-08
| | | | | | | | | | | | | | | | | | | | | | | Removed evil assert from otherwise harmless public accessor gdata::proc(). Very systematically checked that the assert is still done everywhere this function was previously called. New functions bool this_processor_is (int procno); bool lives_on_this_processor (); wrapper messy repeated MPI_Comm_rank code. Had to modify ioascii.cc due to missing assert as explained above, but found no other instances of calls to proc(). This patch contains the following changes: M ./Carpet/CarpetIOASCII/src/ioascii.cc +1 M ./Carpet/CarpetLib/src/data.cc -41 +21 M ./Carpet/CarpetLib/src/gdata.cc -6 +19 M ./Carpet/CarpetLib/src/gdata.hh -1 +3 darcs-hash:20041208171657-32473-11da4a7ec810f97b2c82c322119fa69a9ffb4640.gz
* CarpetLib: Implement using MPI_Waitall for communicatingErik Schnetter2004-11-24
| | | | | | | | Add a parameter CarpetLib::use_waitall that switches from using a series of MPI_Wait statements to using a single MPI_Waitall statement. This might improve performance on many processors. darcs-hash:20041124235118-891bb-034efea054db236a187022b1858e4574da867fa3.gz
* CarpetLib: more consistency checks during communicationErik Schnetter2004-10-20
| | | | darcs-hash:20041020094227-891bb-aac66b506731e35384752ff0fa0df4e91f689a71.gz
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* change indentationschnetter2004-04-19
| | | | darcs-hash:20040419055635-07bb3-6014661f8995a59c900cebf24c03f52b6eed2d21.gz
* Move MPI request and tag handling from data to gdataschnetter2004-03-23
| | | | darcs-hash:20040323111429-07bb3-08a8c6b8d729eb0efc5b757b61401b14712d790c.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
* Follow C++ standard.schnetter2003-11-21
| | | | darcs-hash:20031121125546-07bb3-60338023244fbd0fab9d8663adff4cdf78e31141.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
* 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
generated by cgit v1.2.3 (git 2.39.1) at 2024-06-10 08:25:57 +0000