FIX: Solve axis iteration order problem by sorting axes

Author: bluss

src/impl_owned_array.rs

@@ -252,7 +252,7 @@ impl<A, D> Array<A, D>
     ///            [1., 1., 1., 1.],
     ///            [1., 1., 1., 1.]]);
     /// ```
-    pub fn try_append_array(&mut self, axis: Axis, array: ArrayView<A, D>)
+    pub fn try_append_array(&mut self, axis: Axis, mut array: ArrayView<A, D>)
         -> Result<(), ShapeError>
     where
         A: Clone,
@@ -312,7 +312,7 @@ impl<A, D> Array<A, D>
             // make a raw view with the new row
             // safe because the data was "full"
             let tail_ptr = self.data.as_end_nonnull();
-            let tail_view = RawArrayViewMut::new(tail_ptr, array_shape, strides.clone());
+            let mut tail_view = RawArrayViewMut::new(tail_ptr, array_shape, strides.clone());
 
             struct SetLenOnDrop<'a, A: 'a> {
                 len: usize,
@@ -332,37 +332,86 @@ impl<A, D> Array<A, D>
                 }
             }
 
-            // we have a problem here XXX
-            //
             // To be robust for panics and drop the right elements, we want
             // to fill the tail in-order, so that we can drop the right elements on
-            // panic. Don't know how to achieve that.
+            // panic.
             //
-            // It might be easier to retrace our steps in a scope guard to drop the right
-            // elements.. (PartialArray style).
+            // We have: Zip::from(tail_view).and(array)
+            // Transform tail_view into standard order by inverting and moving its axes.
+            // Keep the Zip traversal unchanged by applying the same axis transformations to
+            // `array`. This ensures the Zip traverses the underlying memory in order.
             //
-            // assign the new elements
+            // XXX It would be possible to skip this transformation if the element
+            // doesn't have drop. However, in the interest of code coverage, all elements
+            // use this code initially.
+
+            if tail_view.ndim() > 1 {
+                for i in 0..tail_view.ndim() {
+                    if tail_view.stride_of(Axis(i)) < 0 {
+                        tail_view.invert_axis(Axis(i));
+                        array.invert_axis(Axis(i));
+                    }
+                }
+                sort_axes_to_standard_order(&mut tail_view, &mut array);
+            } 
             Zip::from(tail_view).and(array)
+                .debug_assert_c_order()
                 .for_each(|to, from| {
                     to.write(from.clone());
                     length_guard.len += 1;
                 });
 
-            //length_guard.len += len_to_append;
-            dbg!(len_to_append);
             drop(length_guard);
 
             // update array dimension
             self.strides = strides;
             self.dim = res_dim;
-            dbg!(&self.dim);
-
         }
         // multiple assertions after pointer & dimension update
         debug_assert_eq!(self.data.len(), self.len());
         debug_assert_eq!(self.len(), new_len);
-        debug_assert!(self.is_standard_layout());
 
         Ok(())
     }
 }
+
+fn sort_axes_to_standard_order<S, S2, D>(a: &mut ArrayBase<S, D>, b: &mut ArrayBase<S2, D>)
+where
+    S: RawData,
+    S2: RawData,
+    D: Dimension,
+{
+    if a.ndim() <= 1 {
+        return;
+    }
+    sort_axes_impl(&mut a.dim, &mut a.strides, &mut b.dim, &mut b.strides);
+    debug_assert!(a.is_standard_layout());
+}
+
+fn sort_axes_impl<D>(adim: &mut D, astrides: &mut D, bdim: &mut D, bstrides: &mut D)
+where
+    D: Dimension,
+{
+    debug_assert!(adim.ndim() > 1);
+    debug_assert_eq!(adim.ndim(), bdim.ndim());
+    // bubble sort axes
+    let mut changed = true;
+    while changed {
+        changed = false;
+        for i in 0..adim.ndim() - 1 {
+            let axis_i = i;
+            let next_axis = i + 1;
+
+            // make sure higher stride axes sort before.
+            debug_assert!(astrides.slice()[axis_i] as isize >= 0);
+            if (astrides.slice()[axis_i] as isize) < astrides.slice()[next_axis] as isize {
+                changed = true;
+                adim.slice_mut().swap(axis_i, next_axis);
+                astrides.slice_mut().swap(axis_i, next_axis);
+                bdim.slice_mut().swap(axis_i, next_axis);
+                bstrides.slice_mut().swap(axis_i, next_axis);
+            }
+        }
+    }
+}
+

src/zip/mod.rs

@@ -673,6 +673,13 @@ macro_rules! map_impl {
                 self.build_and(part)
             }
 
+            #[allow(unused)]
+            #[inline]
+            pub(crate) fn debug_assert_c_order(self) -> Self {
+                debug_assert!(self.layout.is(CORDER) || self.layout_tendency >= 0);
+                self
+            }
+
             fn build_and<P>(self, part: P) -> Zip<($($p,)* P, ), D>
                 where P: NdProducer<Dim=D>,
             {

tests/append.rs

@@ -105,3 +105,44 @@ fn append_array1() {
                [5., 5., 4., 4.],
                [3., 3., 2., 2.]]);
 }
+
+#[test]
+fn append_array_3d() {
+    let mut a = Array::zeros((0, 2, 2));
+    a.try_append_array(Axis(0), array![[[0, 1], [2, 3]]].view()).unwrap();
+    println!("{:?}", a);
+
+    let aa = array![[[51, 52], [53, 54]], [[55, 56], [57, 58]]];
+    let av = aa.view();
+    println!("Send {:?} to append", av);
+    a.try_append_array(Axis(0), av.clone()).unwrap();
+
+    a.swap_axes(0, 1);
+    let aa = array![[[71, 72], [73, 74]], [[75, 76], [77, 78]]];
+    let mut av = aa.view();
+    av.swap_axes(0, 1);
+    println!("Send {:?} to append", av);
+    a.try_append_array(Axis(1), av.clone()).unwrap();
+    println!("{:?}", a);
+    let aa = array![[[81, 82], [83, 84]], [[85, 86], [87, 88]]];
+    let mut av = aa.view();
+    av.swap_axes(0, 1);
+    println!("Send {:?} to append", av);
+    a.try_append_array(Axis(1), av).unwrap();
+    println!("{:?}", a);
+    assert_eq!(a,
+        array![[[0, 1],
+                [51, 52],
+                [55, 56],
+                [71, 72],
+                [75, 76],
+                [81, 82],
+                [85, 86]],
+               [[2, 3],
+                [53, 54],
+                [57, 58],
+                [73, 74],
+                [77, 78],
+                [83, 84],
+                [87, 88]]]);
+}