[C++调试笔记]/* 求解-0.5dt时刻速度 */
生活随笔
收集整理的這篇文章主要介紹了
[C++调试笔记]/* 求解-0.5dt时刻速度 */
小編覺得挺不錯的,現(xiàn)在分享給大家,幫大家做個(gè)參考.
#include "define.h"
#include <stdio.h>
/* 求解-0.5dt時(shí)刻速度 */
void setv(particle* ptr, int numb, double* fiel, double q, double m)
{double vz, az, z, ez;int ip, i;double l1, l2;for (i = 0; i < numb; i++){z = ptr[i].z;vz = ptr[i].vz;ip = (int)(z / dz_plas);l2 = z / dz_plas - ip;l1 = 1.0 - l2;ez = fiel[ip] * l1 + fiel[ip + 1] * l2;az = q * ez / m;vz -= (0.5 * az * dt);ptr[i].vz = vz;}}
總結(jié)
以上是生活随笔為你收集整理的[C++调试笔记]/* 求解-0.5dt时刻速度 */的全部內(nèi)容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 程序员回家过个年,真心不要容易
- 下一篇: [C++调试笔记]推动粒子move.cp