About this role
<p>Objectives<br>1. An optimal implementation of the calculation of the reaction rate convolution integral.<br>• An optimal data structure for the cross sections.<br>• Optimal way of calculating the reaction rate integral using point-wise data where both<br>functions have different grid points. Is it worth creating a “shared” grid to be able to use the<br>trapezoid rule? Should this common grid be stored or regenerateed each time?<br>• Analyse the potential of parallellisation taking into account I/O limitations (each process<br>reading from the same file(s)).<br>2. A strategy to “coarsen” the cross section and neutron flux grid up-front to reduce computation<br>time and memory/storage usage.<br>• Can we, in general or based on a typical application, coarsen both the cross section grid and the<br>energy grid for neutron spectrum binning without loosing accuracy?<br>3. Implement the on-the-fly Doppler broadening<br>• Implement the basic OTFDB approach by calculating the integral above.<br>• Analyze if and how calculations can be cached. Analyze the trade-off between CPU needs, RAM<br>use, I/O burden and storage needs.<br>We expect proof-of-concept implementations that we can integrate in our ALEPH2 code. We will<br>provide the student with the necessary data sets and typical use case examples</p>