# -*- coding: utf-8 -*-
########################################################################
#
#       License: BSD
#       Created: January 11, 2011
#       Author:  Francesc Alted - francesc@blosc.org
#
########################################################################

from __future__ import absolute_import

import numpy as np
from numpy.testing import assert_array_equal
from bcolz.tests.common import (
    MayBeDiskTest, TestCase, unittest)
import bcolz


class constructorTest(MayBeDiskTest):

    open = False

    def test00a(self):
        """Testing `carray` reshape"""
        a = np.arange(16).reshape((2, 2, 4))
        b = bcolz.arange(16, rootdir=self.rootdir).reshape((2, 2, 4))
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test00b(self):
        """Testing `carray` reshape (large shape)"""
        a = np.arange(16000).reshape((20, 20, 40))
        b = bcolz.arange(16000, rootdir=self.rootdir).reshape((20, 20, 40))
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test01a(self):
        """Testing `zeros` constructor (I)"""
        a = np.zeros((2, 2, 4), dtype='i4')
        b = bcolz.zeros((2, 2, 4), dtype='i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test01b(self):
        """Testing `zeros` constructor (II)"""
        a = np.zeros(2, dtype='(2,4)i4')
        b = bcolz.zeros(2, dtype='(2,4)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test01c(self):
        """Testing `zeros` constructor (III)"""
        a = np.zeros((2, 2), dtype='(4,)i4')
        b = bcolz.zeros((2, 2), dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test02(self):
        """Testing `ones` constructor"""
        a = np.ones((2, 2), dtype='(4,)i4')
        b = bcolz.ones((2, 2), dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test03a(self):
        """Testing `fill` constructor (scalar default)"""
        a = np.ones((2, 200), dtype='(4,)i4') * 3
        b = bcolz.fill((2, 200), 3, dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test03b(self):
        """Testing `fill` constructor (array default)"""
        a = np.ones((2, 2), dtype='(4,)i4') * 3
        b = bcolz.fill(
            (2, 2), [3, 3, 3, 3], dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test04(self):
        """Testing `fill` constructor with open and resize (array default)"""
        a = np.ones((3, 200), dtype='(4,)i4') * 3
        b = bcolz.fill(
            (2, 200), [3, 3, 3, 3], dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        c = np.ones((1, 200), dtype='(4,)i4') * 3
        b.append(c)
        # print "b->", `b`, len(b), b[1]
        assert_array_equal(a, b, "Arrays are not equal")

    def test05(self):
        """Testing `fill` constructor with open and resize (nchunks>1)"""
        a = np.ones((3, 2000), dtype='(4,)i4') * 3
        b = bcolz.fill(
            (2, 2000), [3, 3, 3, 3], dtype='(4,)i4', rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        c = np.ones((1, 2000), dtype='(4,)i4') * 3
        b.append(c)
        # print "b->", `b`
        # We need to use the b[:] here to overcome a problem with the
        # assert_array_equal() function
        assert_array_equal(a, b[:], "Arrays are not equal")


class constructorMemoryTest(constructorTest, TestCase):
    disk = False
    open = False


class constructorDiskTest(constructorTest, TestCase):
    disk = True
    open = False


class constructorOpenTest(constructorTest, TestCase):
    disk = True
    open = True


class getitemTest(MayBeDiskTest):

    open = False

    def test00a(self):
        """Testing `__getitem()__` method with only a start (scalar)"""
        a = np.ones((2, 3), dtype="i4") * 3
        b = bcolz.fill((2, 3), 3, dtype="i4", rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = 1
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test00b(self):
        """Testing `__getitem()__` method with only a start (slice)"""
        a = np.ones((27, 2700), dtype="i4") * 3
        b = bcolz.fill((27, 2700), 3, dtype="i4", rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = slice(1)
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test01(self):
        """Testing `__getitem()__` method with a start and a stop"""
        a = np.ones((5, 2), dtype="i4") * 3
        b = bcolz.fill((5, 2), 3, dtype="i4", rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = slice(1, 4)
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test02(self):
        """Testing `__getitem()__` method with a start, stop, step"""
        a = np.ones((10, 2), dtype="i4") * 3
        b = bcolz.fill((10, 2), 3, dtype="i4", rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = slice(1, 9, 2)
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test03a(self):
        """Testing `__getitem()__` method with several slices (I)"""
        a = np.arange(12).reshape((4, 3))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (slice(1, 3, 1), slice(1, 4, 2))
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test03b(self):
        """Testing `__getitem()__` method with several slices (II)"""
        a = np.arange(24 * 1000).reshape((4 * 1000, 3, 2))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (slice(1, 3, 2), slice(1, 4, 2), slice(None))
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test03c(self):
        """Testing `__getitem()__` method with several slices (III)"""
        a = np.arange(120 * 1000).reshape((5 * 1000, 4, 3, 2))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (slice(None, None, 3), slice(1, 3, 2), slice(1, 4, 2))
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test04a(self):
        """Testing `__getitem()__` method with shape reduction (I)"""
        a = np.arange(12000).reshape((40, 300))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (1, 1)
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test04b(self):
        """Testing `__getitem()__` method with shape reduction (II)"""
        a = np.arange(12000).reshape((400, 30))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (1, slice(1, 4, 2))
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test04c(self):
        """Testing `__getitem()__` method with shape reduction (III)"""
        a = np.arange(6000).reshape((50, 40, 3))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (1, slice(1, 4, 2), 2)
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test05a(self):
        """Testing `__getitem()__` method with fancy indexing (I)"""
        a = np.arange(2000).reshape((50, 40))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = [3, 5, 2]
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test05b(self):
        """Testing `__getitem()__` method with fancy indexing (II)"""
        a = np.arange(2000).reshape((50, 40))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = ([0, 2], slice(None))
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test05c(self):
        """Testing `__getitem()__` method with fancy indexing (III)"""
        a = np.arange(2000).reshape((50, 40))
        b = bcolz.carray(a, rootdir=self.rootdir)
        if self.open:
            b = bcolz.open(rootdir=self.rootdir)
        sl = (slice(None), [0, 2])
        # print "b[sl]->", `b[sl]`
        self.assertTrue(a[sl].shape == b[sl].shape, "Shape is not equal")
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")


class getitemMemoryTest(getitemTest, TestCase):
    disk = False
    open = False


class getitemDiskTest(getitemTest, TestCase):
    disk = True
    open = False


class getitemOpenTest(getitemTest, TestCase):
    disk = True
    open = True


class setitemTest(MayBeDiskTest):

    open = False

    def test00a(self):
        """Testing `__setitem()__` method with only a start (scalar)"""
        a = np.ones((2, 3), dtype="i4") * 3
        b = bcolz.fill((2, 3), 3, dtype="i4", rootdir=self.rootdir)
        sl = slice(1)
        a[sl, :] = 0
        b[sl] = 0
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test00b(self):
        """Testing `__setitem()__` method with only a start (vector)"""
        a = np.ones((200, 300), dtype="i4") * 3
        b = bcolz.fill((200, 300), 3, dtype="i4", rootdir=self.rootdir)
        sl = slice(1)
        a[sl, :] = range(300)
        b[sl] = range(300)
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test01a(self):
        """Testing `__setitem()__` method with start,stop (scalar)"""
        a = np.ones((500, 200), dtype="i4") * 3
        b = bcolz.fill((500, 200), 3, dtype="i4", rootdir=self.rootdir,
                       cparams=bcolz.cparams())
        sl = slice(100, 400)
        a[sl, :] = 0
        b[sl] = 0
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")
        # assert_array_equal(a[:], b[:], "Arrays are not equal")

    def test01b(self):
        """Testing `__setitem()__` method with start,stop (vector)"""
        a = np.ones((5, 2), dtype="i4") * 3
        b = bcolz.fill((5, 2), 3, dtype="i4", rootdir=self.rootdir)
        sl = slice(1, 4)
        a[sl, :] = range(2)
        b[sl] = range(2)
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test02a(self):
        """Testing `__setitem()__` method with start,stop,step (scalar)"""
        a = np.ones((1000, 200), dtype="i4") * 3
        b = bcolz.fill((1000, 200), 3, dtype="i4", rootdir=self.rootdir)
        sl = slice(100, 800, 3)
        a[sl, :] = 0
        b[sl] = 0
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test02b(self):
        """Testing `__setitem()__` method with start,stop,step (scalar)"""
        a = np.ones((10, 2), dtype="i4") * 3
        b = bcolz.fill((10, 2), 3, dtype="i4", rootdir=self.rootdir)
        sl = slice(1, 8, 3)
        a[sl, :] = range(2)
        b[sl] = range(2)
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "b[sl]->", `b[sl]`, `b`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test03a(self):
        """Testing `__setitem()__` method with several slices (I)"""
        a = np.arange(12000).reshape((400, 30))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (slice(1, 3, 1), slice(1, None, 2))
        # print "before->", `b[sl]`
        a[sl] = [[1], [2]]
        b[sl] = [[1], [2]]
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[:], b[:], "Arrays are not equal")

    def test03b(self):
        """Testing `__setitem()__` method with several slices (II)"""
        a = np.arange(24000).reshape((400, 3, 20))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (slice(1, 3, 1), slice(1, None, 2), slice(1))
        # print "before->", `b[sl]`
        a[sl] = [[[1]], [[2]]]
        b[sl] = [[[1]], [[2]]]
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[:], b[:], "Arrays are not equal")

    def test03c(self):
        """Testing `__setitem()__` method with several slices (III)"""
        a = np.arange(120).reshape((5, 4, 3, 2))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (slice(1, 3), slice(1, 3, 1), slice(1, None, 2), slice(1))
        # print "before->", `b[sl]`
        a[sl] = [[[[1]], [[2]]]] * 2
        b[sl] = [[[[1]], [[2]]]] * 2
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[:], b[:], "Arrays are not equal")

    def test03d(self):
        """Testing `__setitem()__` method with several slices (IV)"""
        a = np.arange(120).reshape((5, 4, 3, 2))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (slice(1, 3), slice(1, 3, 1), slice(1, None, 2), slice(1))
        # print "before->", `b[sl]`
        a[sl] = 2
        b[sl] = 2
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[:], b[:], "Arrays are not equal")

    def test04a(self):
        """Testing `__setitem()__` method with shape reduction (I)"""
        a = np.arange(12).reshape((4, 3))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (1, 1)
        # print "before->", `b[sl]`
        a[sl] = 2
        b[sl] = 2
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test04b(self):
        """Testing `__setitem()__` method with shape reduction (II)"""
        a = np.arange(12).reshape((4, 3))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (1, slice(1, 4, 2))
        # print "before->", `b[sl]`
        a[sl] = 2
        b[sl] = 2
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")

    def test04c(self):
        """Testing `__setitem()__` method with shape reduction (III)"""
        a = np.arange(24).reshape((4, 3, 2))
        b = bcolz.carray(a, rootdir=self.rootdir)
        sl = (1, 2, slice(None, None, None))
        # print "before->", `b[sl]`
        a[sl] = 2
        b[sl] = 2
        if self.open:
            b.flush()
            b = bcolz.open(rootdir=self.rootdir)
        # print "after->", `b[sl]`
        assert_array_equal(a[sl], b[sl], "Arrays are not equal")


class setitemMemoryTest(setitemTest, TestCase):
    disk = False


class setitemDiskTest(setitemTest, TestCase):
    disk = True


class setitemOpenTest(setitemTest, TestCase):
    disk = True
    open = True


class appendTest(MayBeDiskTest):

    def test00a(self):
        """Testing `append()` method (correct shape)"""
        a = np.ones((2, 300), dtype="i4") * 3
        b = bcolz.fill((1, 300), 3, dtype="i4", rootdir=self.rootdir)
        b.append([(3,) * 300])
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test00b(self):
        """Testing `append()` method (correct shape, single row)"""
        a = np.ones((2, 300), dtype="i4") * 3
        b = bcolz.fill((1, 300), 3, dtype="i4", rootdir=self.rootdir)
        b.append((3,) * 300)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test01(self):
        """Testing `append()` method (incorrect shape)"""
        a = np.ones((2, 3), dtype="i4") * 3
        b = bcolz.fill((1, 3), 3, dtype="i4", rootdir=self.rootdir)
        self.assertRaises(ValueError, b.append, [(3, 3)])

    def test02(self):
        """Testing `append()` method (several rows)"""
        a = np.ones((4, 3), dtype="i4") * 3
        b = bcolz.fill((1, 3), 3, dtype="i4", rootdir=self.rootdir)
        b.append([(3, 3, 3)] * 3)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")


class appendMemoryTest(appendTest, TestCase):
    disk = False


class appendDiskTest(appendTest, TestCase):
    disk = True


class resizeTest(MayBeDiskTest):

    def test00a(self):
        """Testing `resize()` (trim)"""
        a = np.ones((2, 3), dtype="i4")
        b = bcolz.ones((3, 3), dtype="i4", rootdir=self.rootdir)
        b.resize(2)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test00b(self):
        """Testing `resize()` (trim to zero)"""
        a = np.ones((0, 3), dtype="i4")
        b = bcolz.ones((3, 3), dtype="i4", rootdir=self.rootdir)
        b.resize(0)
        # print "b->", `b`
        # The next does not work well for carrays with shape (0,)
        # assert_array_equal(a, b, "Arrays are not equal")
        self.assertTrue("a.dtype.base == b.dtype.base")
        self.assertTrue("a.shape == b.shape+b.dtype.shape")

    def test01(self):
        """Testing `resize()` (enlarge)"""
        a = np.ones((4, 3), dtype="i4")
        b = bcolz.ones((3, 3), dtype="i4", rootdir=self.rootdir)
        b.resize(4)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")


class resizeMemoryTest(resizeTest, TestCase):
    disk = False


class resizeDiskTest(resizeTest, TestCase):
    disk = True


class iterTest(TestCase):

    def test00(self):
        """Testing `iter()` (no start, stop, step)"""
        a = np.ones((3,), dtype="i4")
        b = bcolz.ones((1000, 3), dtype="i4")
        # print "b->", `b`
        for r in b.iter():
            assert_array_equal(a, r, "Arrays are not equal")

    def test01(self):
        """Testing `iter()` (w/ start, stop)"""
        a = np.ones((3,), dtype="i4")
        b = bcolz.ones((1000, 3), dtype="i4")
        # print "b->", `b`
        for r in b.iter(start=10):
            assert_array_equal(a, r, "Arrays are not equal")

    def test02(self):
        """Testing `iter()` (w/ start, stop, step)"""
        a = np.ones((3,), dtype="i4")
        b = bcolz.ones((1000, 3), dtype="i4")
        # print "b->", `b`
        for r in b.iter(15, 100, 3):
            assert_array_equal(a, r, "Arrays are not equal")


class reshapeTest(TestCase):

    def test00a(self):
        """Testing `reshape()` (unidim -> ndim)"""
        a = np.ones((3, 4), dtype="i4")
        b = bcolz.ones(12, dtype="i4").reshape((3, 4))
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test00b(self):
        """Testing `reshape()` (unidim -> ndim, -1 in newshape (I))"""
        a = np.ones((3, 4), dtype="i4")
        b = bcolz.ones(12, dtype="i4").reshape((-1, 4))
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test00c(self):
        """Testing `reshape()` (unidim -> ndim, -1 in newshape (II))"""
        a = np.ones((3, 4), dtype="i4")
        b = bcolz.ones(12, dtype="i4").reshape((3, -1))
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test01(self):
        """Testing `reshape()` (ndim -> unidim)"""
        a = np.ones(12, dtype="i4")
        c = bcolz.ones(12, dtype="i4").reshape((3, 4))
        b = c.reshape(12)
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test02a(self):
        """Testing `reshape()` (ndim -> ndim, I)"""
        a = np.arange(12, dtype="i4").reshape((3, 4))
        c = bcolz.arange(12, dtype="i4").reshape((4, 3))
        b = c.reshape((3, 4))
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test02b(self):
        """Testing `reshape()` (ndim -> ndim, II)"""
        a = np.arange(24, dtype="i4").reshape((2, 3, 4))
        c = bcolz.arange(24, dtype="i4").reshape((4, 3, 2))
        b = c.reshape((2, 3, 4))
        # print "b->", `b`
        assert_array_equal(a, b, "Arrays are not equal")

    def test03(self):
        """Testing `reshape()` (0-dim)"""
        a = np.ones((0, 4), dtype="i4")
        b = bcolz.ones(0, dtype="i4").reshape((0, 4))
        # print "b->", `b`
        # The next does not work well for carrays with shape (0,)
        # assert_array_equal(a, b, "Arrays are not equal")
        self.assertTrue(a.dtype.base == b.dtype.base)
        self.assertTrue(a.shape == b.shape + b.dtype.shape)


class compoundTest():

    def test00(self):
        """Testing compound types (creation)"""
        a = np.ones((300, 4), dtype=self.dtype)
        b = bcolz.ones((300, 4), dtype=self.dtype)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test01(self):
        """Testing compound types (append)"""
        a = np.ones((300, 4), dtype=self.dtype)
        b = bcolz.carray([], dtype=self.dtype).reshape((0, 4))
        b.append(a)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test02(self):
        """Testing compound types (iter)"""
        a = np.ones((3,), dtype=self.dtype)
        b = bcolz.ones((1000, 3), dtype=self.dtype)
        # print "b->", `b`
        for r in b.iter():
            # print "r-->", r
            assert_array_equal(a, r, "Arrays are not equal")


class plainCompoundTest(compoundTest, TestCase):
    dtype = np.dtype("i4,i8")


class nestedCompoundTest(compoundTest, TestCase):
    dtype = np.dtype([('f1', [('f1', 'i2'), ('f2', 'i4')])])


class stringTest(TestCase):

    def test00(self):
        """Testing string types (creation)"""
        a = np.array([["ale", "ene"], ["aco", "ieie"]], dtype="S4")
        b = bcolz.carray(a)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test01(self):
        """Testing string types (append)"""
        a = np.ones((300, 4), dtype="S4")
        b = bcolz.carray([], dtype="S4").reshape((0, 4))
        b.append(a)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test02(self):
        """Testing string types (iter)"""
        a = np.ones((3,), dtype="S40")
        b = bcolz.ones((1000, 3), dtype="S40")
        # print "b->", `b`
        for r in b.iter():
            # print "r-->", r
            assert_array_equal(a, r, "Arrays are not equal")


class unicodeTest(TestCase):

    def test00(self):
        """Testing unicode types (creation)"""
        a = np.array([[u"aŀle", u"eñe"], [u"açò", u"áèâë"]], dtype="U4")
        b = bcolz.carray(a)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test01(self):
        """Testing unicode types (append)"""
        a = np.ones((300, 4), dtype="U4")
        b = bcolz.carray([], dtype="U4").reshape((0, 4))
        b.append(a)
        # print "b.dtype-->", b.dtype
        # print "b->", `b`
        self.assertTrue(a.dtype == b.dtype.base)
        assert_array_equal(a, b[:], "Arrays are not equal")

    def test02(self):
        """Testing unicode types (iter)"""
        a = np.ones((3,), dtype="U40")
        b = bcolz.ones((1000, 3), dtype="U40")
        # print "b->", `b`
        for r in b.iter():
            # print "r-->", r
            assert_array_equal(a, r, "Arrays are not equal")


class evalTest():

    vm = "python"

    def setUp(self):
        self.prev_vm = bcolz.defaults.eval_vm
        if bcolz.numexpr_here:
            bcolz.defaults.eval_vm = self.vm
        else:
            bcolz.defaults.eval_vm = "python"

    def tearDown(self):
        bcolz.defaults.eval_vm = self.prev_vm

    def test00a(self):
        """Testing evaluation of ndcarrays (bool out)"""
        a = np.arange(np.prod(self.shape)).reshape(self.shape)
        b = bcolz.arange(np.prod(self.shape)).reshape(self.shape)
        outa = eval("a>0")
        outb = bcolz.eval("b>0")
        assert_array_equal(outa, outb, "Arrays are not equal")

    def test00b(self):
        """Testing evaluation of ndcarrays (bool out, NumPy)"""
        a = np.arange(np.prod(self.shape)).reshape(self.shape)
        b = bcolz.arange(np.prod(self.shape)).reshape(self.shape)
        outa = eval("a>0")
        outb = bcolz.eval("b>0", out_flavor='numpy')
        assert_array_equal(outa, outb, "Arrays are not equal")

    def test01(self):
        """Testing evaluation of ndcarrays (int out)"""
        a = np.arange(np.prod(self.shape)).reshape(self.shape)
        b = bcolz.arange(np.prod(self.shape)).reshape(self.shape)
        outa = eval("a*2.+1")
        outb = bcolz.eval("b*2.+1")
        assert_array_equal(outa, outb, "Arrays are not equal")

    def test02(self):
        """Testing evaluation of ndcarrays (reduction, no axis)"""
        a = np.arange(np.prod(self.shape)).reshape(self.shape)
        b = bcolz.arange(np.prod(self.shape)).reshape(self.shape)
        if bcolz.defaults.eval_vm == "python":
            assert_array_equal(sum(a), bcolz.eval("sum(b)"),
                               "Arrays are not equal")
        else:
            self.assertEqual(a.sum(), bcolz.eval("sum(b)"))

    def test02b(self):
        """Testing evaluation of ndcarrays (reduction, with axis)"""
        a = np.arange(np.prod(self.shape)).reshape(self.shape)
        b = bcolz.arange(np.prod(self.shape)).reshape(self.shape)
        if bcolz.defaults.eval_vm == "python":
            # The Python VM does not have support for `axis` param
            assert_array_equal(sum(a), bcolz.eval("sum(b)"),
                               "Arrays are not equal")
        else:
            assert_array_equal(a.sum(axis=1), bcolz.eval("sum(b, axis=1)"),
                               "Arrays are not equal")


class d2eval_python(evalTest, TestCase):
    shape = (3, 4)


class d2eval_ne(evalTest, TestCase):
    shape = (3, 4)
    vm = "numexpr"


class d3eval_python(evalTest, TestCase):
    shape = (3, 4, 5)


class d3eval_ne(evalTest, TestCase):
    shape = (3, 4, 5)
    vm = "numexpr"


class d4eval_python(evalTest, TestCase):
    shape = (3, 40, 50, 2)


class d4eval_ne(evalTest, TestCase):
    shape = (3, 40, 50, 2)
    vm = "numexpr"


class computeMethodsTest(TestCase):

    def test00(self):
        """Testing sum()."""
        a = np.arange(int(1e5)).reshape(10, int(1e4))
        sa = a.sum()
        ac = bcolz.carray(a)
        sac = ac.sum()
        # print "numpy sum-->", sa
        # print "carray sum-->", sac
        self.assertEqual(
            sa.dtype, sac.dtype, "sum() is not working correctly.")
        self.assertEqual(sa, sac, "sum() is not working correctly.")


if __name__ == '__main__':
    unittest.main(verbosity=2)


# Local Variables:
# mode: python
# tab-width: 4
# fill-column: 72
# End:
