aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetReduce
Commit message (Collapse)AuthorAge
* 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
* typecase_INCLUDE_HEADER_INswhite2004-12-01
| | | | | | | | 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
* Do not use CVS $Date$ for latex documentation any moreErik Schnetter2004-11-18
| | | | darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
* CarpetRegrid: clean up time interpolationErik Schnetter2004-10-08
| | | | darcs-hash:20041008214326-891bb-92f62ca9e16ba5acb5d4b79b1865057fb345711d.gz
* Small cleanup in CarpetReduceErik Schnetter2004-09-22
| | | | darcs-hash:20040922204131-891bb-e4070a4c2cee6417e961ca1aa283aeb98bc695d8.gz
* Do not set weight to zero on ghost zonesErik Schnetter2004-09-22
| | | | | | Do not set weight to zero on ghost zones. Ghost zones are ignored anyway. darcs-hash:20040922204017-891bb-26f1608179b6415cebf161275751ea767cf8d8f8.gz
* Correct error in weight calculation on multiple processorsErik Schnetter2004-09-22
| | | | darcs-hash:20040922203956-891bb-9d66f72ca19242dc7eda1208732a4acd6c5ee923.gz
* Say that CarpetRegrid::weight has only one time levelErik Schnetter2004-09-22
| | | | darcs-hash:20040922203920-891bb-59f30143f8d685db70da4626c5129e1908484eae.gz
* Relax consistency test in reduction operationsErik Schnetter2004-09-21
| | | | | | | Use relative instead of absolute error when checking that the times are aligned during reduction. darcs-hash:20040921161026-891bb-361d8f344ec91f7d1392101267fa94114d1d18e5.gz
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* Do not loop over empty bboxes, because the imin/imax calculation wouldschnetter2004-08-05
| | | | | | | Do not loop over empty bboxes, because the imin/imax calculation would become confused. darcs-hash:20040805091549-07bb3-4523e04a4f95a2acdd80816e794215ae1f281bcf.gz
* Correct error in setting up the edges and corners of the weightschnetter2004-08-04
| | | | | | | Correct error in setting up the edges and corners of the weight function. darcs-hash:20040804110328-07bb3-79be97be201e736d66480f0774e877fa9006c39d.gz
* Correct output text.schnetter2004-08-04
| | | | darcs-hash:20040804110309-07bb3-258131fe3ce936fd74cf27fc043c3342505d5d99.gz
* Calculate the weight function correctly.schnetter2004-08-02
| | | | darcs-hash:20040802094335-07bb3-8cf89dfef6eab7c25a50842c68b62d8a7c230d57.gz
* Support multi-patch.schnetter2004-07-06
| | | | darcs-hash:20040706124432-07bb3-194cca535ccc5d222c4abc99c6a6dd678e61732b.gz
* Add "if weight > 0" statements to prevent reducing nans on coarseschnetter2004-07-02
| | | | | | | Add "if weight > 0" statements to prevent reducing nans on coarse grids, which lead to an unjustified nan result. darcs-hash:20040702154738-07bb3-b4078aeca351666affec53ae2b6e0a5f4747b897.gz
* Implement true global reduction operations, i.e. time interpolation.schnetter2004-06-27
| | | | darcs-hash:20040627191808-07bb3-6cfb448d84dba6049ee25ba1fccc476cf50be746.gz
* Remove warnings about group and variable names being distinct.hawke2004-06-21
| | | | darcs-hash:20040621103707-58737-ae46e243d43f093ca85965fbf4c0d3f11e710a08.gz
* Don't call MPI reduction operators for complex numbers; call the realschnetter2004-06-21
| | | | | | | Don't call MPI reduction operators for complex numbers; call the real operators twice instead. darcs-hash:20040621102753-07bb3-2391f206857569e3d6b38bef8dc7e66383446433.gz
* Don't require the local boxes to be non-empty.schnetter2004-06-15
| | | | darcs-hash:20040615144529-07bb3-781ff26e87e690804ad8f2faf5294df94847be25.gz
* Implement global reduction operators:schnetter2004-06-14
| | | | | | | | | | | | Implement global reduction operators: Add a weight grid function. Initialise it. Set it according to the symmetries. Set it according to the refinement structure. Take the weight into account when reducing. darcs-hash:20040614050121-07bb3-39cf335a52bc2c1071e56375bbedf771a5f118d4.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
* no changeschnetter2004-05-21
| | | | darcs-hash:20040521161353-07bb3-88f01ca8c5fb75607ec40b9229f8ea170c1cd144.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
* Enable multiple mapsschnetter2004-05-04
| | | | darcs-hash:20040504201149-07bb3-cf7f9bfd2004d97467cac41637826e0fe00193b7.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
* Explicitly cast constness awayschnetter2004-03-31
| | | | darcs-hash:20040331143832-07bb3-7a1e263661a7034f7ad6d38b048ce79015d4f20d.gz
* Fix bug in looping over refinement levels.schnetter2004-03-31
| | | | darcs-hash:20040331105117-07bb3-75627a90c4e513c813c224f68dc22892e39aac35.gz
* Restructure the conversion between CCTK_COMPLEX* andschnetter2004-03-01
| | | | | | | | Restructure the conversion between CCTK_COMPLEX* and complex<CCTK_REAL*> types. Cleaner, more readable, more structured, more comments, etc. darcs-hash:20040301203639-07bb3-b5eca730ad2f601ee03b83cef3f5dacd00488941.gz
* Fix definition of mymin and mymax for complex values.schnetter2004-02-01
| | | | darcs-hash:20040201135456-07bb3-9c7b217df919037f256454b49062f7231a68fd00.gz
* Fixed some namespace problems by explicitly specifying namespace stdswhite2004-01-27
| | | | | | | | | | | | | | | | | | 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: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
* Fix bug in the maximum reduction operator.schnetter2003-11-12
| | | | darcs-hash:20031112162930-07bb3-eb9a8b4a6c167c98bfeca86bc0c7b10baed922b5.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
* 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
* Make fit for complex values.schnetter2003-07-30
| | | | darcs-hash:20030730133311-07bb3-a84302f705a727b41c92d3166913b392c9799910.gz
* Adapt local array reduction to PUGH API.schnetter2003-07-23
| | | | darcs-hash:20030723124610-07bb3-0e89ee3460955e33f8ac3c2b90209e52b0093295.gz
* Remove leftover debugging statements.schnetter2003-07-23
| | | | darcs-hash:20030723123044-07bb3-d56e7772784e192029203b5d7cecfde4eb4b7ba9.gz
* Change API of grid variable reductions to match PUGH.schnetter2003-07-23
| | | | darcs-hash:20030723074706-07bb3-188ace74c5d748c8c19437ef72bbfa525450e1dd.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
* Allow reducing grid arrays.schnetter2003-05-21
| | | | darcs-hash:20030521123103-07bb3-b55fc535375111abe37293b2e9e4b907ed8f5e30.gz
* Make C++ arrays sizes known at compile time.schnetter2003-05-14
| | | | darcs-hash:20030514061252-07bb3-5f3d15a6be7b17fc7599544c4bb193c1a1546bf1.gz
* Make code const-correct. Well, sort of.schnetter2003-05-13
| | | | darcs-hash:20030513101942-07bb3-2a9b4aa40aafc2566a13522ff7704e480e9817b6.gz
* Handle empty grids.schnetter2003-05-13
| | | | | | | | | | Handle empty grids. Fix global and level mode handling of the interpolator. Allow local mode calls for the interpolator, reduction operator, and hyperslabber, assuming that people know what they are doing. darcs-hash:20030513101400-07bb3-42af7544c2ea98dbef5a7ad6d5112bea2dcc404f.gz
* Abort when in global mode, because this is not yet handled.schnetter2003-05-12
| | | | darcs-hash:20030512142528-07bb3-f6924002fd453379da2af0a15130f772cb7a05b0.gz
* Remove the special case of "extended local mode" which is now handledschnetter2003-05-02
| | | | | | | Remove the special case of "extended local mode" which is now handled automatically. darcs-hash:20030502122344-07bb3-a0b56e30126a8625374118c026ca090a1082b451.gz
* Use vector<> instead of pointers.schnetter2003-04-30
| | | | darcs-hash:20030430104102-07bb3-d58d092da89450a170cdd2cb6c6fc0c73b81f1fb.gz
* Fix serious bugs with counting dimensions.schnetter2003-04-29
| | | | | | | | Fix serious bugs with counting dimensions. Add explicit constructors for some classes. Add documentation template. darcs-hash:20030429120152-07bb3-61d82a853ea12abcf94d6d300217fd093f5ef0b9.gz