Logo PTI
Polish Information Processing Society
Logo FedCSIS

Annals of Computer Science and Information Systems, Volume 17

Communication Papers of the 2018 Federated Conference on Computer Science and Information Systems

Parallelizing the code of the Fokker-Planck equation solution by stochastic approach in Julia programming language

DOI: http://dx.doi.org/10.15439/2018F253

Citation: Communication Papers of the 2018 Federated Conference on Computer Science and Information Systems, M. Ganzha, L. Maciaszek, M. Paprzycki (eds). ACSIS, Vol. 17, pages 115120 ()

Full text

Abstract. Presenting a reliable physical simulation requires very often use of the supercomputers and models run for many days or weeks. The numerical computing is divided into two groups. One uses highly efficient low-level languages like Fortran, C, and C++. The second applies high-level languages like Python or Matlab, being usually quite slow in HPC applications. This paper presents the application of the relatively new programming language Julia, advertised as the as``a high-level, high-performance dynamic programming language for numerical computing''. We employ Julia is to solve the Fokker-Planck equation by the stochastic approach with the use of the corresponding set of ordinary differential equations. We propose the method of parallelizing the algorithm with use of the distributed arrays. We test the speedup and efficiency of the given code with use of the cluster set at the {\'S}wierk Computing Centre and show that Julia is capable of achieving a good performance.


  1. Rei, L., Carvalho, S., Alves, M., Brito, J., A look at dynamic languages, Tech. report, 2007, Faculty of Engineering University of Porto.
  2. The Julia Language https://julialang.org/
  3. Bezanson J. et. al., Julia: A Fresh Approach to Numerical Computing, SIAM Review, vol. 59, 1 (2017) 65-98
  4. GitHub https://github.com/JuliaLang/julia
  5. The Jupyter Project, http://jupyter.org/
  6. J. Bezanson, Abstraction in Technical Computing, Ph.D. thesis, Massachusetts Institute of Technology, MA, 2015.
  7. Bezanson J., S. Karpinski, V. B. Shah, and A. Edelman, Julia: A Fast Dynamic Language for Technical Computing, preprint, https://arxiv.org/abs/1209.5145 [cs.PL], 2012.
  8. Lattner C., Adve V., LLVM: A compilation framework for lifelong program analysis & transformation. In: Proceedings of the international symposium on Code generation and optimization: feedback-directed and runtime optimization. IEEE Computer Society, 2004, 75
  9. Distributed Arrays Packadge https://github.com/JuliaParallel/DistributedArrays.jl
  10. Risken,H., The Fokker–Planck Equation Method of Solution and Applications, SpringerVerlag,Berlin,Heidelberg, 1989
  11. Moraal H , Cosmic-Ray Modulation Equations, Space Science Reviews,2013, vol.176, 299–319 http://dx.doi.org/10.1007/s11214-011-9819-3
  12. Parker E. The passage of energetic charged particles through interplanetary space, Planetary and Space Science, 1965, vol.13, 9-49
  13. Wawrzynczak A., Alania M.V., Numerical Solution of the Time and Rigidity Dependent Three Dimensional Second Order Partial Differential Equation, Lecture Notes in Computer Science, 2010, vol.6067, pp. 105–114, http://dx.doi.org/10.1007/9783−642−14390−81_2
  14. Wawrzynczak A, Modzelewska R and Gil A, Stochastic approach to the numerical solution of the non-stationary Parker’s transport equation, Journal of Physics: Conference Series, 2015, vol. 574, 012078,http://dx.doi.org/10.1088/1742-6596/574/1/012078
  15. Wawrzynczak A, Modzelewska R and Gil A, The algorithms for forward and backward solution of the Fokker-Planck equation in the heliospheric transport of cosmic rays, Lecture Notes in Computer Science,2018, vol. 10777, 14–23, http://dx.doi.org/10.1007/978-3-319-78024-5-2
  16. Gardiner C.W., Handbook of stochastic methods. For physics, chemistry and the natural sciences, Springer Series in Synergetics, 2009
  17. Kloeden P E, Platen E, Schurz H , Numerical solution of SDE through computer experiments, Springer–Verlag Berlin Heidelberg, 1992
  18. Alania M. V., Stochastic Variations of Galactic Cosmic Rays, Acta Physica Pol. B, 2002, vol.33(4), 1149–1166
  19. Cluster Manager Packadge https://github.com/JuliaParallel/ClusterManagers.jl
  20. Amdahl, G.M., Validity of the single-processor approach to achieving large scale computing capabilities, in: Proc. Am. Federation of Information Processing Societies Conf., AFIPS Press, (1967) 483-485.