Typo in the DEFAULT_SPH case

src/hydro/Default/hydro.h

@@ -194,12 +194,12 @@ __attribute__((always_inline)) INLINE static void hydro_end_density(
   p->rho_dh = 1.f / (1.f + hydro_dimension_inv * p->h * p->rho_dh * irho);
 
   /* Finish calculation of the velocity curl components */
-  p->density.rot_v[0] *= hydro_dimension_inv * irho;
-  p->density.rot_v[1] *= hydro_dimension_inv * irho;
-  p->density.rot_v[2] *= hydro_dimension_inv * irho;
+  p->density.rot_v[0] *= h_inv_dim_plus_one * irho;
+  p->density.rot_v[1] *= h_inv_dim_plus_one * irho;
+  p->density.rot_v[2] *= h_inv_dim_plus_one * irho;
 
   /* Finish calculation of the velocity divergence */
-  p->density.div_v *= hydro_dimension_inv * irho;
+  p->density.div_v *= h_inv_dim_plus_one * irho;
 }
 
 /**

tests/test27cells.c

@@ -321,7 +321,7 @@ int main(int argc, char *argv[]) {
   message("Hydro implementation: %s", SPH_IMPLEMENTATION);
   message("Smoothing length: h = %f", h * size);
   message("Kernel:               %s", kernel_name);
-  message("Neighbour target: N = %f", h * h * h * kernel_norm);
+  message("Neighbour target: N = %f", pow_dimension(h) * kernel_norm);
   message("Density target: rho = %f", rho);
   message("div_v target:   div = %f", vel == 2 ? 3.f : 0.f);
   message("curl_v target: curl = [0., 0., %f]", vel == 3 ? -2.f : 0.f);