reshape.rs

use ndarray::prelude::*;

use itertools::enumerate;

use ndarray::Order;

#[test]
fn reshape()
{
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.into_shape_with_order((3, 3));
    assert!(u.is_err());
    let u = v.into_shape_with_order((2, 2, 2));
    assert!(u.is_ok());
    let u = u.unwrap();
    assert_eq!(u.shape(), &[2, 2, 2]);
    let s = u.into_shape_with_order((4, 2)).unwrap();
    assert_eq!(s.shape(), &[4, 2]);
    assert_eq!(s, aview2(&[[1, 2], [3, 4], [5, 6], [7, 8]]));
}

#[test]
#[should_panic(expected = "IncompatibleShape")]
fn reshape_error1()
{
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let _u = v.into_shape_with_order((2, 5)).unwrap();
}

#[test]
#[should_panic(expected = "IncompatibleLayout")]
fn reshape_error2()
{
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let mut u = v.into_shape_with_order((2, 2, 2)).unwrap();
    u.swap_axes(0, 1);
    let _s = u.into_shape_with_order((2, 4)).unwrap();
}

#[test]
fn reshape_f()
{
    let mut u = Array::zeros((3, 4).f());
    for (i, elt) in enumerate(u.as_slice_memory_order_mut().unwrap()) {
        *elt = i as i32;
    }
    let v = u.view();
    println!("{:?}", v);

    // noop ok
    let v2 = v.into_shape_with_order(((3, 4), Order::F));
    assert!(v2.is_ok());
    assert_eq!(v, v2.unwrap());

    let u = v.into_shape_with_order(((3, 2, 2), Order::F));
    assert!(u.is_ok());
    let u = u.unwrap();
    println!("{:?}", u);
    assert_eq!(u.shape(), &[3, 2, 2]);
    let s = u.into_shape_with_order(((4, 3), Order::F)).unwrap();
    println!("{:?}", s);
    assert_eq!(s.shape(), &[4, 3]);
    assert_eq!(s, aview2(&[[0, 4, 8], [1, 5, 9], [2, 6, 10], [3, 7, 11]]));
}

#[test]
fn to_shape_easy()
{
    // 1D -> C -> C
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.to_shape(((3, 3), Order::RowMajor));
    assert!(u.is_err());

    let u = v.to_shape(((2, 2, 2), Order::C));
    assert!(u.is_ok());

    let u = u.unwrap();
    assert!(u.is_view());
    assert_eq!(u.shape(), &[2, 2, 2]);
    assert_eq!(u, array![[[1, 2], [3, 4]], [[5, 6], [7, 8]]]);

    let s = u.to_shape((4, 2)).unwrap();
    assert_eq!(s.shape(), &[4, 2]);
    assert_eq!(s, aview2(&[[1, 2], [3, 4], [5, 6], [7, 8]]));

    // 1D -> F -> F
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.to_shape(((3, 3), Order::ColumnMajor));
    assert!(u.is_err());

    let u = v.to_shape(((2, 2, 2), Order::ColumnMajor));
    assert!(u.is_ok());

    let u = u.unwrap();
    assert!(u.is_view());
    assert_eq!(u.shape(), &[2, 2, 2]);
    assert_eq!(u, array![[[1, 5], [3, 7]], [[2, 6], [4, 8]]]);

    let s = u.to_shape(((4, 2), Order::ColumnMajor)).unwrap();
    assert_eq!(s.shape(), &[4, 2]);
    assert_eq!(s, array![[1, 5], [2, 6], [3, 7], [4, 8]]);
}

#[test]
fn to_shape_copy()
{
    // 1D -> C -> F
    let v = ArrayView::from(&[1, 2, 3, 4, 5, 6, 7, 8]);
    let u = v.to_shape(((4, 2), Order::RowMajor)).unwrap();
    assert_eq!(u.shape(), &[4, 2]);
    assert_eq!(u, array![[1, 2], [3, 4], [5, 6], [7, 8]]);

    let u = u.to_shape(((2, 4), Order::ColumnMajor)).unwrap();
    assert_eq!(u.shape(), &[2, 4]);
    assert_eq!(u, array![[1, 5, 2, 6], [3, 7, 4, 8]]);

    // 1D -> F -> C
    let v = ArrayView::from(&[1, 2, 3, 4, 5, 6, 7, 8]);
    let u = v.to_shape(((4, 2), Order::ColumnMajor)).unwrap();
    assert_eq!(u.shape(), &[4, 2]);
    assert_eq!(u, array![[1, 5], [2, 6], [3, 7], [4, 8]]);

    let u = u.to_shape((2, 4)).unwrap();
    assert_eq!(u.shape(), &[2, 4]);
    assert_eq!(u, array![[1, 5, 2, 6], [3, 7, 4, 8]]);
}

#[test]
fn to_shape_add_axis()
{
    // 1D -> C -> C
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.to_shape(((4, 2), Order::RowMajor)).unwrap();

    assert!(u.to_shape(((1, 4, 2), Order::RowMajor)).unwrap().is_view());
    assert!(u.to_shape(((1, 4, 2), Order::ColumnMajor)).unwrap().is_view());
}

#[test]
fn to_shape_copy_stride()
{
    let v = array![[1, 2, 3, 4], [5, 6, 7, 8]];
    let vs = v.slice(s![.., ..3]);
    let lin1 = vs.to_shape(6).unwrap();
    assert_eq!(lin1, array![1, 2, 3, 5, 6, 7]);
    assert!(lin1.is_owned());

    let lin2 = vs.to_shape((6, Order::ColumnMajor)).unwrap();
    assert_eq!(lin2, array![1, 5, 2, 6, 3, 7]);
    assert!(lin2.is_owned());
}

#[test]
fn to_shape_zero_len()
{
    let v = array![[1, 2, 3, 4], [5, 6, 7, 8]];
    let vs = v.slice(s![.., ..0]);
    let lin1 = vs.to_shape(0).unwrap();
    assert_eq!(lin1, array![]);
    assert!(lin1.is_view());
}

#[test]
#[should_panic(expected = "IncompatibleShape")]
fn to_shape_error1()
{
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let _u = v.to_shape((2, 5)).unwrap();
}

#[test]
#[should_panic(expected = "IncompatibleShape")]
fn to_shape_error2()
{
    // overflow
    let data = [3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let _u = v.to_shape((2, usize::MAX)).unwrap();
}

#[test]
fn to_shape_discontig()
{
    for &create_order in &[Order::C, Order::F] {
        let a = Array::from_iter(0..64);
        let mut a1 = a.to_shape(((4, 4, 4), create_order)).unwrap();
        a1.slice_collapse(s![.., ..;2, ..]); // now shape (4, 2, 4)
        assert!(a1.as_slice_memory_order().is_none());

        for &order in &[Order::C, Order::F] {
            let v1 = a1.to_shape(((2, 2, 2, 2, 2), order)).unwrap();
            assert!(v1.is_view());
            let v1 = a1.to_shape(((4, 1, 2, 1, 2, 2), order)).unwrap();
            assert!(v1.is_view());
            let v1 = a1.to_shape(((4, 2, 4), order)).unwrap();
            assert!(v1.is_view());
            let v1 = a1.to_shape(((8, 4), order)).unwrap();
            assert_eq!(v1.is_view(), order == create_order && create_order == Order::C,
                       "failed for {:?}, {:?}", create_order, order);
            let v1 = a1.to_shape(((4, 8), order)).unwrap();
            assert_eq!(v1.is_view(), order == create_order && create_order == Order::F,
                       "failed for {:?}, {:?}", create_order, order);
            let v1 = a1.to_shape((32, order)).unwrap();
            assert!(!v1.is_view());
        }
    }
}

#[test]
fn to_shape_broadcast()
{
    for &create_order in &[Order::C, Order::F] {
        let a = Array::from_iter(0..64);
        let mut a1 = a.to_shape(((4, 4, 4), create_order)).unwrap();
        a1.slice_collapse(s![.., ..1, ..]); // now shape (4, 1, 4)
        let v1 = a1.broadcast((4, 4, 4)).unwrap(); // Now shape (4, 4, 4)
        assert!(v1.as_slice_memory_order().is_none());

        for &order in &[Order::C, Order::F] {
            let v2 = v1.to_shape(((2, 2, 2, 2, 2, 2), order)).unwrap();
            assert_eq!(v2.strides(), match (create_order, order) {
                (Order::C, Order::C) => { &[32, 16, 0, 0, 2, 1] }
                (Order::C, Order::F) => { &[16, 32, 0, 0, 1, 2] }
                (Order::F, Order::C) => { &[2, 1, 0, 0, 32, 16] }
                (Order::F, Order::F) => { &[1, 2, 0, 0, 16, 32] }
                _other => unreachable!()
            });

            let v2 = v1.to_shape(((4, 4, 4), order)).unwrap();
            assert!(v2.is_view());
            let v2 = v1.to_shape(((8, 8), order)).unwrap();
            assert!(v2.is_owned());
        }
    }
}

#[test]
fn into_shape_with_order()
{
    // 1D -> C -> C
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.into_shape_with_order(((3, 3), Order::RowMajor));
    assert!(u.is_err());

    let u = v.into_shape_with_order(((2, 2, 2), Order::C));
    assert!(u.is_ok());

    let u = u.unwrap();
    assert_eq!(u.shape(), &[2, 2, 2]);
    assert_eq!(u, array![[[1, 2], [3, 4]], [[5, 6], [7, 8]]]);

    let s = u.into_shape_with_order((4, 2)).unwrap();
    assert_eq!(s.shape(), &[4, 2]);
    assert_eq!(s, aview2(&[[1, 2], [3, 4], [5, 6], [7, 8]]));

    // 1D -> F -> F
    let data = [1, 2, 3, 4, 5, 6, 7, 8];
    let v = aview1(&data);
    let u = v.into_shape_with_order(((3, 3), Order::ColumnMajor));
    assert!(u.is_err());

    let u = v.into_shape_with_order(((2, 2, 2), Order::ColumnMajor));
    assert!(u.is_ok());

    let u = u.unwrap();
    assert_eq!(u.shape(), &[2, 2, 2]);
    assert_eq!(u, array![[[1, 5], [3, 7]], [[2, 6], [4, 8]]]);

    let s = u
        .into_shape_with_order(((4, 2), Order::ColumnMajor))
        .unwrap();
    assert_eq!(s.shape(), &[4, 2]);
    assert_eq!(s, array![[1, 5], [2, 6], [3, 7], [4, 8]]);
}

#[test]
fn into_shape_clone()
{
    // 1D -> C -> C
    {
        let data = [1, 2, 3, 4, 5, 6, 7, 8];
        let v = Array::from(data.to_vec());
        let u = v.clone().into_shape_clone(((3, 3), Order::RowMajor));
        assert!(u.is_err());

        let u = v.clone().into_shape_clone(((2, 2, 2), Order::C));
        assert!(u.is_ok());

        let u = u.unwrap();
        assert_eq!(u.shape(), &[2, 2, 2]);
        assert_eq!(u, array![[[1, 2], [3, 4]], [[5, 6], [7, 8]]]);

        let s = u.into_shape_clone((4, 2)).unwrap();
        assert_eq!(s.shape(), &[4, 2]);
        assert_eq!(s, aview2(&[[1, 2], [3, 4], [5, 6], [7, 8]]));

        let u = v.clone().into_shape_clone(((2, 2, 2), Order::F));
        assert!(u.is_ok());

        let u = u.unwrap();
        assert_eq!(u.shape(), &[2, 2, 2]);
        assert_eq!(u, array![[[1, 5], [3, 7]], [[2, 6], [4, 8]]]);
    }

    // 1D -> F -> F
    {
        let data = [1, 2, 3, 4, 5, 6, 7, 8];
        let v = Array::from(data.to_vec());
        let u = v.clone().into_shape_clone(((3, 3), Order::ColumnMajor));
        assert!(u.is_err());

        let u = v.into_shape_clone(((2, 2, 2), Order::ColumnMajor));
        assert!(u.is_ok());

        let u = u.unwrap();
        assert_eq!(u.shape(), &[2, 2, 2]);
        assert_eq!(u, array![[[1, 5], [3, 7]], [[2, 6], [4, 8]]]);

        let s = u.into_shape_clone(((4, 2), Order::ColumnMajor)).unwrap();
        assert_eq!(s.shape(), &[4, 2]);
        assert_eq!(s, array![[1, 5], [2, 6], [3, 7], [4, 8]]);
    }
}