1
0
Fork 0
mirror of https://github.com/cosmo-sims/MUSIC.git synced 2024-09-19 17:03:46 +02:00
Commit graph

25 commits

Author SHA1 Message Date
Oliver Hahn
d6da7b16fc More 64 bit compliance changes.
New convergence criterium for multi-grid (No longer reduction of initial residual by X, but rel. error of X).
2011-06-01 18:05:12 -07:00
Oliver Hahn
aac5720a41 more fixes to ensure 64bit compliance throughout the code 2011-06-01 11:25:53 -07:00
Oliver Hahn
b7756e22a0 full64bitcompliance 2011-06-01 11:01:24 -07:00
Oliver Hahn
61a460b3b6 * Fixed a compiler error with hybrid poisson solve when using single precision
* Some minor code cleanup
2011-02-22 14:17:34 -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
3f332540dc Disabled that FFTW3 is used by default in the Makefile 2010-10-26 22:48:34 -07:00
Oliver Hahn
61caef1575 2LPT now reuses velocity data to compute displacements in pure DM case (instead of recomputing everything).
Random number generator now creates 'dummy seeds' if seeds are not specified
2010-10-26 22:46:33 -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
c09f2adbe4 Random numbers are now only generated once. Can choose between memory caching and disk caching of random numbers
in conf file (option '[random]/disk_cache={yes|no}').
2010-10-04 18:57:24 -07:00
Oliver Hahn
7eeba48408 Changed T(r=0) determination to 3D adaptive Gauss-Kronrod integration rather than spherical approximation.
Seems very accurate in very first tests.
2010-10-04 16:34:01 -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
745f9dc543 Added a strictly conservative cubic spline interpolation.
This reduces errors further below 1e-4 in the entire interior.
2010-09-30 16:05:59 -07:00
Oliver Hahn
ad709d7400 Various minor fixes.
* Fixed a bug in the coarse grid density correction that lead to a shift in the correction in some rare cases
* Fixed a bug in restricting fine grid noise when the double size fine grid (for FFT) is actually larger than the domain.
* Added new log file routines that use printf syntax
* Added log file output to most parts of the code
* Removed the stats file output, it is now appended to the log file
2010-09-29 15:42:07 -07:00
Oliver Hahn
e424f3c205 Some fixes
* Fixed a normalisation bug in the hybrid Poisson solver
* Added calculation of the exact T(r=0) by summing modes rather than by analytical integration
* Fixed a bug in the deconvolution that caused the Nyquist mode to be excluded in the deconv.
* Fixed a bug in the random number generator when using single thereaded FFTW.
2010-09-28 15:58:41 -07:00
Oliver Hahn
876740e8d9 Fixed a bug in the random number generator that could lead to segfaults in some situations. 2010-09-27 15:06:44 -07:00
Oliver Hahn
fcc53f6edf Massive code cleanup
New conservative transfer function scheme.
New random number generator that preserves coarse Fourier modes instead of applying Hoffman-Ribak
Some minor bugs fixed.
2010-09-24 16:14:09 -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
785601d4ae Speeded up transfer function calculation significantly.
Added an approximate log10 algorithm.
Use spherical symmetry of problem to avoid computing duplicate values.
2010-08-14 16:52:25 -07:00
hahn
506b882893 Fixed some stuff that caused intel c++ compiler warnings 2010-08-11 05:35:29 +02:00
Oliver Hahn
b5bd5fe15a Fixed a bug when reading an external seed file 2010-08-10 15:50:02 -07:00
Oliver Hahn
fbebfabdfb Fixed a minus in 2LPT. Fixed an issue with a deconvolution filter.
Code needs cleanup! Added (undocumented) k-space 2LPT.
2010-07-27 15:21:50 -07:00
Oliver Hahn
002763346f major mod: added subgrid deconvolution for DM
* The DM transfer function is now deconvolved with the cell assignment
  function in order to restore power near the Nyquist wavenumber.
* Added an option to change the random number sample cube size
  [random]/cubesize. This should be set >128 to generate un-
  correlated numbers. True fix pending.
* Removed some unused code
* Increased Buffer size for Gadget-2 output plug-in
2010-07-21 01:10:29 -07:00
Oliver Hahn
c48c73488a Initial commit, beta version release candidate 2010-07-02 11:49:30 -07:00