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
0f5cbaa8a8
Added tag beta_candidate for changeset 36271198aaf9
2010-10-26 21:54:18 -07:00
Oliver Hahn
6735894892
More 64bit fixes and some more OpenMP parallelizations
2010-10-26 21:54:15 -07:00
Oliver Hahn
32e3ef185f
Added tag beta_candidate for changeset 1e095fb6c31b
2010-10-26 21:44:19 -07:00
Oliver Hahn
85ab2aa1ce
Hybrid poisson solver is now default
2010-10-26 21:44:17 -07:00
Oliver Hahn
032031c078
Added tag beta_candidate for changeset ed36154e55dd
2010-10-26 21:43:02 -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
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
eff7846ec9
Changed to 'unsigned long long' in exact shotnoise sum to avoid overflow.
2010-10-01 19:20:49 -07:00
Oliver Hahn
a25b9ffe3c
Added explicit inclusion of <cstdarg> in log.hh, otherwise got compile errors on some systems
2010-10-01 15:26:09 -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
db164ddc83
Added tag beta_candidate for changeset 4bf26f93896b
2010-09-30 16:25:30 -07:00
Oliver Hahn
a07e5e59ba
Added some more transparent infrastructure for default parameter values
2010-09-30 16:19:10 -07:00
Oliver Hahn
bf3a3ee104
Added tag beta_candidate for changeset ef17480b490a
2010-09-30 16:06:03 -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
d94763e0d3
Added tag beta_candidate for changeset 960b39975d06
2010-09-29 15:42:41 -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
e1af86f9cc
Added OpenMP parallelization for exact T(r=0) calculation.
2010-09-28 16:05:50 -07:00
Oliver Hahn
a740e65a20
Added tag beta_candidate for changeset df8d9bf3d822
2010-09-28 16:05:15 -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
03d6ca8a71
Added tag beta_candidate for changeset 2f262781b4fb
2010-09-28 15:58:03 -07:00
Oliver Hahn
9ac31e4ffd
Added tag beta_candidate for changeset 6219a358469e
2010-09-28 15:51:17 -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
3ba59a04f9
Added tag beta_candidate for changeset c142f8da13da
2010-09-27 15:06:21 -07:00
Oliver Hahn
e237cd80c3
Deconvolution of transfer function is now switchable
...
with [setup]/deconvolve=[yes/no]
2010-09-24 16:16:31 -07:00
Oliver Hahn
b03b1c721c
Added tag beta_candidate for changeset 78bd61d6b7b2
2010-09-24 16:14:18 -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
aa9e88310e
Fixed a normalization issue for they hybrid Poisson solver.
...
Option 'avg_fine' now defaults to 'true'.
2010-09-08 01:06:18 -07:00
Oliver Hahn
3ec25ad49a
Added tag beta_candidate for changeset 9040d899b80f
2010-09-08 01:04:29 -07:00
Oliver Hahn
f973bf32ed
Added tag beta_candidate for changeset 6f5b12c825c8
2010-08-31 21:59:44 -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
Oliver Hahn
da227a0a04
Added tag beta_candidate for changeset f621982eb4eb
2010-08-15 18:18:00 -07:00
Oliver Hahn
b7083b0217
Important: Fixed a missing logarithm
...
Last changeset used a wrong transfer function as a logarithm
had been removed and not put back. This is fixed now!
2010-08-15 18:16:08 -07:00
Oliver Hahn
72739b8ea7
Updated Makefile for new logfile infrastructure
2010-08-14 20:01:43 -07:00
Oliver Hahn
c815c2db1f
Added tag beta_candidate for changeset 135e5f558833
2010-08-14 20:01:26 -07:00
Oliver Hahn
965bc5c1c8
Added tag beta_candidate for changeset eabada1ee28b
2010-08-14 18:20:50 -07:00
Oliver Hahn
28701c3d8e
Added some logging infra structure.
2010-08-14 18:19:48 -07:00
Oliver Hahn
b5d5cd0d90
Added tag beta_candidate for changeset 9473d7ef2abd
2010-08-14 18:01:54 -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
Oliver Hahn
ae0d833155
Fixed a bug that caused grid shifts to be wrong in some cases when using
...
levelmin_TF>levelmin
2010-08-14 16:51:17 -07:00
hahn
506b882893
Fixed some stuff that caused intel c++ compiler warnings
2010-08-11 05:35:29 +02:00
Oliver Hahn
5d3b66e96a
Added new interpolation scheme for long-range component.
...
The LR part is now interpolated directly to fine grids from coarse grids not
involving intermediate grids. This improves accuracy for >=2 levels.
2010-08-10 15:51:49 -07:00
Oliver Hahn
b5bd5fe15a
Fixed a bug when reading an external seed file
2010-08-10 15:50:02 -07:00
Oliver Hahn
67d3038e25
Fixed a bug that caused an error message when reading random values from
...
file and only using one unigrid level.
2010-08-06 10:23:45 -07:00
Oliver Hahn
587bf37802
Added tag beta_candidate for changeset 696daefc860e
2010-08-05 16:25:42 -07:00