banner

The site's hierarchy



Home page :: Private area :: Site Map :: Credits
Software / Degraf

Navigation

Mar 2009
MTWTFSS
2324252627281
2345678
9101112131415
16171819202122
23242526272829
303112345
News
  • Les premiers zoom de galaxies sont accessibles
  • Le halo FOF 6133 de la boite Horizon L et le halo FOF 544 de Horizon S ont été resimulés avec plusieurs techniques de zoom. Il est important que chaque "zoomer" valide sa méthode, avant de lancer une campagne de zoom sur un plus grand cataloque.
  • Les données de la simulation Mare Nostrum sont disponibles
  • 34 snapshots jusqu’à z=4 sont accessibles sur horizon3 et sur le serveur de fichiers de l’IDRIS à la collaboration Horizon.
  • Méso Machine HPC1 opérationnelle
  • Depuis le 23 Octobre 2005, la Méso machine du site horizon est operationnelle. Elle correspond à 3 quadriprocesseurs avec chacun 64 Giga de RAM reliés par infiniband, ainsi qu’un access conséquent (sur une base de projet dédié) au reste de la ferme). Son acces est ouvert a toute personne de la collaboration ayant acces à la minigrille et qui en fait la demande a admin-minigrille
  • http://
  • Workshop Horizon le 14 et 15 novembre 2005
  • Il aura lieu à Paris les 14 et 15 novembre 2005 (prévoir une nuit sur place). L’enregistrement est ouvert dans la rubrique "meeting!"
  • http://

Degraf


by Pichon Christophe (Tuesday 18 April 2006)

Degraf 0.1

Degraf is a parallel MPI code to degrade in resolution large cosmological initial conditions on a cluster without requiring shared memory, while possibly shifting the box for the purpose of resimulations. The real Fourier transforms are carried in place using fftw, while minimizing the amount of used memory (at the expense of performance), in the spirit of grafic-1. The output of the result(s) is also carried in parallel.

It runs interactively using the command:

mpirun -np [number of processors] degraf [namelistfile]

where namelistfile is a fortran namelist parameter file. The structure of the parameter file is the following:


&parameters

nfiles = 1

inputfiles = 'ic_deltab'

degramax = 2

recurse = .false.

shift = 0 0 0
/

* nfiles is the number of input files to be processed

* inputfiles is the list of input file names, in simple quotes, separated by white spaces

* degramax is a power of 2 (greater or equal to 2) indicating by how much the resolution of the input box(es) is degraded

* if recurse is set to .false., the output resolution is uniquely set by degramax. If set to .true., then all resolutions corresponding to degradation factors from to 2 up to degramax (in powers of 2) are computed

* if shift is present the whole hierarchy will be shifted according to shift (in relative units of the boxsize).

* The output file names are constructed from the input file names, adding at the end of the file names the output linear size of the box in the x direction (note that the code should work also on non-cubic boxes)

For velocity fields (or in general for fields where the grid is offset by a fraction of a mesh length), care is taken during the degradation to first shift the grid points back to the original (zero offset) position at full resolution. After degradation, the shift is reapplied, with a value equal to the same offset fraction, but computed in units of the new (bigger) mesh size.