void FluidField::update() { get_from_UI ( dens_prev, u_prev, v_prev ); /*u_prev[IX(50,20)] = -100.0f; v_prev[IX(50,20)] = -100.0f; u_prev[IX(50,21)] = 100.0f; v_prev[IX(50,21)] = -100.0f; u_prev[IX(51,21)] = 100.0f; v_prev[IX(51,21)] = 100.0f; u_prev[IX(51,20)] = -100.0f; v_prev[IX(51,20)] = 100.0f; */ densityStep( dens, dens_prev, u, v, diff, dt ); velocityStep( u, v, u_prev, v_prev, visc, dt); }
void leapFrogStep(particle * particles, int n, double h){ positionStep(particles, n, h/2); velocityStep(particles, n, h); positionStep(particles, n, h/2); }