fn kernel(
    dt: f64,
    nitems: usize,
    poss: &mut Vec<f64>,
    vels: &mut Vec<f64>
) -> bool
Expand description

Core function of the integrator.

This function solves the discrete governing equations. The main purpose is to iteratively update the difference of the angular velocities dv until converged.