Add Grid::repair test

Author: felixhekhorn

pineappl/src/grid.rs

@@ -1609,6 +1609,39 @@ mod tests {
         assert_eq!(grid.orders().len(), 1);
     }
 
+    #[test]
+    fn grid_repair() {
+        use super::super::packed_array::PackedArray;
+        // create emtpy grid
+        let mut grid = Grid::new(
+            BinsWithFillLimits::from_fill_limits([0.0, 0.25, 0.5, 0.75, 1.0].to_vec()).unwrap(),
+            vec![Order::new(0, 2, 0, 0, 0)],
+            vec![channel![1.0 * (2, 2) + 1.0 * (4, 4)]],
+            PidBasis::Pdg,
+            vec![Conv::new(ConvType::UnpolPDF, 2212); 2],
+            v0::default_interps(false, 2),
+            vec![Kinematics::Scale(0), Kinematics::X(0), Kinematics::X(1)],
+            Scales {
+                ren: ScaleFuncForm::Scale(0),
+                fac: ScaleFuncForm::Scale(0),
+                frg: ScaleFuncForm::NoScale,
+            },
+        );
+        let was_reparied = grid.repair();
+        assert!(!was_reparied);
+        // insert nothing
+        let x = vec![
+            0.015625, 0.03125, 0.0625, 0.125, 0.1875, 0.25, 0.375, 0.5, 0.75, 1.0,
+        ];
+        let mut ar = PackedArray::new(vec![1, 10, 10]);
+        ar[[0, 0, 0]] = 0.;
+        let sg: SubgridEnum =
+            ImportSubgridV1::new(ar, vec![vec![0.0], x.clone(), x.clone()]).into();
+        grid.subgrids_mut()[[0, 0, 0]] = sg;
+        let was_reparied = grid.repair();
+        assert!(was_reparied);
+    }
+
     #[test]
     fn grid_merge_orders() {
         let mut grid = Grid::new(