# ALYA ## Summary Version 1.0 ## Purpose of Benchmark The Alya System is a Computational Mechanics code capable of solving different physics, each one with its own modelization characteristics, in a coupled way. Among the problems it solves are: convection-diffusion reactions, incompressible flows, compressible flows, turbulence, bi-phasic flows and free surface, excitable media, acoustics, thermal flow, quantum mechanics (DFT) and solid mechanics (large strain). ALYA is written in Fortran 90/95 and parallelized using MPI and OpenMP. * Web site: https://www.bsc.es/computer-applications/alya-system * Code download: https://repository.prace-ri.eu/ueabs/ALYA/2.1/Alya.tar.gz * Test Case A: https://repository.prace-ri.eu/ueabs/ALYA/2.1/TestCaseA.tar.gz * Test Case B: https://repository.prace-ri.eu/ueabs/ALYA/2.1/TestCaseB.tar.gz ## Mechanics of Building Benchmark Alya builds the makefile from the compilation options defined in config.in. In order to build ALYA (Alya.x), please follow these steps: Goto to directory: Executables/unix ``` cd Executables/unix ``` Edit config.in (some default config.in files can be found in directory configure.in): * Select your own MPI wrappers and paths * Select size of integers. Default is 4 bytes, For 8 bytes, select -DI8 * Choose your metis version, metis-4.0 or metis-5.1.0_i8 for 8-bytes integers Configure Alya: ./configure -x nastin parall Compile metis: make metis4 or make metis5 ## Mechanics of Running Benchmark ### Datasets The parameters used in the datasets try to represent at best typical industrial runs in order to obtain representative speedups. For example, the iterative solvers are never converged to machine accuracy, but only as a percentage of the initial residual. The different datasets are: SPHERE_16.7M ... 16.7M sphere mesh SPHERE_132M .... 132M sphere mesh ### How to execute Alya with a given dataset In order to run ALYA, you need at least the following input files per execution: X.dom.dat X.ker.dat X.nsi.dat X.dat In our case X=sphere To execute a simulation, you must be inside the input directory and you should submit a job like: mpirun Alya.x sphere How to measure the speedup -------------------------- There are many ways to compute the scalability of Nastin module. 1. **For the complete cycle including: element assembly + boundary assembly + subgrid scale assembly + solvers, etc.** > In *.nsi.cvg file, column "30. Elapsed CPU time" 2. **For single kernels: element assembly, boundary assembly, subgrid scale assembly, solvers**. Single kernels. Here, average and maximum times are indicated in *.nsi.cvg at each iteration of each time step: > Element assembly: 19. Ass. ave cpu time 20. Ass. max cpu time > > Boundary assembly: 33. Bou. ave cpu time 34. Bou. max cpu time > > Subgrid scale assembly: 31. SGS ave cpu time 32. SGS max cpu time > > Iterative solvers: 21. Sol. ave cpu time 22. Sol. max cpu time > > Note that in the case of using Runge-Kutta time integration (the case > of the sphere), the element and boundary assembly times are this of > the last assembly of current time step (out of three for third order). 3. **Using overall times**. > At the end of *.log file, total timings are shown for all modules. In > this case we use the first value of the NASTIN MODULE. Contact ------- If you have any question regarding the runs, please feel free to contact Guillaume Houzeaux: guillaume.houzeaux@bsc.es