void applyTorque() {
		if (torque == null)
			return;
		if (Math.abs(torque.getForce()) < Particle.ZERO)
			return;
		computeCenter();
		int[] exclusion = torque.getExclusion();
		double k, sinx, cosx;
		outerloop2: for (Atom a : atoms) {
			if (exclusion != null && exclusion.length > 0) {
				for (int i : exclusion) {
					if (i == a.index)
						continue outerloop2;
				}
			}
			k = (xCenter - a.rx) / (a.ry - yCenter);
			cosx = (a.ry > yCenter ? -torque.getForce() : torque.getForce()) / Math.sqrt(1.0 + k * k);
			sinx = k * cosx;
			a.fx += cosx;
			a.fy += sinx;
		}
	}