代碼initial.m

function initial(prm, hdiag)global q mass rho0global npt % total number perticlesglobal slx % system length xglobal bx0 by0global ex ey ezglobal bx by bzglobal ajx ajy ajzglobal rhoglobal vx vy vzglobal xglobal nxp1 nxp2 % nx+1, nx+2global X1 X2 X3global cs tcs % cv^2, 2*cv^2global ifdiag%---diagnostics---%global field % field data for wk plotglobal kspec % k-spectrumglobal eng % energyslx = prm.nx;% nx=128，system length xnpt = sum(prm.np(1:prm.ns));%總粒子數(shù)nxp1 = prm.nx+1;nxp2 = prm.nx+2;X1 = 1:prm.nx;X2 = 2:(prm.nx+1);X3 = 3:(prm.nx+2);cs = prm.cv*prm.cv;%光速的平方tcs = 2.0*cs;end

總結(jié)

以上是生活随笔為你收集整理的[一维粒子模拟 version3.6]实现initial.m函数的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。