JFEngels
48e4f7004d
we now enforce the blocking factor...this time for real
2013-05-10 12:24:50 +02:00
JFEngels
024e81dc8d
some merge...
2013-04-23 18:02:16 +02:00
Oliver Hahn
ad67f23f30
reduced the padding around ellipsoids when not using RAMSES output
2013-04-21 11:55:36 +02:00
Oliver Hahn
ead40fd5b6
added support for ellipsoid region files that include velocity information
...
to correctly deal with the initial displacement of the particles from their
Lagrangian coordinates
2013-04-19 16:45:48 +02:00
JFEngels
2c8023c1f2
Most of the stuff needed to produce initial conditions for Nyx.
...
Still a little fragile...
2013-02-07 15:02:07 +01:00
Oliver Hahn
549c1dbd47
fixed wrong sign in new velocity frame calculation
2013-01-10 11:33:32 +01:00
Oliver Hahn
2705f78e69
fixed correction for previous grid shift when centering velocities via explicit shift
2013-01-08 17:14:02 +01:00
Oliver Hahn
2f4794e408
initial support for explicit shift compensation with center_vel, center_vel_posfinal...
2013-01-08 17:11:32 +01:00
Oliver Hahn
83e381f5cb
some code cleanup
2012-12-17 10:33:18 +01:00
Oliver Hahn
7cbfb01605
updates and fixes to refinement region infrastructure
...
support for gadget preliminary tests ok
2012-12-11 02:11:04 +01:00
Oliver Hahn
966e670a9d
fixed a bug in the grid hierarchy generation that could lead to inconsistent grid placement when using also the option 'force_equal_extent=yes'
2012-11-23 16:41:19 +01:00
Oliver Hahn
7a72cd4f0b
updated code version to 1.2b
2012-11-20 16:07:50 +01:00
Oliver Hahn
cd902e787b
added CIC deconvolution step to hybrid Poisson solver. This is necessary when glass or resampling is employed via CIC particle interpolation
2012-11-20 16:01:45 +01:00
Oliver Hahn
d9fbe83bbc
added new parameter [setup]/ref_dims=NX,NY,NZ to explicitly specify the resoluti
...
on of the high resolution region.
added output of comoving size of high-res region to screen
added output of densest point also when using the generic output plugin
2012-09-28 15:17:25 +02:00
Oliver Hahn
5589695ef6
updated code version to 1.1b, first version with significantly reduced memory usage for very deep zooms.
2012-09-07 16:53:43 +02:00
Oliver Hahn
3802e03da2
merge
2012-07-11 17:18:20 -07:00
Oliver Hahn
f10f8f9a6c
version freeze for Santa Cruz Code Comparison Project
2012-07-11 17:15:45 -07:00
Oliver Hahn
0754897c60
added possibility to restart using previously generated white noise fields
...
via option [random]/restart = yes/no
if set to yes, it will try to load the data, but all files have to exist,
and the refinement hierarchy cannot be changed.
2012-03-28 14:55:16 -07:00
Oliver Hahn
19cc12a603
moved all omp function calls into #ifndef SINGLETHREAD_FFTW so that the code also compiles on non-openmp compilers
2012-03-28 12:38:31 -07:00
Oliver Hahn
924092d96c
Fixed a missing velocity unit factor when using "linger++" transfer functions, which
...
have different velocity units than the other transfer functions, and 2LPT.
2011-09-19 18:49:43 -07:00
Oliver Hahn
826071634a
updated version to 1.01
2011-08-29 16:12:40 -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
e5f663ab0a
Added a missing pragma that leads to errors when compiling for single-threaded FFTW3
2011-07-21 10:43:00 -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
0f012ee9f9
a 64bit compliance fix in convolution_kernel
...
fixed a bug introduced in this branch in random.cc
2011-06-02 15:10:05 -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
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
088fb32c23
Fixed some compiler warnings (and possibly errors with GCC 4.3)
2011-03-16 15:33:03 -07: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
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
ccf7c3dab1
changed version number to 0.9b, this is the beta version.
...
no more features will be added.
2010-11-16 21:01:15 -08:00
Oliver Hahn
3ad0ac6a10
Fixed a memory leak.
...
Slightly better memory conservation.
Fixed a bug introduced in the previous revision that led to code termination.
2010-11-16 20:56:08 -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
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
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
6735894892
More 64bit fixes and some more OpenMP parallelizations
2010-10-26 21:54:15 -07:00
Oliver Hahn
85ab2aa1ce
Hybrid poisson solver is now default
2010-10-26 21:44:17 -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
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
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
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
28701c3d8e
Added some logging infra structure.
2010-08-14 18:19:48 -07:00