#   Copyright (c) 2018 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 decorator_helper import prog_scope
from op_test import get_device_place
from utils import dygraph_guard

import paddle
from paddle import base


class TestMathOpPatches(unittest.TestCase):
    @classmethod
    def setUp(self):
        np.random.seed(1024)
        paddle.enable_static()

    @prog_scope()
    def test_add_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = a + 10
        ab = paddle.concat([a, b], axis=1)
        c = ab + 10
        d = ab + a
        # e = a + ab
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        b_np, c_np, d_np = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b, c, d]
        )
        np.testing.assert_allclose(a_np + 10, b_np, rtol=1e-05)
        ab_np = np.concatenate([a_np, b_np], axis=1)
        np.testing.assert_allclose(ab_np + 10, c_np, rtol=1e-05)
        d_expected = ab_np + np.concatenate([a_np, a_np], axis=1)
        np.testing.assert_allclose(d_expected, d_np, rtol=1e-05)

    @prog_scope()
    def test_radd_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = 10 + a
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(a_np + 10, b_np, rtol=1e-05)

    @prog_scope()
    def test_sub_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = a - 10
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(a_np - 10, b_np, rtol=1e-05)

    @prog_scope()
    def test_rsub_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = 10 - a
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(10 - a_np, b_np, rtol=1e-05)

    @prog_scope()
    def test_mul_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = a * 10
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(a_np * 10, b_np, rtol=1e-05)

    @prog_scope()
    def test_rmul_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = 10 * a
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(10 * a_np, b_np, rtol=1e-05)

    @prog_scope()
    def test_div_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = a / 10
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(a_np / 10, b_np, rtol=1e-05)

    @prog_scope()
    def test_rdiv_scalar(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = 10 / a
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32') + 1e-2

        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(10 / a_np, b_np, rtol=1e-05)

    @prog_scope()
    def test_div_two_tensor(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = paddle.static.data(name="b", shape=[-1, 1])
        c = a / b
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        b_np = np.random.random(size=[10, 1]).astype('float32') + 1e-2
        (c_np,) = exe.run(
            base.default_main_program(),
            feed={"a": a_np, 'b': b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np / b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_mul_two_tensor(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = paddle.static.data(name="b", shape=[-1, 1])
        c = a * b
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        b_np = np.random.random(size=[10, 1]).astype('float32')
        (c_np,) = exe.run(
            base.default_main_program(),
            feed={"a": a_np, 'b': b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np * b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_add_two_tensor(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = paddle.static.data(name="b", shape=[-1, 1])
        c = a + b
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        b_np = np.random.random(size=[10, 1]).astype('float32')
        (c_np,) = exe.run(
            base.default_main_program(),
            feed={"a": a_np, 'b': b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np + b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_sub_two_tensor(self):
        a = paddle.static.data(name="a", shape=[-1, 1])
        b = paddle.static.data(name="b", shape=[-1, 1])
        c = a - b
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.random(size=[10, 1]).astype('float32')
        b_np = np.random.random(size=[10, 1]).astype('float32')
        (c_np,) = exe.run(
            base.default_main_program(),
            feed={"a": a_np, 'b': b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np - b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_integer_div(self):
        a = paddle.static.data(name="a", shape=[-1, 1], dtype='int64')
        b = a / 7
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.array([3, 4, 10, 14, 9, 18]).astype('int64')
        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )

        b_np_actual = (a_np / 7).astype('float32')
        np.testing.assert_allclose(b_np, b_np_actual, rtol=1e-05)

    @prog_scope()
    def test_equal(self):
        a = paddle.static.data(name="a", shape=[-1, 1], dtype='float32')
        b = paddle.static.data(name="b", shape=[-1, 1], dtype='float32')
        c = a == b

        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.array([3, 4, 10, 14, 9, 18]).astype('float32')
        b_np = np.array([3, 4, 11, 15, 8, 18]).astype('float32')

        (c_np,) = exe.run(
            base.default_main_program(),
            feed={"a": a_np, "b": b_np},
            fetch_list=[c],
        )

        np.testing.assert_array_equal(c_np, a_np == b_np)
        self.assertEqual(c.dtype, paddle.bool)

    @prog_scope()
    def test_neg(self):
        a = paddle.static.data(name="a", shape=[-1, 10, 1], dtype='float32')
        if not paddle.framework.use_pir_api():
            a.desc.set_need_check_feed(False)
        b = -a
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.uniform(-1, 1, size=[10, 1]).astype('float32')

        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(-a_np, b_np, rtol=1e-05)

    @prog_scope()
    def test_abs(self):
        # test for real number
        a = paddle.static.data(name="a", shape=[-1, 10, 1], dtype='float32')
        if not paddle.framework.use_pir_api():
            a.desc.set_need_check_feed(False)
        b = abs(a)  # call __abs__
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.uniform(-1, 1, size=[10, 1]).astype('float32')

        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(np.abs(a_np), b_np, rtol=1e-05)

    @prog_scope()
    def test_abs_complex(self):
        # test for complex number
        a = paddle.static.data(name="a", shape=[-1, 10, 1], dtype='complex64')
        if not paddle.framework.use_pir_api():
            a.desc.set_need_check_feed(False)
        b = abs(a)  # call __abs__
        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.uniform(-1, 1, size=[10, 1]).astype(
            'float32'
        ) + 1j * np.random.uniform(-1, 1, size=[10, 1]).astype('float32')

        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(np.abs(a_np), b_np, rtol=1e-05)

    @prog_scope()
    def test_astype(self):
        a = paddle.static.data(name="a", shape=[-1, 10, 1])
        if not paddle.framework.use_pir_api():
            a.desc.set_need_check_feed(False)
        b = a.astype('float64')

        c = a.astype(a.dtype)
        self.assertTrue(c.is_same(a))

        place = base.CPUPlace()
        exe = base.Executor(place)
        a_np = np.random.uniform(-1, 1, size=[10, 1]).astype('float64')

        (b_np,) = exe.run(
            base.default_main_program(), feed={"a": a_np}, fetch_list=[b]
        )
        np.testing.assert_allclose(a_np.astype('float32'), b_np, rtol=1e-05)

    def test_bitwise_and(self):
        temp = 2
        x_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        y_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        out_np = x_np & y_np
        e_np = temp & x_np
        x = paddle.static.data(name="x", shape=[2, 3, 5], dtype="int32")
        y = paddle.static.data(name="y", shape=[2, 3, 5], dtype="int32")
        z = x & y
        e = temp & x
        exe = base.Executor()
        (out, e_out) = exe.run(
            base.default_main_program(),
            feed={"x": x_np, "y": y_np},
            fetch_list=[z, e],
        )
        np.testing.assert_array_equal(out, out_np)
        np.testing.assert_array_equal(e_out, e_np)

    @prog_scope()
    def test_bitwise_or(self):
        x_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        y_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        out_np = x_np | y_np

        x = paddle.static.data(name="x", shape=[2, 3, 5], dtype="int32")
        y = paddle.static.data(name="y", shape=[2, 3, 5], dtype="int32")
        z = x | y

        exe = base.Executor()
        out = exe.run(
            base.default_main_program(),
            feed={"x": x_np, "y": y_np},
            fetch_list=[z],
        )
        np.testing.assert_array_equal(out[0], out_np)

    @prog_scope()
    def test_ror(self):
        place = get_device_place()
        x_int = 5
        y_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        y = paddle.static.data("y", y_np.shape, dtype=y_np.dtype)
        z = x_int | y
        exe = paddle.static.Executor(place)
        out = exe.run(
            feed={'y': y_np},
            fetch_list=[z],
        )
        out_ref = x_int | y_np
        np.testing.assert_array_equal(out[0], out_ref)
        x_bool = True
        res_ror_bool = x_bool | y
        out_bool = exe.run(
            feed={'y': y_np},
            fetch_list=[res_ror_bool],
        )
        res_py_bool = x_bool | y_np
        np.testing.assert_array_equal(out_bool[0], res_py_bool)

        for x_invalid in (
            np.float32(5.0),
            np.float64(5.0),
            np.complex64(5),
            np.complex128(5.0 + 2j),
        ):
            with self.assertRaises(TypeError):
                x_invalid | y

    @prog_scope()
    def test_bitwise_xor(self):
        x_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        y_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        out_np = x_np ^ y_np

        x = paddle.static.data(name="x", shape=[2, 3, 5], dtype="int32")
        y = paddle.static.data(name="y", shape=[2, 3, 5], dtype="int32")
        z = x ^ y

        exe = base.Executor()
        out = exe.run(
            base.default_main_program(),
            feed={"x": x_np, "y": y_np},
            fetch_list=[z],
        )
        np.testing.assert_array_equal(out[0], out_np)

    @prog_scope()
    def test_rxor(self):
        place = get_device_place()
        x_int = 5
        y_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        y = paddle.static.data("y", y_np.shape, dtype=y_np.dtype)
        z = x_int ^ y
        exe = paddle.static.Executor(place)
        out = exe.run(
            feed={'y': y_np},
            fetch_list=[z],
        )
        out_ref = x_int ^ y_np
        np.testing.assert_array_equal(out[0], out_ref)
        x_bool = True
        res_rxor_bool = x_bool ^ y
        out_bool = exe.run(
            feed={'y': y_np},
            fetch_list=[res_rxor_bool],
        )
        res_py_bool = x_bool ^ y_np
        np.testing.assert_array_equal(out_bool[0], res_py_bool)

        for x_invalid in (
            np.float32(5.0),
            np.float64(5.0),
            np.complex64(5),
            np.complex128(5.0 + 2j),
        ):
            with self.assertRaises(TypeError):
                x_invalid ^ y

    @prog_scope()
    def test_bitwise_not(self):
        x_np = np.random.randint(-100, 100, [2, 3, 5]).astype("int32")
        out_np = ~x_np

        x = paddle.static.data(name="x", shape=[2, 3, 5], dtype="int32")
        z = ~x

        exe = base.Executor()
        out = exe.run(
            base.default_main_program(), feed={"x": x_np}, fetch_list=[z]
        )
        np.testing.assert_array_equal(out[0], out_np)

    @prog_scope()
    def test_T(self):
        x_np = np.random.randint(-100, 100, [2, 8, 5, 3]).astype("int32")
        out_np = x_np.T

        x = paddle.static.data(name="x", shape=[2, 8, 5, 3], dtype="int32")
        z = x.T

        exe = base.Executor()
        out = exe.run(
            base.default_main_program(), feed={"x": x_np}, fetch_list=[z]
        )
        np.testing.assert_array_equal(out[0], out_np)

    @prog_scope()
    def test_ndim(self):
        a = paddle.static.data(name="a", shape=[10, 1])
        self.assertEqual(a.dim(), 2)
        self.assertEqual(a.ndimension(), 2)
        self.assertEqual(a.ndim, 2)

    @prog_scope()
    def test_matmul(self):
        a = paddle.static.data(name='a', shape=[2, 3], dtype='float32')
        b = paddle.static.data(name='b', shape=[3, 5], dtype='float32')
        c = a @ b  # __matmul__
        a_np = np.random.uniform(-1, 1, size=[2, 3]).astype('float32')
        b_np = np.random.uniform(-1, 1, size=[3, 5]).astype('float32')
        place = paddle.CPUPlace()
        exe = paddle.static.Executor(place)
        (c_np,) = exe.run(
            paddle.static.default_main_program(),
            feed={"a": a_np, "b": b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np @ b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_rmatmul(self):
        a = paddle.static.data(name='a', shape=[2, 3], dtype='float32')
        b = paddle.static.data(name='b', shape=[3, 5], dtype='float32')
        c = b.__rmatmul__(a)
        a_np = np.random.uniform(-1, 1, size=[2, 3]).astype('float32')
        b_np = np.random.uniform(-1, 1, size=[3, 5]).astype('float32')
        place = paddle.CPUPlace()
        exe = paddle.static.Executor(place)
        (c_np,) = exe.run(
            paddle.static.default_main_program(),
            feed={"a": a_np, "b": b_np},
            fetch_list=[c],
        )
        np.testing.assert_allclose(a_np @ b_np, c_np, rtol=1e-05)

    @prog_scope()
    def test_builtin_type_conversion(self):
        a = paddle.static.data(name="a", shape=[])
        with self.assertRaises(TypeError):
            int(a)
        with self.assertRaises(TypeError):
            float(a)
        with self.assertRaises(TypeError):
            complex(a)


class CustomTensor:
    def __init__(self, value):
        assert isinstance(value, paddle.Tensor)

        self._value = value

    def __radd__(self, other):
        return other + self._value

    def __rsub__(self, other):
        return other - self._value

    def __rmul__(self, other):
        return other * self._value

    def __rdiv__(self, other):
        return other / self._value

    def __rfloordiv__(self, other):
        return other.__floordiv__(self._value)

    def __rtruediv__(self, other):
        return other.__truediv__(self._value)

    def __rpow__(self, other):
        return other.__pow__(self._value)

    def __rmatmul__(self, other):
        return other.__matmul__(self._value)

    def __rmod__(self, other):
        return other.__mod__(self._value)


class TestDygraphMathOpPatches(unittest.TestCase):
    def init_data(self):
        self.np_a = np.random.random((2, 3, 4)).astype(np.float32)
        self.np_b = np.random.random((2, 3, 4)).astype(np.float32)
        self.np_a[np.abs(self.np_a) < 0.0005] = 0.002
        self.np_b[np.abs(self.np_b) < 0.0005] = 0.002

        self.tensor_a = paddle.to_tensor(self.np_a, dtype="float32")
        self.tensor_b = paddle.to_tensor(self.np_b, dtype="float32")

    def test_dygraph_greater_than(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor > nparray
            expect_out = self.np_a > self.np_b
            actual_out = self.tensor_a > self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_greater_equal(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor >= nparray
            expect_out = self.np_a >= self.np_b
            actual_out = self.tensor_a >= self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_remainder(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor % nparray
            expect_out = self.np_a % self.np_b
            actual_out = self.tensor_a % self.np_b
            np.testing.assert_allclose(
                actual_out, expect_out, rtol=1e-7, atol=1e-7
            )

    def test_dygraph_rmod(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor % nparray
            expect_out = self.np_a % self.np_b
            actual_out = self.tensor_b.__rmod__(self.tensor_a)
            np.testing.assert_allclose(
                actual_out, expect_out, rtol=1e-7, atol=1e-7
            )

    def test_dygraph_less_than(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor < nparray
            expect_out = self.np_a < self.np_b
            actual_out = self.tensor_a < self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_less_equal(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor <= nparray
            expect_out = self.np_a <= self.np_b
            actual_out = self.tensor_a <= self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_floor_divide(self):
        with dygraph_guard():
            np_a = np.random.random((2, 3, 4)).astype(np.int32)
            np_b = np.random.random((2, 3, 4)).astype(np.int32)
            np_b[np.abs(np_b) < 1] = 2
            # normal case: tenor // nparray
            tensor_a = paddle.to_tensor(np_a, dtype="int32")
            tensor_b = paddle.to_tensor(np_b, dtype="int32")
            expect_out = np_a // np_b
            actual_out = tensor_a // np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_rfloordiv(self):
        with dygraph_guard():
            np_a = np.random.random((2, 3, 4)).astype(np.int32)
            np_b = np.random.random((2, 3, 4)).astype(np.int32)
            np_b[np.abs(np_b) < 1] = 2
            # normal case: nparray // tensor
            tensor_a = paddle.to_tensor(np_a, dtype="int32")
            tensor_b = paddle.to_tensor(np_b, dtype="int32")
            expect_out = np_b // np_a
            actual_out = tensor_b.__rfloordiv__(np_a)
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_elementwise_pow(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor ** nparray
            expect_out = self.np_a**self.np_b
            actual_out = self.tensor_a**self.np_b
            np.testing.assert_allclose(
                actual_out, expect_out, rtol=1e-7, atol=1e-7
            )

            # normal case: nparray ** tensor
            expect_out = self.np_a**self.np_b
            actual_out = self.np_a**self.tensor_b
            np.testing.assert_allclose(
                actual_out, expect_out, rtol=1e-7, atol=1e-7
            )

    def test_dygraph_not_equal(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor != nparray
            expect_out = self.np_a != self.np_b
            actual_out = self.tensor_a != self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_equal(self):
        with dygraph_guard():
            self.init_data()
            # normal case: tenor == nparray
            expect_out = self.np_a == self.np_b
            actual_out = self.tensor_a == self.np_b
            np.testing.assert_equal(actual_out, expect_out)

    def test_dygraph_rmatmul(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((3, 5)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = b.__rmatmul__(a)
            np.testing.assert_allclose(a @ b, c.numpy(), rtol=1e-5, atol=1e-5)

    # @unittest.skipIf(
    #     paddle.device.is_compiled_with_xpu(),
    #     reason="XPU not support custom tensor.",
    # )
    def test_dygraph_custom_add(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a + CustomTensor(b)
            np.testing.assert_allclose((a + b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_add(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__add__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a + b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a + a.T

    # @unittest.skipIf(
    #     paddle.device.is_compiled_with_xpu(),
    #     reason="XPU not support custom tensor.",
    # )
    def test_dygraph_custom_mul(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a * CustomTensor(b)
            np.testing.assert_allclose((a * b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_mul(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__mul__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a * b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a * a.T

    def test_dygraph_custom_sub(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a - CustomTensor(b)
            np.testing.assert_allclose((a - b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_sub(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__sub__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a - b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a - a.T

    def test_dygraph_custom_pow(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a ** CustomTensor(b)
            np.testing.assert_allclose((a**b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_pow(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__pow__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a**b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a**a.T

    def test_dygraph_custom_mod(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a % CustomTensor(b)
            np.testing.assert_allclose((a % b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_mod(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__mod__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a % b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a % a.T

    def test_dygraph_custom_matmul(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((3, 2)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a @ CustomTensor(b)
            np.testing.assert_allclose((a @ b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_matmul(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__matmul__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a @ b_f

            with self.assertRaises(ValueError):
                _ = a @ a

    def test_dygraph_custom_div(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a / CustomTensor(b)
            np.testing.assert_allclose((a / b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_div(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__div__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a / b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a / a.T

    def test_dygraph_custom_truediv(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a / CustomTensor(b)
            np.testing.assert_allclose((a / b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_truediv(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__div__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a / b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a / a.T

    def test_dygraph_custom_floordiv(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_np = np.random.random((2, 3)).astype(np.float32) * 100
            a = paddle.to_tensor(a_np)
            b = paddle.to_tensor(b_np)
            c = a // CustomTensor(b)
            np.testing.assert_allclose((a // b).numpy(), c.numpy(), atol=0)

    def test_dygraph_vanilla_floordiv(self):
        with dygraph_guard():
            a_np = np.random.random((2, 3)).astype(np.float32) * 100
            b_f = lambda y: y + 1
            a = paddle.to_tensor(a_np)
            with self.assertRaisesRegex(
                TypeError,
                r"__floordiv__\(\): argument \(position 1\) must be int, float, bool or Tensor, but got ",
            ):
                _ = a // b_f

            with self.assertRaisesRegex(
                ValueError, "Broadcast dimension mismatch"
            ):
                _ = a // a.T


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