Overdamped#

Forward Euler#

// position and velocity
xt::xtensor<double,2> u = xt::zeros<double>({nnode, ndim});
xt::xtensor<double,2> v = xt::zeros<double>({nnode, ndim});

// time increments
for ( ... )
{
  // new displacement
  xt::noalias(u) = u + dt * v;

  // new velocity based on residual force
  ...
}