<< WebHome

LATfield2 + Gevolution tutorial

For files, see https://bitbucket.org/m-ahsan-nazer/cosmotorun19 .

-- BoudRoukema - 15 Jul 2019

Using the tutorial examples

Step 0 (optional)

For those who have never coded in c++ a basic intro using three examples. Use cppBasics.cpp to gain enough familiarity with the C++ language to be able to follow the gevolutions session discussions. In the main function uncomment the relevant lines to see output of the example.

a) The header replaces <stdio.h> for input/output b) c++ template. A simple template function that can receive two integers or floats and returns the maximum value of the correct data type. (Gevolution uses function and class templates to define fields such as gravitational potential, tensor fields etc on lattices, the intention is that the workshop participant recognises template notation.) c) A simple c++ class that has two member functions. Classes are a foreign concept for those who have only coded in C or Fortran 90. d) Over loaded operators. (Required if the participant wants to decode physics equations in gevolution. Only a superficial understanding required)

Step 1 (Install required software)

Install gevolution, LATfield2, hdf5 and fftw3.

Step 2 (Run potForce.cpp, dft_forwards.cpp)

Purpose: Be able to use "Lattice", "Site", "Field", "PlanFFT" classes. (Gevolution takes Fourier transforms (computes discrete Fourier transforms) when setting up the simulation initial conditions. It uses "Lattice", "Site", "Field", "PlanFFT" classes. These classes are inherited from LATfield2 library. To be able to decipher equations coded in gevolution or modify gevolution a working knowledge of these classes are required).

a) A simple example that differentiates phi=1/r (Newtonian potential for a point source placed at the centre of the lattice, without units) to obtain Force=1/r^2 (without units). Requires defining "Lattice", "Site", "Fields", and taking finite differences the LATfield2 way to compute the force from potential.

b) Use LATfield2 to compute FFT of cos(2pi/5 x)+ 3*cos(2pi/10 x) (a made up example) on the lattice. Plot the FFT coefficients using python and identity the two corresponding frequency peaks.

Step 3 (Run gevolution with default settings)

Purpose: Run a simulation using gevolution. (This step can be performed without going through steps 1 to 2)

Step 4 (Decipher gevolution code to recover the equations in arXiv 1604.06065)

Purpose: Recover underlying equations from gevolution code. (In ic_basic.hpp the initial conditions for the simulation are set. The initial particle displacements are obtained from the gradient of a scalar field. This equation is provided in arXiv 1604.06065 i.e the main gevolution paper. The workshop participant should be able to identity this equation from the ic_basic.hpp and perhaps modify it.)
Edit | Attach | Print version | History: r4 < r3 < r2 < r1 | Backlinks | View wiki text | Edit WikiText | More topic actions...
Topic revision: r2 - 15 Jul 2019, AhsanNazer
This site is powered by FoswikiCopyright © CC-BY-SA by the contributing authors. All material on this collaboration platform is copyrighted under CC-BY-SA by the contributing authors unless otherwise noted.
Ideas, requests, problems regarding Foswiki? Send feedback