paddlepaddle--paddle/test/fft/test_spectral_op.py

# Copyright (c) 2021 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 re
import sys
import unittest

import numpy as np
from op_test import OpTest

sys.path.append("../../fft")
from spectral_op_np import fft_c2c, fft_c2r, fft_r2c

import paddle
from paddle import _C_ops

paddle.enable_static()

TEST_CASE_NAME = 'test_case'


def parameterize(attrs, input_values=None):
    if isinstance(attrs, str):
        attrs = [attrs]
    input_dicts = (
        attrs
        if input_values is None
        else [dict(zip(attrs, vals)) for vals in input_values]
    )

    def decorator(base_class):
        test_class_module = sys.modules[base_class.__module__].__dict__
        for idx, input_dict in enumerate(input_dicts):
            test_class_dict = dict(base_class.__dict__)
            test_class_dict.update(input_dict)

            name = class_name(base_class, idx, input_dict)

            test_class_module[name] = type(name, (base_class,), test_class_dict)

        for method_name in list(base_class.__dict__):
            if method_name.startswith("test"):
                delattr(base_class, method_name)
        return base_class

    return decorator


def to_safe_name(s):
    return str(re.sub("[^a-zA-Z0-9_]+", "_", s))


def class_name(cls, num, params_dict):
    suffix = to_safe_name(
        next((v for v in params_dict.values() if isinstance(v, str)), "")
    )
    if TEST_CASE_NAME in params_dict:
        suffix = to_safe_name(params_dict["test_case"])
    return "{}_{}{}".format(cls.__name__, num, suffix and "_" + suffix)


def fft_c2c_python_api(x, axes, norm, forward):
    return _C_ops.fft_c2c(x, axes, norm, forward)


def fft_r2c_python_api(x, axes, norm, forward, onesided):
    return _C_ops.fft_r2c(x, axes, norm, forward, onesided)


def fft_c2r_python_api(x, axes, norm, forward, last_dim_size=0):
    return _C_ops.fft_c2r(x, axes, norm, forward, last_dim_size)


@parameterize(
    (TEST_CASE_NAME, 'x', 'axes', 'norm', 'forward'),
    [
        (
            'test_axes_is_sqe_type',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            [0, 1],
            'forward',
            True,
        ),
        (
            'test_axis_not_last',
            (
                np.random.random((4, 8, 4)) + 1j * np.random.random((4, 8, 4))
            ).astype(np.complex128),
            (0, 1),
            "backward",
            False,
        ),
        (
            'test_norm_forward',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (0,),
            "forward",
            False,
        ),
        (
            'test_norm_backward',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (0,),
            "backward",
            True,
        ),
        (
            'test_norm_ortho',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (1,),
            "ortho",
            True,
        ),
    ],
)
class TestFFTC2COp(OpTest):
    def setUp(self):
        self.op_type = "fft_c2c"
        self.dtype = self.x.dtype
        self.python_api = fft_c2c_python_api

        out = fft_c2c(self.x, self.axes, self.norm, self.forward)

        self.inputs = {'X': self.x}
        self.attrs = {
            'axes': self.axes,
            'normalization': self.norm,
            "forward": self.forward,
        }
        self.outputs = {'Out': out}

    def test_check_output(self):
        self.check_output()

    def test_check_grad(self):
        self.check_grad(
            ["X"],
            "Out",
        )


@parameterize(
    (TEST_CASE_NAME, 'x', 'axes', 'norm', 'forward', 'last_dim_size'),
    [
        (
            'test_axes_is_sqe_type',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            [0, 1],
            'forward',
            False,
            26,
        ),
        (
            'test_axis_not_last',
            (
                np.random.random((4, 7, 4)) + 1j * np.random.random((4, 7, 4))
            ).astype(np.complex128),
            (0, 1),
            "backward",
            False,
            None,
        ),
        (
            'test_norm_forward',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (0,),
            "forward",
            False,
            22,
        ),
        (
            'test_norm_backward',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (0,),
            "backward",
            False,
            22,
        ),
        (
            'test_norm_ortho',
            (
                np.random.random((12, 14)) + 1j * np.random.random((12, 14))
            ).astype(np.complex128),
            (1,),
            "ortho",
            True,
            26,
        ),
    ],
)
class TestFFTC2ROp(OpTest):
    def setUp(self):
        self.op_type = "fft_c2r"
        self.dtype = self.x.dtype
        self.python_api = fft_c2r_python_api

        out = fft_c2r(
            self.x, self.axes, self.norm, self.forward, self.last_dim_size
        )

        self.inputs = {'X': self.x}
        self.attrs = {
            "axes": self.axes,
            "normalization": self.norm,
            "forward": self.forward,
            "last_dim_size": self.last_dim_size,
        }
        self.outputs = {'Out': out}

    def test_check_output(self):
        self.check_output()

    def test_check_grad(self):
        self.check_grad(
            ["X"],
            "Out",
        )


@parameterize(
    (TEST_CASE_NAME, 'x', 'axes', 'norm', 'forward', 'onesided'),
    [
        (
            'test_axes_is_sqe_type',
            np.random.randn(12, 18).astype(np.float64),
            (0, 1),
            'forward',
            True,
            True,
        ),
        (
            'test_axis_not_last',
            np.random.randn(4, 8, 4).astype(np.float64),
            (0, 1),
            "backward",
            False,
            False,
        ),
        (
            'test_norm_forward',
            np.random.randn(12, 18).astype(np.float64),
            (0, 1),
            "forward",
            False,
            False,
        ),
        (
            'test_norm_backward',
            np.random.randn(12, 18).astype(np.float64),
            (0,),
            "backward",
            True,
            False,
        ),
        (
            'test_norm_ortho',
            np.random.randn(12, 18).astype(np.float64),
            (1,),
            "ortho",
            True,
            False,
        ),
    ],
)
class TestFFTR2COp(OpTest):
    def setUp(self):
        self.op_type = "fft_r2c"
        self.dtype = self.x.dtype
        self.python_api = fft_r2c_python_api

        out = fft_r2c(self.x, self.axes, self.norm, self.forward, self.onesided)

        self.inputs = {'X': self.x}
        self.attrs = {
            'axes': self.axes,
            'normalization': self.norm,
            "forward": self.forward,
            'onesided': self.onesided,
        }
        self.outputs = {'Out': out}

    def test_check_output(self):
        self.check_output()

    def test_check_grad(self):
        self.check_grad(
            ["X"],
            "Out",
        )


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