ParticleInCell.jl:Julia中的“单元格中的粒子”代码 源码

ParticleInCell.jl Julia中的粒子在单元格代码 正在进行中,您可以在文档中找到第一个测试。 我正在研究这个Julia代码,以比较性能与Fortran代码的性能,该代码可以使用相同的参数和相同的数值方法来求解相同的方程式。 你可以在找到它。 这是2005年编写的旧代码,没有经过优化,但是使用204800个粒子需要6秒钟。 我在这里注意到了这段Julia代码的时间以及我为加快处理速度所做的工作。 323秒:第一版 303秒:使用julia -O3 --check-bounds = no 302秒:更改粒子的位置和速度数组的形状 156秒:将ex,ey和bz重新排列在同一数组eb(3,nx,ny)中 176秒:将数组ebj放入fdtd类型 140秒:将粒子数据放在一个数组中,粒子push_v现在非常快。 077秒:通过调用fortran子例程替换julia插值函数。