aboutsummaryrefslogtreecommitdiff
path: root/Carpet/CarpetIOHDF5
Commit message (Collapse)AuthorAge
...
* CarpetIOHDF5: avoid accessing non existend refinement levels for grid arraysroland haas2011-12-14
|
* CarperIOHDF5: Correct stride in recoveryErik Schnetter2011-12-14
|
* CarpetIOHDF5: Correct stride for grid arraysErik Schnetter2011-12-14
|
* CarpetIOHDF5: Use hg::baseextent instead of Carpet::maxspacereflevelfactErik Schnetter2011-12-14
| | | | | | Use hg::baseextent instead of Carpet::maxspacereflevelfact to determine the stride of a refinement level, because this works independent of the stride on the finest level.
* * encapsulate read of potentially not present attributecott2011-12-14
| | | | | in an H5E_BEGIN_TRY H5E_END_TRY environment. This prevents hate crimes.
* CarpetIOHDF5: Store cell centering offset with grid function attributesErik Schnetter2011-12-14
|
* CarpetIOHDF5: make sliced 3D output depend on out3D_XXXRoland Haas2011-12-14
| | | | rather than always ocuring together with the old-style 3D output.
* CarpetIOHDF5: Checkpoint and recover delta times (time step sizes)Erik Schnetter2011-12-14
|
* CarpetIOHDF5: Check that MPI datatypes are defined before using themErik Schnetter2011-12-14
|
* CarpetIOHDF5: Remove grid hierarchy consistency tests that don't apply hereErik Schnetter2011-12-14
|
* CarpetIOHDF5: add example parameter file for output_buffersRoland Haas2011-12-14
| | | | | | right now there is no support for hdf5 files in testsuites so users will have to manually use 'h5diff -v' and 'h5dump -r' to compare and test the output when changes are made to the code
* CarpetIOHDF5: fix way hyperslab is computedRoland Haas2011-12-14
| | | | | to compute the hyperslab one needs to base its location on the extent of the data in memory not in the extent that we want to output.
* CarpetIOHDF5: make sure metadata matches old codeRoland Haas2011-12-14
| | | | | * add some attributes to metadata to match what the old code wrote * only tag object names with components of there are more than one
* CarpetIOHDF5: add option to not output the buffer zonesRoland Haas2011-12-14
| | | | | | | This is implemented by re-computing the allactive bboxset that dh::regrid computes and outputting only the part of each component that intersects allactive. This changes the number of components since the intersection might be non-rectangular-box.
* CarpetIOHDF5: add support for output on multiple processorsRoland Haas2011-12-14
| | | | | add code to honor IO::ioproc_every to output data unchunked (for 3d output only)
* CarpetIOHDF5: turn on 3d output in OutputSlice.ccRoland Haas2011-12-14
| | | | | | | | Register new-style output code as an IO method with IOUtil for 3d output. The new code is not quite as capable as the old code, since it does not include Ian Hinder's indexing facility and so far outputs all data on the root processor. It does support the new-style output_symmetry_points etc. options however.
* CarpetIOHDF5: Set times correctly after recoveringErik Schnetter2011-12-14
| | | | | | Recovery happens level-by-level in Cactus. When recovering the refinement level times and the global time, set them correctly according to the current refinement level.
* CarpetIOHDF5: Correct iorigin attributeErik Schnetter2011-12-14
|
* CarpetIOHDF5: Index file supportIan Hinder2011-12-14
| | | | | | | | | | | | | Scanning the attributes of a large CarpetIOHDF5 output file, as is necessary in the visitCarpetHDF5 plugin, can be very time consuming. This commit adds support for writing an "index" HDF5 file at the same time as the data file, conditional on a parameter "CarpetIOHDF5::output_index". The index file is the same as the data file except it contains null datasets, and hence is very small. The attributes can be read from this index file instead of the data file, greatly increasing performance. The datasets will have size 1 in the index file, so an additional attribute (h5space) is added to the dataset to specify the correct dataset dimensions.
* CarpetIOHDF5: Correct error in outputting distributed grid arraysErik Schnetter2011-12-14
|
* CarpetIOHDF5: Comment out unused local variableErik Schnetter2011-12-14
|
* CarpetIOHDF5: Update to new dh classesErik Schnetter2011-12-14
|
* CarpetIOHDF5: Remove explicit bboxset normalizationsErik Schnetter2011-12-14
|
* CarpetIOHDF5: Initialise times on all time levels while recoveringErik Schnetter2011-12-14
| | | | | Initialise the times of all time levels of grid arrays while recovering.
* CarpetIOHDF5: Improve handling termination checkpointsErik Schnetter2011-12-14
| | | | | When checkpoints of initial data are disabled, but termination checkpoints are enabled, then do checkpoint the initial data.
* CarpetIOHDF5: Add comments and indent one lineErik Schnetter2011-12-14
|
* Re-organise time level handlingErik Schnetter2011-12-14
| | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | | Store the current Cactus time (and not a fake Carpet time) in the th "time hiearchy". This removes the now redundant "leveltimes" data structure in Carpet. Add past time levels to th, so that it can store the time for past time levels instead of assuming the time step size is constant. This allows changing the time step size during evolution. Share the time hierarchy between all maps, instead of having one time hierarchy per map. Simplify the time level cycling and time stepping code used during evolution. Improve structure of the code that loops over time levels for certain schedule bins. Introduce a new Carpet variable "timelevel", similar to "reflevel". This also makes it possible to avoid time interpolation for the past time levels during regridding. The past time levels of the fine grid then remain aligned (in time) with the past time levels of the coarse grid. This is controlled by a new parameter "time_interpolation_during_regridding", which defaults to "yes" for backwards compatibility. Simplify the three time level initialisation. Instead of initialising all three time levels by taking altogether three time steps (forwards and backwards), initialise only one past time level by taking one time step backwards. The remaining time level is initialised during the first time step of the evolution, which begins by cycling time levels, which drops the non-initialised last time level anyway. Update Carpet and the mode handling correspondingly. Update the CarpetIOHDF5 checkpoint format correspondingly. Update CarpetInterp, CarpetReduce, and CarpetRegrid2 correspondingly. Update CarpetJacobi and CarpetMG correspondingly.
* CarpetIOHDF5: Correct more attribute namesErik Schnetter2011-12-14
|
* CarpetIOHDF5: Correct attribute nameErik Schnetter2011-12-14
|
* CarpetIOHDF5: Add commentErik Schnetter2011-12-14
|
* CarpetIOHDF5: Use parameter IO::out_timesteps_per_fileErik Schnetter2011-12-14
|
* CarpetLib: Introduce new virtual function new_typed_dataErik Schnetter2011-12-14
|
* Combine CarpetLib's INSTANTIATE and Carpet's TYPECASE mechanism into aErik Schnetter2011-12-14
| | | | | | single mechanism provided by CarpetLib. Use this mechanism everywhere.
* CarpetIOHDF5: Handle both old and new IOUtil APIErik Schnetter2011-12-14
|
* CarpetIOHDF5: Allow different numbers of ghost zones on different levelsErik Schnetter2011-12-14
| | | | | | | Allow different numbers of ghost zones and different spatial prolongation orders on different refinement levels. This causes incompatible changes to the checkpoint file format.
* CarpetIOHDF5: Increase output accuracy for time hierarchyErik Schnetter2011-12-14
| | | | | | Serialise the time hierarchy with 17 digits of accuracy, instead of the standard 6 digits. This is required for correctly recovering the time hierarchy.
* CarpetIOHDF5: Rename some parametersErik Schnetter2011-12-14
| | | | | | Rename out3D_ghosts to output_ghost_points. Rename out3D_outer_ghosts to output_boundary_points. Keep the old parameter name for compatibility.
* CarpetIOHDF5: Output a warning if multiple input files need to be readErik Schnetter2011-12-14
| | | | | | | Output a warning message if multiple input files need to be read from one MPI process, since this is usually very slow. When reading files from the same number of processes that wrote them, each process is only supposed to need to open one file.
* CarpetIOHDF5: #include <cstddef>Erik Schnetter2011-12-14
| | | | Ignore-this: e7d15a216ec4fbb524f1ebe1fdeff905
* CarpetIOHDF5: Correct small error in commentErik Schnetter2011-12-14
| | | | Ignore-this: 9e10c4fe3b231a3913fca9572d319474
* Import CarpetErik Schnetter2011-12-14
| | | | Ignore-this: 309b4dd613f4af2b84aa5d6743fdb6b3
* Merge branch 'master' of /Users/eschnett/Cbeta/carpet/Erik Schnetter2009-01-25
|\
| * CarpetIOHDF5: fix a bug in the recovery of multipatch simulationsThomas Radke2009-01-19
| | | | | | | | Due to a wrong upper range in the time hierarchy initialisation loop, only maps on the coarsest refinement level were initialised. This caused an assertion failure when recovering multiple refinement levels which weren't aligned.
* | Merge branch 'master' of /Users/eschnett/Cbeta/carpet/Erik Schnetter2008-12-08
|\|
| * CarpetIOHDF5: in case the grid structure cannot be recovered properly, abort ↵Thomas Radke2008-12-08
| | | | | | | | simulation with a proper error message (rather than just an assertion failure)
* | Merge branch 'master' of /Users/eschnett/Cbeta/carpet/Erik Schnetter2008-10-30
|\|
| * CarpetIOHDF5: also activate LoopControl (as required by CarpetReduce) in ↵Thomas Radke2008-10-30
| | | | | | | | example and testsuite parfiles
| * CarpetIOHDF5: let the hdf5_slicer write out a few more attributes which are ↵Thomas Radke2008-10-30
| | | | | | | | needed by OpenDX and VisIt
* | Merge branch 'master' of /Users/eschnett/Cbeta/carpet/Erik Schnetter2008-10-23
|\|
| * CarpetIOHDF5: added hdf5_slicer utility to extract 1D/2D slices from 3D ↵Thomas Radke2008-10-17
| | | | | | | | Carpet HDF5 output data. The slice data are output in HDF5 again.