aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetInterp
Commit message (Collapse)AuthorAge
...
* CarpetInterp: Correct error in allocating data objectsErik Schnetter2005-04-11
| | | | | | | | | data<> objects cannot be implicitly copied. The standard template library containers copy objects arbitrarily. That means that one has to store pointers to data<> objects instead of the objects themselves, and has to allocate and free them manually. darcs-hash:20050411170907-891bb-406b9b6bb6b97df2f47c349f32d91398338439ae.gz
* CarpetInterp: Correct check whether a processor owns a componentErik Schnetter2005-03-07
| | | | darcs-hash:20050307170026-891bb-03754477692ad245563fda22ecdd4510da4549ab.gz
* CarpetInterp: Handle groups that have not all time levels allocatedErik Schnetter2005-03-03
| | | | darcs-hash:20050303102506-891bb-118df8e38d1ca7e981bb4f97b14d5a747c1b3928.gz
* global: Change the way in which the grid hierarchy is storedErik Schnetter2005-02-01
| | | | | | | | | | | | | | | | | | | | 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
* 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
* global: Remove $Header$ keywords from Carpet arrangementErik Schnetter2005-01-01
| | | | darcs-hash:20050101162121-891bb-ac9d070faecc19f91b4b57389d3507bfc6c6e5ee.gz
* CarpetLib:gh exposed data clamp-downswhite2004-12-10
| | | | | | | | | | | | | | | | | | | | | | | | | | | 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
* Do not use CVS $Date$ for latex documentation any moreErik Schnetter2004-11-18
| | | | darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
* CarpetInterp: return local_interpolator_statusErik Schnetter2004-10-20
| | | | | | | Return the local interpolators' statuses in a table entry local_interpolator_status. This is e.g. required by AHFinderDirect. darcs-hash:20041020092737-891bb-495add6874e22d32e23eacbc340493a1e4656da3.gz
* CarpetInterp: clean up time interpolationErik Schnetter2004-10-08
| | | | darcs-hash:20041008214254-891bb-222fb022121e1d3b6f18ca8176640d400039a3d4.gz
* CarpetInterp: Fix segfault when determining interpolation timesErik Schnetter2004-10-02
| | | | darcs-hash:20041002183708-891bb-96a5a3227becd0e40834fb9ab25ef5f9d75c1df5.gz
* Whitespace change in CarpetInterpErik Schnetter2004-09-22
| | | | darcs-hash:20040922203859-891bb-7e628f331129d9afeded4ea54732c37ca48cf26c.gz
* Implement time derivatives during interpolationErik Schnetter2004-09-21
| | | | | | | Add interpolation option key "time_interpolation_order" that allows taking time derivatives while interpolating. darcs-hash:20040921160917-891bb-236dc322c4bacdd2fa8d13b5c535541640e02537.gz
* Add true interpolation in timeErik Schnetter2004-09-18
| | | | | | | | | Add an interpolation option "interpolation_times" that specifies for each grid point at what time the value should be calculated. This allows for some interpolation in time when there are sufficient time levels. darcs-hash:20040918132542-891bb-887b348b0e510d7f17ac32b41cb78176079ff1fb.gz
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* Map output indices to input indices when interpolating in time.schnetter2004-09-13
| | | | darcs-hash:20040913123613-07bb3-060a47d1900cb82aa8260be38049e1e4b70f4e50.gz
* Add a new tags table entry to check: InterpNumTimelevels.hawke2004-08-20
| | | | | | | | | | | | | | | | | | | | | | | | Add a new tags table entry to check: InterpNumTimelevels. When this is set then the interpolator will take the given number of timelevels as the number to use when interpolating in time. So if you have grid functions with fewer than 3 timelevels mixed with grid functions with 3 timelevels you can use this tag, set it to the number of timelevels that you have available, and the interpolator will "do the right thing". Note that this means that GFs with only one timelevel are implicitly treated as being constant in time. The obvious examples, soon to be committed (I hope) are the static conformal factor and the coordinates (although why you'd want to interpolate the coordinates...). E.g., REAL confac TYPE = GF timelevels = 1 tags='tensortypealias="Scalar" Prolongation="None" InterpNumTimelevels=1' ... darcs-hash:20040820072138-58737-5ff68c487af14a1d72228aa6697cc9eea1ae3a89.gz
* Correct error in calculating the components' coordinate extents whenschnetter2004-08-02
| | | | | | | Correct error in calculating the components' coordinate extents when convergence_level is not zero. darcs-hash:20040802163927-07bb3-639b6b816f790ab12935ccb4052c7c8808b7a6ff.gz
* Ensure that the array pointers, if given, are not null.schnetter2004-08-02
| | | | | | | Ensure that the array pointers, if given, are not null. Allow zero interpolation points. darcs-hash:20040802093953-07bb3-2af122bfdc454e59d31b4a7f73cc44f59c9bb9d7.gz
* A more useful error message if time interpolation can't be done.hawke2004-07-02
| | | | darcs-hash:20040702131619-58737-b02419fc8b4f4514ab1bb31f341054785bef0a21.gz
* Correct a bug in the time interpolation. We probably neverschnetter2004-06-27
| | | | | | | | Correct a bug in the time interpolation. We probably never interpolated in time after all, and if so, then with a sign error. This correction is by analogy from CarpetReduce and not tested. darcs-hash:20040627191732-07bb3-326b9ecadca08013a8a463f5ce1a7883b0e49bac.gz
* Use "requires thorns" to ensure that the thorns are activated in theschnetter2004-06-08
| | | | | | | Use "requires thorns" to ensure that the thorns are activated in the correct order. darcs-hash:20040608190925-07bb3-70b0ac8ee8e570c0b7fff50b02395a288482257d.gz
* Remove restrict qualifier.schnetter2004-05-28
| | | | darcs-hash:20040528115004-07bb3-dfa6908bd400b1096f8406de3b7e4cdd1283ecfc.gz
* Add support for symmetry interpolation.schnetter2004-05-21
| | | | darcs-hash:20040521161223-07bb3-6cc7ddc3777345de379eb54daf2cf6d532925c23.gz
* Undo interface/capability update, because now Carpet conflicts withschnetter2004-05-04
| | | | | | | | Undo interface/capability update, because now Carpet conflicts with PUGH. Need to sort things out, maybe switch temporarily back to the "requires thorns" mechanism. darcs-hash:20040504205239-07bb3-771c6ec96a2b454d63c87a6dc7a1f44f62cd8fd7.gz
* Use the local instead of a registered coordinate system forschnetter2004-05-04
| | | | | | | | Use the local instead of a registered coordinate system for interpolation. This will enable patch-local multi-patch interpolation. darcs-hash:20040504201254-07bb3-2c00ac4fb8a00d211289814abc892908d4c78a42.gz
* Begin to clean up inheritance and requires relationshipsschnetter2004-05-04
| | | | darcs-hash:20040504200954-07bb3-cd9f17c981ce37247fa2122858b10b46425f678a.gz
* Use "requires" instead of "requires thorns".schnetter2004-04-18
| | | | darcs-hash:20040418120959-07bb3-05050a71e475171cac6b293de35118e4af8df5ba.gz
* When a grid point is not in any component, map it to an arbitraryschnetter2004-04-16
| | | | | | | When a grid point is not in any component, map it to an arbitrary component instead of aborting. darcs-hash:20040416185756-07bb3-9faf77dee8876c69e83780620702df6d2c660bec.gz
* Explicitly cast constness awayschnetter2004-03-31
| | | | darcs-hash:20040331143832-07bb3-7a1e263661a7034f7ad6d38b048ce79015d4f20d.gz
* Handle looping over refinement levels more elegantly.schnetter2004-03-31
| | | | darcs-hash:20040331105058-07bb3-0c0495c75cd248152699703e2e6d1520c9c80052.gz
* Change operator[] to at()schnetter2004-03-08
| | | | darcs-hash:20040308081244-07bb3-353fed99bc85bc4ba6c8aa6f6c9739a31831ffa3.gz
* Use the new vect::floor instead of vect::map(floor,)schnetter2004-02-27
| | | | darcs-hash:20040227152205-07bb3-33135e23f97ccb01f347189ca57b01003732863d.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
* 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
* 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
* Return the local interpolator's error code if something goes wrong.schnetter2003-11-13
| | | | darcs-hash:20031113162947-07bb3-267cc4954cf1498f2a1f058eb4c23904976e92e9.gz
* Use vector<> instead of dynamic [] arrays.schnetter2003-11-07
| | | | darcs-hash:20031107122106-07bb3-a21c30bb42714054caed0490f79121cfef385697.gz
* Return local interpolation errors to the caller instead of aborting.schnetter2003-11-06
| | | | darcs-hash:20031106133319-07bb3-f59f24cdc8a5f5b170693e56b1b8f2a14b1adf0e.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
* 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:20031105141838-07bb3-d348d033307d82905893e241be485c4db2ce8ce3.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
* Carpet infrastructure thorns now require the thorns they depend upontradke2003-09-04
| | | | | | | | (CarpetLib, Carpet, FlexIO) to be compiled in and activated via the "Requires thorns: ..." attribute in the configuration.ccl files. No need to inherit from those thorns/implementations anymore. darcs-hash:20030904142321-1d9bf-309be01f0950d97b8cb45d041d44edce01aff2a5.gz
* Print an error message instead of an assertion failure whenschnetter2003-08-15
| | | | | | | Print an error message instead of an assertion failure when CarpetInterp is asked to interpolate in time. darcs-hash:20030815094225-07bb3-ca72656b5a70fca45e186e9fd299fffbe5a08990.gz
* Do not include header files from other Carpet thorns directly.schnetter2003-06-18
| | | | | | | | Do not include header files from other Carpet thorns directly. Instead, use the "INCLUDES HEADER" and "USES INCLUDE HEADER" mechanism. darcs-hash:20030618162807-07bb3-a81444cde6c76e6a24516d108861fc1b5541c643.gz
* Major update after a quiet time.schnetter2003-06-18
| | | | | | | | | | | | | | | | | | | | Major update after a quiet time. Carpet: The flesh now has new cGH fields cctk_levoff[], cctk_levoffdenom[], and cctk_timefac that describe the spatial offset and temporal refinement factor between the base and the current refinement level. These fields are now set and used; they change how coordinates are handled. CarpetIOASCII: Fix bugs regarding choosing the output hyperslab and the output coordinates. ID*, *Toy*: New WaveToy examples with various formulations and different integrations methods. Currently, none of them converge to second order except the standard WaveToy formulation. These updates require the recent flesh, base thorn (and MoL) updates. darcs-hash:20030618162427-07bb3-70761f74bce6ae246b5a2943a385647657d46d34.gz
* Do not use "restrict".schnetter2003-05-21
| | | | darcs-hash:20030521140332-07bb3-e588734900063e2139bc9a62e0f697c1faafa68d.gz
generated by cgit v1.2.3 (git 2.39.1) at 2024-07-02 14:47:06 +0000