aboutsummaryrefslogtreecommitdiff
path: root/CarpetDev
Commit message (Collapse)AuthorAge
...
* CarpetJacobi: Correct time level cyclingErik Schnetter2005-12-22
| | | | | | | | | | | | When the solution has multiple time levels, the solution needs to be copyied from the past to the current time level after time level cycling. Add some debugging code, disabled by default. Correct handling of restrict qualifier. darcs-hash:20051222195540-dae7b-5c276af06aed498dbbea3f23f589641377c2bfa7.gz
* CarpetJacobi: Update parameter filesErik Schnetter2005-12-22
| | | | darcs-hash:20051222195438-dae7b-f54bd47a19eaa525a6c915aaacb77aa3ea83a900.gz
* CarpetJacobie: Remove $Header$ tagsErik Schnetter2005-12-22
| | | | darcs-hash:20051222195356-dae7b-dd6726ec26eb58f34d649c39e0ad2bbfc100291d.gz
* CarpetMG: Correct #define restrictErik Schnetter2005-11-22
| | | | darcs-hash:20051122010740-dae7b-8cf23d20a225ea887eca3bf8b9833accd344682c.gz
* CarpetJacobi: Correct #define restrictErik Schnetter2005-11-22
| | | | darcs-hash:20051122010715-dae7b-5cc113b99f976336518c75076ad0ec0761e1deeb.gz
* CarpetCG: Correct #define restrictErik Schnetter2005-11-22
| | | | darcs-hash:20051122010636-dae7b-b8eb2c0c566e724e655010ca3eae0b5d5ecd81f0.gz
* CarpetMG: #define restrict according to CCTK_CXX_RESTRICTErik Schnetter2005-11-20
| | | | darcs-hash:20051120212633-dae7b-06eb5c442f31b1c06c9237dbd37d1d46404274fc.gz
* CarpetJacobi: #define restrict according to CCTK_CXX_RESTRICTErik Schnetter2005-11-20
| | | | darcs-hash:20051120211722-dae7b-a69341f50c2f8d4f0d194c75bdbf63f7dc878aad.gz
* CarpetCG: Add some empty linesErik Schnetter2005-11-20
| | | | darcs-hash:20051120211704-dae7b-cb05821a52392f2210a868360020d4b9a1ab75c2.gz
* CarpetCG: #define restrict according to CCTK_CXX_RESTRICTErik Schnetter2005-11-20
| | | | darcs-hash:20051120211116-dae7b-51488736c2244728c943db954c06eb4c2e4d880e.gz
* global: Add varying refinement factorsErik Schnetter2005-05-01
| | | | | | | | | | 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
* CarpetMG: Turn CarpetLib templates into classesErik Schnetter2005-01-01
| | | | darcs-hash:20050101182657-891bb-ae26379b5a248728005855c6db4dcdee6264ed62.gz
* Use ipow instead of pow in various placesErik Schnetter2004-12-25
| | | | darcs-hash:20041225201943-891bb-0c5fe9240c54991c38b7710ab54a0de71d69e221.gz
* CarpetMG: Snapshot of current stateErik Schnetter2004-12-08
| | | | darcs-hash:20041208222153-891bb-d9ac7d30cfb9ed9eda22377bcb72317b349cbc73.gz
* CarpetMG: A multigrid solver for CarpetErik Schnetter2004-11-24
| | | | darcs-hash:20041124061259-891bb-566d93667d53c7746838a83c4171f6874dd496f2.gz
* CarpetAdaptiveRegrid:Really fix the inheriting child errors this time, I hope.Ian Hawke2005-03-04
| | | | darcs-hash:20050304163739-58c7f-6e0236559166e77d7562b96370e4ffa95068eee1.gz
* CarpetAdaptiveRegrid: Fix checking the child level for errorsIan Hawke2005-03-04
| | | | | | | | | Firstly only check the errors of a child level if a grandchild level exists. Secondly take the refinement factor into account when calculating where this error should be placed! darcs-hash:20050304134516-58c7f-484d0c7d61a089b4eb7f294c427c76bb529d8dde.gz
* CarpetAdaptiveRegrid: ASCII pictures of the maskIan Hawke2005-03-04
| | | | | | | Print to screen (using cout; not good on multiprocs) an ASCII art picture of the mask for debugging purposes. darcs-hash:20050304120741-58c7f-9cc7fbbb8462dafc558b974dcc1a04a28b5e7987.gz
* CarpetAdaptiveRegrid: Change implementation name to CarpetAdaptiveRegridErik Schnetter2005-03-03
| | | | | | | | | Change implementation name from CarpetRegrid to CarpetAdaptiveRegrid. This is necessary because the thorns CarpetRegrid and CarpetAdaptiveRegrid export different parameters -- CarpetRegrid supports moving boxes, and CarpetAdaptiveRegrid does not. darcs-hash:20050303102220-891bb-fb35e947b47f660ceb384a8ecf175974a1a60c5c.gz
* CarpetAdaptiveRegrid: Code documentationih2005-02-25
| | | | | | | | As well as documenting the code a bit, add some fixme's (including noting that the code probably won't work correctly with checkpoint/restart yet!). darcs-hash:20050225203643-0ff1f-ae45e424aa69d52e9381a9a011d2242b58944fe7.gz
* CarpetAdaptiveRegrid: Documentationih2005-02-25
| | | | | | | | | | | | | Add some documentation for the AMR, and also an example parameter file (shock tube, 6 levels). Note that my last message about parallel AMR was incorrect; the difference between the 1 and 2 processor runs was entirely due to a bug in Whisky, not in the AMR. Once that was fixed, the results now agree to output precision. Therefore I pronounce parallel AMR "working". darcs-hash:20050225194759-0ff1f-967f0360b73244634d422198c9446f812302c916.gz
* CarpetAdaptiveRegrid: parallel supportih2005-02-24
| | | | | | | | | | | | Make CAR work on more than one processor (that is, I've tested it on my laptop using mpirun -np 2...). Note that at present whilst it runs and gives reasonable answers you do _not_ get the same answers as you would on 1 proc. This is a bug, as yet unfound. But, there probably are lots of other bugs in this code still to be found... darcs-hash:20050224223702-0ff1f-f766381f5309c747ea438dbd10976a99c86abd48.gz
* CarpetAdaptiveRegrid: Allow non-trivial initial gridsih2005-02-23
| | | | | | | Use the CarpetRegrid::manual-coordinate-list notation to allow the setting up of initially refined grids. Literal cut 'n paste in many places. darcs-hash:20050223235501-0ff1f-4be0e5c479aadabc7e1c5682476f95a23c00ec71.gz
* CarpetAdaptiveRegrid: Correct striding correctionIan Hawke2005-02-23
| | | | darcs-hash:20050223105949-58c7f-e087970950767cb32adb21fc617d7789a6b988ec.gz
* CarpetAdaptiveRegrid: Only regrid once per iterationih2005-02-23
| | | | | | | Because the routine now regrids all finer levels than the one on which it is called it is not necessary to regrid more than once per iteration. darcs-hash:20050223143548-0ff1f-d50f953738e11a8653e5eb98a13491b7810cf0dd.gz
* CarpetAdaptiveRegrid: Enforce proper nestingIan Hawke2005-02-22
| | | | | | | | Enforce proper nesting both for parents and children by looping over all finer levels regridding them all at once. This completely destroys efficiency... darcs-hash:20050222174437-58c7f-d280bacafa2fb4c634613381b4c7134733348ec5.gz
* Correct paddingih2005-02-21
| | | | | | | Correct the way the errors are padded (an off-by-one error). Also be a bit more careful with F90 intent statements. darcs-hash:20050221195416-0ff1f-9a41dbe1c6c26563fb6b69513fdcabfa887f29d7.gz
* Fix adaptive patch hitting boundaryIan Hawke2005-02-18
| | | | darcs-hash:20050218164449-58c7f-a63ebeaedb7ea80c398986e7267af95dd0dd4679.gz
* Multilevel AMRIan Hawke2005-02-18
| | | | | | | | | Basic support for more than 2 refinement levels with AMR. Basic tests done with shock tubes seem to show that it works. Note that proper nesting is not ensured as yet as we don't check for orphaned grandchildren. darcs-hash:20050218120448-58c7f-5d2e6605ac05f1363596c0e61f6026e8fc057ff8.gz
* CarpetAdaptiveRegrid: Use CarpetIOScalar instead of IOBasic for scalar outputErik Schnetter2005-02-11
| | | | darcs-hash:20050211161343-891bb-819221d7878f630867fc2cf7a6d94c2893e5a37b.gz
* AMR boundariesIan Hawke2005-02-11
| | | | | | | | Add support for using boundaries with AMR. Seems to work with symmetry boundaries but gets a bit confused if you try refining outer boundaries. Requires that you setup the domain using CoordBase. darcs-hash:20050211135704-58c7f-6ed6ce1449de17599dd8b23abc909f5ba313d048.gz
* Correct outer boundary checkIan Hawke2005-02-10
| | | | | | | | Correct the outer boundary check so that there's no more assertion failures. However, the outer boundary is of course not dealt with correctly yet. darcs-hash:20050210134049-58c7f-e1e31185738d7ef7cd4f2ae7517b3722e0103cfa.gz
* Comment out debugging statementsIan Hawke2005-02-09
| | | | darcs-hash:20050209173245-58c7f-862ef5eb1de5c153c629079ebe2e866888348aa6.gz
* Adaptive regridding symmetry / indexing partial fixIan Hawke2005-02-09
| | | | | | | | | | Fixed a bug where the array indexing in C++ was wrong. Now something that should be symmetric almost is. There is still an assertion failure when it tries to derefine - for the moment this appears to be a Carpet internal problem. darcs-hash:20050209161923-58c7f-7d13dd7c72c44d935bff9050e3d9281f5f39ec90.gz
* Adaptive Regridding might workIan Hawke2005-02-08
| | | | | | | | Changes to the adaptive regridding routine such that it appears to work for short times on one processor with various changes to, e.g., WaveMoL, before hitting internal Carpet problems. darcs-hash:20050208180418-58c7f-22dc56373df5a208b82d45604267f398bfd3ca0f.gz
* CarpetJacobi: Apply boundary conditions in level modeErik Schnetter2004-12-08
| | | | darcs-hash:20041208222256-891bb-ae079ce98bccbd631cc8bf41e32902bec31edcf9.gz
* CarpetAdaptiveRegridhawke2004-11-13
| | | | | | | | The first part of CarpetAdaptiveRegrid. At the moment the only piece of "working" code is the F90 that does the box manipulation. This can be tested using the standalone makefile in src/test. darcs-hash:20041113152820-5d2e4-7bbcb3804e5643c7753c9ab663df4dc7d357d7e9.gz
* Do not use CVS $Date$ for latex documentation any moreErik Schnetter2004-11-18
| | | | darcs-hash:20041118233344-891bb-44b65f16949b0d78bee0a66ad6f65c262652c671.gz
* Remove CVS header tagsErik Schnetter2004-09-18
| | | | | | Replace all CVS header tags with the standard "$Header:$". darcs-hash:20040918132147-891bb-dea889bdd94a479ec412d14d08e9efca63e5c24d.gz
* Activate thorn SymBaseschnetter2004-03-23
| | | | darcs-hash:20040323095920-07bb3-b9e755b6b0a64287b8e569d1f70b16633e473b1e.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
* Fix the memory corruption.hawke2003-11-20
| | | | | | | Fix the memory corruption. Doesn't converge with refinement, though. darcs-hash:20031120073403-58737-0a5f287cad04535eac0f5403862605f33e020a1f.gz
* Only switch to global mode if not already in global mode.schnetter2003-11-19
| | | | darcs-hash:20031119130536-07bb3-7d6722ae739c4f3315fdfe85b3e0ac3dc899454e.gz
* Activate CoordBase.schnetter2003-11-19
| | | | darcs-hash:20031119130523-07bb3-613b78a757bf8723f7dbfd8d679a70125674bc87.gz
* Initial revisionhawke2003-11-19
| | | | darcs-hash:20031119094047-58737-ef9911a12e56712a462cb6e555f79d0fa4a7dd6a.gz
* Initial revisionhawke2003-11-19
| | | | darcs-hash:20031119084047-58737-b7481050914bb163f0867f194597c7a50b26e64e.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
* Add elliptic solver thorn that works with Carpet.schnetter2003-09-02
| | | | darcs-hash:20030902123557-07bb3-67dd485b689f921257309f4eeb7389821fd5fa7a.gz
* Initial revisioneschnett2001-03-01
darcs-hash:20010301114010-f6438-12fb8a9ffcc80e86c0a97e37b5b0dae0dbc59b79.gz