void propagate_trajectories_block(double* variables, double nu, double chi, double xi_rot, double xi_alpha,
  double dt, int num_steps);
void velocity_field(double* defect_config, int n, double* lattice_thetas, double* lattice_phis,
  double rho, double* velocities);