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
Oliver Hahn
8e1ed3d284
Added option '[setup]/center_velocities' that centers the box rest-frame to
...
the mean velocity in the finest refinement grid.
2010-08-05 16:22:40 -07:00