jonorbeb
f45bed11a9
Dark matter only art module tested and working ok. Units checked and ok.
...
All previous questions solved.
Still waiting for one last zoom test (moderate resolution 512^3) to compare it with
gadget results but full box and low res zoom look ok.
We should start thinking in the baryon format output.
2012-09-04 20:28:22 -07:00
jonorbe
bf38faaf23
Swap endianness option fully enabled for ART module: header, particle data and pt files.
2012-08-19 16:43:31 -07:00
Oliver Hahn
9b34875aa5
added change of endianness to data arrays in ART plugin, header etc. still need to be adjusted
2012-08-17 17:38:35 -07:00
jonorbe
ab24981135
- Fixed velocity units conversion. Still final test needed.
...
- Header ascii variable can be set in the config file.
- astep variable for ART has to be given in the config file.
2012-07-31 01:36:58 -07:00
jonorbe
c7e6cd5825
ART module for dm-only operative. Right now the module generates all the files necessary with the correct format and (re-check mandatory) units. Therefore, ICs generated should run!
...
Here I write a list of to-do things:
- Find a clever way to set the char array in header. Should be set in the config file and then probably use some getValue?
- Re-check velocity units
- How important is the EKIN parameter of the header to start ART. We are not calculating it. It could be done when writing the files but seems to much work taking into account that it is probably calculated at each timestep.
- astep parameter: clarification on how is set/chosen for the initial conditions. It is important for subsequent steps or ART auto-regulates the next steps?
- Need some clarifications for zoom-in simulations:
+ which NROWC and NGRIDC is used (now using levelmin, but levelmax?)
+ Number of files. Not clear if all together in one or one for each level (in the last case we will need to do some important changes in the module)
- Confirm our current format is standard for all ART versions.
2012-07-27 20:17:12 -07:00
Oliver Hahn
3806e940c6
added code outline to read DM temporary particle files and assemble them
2012-07-27 12:27:27 -07:00
jonorbe
39025f7173
Assigning correct values to the parameters in header file.
2012-07-27 11:43:58 -07:00
Oliver Hahn
2ed1414fba
added correct averaging of SPH+DM particles in coarse region
2012-07-23 12:16:05 -07:00
Oliver Hahn
448796873d
added temp file output for DM data in ART plugin
2012-07-20 19:36:59 -07:00
Oliver Hahn
c3f5a14438
made art plugin compilable
2012-07-20 19:27:33 -07:00
Oliver Hahn
b7dd1cd8e3
added baryons to tipsy file output plugin
2012-07-20 19:17:18 -07:00
jonorbe
6d43639fc6
changes output_art.cc
2012-07-19 20:26:19 +02:00
jonorbe
aa80db066a
Changes is output_art.cc
2012-07-19 19:41:57 +02:00
jonorbe
e121d9547c
Some changes in output_art.cc
2012-07-19 19:14:34 +02:00
jonorbe
bb700ebd7c
Added output_art.cc file.
2012-07-19 18:31:29 +02:00
Oliver Hahn
f10f8f9a6c
version freeze for Santa Cruz Code Comparison Project
2012-07-11 17:15:45 -07:00
Oliver Hahn
3c12876d5e
Added a temporary workaround for the 'music' transfer function plugin that
...
catches cases where the transfer functions become negative at very large k.
2011-09-07 22:28:39 -07:00
Oliver Hahn
06d1f88aa2
Gadget plugin: Fixed an overflow bug when using more than 2048**3 particles and baryons.
2011-09-06 17:08:36 -07:00
Oliver Hahn
cfc16f509d
Added new option cosmology/z0norm to 'linger++' plugin. This allows to choose
...
the normalization at z=0 or the linearly scaled from z_ini.
2011-09-06 16:39:41 -07:00
Oliver Hahn
e7f8c2e1ec
Added a first DM only version of a tipsy output plugin.
...
Credits go to Kyle Stewart and Shea Garrison-Kimmel.
2011-09-02 15:21:55 -07:00
Oliver Hahn
e6137dcf07
Fixed a bug that prevented the 'MUSIC' transfer function plugin to be used when
...
also the 'LINGER++' one was present (they used the same class name, so which one
was actually used depended on the compile order rather than the user choice).
2011-08-29 16:12:14 -07:00
Oliver Hahn
29822726ec
Removed option "shift_back" for Gadget2 output as it can cause problems with zoom sims overlapping the periodic boundaries.
2011-08-29 14:00:47 -07:00
Oliver Hahn
ebd300f2fa
Some more bugfixes for Gadget-2 output with very large grids.
...
Verified to work with 2048**3.
2011-06-09 12:08:24 -07:00
Oliver Hahn
6bc18e1aaf
More 64 bit compliance fixes for output plugins.
...
Gadget-2 output plugin now can write multiple output files.
Saved a bit of memory in main loop
2011-06-08 20:58:42 -07:00
Oliver Hahn
03fb47cc3d
Added support for very large arrays to enzo output plugin.
...
Enzo plugin now writes in slabs rather than in one monoloithic block (saving memory).
2011-06-03 18:41:39 -07:00
Oliver Hahn
e470be5da7
added support for new linger++ transfer function file format
...
cleaned up a bit of code
2011-06-02 11:12:44 -07:00
Oliver Hahn
b7756e22a0
full64bitcompliance
2011-06-01 11:01:24 -07:00
Oliver Hahn
df3ee81111
Fixed a bug in the ENZO output plugin for non-cubical subgrids when used in conjunction with PartitionNestedGrids=1
2011-04-27 10:46:44 -07:00
Oliver Hahn
af2f45763d
Fixed a bug for Gadget2 output with baryons that caused low-res DM particles to have the wrong mass
...
when using more than two levels
2011-04-15 10:04:40 -07:00
Oliver Hahn
731641524c
Added output of particle displacement file names to template parameter file in enzo output plugin.
2011-04-12 09:55:34 -07:00
Oliver Hahn
7cef5e69cc
Fixed a bug when the transfer function tables contain negative values by clipping them off
2011-04-07 16:19:16 -07:00
Oliver Hahn
76a7a0559e
Fixed a bug introduced with the last "fix", where a conversion to size_t fails.
2011-03-16 17:20:06 -07:00
Oliver Hahn
088fb32c23
Fixed some compiler warnings (and possibly errors with GCC 4.3)
2011-03-16 15:33:03 -07:00
Oliver Hahn
8afc1f9b4a
* fixed a bug in 'music' transfer function, where spline tables where not allocated correctly
...
* fixed a warning message of an unsigned loop variable with openmp in poisson.cc
2011-02-23 02:20:04 -08:00
Oliver Hahn
fdd12c3d9c
minor fixes (mainly constraints)
2010-12-16 12:05:03 +01:00
Oliver Hahn
ff5d87273c
Some minor fixes. No longer uses FFT to compute 2LPT terms when in unigrid mode and
...
when using deconolve=true (default=true).
2010-12-13 08:47:01 -08:00
Oliver Hahn
fa8cde6d29
Fixed a bug that cause negative grid offsets in some cases when overlapping with box edges.
2010-11-24 14:37:19 -08:00
Oliver Hahn
c927462924
Fixed a bug that caused the baryon velocity field to be not shifted in the SPH
...
case when a transfer function with no velocity info was used.
2010-11-18 22:30:16 -08:00
Oliver Hahn
2ea869e85a
Fixed a bug where memory was deallocated that shouldn't be
...
Re-enabled the staggering for SPH particles.
2010-11-17 21:24:21 -08:00
Oliver Hahn
8d9068852e
* IMPORTANT: removed code left over from testing that caused suppression of small scale perturbations
...
* added baryon staggering for SPH
* added 2LPT for baryons/DM two-component fluid
* cleaned up some code
2010-11-16 16:41:40 -08:00
Oliver Hahn
27fe21e3d4
Migrated several implementations to .cc files.
...
Some minor bugfixes.
No longer use Akima splines, use cubic splines in tabulated transfer functions
Added more ocnstraint infrastructure.
2010-11-12 15:05:45 -08:00
Oliver Hahn
485d6f2a4e
Gadget2 plugin: Important Bugfix. Bug introduced two versions ago that caused the first particle
...
to have zero mass, which caused Gadget to choose a zero time step.
2010-11-10 22:56:46 -08:00
Oliver Hahn
f3364e9877
* Added support for velocity potentials for baryons and CDM
...
* Added new transfer function type 'music' which reads in a file with
density and velocity transfer functions for baryons and CDM.
* Powerspectrum evolution in two-fluid sim appears to scale as
predicted now at the per cent level
* Tool to compute velocity transfer functions is not included as of
now. Let's see..
* Velocity potentials are not yet supported with 2LPT.
2010-11-02 23:35:57 -07:00
Oliver Hahn
dec2181626
Some more 64bit related fixes.
...
Fixed a bug in the Gadget-2 output that was left from an unfinished modification.
2010-10-26 21:42:49 -07:00
Oliver Hahn
c98e030c7c
Updated MUSIC internal array access to 64bit compatibility.
...
Added support for FFTW3 (can be switched on in the Makefile), needed for
64bit compatible FFTs.
2010-10-26 11:37:31 -07:00
Oliver Hahn
332362ce0c
Added buffered writing for Gadget-2 I/O plugin that prevents overflows of STL ofstream operations.
...
Buffer block size can be set by [output]/gadget_blksize, default is 2**20=1048576 so that the
buffer contains 8*2**20 bytes for doubles = 8MB for one write/read operation
2010-10-15 23:11:10 -07:00
Oliver Hahn
62b765002a
fixed a bug in the RAMSES output plugin when grids are not cubical
2010-10-05 11:22:14 -07:00
Oliver Hahn
1dee46e390
Fixed a bug that caused k-mode averaging not to be used when it should be in some cases.
...
Disabled staggered particles for Gadget temporarily.
2010-10-01 15:18:06 -07:00
Oliver Hahn
275753ede2
Many minor changes
...
* Moved hybrid poisson solver from convolution_kernel.cc to poisson.cc
* Fine-tuned deconvolution schemes for various running modes
* Added lots of doxygen documentation
2010-08-31 21:59:31 -07:00
Oliver Hahn
af02e4ef4e
Added new fine-grid hybrid poisson solver. This can be turned on with '[poisson]/fft_fine=yes'. Needs more testing though.
...
Also added a cut-off for baryons when using CAMB tabulated transfer functions to reduce ringing.
2010-08-21 19:43:21 -07:00