Пример #1
0
  void Particle::integrate(fseconds duration, int type) { // Defaults to Euler

    // don't move objects that have "infinite mass."
    if (invMass <= 0.0f) {
      return;
    }

	// Universal properties
	//lastPos = pos; // Only necessary for position Verlet
	vec3 resultantAcc = acc + (invMass * accumulatedForces);

    auto dt = duration.count();

	if ( type == 1 ) { // Euler
		// update position using Euler integration
		pos = pos + dt * vel;

		// update velocity using Euler integration
		vel = vel + dt * resultantAcc;

		// incorporate damping
		vel = vel * std::pow(damping, dt);
	} else if ( type == 2 ) { // RK
		// Pop RK4 in here at some point because... why not
	} else { // Velocity Verlet
		pos = pos + ( vel * dt ) + ( 0.5f * resultantAcc * dt * dt );

		// Assumes acc doesn't change between steps - still need to solve for acc+1 to be accurate if a changes
		vel = vel + ( 0.5f * 2 * resultantAcc ) * dt;	
	}   

    clearForces();

  }
Пример #2
0
  void Particle::integrate(fseconds duration) {

    // don't move objects that have "infinite mass."
    if (invMass <= 0.0f) {
      return;
    }

    //update velocity using Verlet
    vel = vel + 1.0f/2.0f*acc*duration.count();
    // update position using Velocity Verlet
    pos = pos + vel * duration.count();
    vec3 resultantAcc  = (pos+duration.count())-(2.0f*pos+pos)-(duration.count());
    // update velocity using Velocity Verlet
    vel = vel + 1.0f/2.0f*resultantAcc*duration.count();
    // incorporate damping
    vel = vel * damping;

    clearForces();

  }
Пример #3
0
  void Particle::integrate(fseconds duration) {

    // don't move objects that have "infinite mass."
    if (invMass <= 0.0f) {
      return;
    }

    // update position using Euler integration
    pos = pos + duration.count() * vel;

    vec3 resultantAcc = acc;
    resultantAcc = resultantAcc + (invMass * accumulatedForces);

    // update velocity using Euler integration
    vel = vel + duration.count() * resultantAcc;

    // incorporate damping
    vel = vel * damping;

    clearForces();

  }