# Copyright (c) 2020 PaddlePaddle Authors. All Rights Reserved.
#
# Licensed under the Apache License, Version 2.0 (the "License");
# you may not use this file except in compliance with the License.
# You may obtain a copy of the License at
#
#     http://www.apache.org/licenses/LICENSE-2.0
#
# Unless required by applicable law or agreed to in writing, software
# distributed under the License is distributed on an "AS IS" BASIS,
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
# See the License for the specific language governing permissions and
# limitations under the License.

import unittest

import numpy as np
from op_test import (
    OpTest,
    convert_float_to_uint16,
    get_device_place,
    is_custom_device,
)

import paddle
from paddle import base
from paddle.base import core


class TestCrossOp(OpTest):
    def setUp(self):
        self.op_type = "cross"
        self.python_api = paddle.cross
        self.initTestCase()
        self.inputs = {
            'X': np.random.random(self.shape).astype(self.dtype),
            'Y': np.random.random(self.shape).astype(self.dtype),
        }
        if self.dtype is np.complex64 or self.dtype is np.complex128:
            self.inputs = {
                'X': (
                    np.random.random(self.shape)
                    + 1j * np.random.random(self.shape)
                ).astype(self.dtype),
                'Y': (
                    np.random.random(self.shape)
                    + 1j * np.random.random(self.shape)
                ).astype(self.dtype),
            }
        self.init_output()

    def initTestCase(self):
        self.attrs = {'dim': -2}
        self.dtype = np.float64
        self.shape = (1024, 3, 1)

    def init_output(self):
        x = np.squeeze(self.inputs['X'], 2)
        y = np.squeeze(self.inputs['Y'], 2)
        z_list = []
        for i in range(1024):
            z_list.append(np.cross(x[i], y[i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}

    def test_check_output(self):
        self.check_output(check_pir=True)

    def test_check_grad_normal(self):
        self.check_grad(['X', 'Y'], 'Out', check_pir=True)


class TestCrossOpCase1(TestCrossOp):
    def initTestCase(self):
        self.shape = (2048, 3)
        self.dtype = np.float32

    def init_output(self):
        z_list = []
        for i in range(2048):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


@unittest.skipIf(
    not (core.is_compiled_with_cuda() or is_custom_device()),
    "core is not compiled with CUDA",
)
class TestCrossFP16Op(TestCrossOp):
    def initTestCase(self):
        self.shape = (2048, 3)
        self.dtype = np.float16

    def init_output(self):
        z_list = []
        for i in range(2048):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossComplex64Op(TestCrossOp):
    def initTestCase(self):
        self.shape = (2048, 3)
        self.dtype = np.complex64

    def init_output(self):
        z_list = []
        for i in range(2048):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossComplex128Op(TestCrossOp):
    def initTestCase(self):
        self.shape = (2048, 3)
        self.dtype = np.complex128

    def init_output(self):
        z_list = []
        for i in range(2048):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


@unittest.skipIf(
    not (core.is_compiled_with_cuda() or is_custom_device())
    or not core.is_bfloat16_supported(get_device_place()),
    "core is not compiled with CUDA and not support the bfloat16",
)
class TestCrossBF16Op(OpTest):
    def setUp(self):
        self.op_type = "cross"
        self.python_api = paddle.cross
        self.initTestCase()
        self.x = np.random.random(self.shape).astype(np.float32)
        self.y = np.random.random(self.shape).astype(np.float32)
        self.inputs = {
            'X': convert_float_to_uint16(self.x),
            'Y': convert_float_to_uint16(self.y),
        }
        self.init_output()

    def initTestCase(self):
        self.attrs = {'dim': -2}
        self.dtype = np.uint16
        self.shape = (1024, 3, 1)

    def init_output(self):
        x = np.squeeze(self.x, 2)
        y = np.squeeze(self.y, 2)
        z_list = []
        for i in range(1024):
            z_list.append(np.cross(x[i], y[i]))
        out = np.array(z_list).astype(np.float32).reshape(self.shape)
        self.outputs = {'Out': convert_float_to_uint16(out)}

    def test_check_output(self):
        if core.is_compiled_with_cuda() or is_custom_device():
            place = get_device_place()
            if core.is_bfloat16_supported(place):
                self.check_output_with_place(place, check_pir=True)

    def test_check_grad_normal(self):
        if core.is_compiled_with_cuda() or is_custom_device():
            place = get_device_place()
            if core.is_bfloat16_supported(place):
                self.check_grad_with_place(
                    place, ['X', 'Y'], 'Out', check_pir=True
                )


class TestCrossAPI(unittest.TestCase):
    def input_data(self):
        self.data_x = np.array(
            [[1.0, 1.0, 1.0], [2.0, 2.0, 2.0], [3.0, 3.0, 3.0]]
        ).astype('float32')
        self.data_y = np.array(
            [[1.0, 1.0, 1.0], [1.0, 1.0, 1.0], [1.0, 1.0, 1.0]]
        ).astype('float32')
        self.data_x_zero = np.array([]).reshape(0, 3).astype('float32')
        self.data_y_zero = np.array([]).reshape(0, 3).astype('float32')

    def test_cross_api(self):
        self.input_data()

        main = paddle.static.Program()
        startup = paddle.static.Program()
        # case 1:
        with paddle.static.program_guard(main, startup):
            x = paddle.static.data(name='x', shape=[-1, 3], dtype="float32")
            y = paddle.static.data(name='y', shape=[-1, 3], dtype="float32")
            z = paddle.cross(x, y, axis=1)
            exe = base.Executor(base.CPUPlace())
            (res,) = exe.run(
                main,
                feed={'x': self.data_x, 'y': self.data_y},
                fetch_list=[z],
                return_numpy=False,
            )
        expect_out = np.array(
            [[0.0, 0.0, 0.0], [0.0, 0.0, 0.0], [0.0, 0.0, 0.0]]
        )
        np.testing.assert_allclose(expect_out, np.array(res), rtol=1e-05)

        main = paddle.static.Program()
        startup = paddle.static.Program()
        # case 2:
        with paddle.static.program_guard(main, startup):
            x = paddle.static.data(name='x', shape=[-1, 3], dtype="float32")
            y = paddle.static.data(name='y', shape=[-1, 3], dtype="float32")
            z = paddle.cross(x, y)
            exe = base.Executor(base.CPUPlace())
            (res,) = exe.run(
                main,
                feed={'x': self.data_x, 'y': self.data_y},
                fetch_list=[z],
                return_numpy=False,
            )
        expect_out = np.array(
            [[-1.0, -1.0, -1.0], [2.0, 2.0, 2.0], [-1.0, -1.0, -1.0]]
        )
        np.testing.assert_allclose(expect_out, np.array(res), rtol=1e-05)

        main = paddle.static.Program()
        startup = paddle.static.Program()
        # case 3:
        with paddle.static.program_guard(main, startup):
            x = paddle.static.data(name='x', shape=[0, 3], dtype="float32")
            y = paddle.static.data(name='y', shape=[0, 3], dtype="float32")
            z = paddle.cross(x, y, axis=1)
            exe = base.Executor(base.CPUPlace())
            (res,) = exe.run(
                main,
                feed={'x': self.data_x_zero, 'y': self.data_y_zero},
                fetch_list=[z],
                return_numpy=False,
            )
        expect_out = np.empty((0, 3))
        np.testing.assert_allclose(expect_out, np.array(res), rtol=1e-05)

        main = paddle.static.Program()
        startup = paddle.static.Program()

    def test_dygraph_api(self):
        self.input_data()
        # case 1:
        # with base.dygraph.guard():
        #     x = paddle.to_tensor(self.data_x)
        #     y = paddle.to_tensor(self.data_y)
        #     z = paddle.cross(x, y)
        #     np_z = z.numpy()
        # expect_out = np.array([[-1.0, -1.0, -1.0], [2.0, 2.0, 2.0],
        #                        [-1.0, -1.0, -1.0]])
        # np.testing.assert_allclose(expect_out, np_z, rtol=1e-05)

        # case 2:
        with base.dygraph.guard():
            x = paddle.to_tensor(self.data_x)
            y = paddle.to_tensor(self.data_y)
            z = paddle.cross(x, y, axis=1)
            np_z = z.numpy()
        expect_out = np.array(
            [[0.0, 0.0, 0.0], [0.0, 0.0, 0.0], [0.0, 0.0, 0.0]]
        )
        np.testing.assert_allclose(expect_out, np_z, rtol=1e-05)

        # case 3:
        with base.dygraph.guard():
            x = paddle.to_tensor(self.data_x_zero)
            y = paddle.to_tensor(self.data_y_zero)
            z = paddle.cross(x, y, axis=1)
            np_z = z.numpy()
        expect_out = np.empty((0, 3))
        np.testing.assert_allclose(expect_out, np_z, rtol=1e-05)


class TestCrossOpZeroSizeTest(TestCrossOp):
    def initTestCase(self):
        self.shape = (0, 3, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(0):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossOpZeroSizeTest1(TestCrossOp):
    def initTestCase(self):
        self.shape = (3, 0, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(3):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossOpZeroSizeTest2(TestCrossOp):
    def initTestCase(self):
        self.shape = (0, 0, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(0):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossOpZeroSizeCPUTest(TestCrossOp):
    def initTestCase(self):
        self.shape = (0, 0, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(0):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}

    def test_check_output(self):
        place = paddle.CPUPlace()
        self.check_output_with_place(place, check_pir=True)

    def test_check_grad_normal(self):
        place = paddle.CPUPlace()
        self.check_grad_with_place(place, ['X', 'Y'], 'Out', check_pir=True)


class TestCrossOpZeroSizeCPUTest1(TestCrossOpZeroSizeCPUTest):
    def initTestCase(self):
        self.shape = (3, 0, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(3):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestCrossOpZeroSizeCPUTest2(TestCrossOpZeroSizeCPUTest):
    def initTestCase(self):
        self.shape = (0, 0, 3)
        self.dtype = np.float64
        self.attr = {'dim': -1}

    def init_output(self):
        z_list = []
        for i in range(0):
            z_list.append(np.cross(self.inputs['X'][i], self.inputs['Y'][i]))
        self.outputs = {'Out': np.array(z_list).reshape(self.shape)}


class TestLinalgCrossDefaultDim(unittest.TestCase):
    def test_linalg_cross_default_dim(self):
        # Test that paddle.linalg.cross defaults to dim=-1, not axis=9 auto
        # Using shape [3, 2, 3] where auto-axis picks dim 0, but dim=-1 picks dim 2
        paddle.disable_static()
        np_x = np.random.randn(3, 2, 3).astype('float32')
        np_y = np.random.randn(3, 2, 3).astype('float32')

        x = paddle.to_tensor(np_x)
        y = paddle.to_tensor(np_y)

        # linalg.cross with default (should use dim=-1)
        out_default = paddle.linalg.cross(x, y)
        # linalg.cross with explicit dim=-1
        out_neg1 = paddle.linalg.cross(x, y, dim=-1)
        # linalg.cross with explicit dim=2
        out_dim2 = paddle.linalg.cross(x, y, dim=2)
        # linalg.cross with explicit dim=0
        out_dim0 = paddle.linalg.cross(x, y, dim=0)

        np.testing.assert_allclose(
            out_default.numpy(), out_neg1.numpy(), rtol=1e-5
        )
        np.testing.assert_allclose(
            out_default.numpy(), out_dim2.numpy(), rtol=1e-5
        )
        # dim=0 should give different result when shape is [3, 2, 3]
        with self.assertRaises(AssertionError):
            np.testing.assert_allclose(
                out_default.numpy(), out_dim0.numpy(), rtol=1e-5
            )

        paddle.enable_static()


if __name__ == '__main__':
    paddle.enable_static()
    unittest.main()