From 64d6cd8a7ed46be4a34b4818705186f8e0998ab1 Mon Sep 17 00:00:00 2001 From: Jakob Odersky Date: Wed, 11 May 2011 14:59:00 +0000 Subject: Rajoute getHorizontalAt() pour tous les elements, en vue des faisceaux --- src/main/Accelerator.cc | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) (limited to 'src/main/Accelerator.cc') diff --git a/src/main/Accelerator.cc b/src/main/Accelerator.cc index 37f8dbb..ca5de8a 100644 --- a/src/main/Accelerator.cc +++ b/src/main/Accelerator.cc @@ -144,13 +144,15 @@ void Accelerator::step(double dt) { for (ParticleIterator i = particleCollec.begin(); i != particleCollec.end(); ++i) { Particle& particle = **i; + particle.setForce(Vector3D::Null); + particle.applyMagneticForce(particle.getElement()->magneticFieldAt(particle.getPosition()), dt); Vector3D a = particle.getForce() / (particle.getGamma() * particle.getMass()); particle.setVelocity(particle.getVelocity() + a * dt); particle.translate(particle.getVelocity() * dt); - particle.setForce(Vector3D::Null); + } -- cgit v1.2.3