README.md

HeLaZ (Hermite-Laguerre Z-pinch solver, 2021)

To compile it check INSTALLATION.txt

How to run it

1. Be sure to have correct paths in local/dirs.inc for the different libraries
2. Compile from HeLaZ/ using make
3. To run the code, use HeLaZ/wk/local_run.m and run it to set up the parameters and the results folder
4. Then go to the results folder and launch HeLaZ using mpirun -np num_procs ./../../../bin/helaz num_p num_kr
5. You can obtain various plots and gifs using HeLaZ/wk/analysis_2D.m once the simulation is done. To select the correct output file, run parameters*.m with the corresponding simulation parameters and then run analysis_2D.m (everything with matlab from wk/)

// Comment : For some collision operators (Sugama and Full Coulomb) you have to run COSOlver from B.J.Frei first in order to generate the required matrices in HeLaZ/iCa folder.

Changelog

-2. MPI parallel version

<<<<<<< HEAD 2.5 //GK COSOlver collision operator//

-2.5 //GK COSOlver collision operator//

d7a9164c

-2.4 2D cartesian parallel (along p and kr)

-2.3 GK Dougherty operator

-2.2 Allow restart with different P,J values (results are not concluents)

-2.1 First compilable parallel version (1D parallel along kr)

-1. Implementation of the non linear Poisson brackets term

-1.4 Quantitative study with stationary average particle flux \Gamma_\infty

-1.3 Linear analysis showed that a certain amount of PJ are recquired to trigger mode

-1.2 Zonal flows are observed in a similar way to Ricci Rogers 2006 with GS2

-1.1 Qualitative test : find similar turbulences as Hasegawa Wakatani system with few moments

-1.1 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system

-1.1 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system

-1.0 FFTW3 has been used to treat the convolution as a product and discrete fourier transform

-0. Write MOLI matlab solver in Fortran using Monli1D as starting point

-0.6 Benchmarks now include Dougherty, Lenard-Bernstein and Full Coulomb collision operators

-0.5 Load COSOlver matrices

-0.4 Benchmark with MOLI matlab results for Z-pinch (cf. kz_linear script)

-0.3 RK4 time solver

-0.2 implement moment hierarchy linear terms

-0.1 implement linear Poisson equation in fourier space

-0.0 go from 1D space to 2D fourier and from Hermite basis to Hermite-Laguerre basis

Roadmap

-2. MPI parallel version

-2.6 GPU accelerated version

-2.7 GK Full Coulomb collision operator

-3. GK 3D version, kr,kz,kpar for linear device

-4. DK 3D version, kr,kz,kpar for linear device

-5. DK+GK 3D version, kr,kz,kpar for linear device

-6. 3D version with curvature