Files
2026-07-13 12:40:42 +08:00

734 lines
22 KiB
Python

# 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 op_test import (
OpTest,
convert_float_to_uint16,
convert_uint16_to_float,
get_device_place,
get_numeric_gradient,
is_custom_device,
)
from testsuite import create_op
import paddle
from paddle.base import core
def adaptive_start_index(index, input_size, output_size):
return int(np.floor(index * input_size / output_size))
def adaptive_end_index(index, input_size, output_size):
return int(np.ceil((index + 1) * input_size / output_size))
def max_pool3D_forward_naive(
x, ksize, strides, paddings, global_pool=False, adaptive=False
):
N, C, D, H, W = x.shape
if global_pool:
ksize = [D, H, W]
paddings = [0, 0, 0]
if adaptive:
D_out, H_out, W_out = ksize
else:
D_out = (D - ksize[0] + 2 * paddings[0]) // strides[0] + 1
H_out = (H - ksize[1] + 2 * paddings[1]) // strides[1] + 1
W_out = (W - ksize[2] + 2 * paddings[2]) // strides[2] + 1
out = np.zeros((N, C, D_out, H_out, W_out))
mask = np.zeros((N, C, D_out, H_out, W_out))
for k in range(D_out):
if adaptive:
d_start = adaptive_start_index(k, D, ksize[0])
d_end = adaptive_end_index(k, D, ksize[0])
else:
d_start = np.max((k * strides[0] - paddings[0], 0))
d_end = np.min((k * strides[0] + ksize[0] - paddings[0], D))
for i in range(H_out):
if adaptive:
h_start = adaptive_start_index(i, H, ksize[1])
h_end = adaptive_end_index(i, H, ksize[1])
else:
h_start = np.max((i * strides[1] - paddings[1], 0))
h_end = np.min((i * strides[1] + ksize[1] - paddings[1], H))
for j in range(W_out):
if adaptive:
w_start = adaptive_start_index(j, W, ksize[2])
w_end = adaptive_end_index(j, W, ksize[2])
else:
w_start = np.max((j * strides[2] - paddings[2], 0))
w_end = np.min((j * strides[2] + ksize[2] - paddings[2], W))
x_masked = x[:, :, d_start:d_end, h_start:h_end, w_start:w_end]
out[:, :, k, i, j] = np.max(x_masked, axis=(2, 3, 4))
for n in range(N):
for c in range(C):
arr = x_masked[n, c, :, :, :]
index = np.where(arr == np.max(arr))
sub_deep = index[0][0]
sub_row = index[1][0]
sub_col = index[2][0]
index = (
((d_start + sub_deep) * H + (h_start + sub_row)) * W
+ w_start
+ sub_col
)
mask[n, c, k, i, j] = index
return out, mask
def max_pool2D_forward_naive(
x, ksize, strides, paddings, global_pool=False, adaptive=False
):
N, C, H, W = x.shape
if global_pool:
ksize = [H, W]
paddings = [0, 0]
if adaptive:
H_out, W_out = ksize
else:
H_out = (H - ksize[0] + 2 * paddings[0]) // strides[0] + 1
W_out = (W - ksize[1] + 2 * paddings[1]) // strides[1] + 1
out = np.zeros((N, C, H_out, W_out))
mask = np.zeros((N, C, H_out, W_out))
for i in range(H_out):
for j in range(W_out):
if adaptive:
r_start = adaptive_start_index(i, H, ksize[0])
r_end = adaptive_end_index(i, H, ksize[0])
c_start = adaptive_start_index(j, W, ksize[1])
c_end = adaptive_end_index(j, W, ksize[1])
else:
r_start = np.max((i * strides[0] - paddings[0], 0))
r_end = np.min((i * strides[0] + ksize[0] - paddings[0], H))
c_start = np.max((j * strides[1] - paddings[1], 0))
c_end = np.min((j * strides[1] + ksize[1] - paddings[1], W))
x_masked = x[:, :, r_start:r_end, c_start:c_end]
out[:, :, i, j] = np.max(x_masked, axis=(2, 3))
for n in range(N):
for c in range(C):
arr = x_masked[n, c, :, :]
index = np.where(arr == np.max(arr))
sub_row = index[0][0]
sub_col = index[1][0]
index = (r_start + sub_row) * W + c_start + sub_col
mask[n, c, i, j] = index
return out, mask
def max_pool3d_with_index_wrapper(
x,
kernel_size=[],
strides=[],
paddings=[],
global_pooling=False,
adaptive=False,
ceil_mode=False,
):
dilations = [1, 1, 1]
return paddle._C_ops.max_pool3d_with_index(
x,
kernel_size,
strides,
paddings,
dilations,
global_pooling,
adaptive,
ceil_mode,
)
class TestMaxPoolWithIndex_Op(OpTest):
def setUp(self):
self.init_test_case()
self.init_global()
self.init_adaptive()
self.init_dtype()
if self.is_bfloat16_op():
input = np.random.random(self.shape).astype(np.float32)
input = convert_uint16_to_float(
convert_float_to_uint16(np.round(input * 100.0, 2))
)
else:
input = np.random.random(self.shape).astype(self.dtype)
input = np.round(input * 100.0, 2)
output, mask = self.pool_forward_naive(
input,
self.ksize,
self.strides,
self.paddings,
self.global_pool,
self.adaptive,
)
mask = mask.astype("int32")
if self.is_bfloat16_op():
output = output.astype(np.float32)
else:
output = output.astype(self.dtype)
self.attrs = {
'strides': self.strides,
'paddings': self.paddings,
'ksize': self.ksize,
'global_pooling': self.global_pool,
'adaptive': self.adaptive,
}
if self.is_bfloat16_op():
self.inputs = {'X': convert_float_to_uint16(input)}
self.outputs = {
'Out': convert_float_to_uint16(output),
"Mask": mask,
}
self.inputs_fp32 = {'X': input}
else:
self.inputs = {'X': input}
self.outputs = {'Out': output, "Mask": mask}
def init_dtype(self):
self.dtype = np.float64
def test_check_output(self):
self.check_output()
def test_check_grad(self):
self.check_grad({'X'}, ['Out'])
def init_test_case(self):
self.op_type = "max_pool3d_with_index"
self.python_api = max_pool3d_with_index_wrapper
self.pool_forward_naive = max_pool3D_forward_naive
self.shape = [2, 3, 7, 7, 7]
self.ksize = [3, 3, 3]
self.strides = [2, 2, 2]
self.paddings = [1, 1, 1]
def init_global(self):
self.global_pool = False
def init_adaptive(self):
self.adaptive = False
class TestCase1(TestMaxPoolWithIndex_Op):
def init_global(self):
self.global_pool = True
class TestCase2(TestMaxPoolWithIndex_Op):
def init_test_case(self):
self.op_type = "max_pool3d_with_index"
self.python_api = max_pool3d_with_index_wrapper
self.pool_forward_naive = max_pool3D_forward_naive
self.shape = [2, 3, 7, 7, 7]
self.ksize = [3, 3, 3]
self.strides = [2, 2, 2]
self.paddings = [0, 0, 0]
def init_global(self):
self.global_pool = True
class TestCase3(TestCase2):
def init_global(self):
self.global_pool = False
class TestCastAdaptive3d(TestMaxPoolWithIndex_Op):
def init_adaptive(self):
self.adaptive = True
# ----------------max_pool3d_with_index_fp16----------------
def create_test_fp16_class(parent):
@unittest.skipIf(
not (core.is_compiled_with_cuda() or is_custom_device()),
"core is not compiled with CUDA",
)
class TestMaxPool3dFP16(parent):
def init_dtype(self):
self.dtype = np.float16
def test_check_output(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
if core.is_float16_supported(place):
self.check_output_with_place(place)
def test_check_grad(self):
place = get_device_place()
if core.is_float16_supported(place):
self.check_grad_with_place(place, {'X'}, ['Out'])
cls_name = "{}_{}".format(parent.__name__, "FP16OP")
TestMaxPool3dFP16.__name__ = cls_name
globals()[cls_name] = TestMaxPool3dFP16
create_test_fp16_class(TestMaxPoolWithIndex_Op)
create_test_fp16_class(TestCase1)
create_test_fp16_class(TestCase2)
create_test_fp16_class(TestCase3)
create_test_fp16_class(TestCastAdaptive3d)
# ----------------max_pool3d_with_index_bf16----------------
def create_test_bf16_class(parent):
@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 do not support bfloat16",
)
class TestMaxPool3dBF16(parent):
def init_dtype(self):
self.dtype = np.uint16
def get_numeric_grad(self, place, check_name):
scope = core.Scope()
self._check_grad_helper()
op = create_op(
scope, self.op_type, self.inputs, self.outputs, self.attrs
)
return get_numeric_gradient(
place, scope, op, self.inputs_fp32, check_name, ['Out']
)
def test_check_output(self):
place = get_device_place()
if core.is_bfloat16_supported(place):
self.check_output_with_place(place)
def test_check_grad(self):
place = get_device_place()
numeric_grads = self.get_numeric_grad(place, 'X')
if core.is_bfloat16_supported(place):
self.check_grad_with_place(
place,
{'X'},
['Out'],
)
cls_name = "{}_{}".format(parent.__name__, "BF16OP")
TestMaxPool3dBF16.__name__ = cls_name
globals()[cls_name] = TestMaxPool3dBF16
create_test_bf16_class(TestMaxPoolWithIndex_Op)
create_test_bf16_class(TestCase1)
create_test_bf16_class(TestCase2)
create_test_bf16_class(TestCase3)
create_test_bf16_class(TestCastAdaptive3d)
# ----------------max_pool2d_with_index----------------
def max_pool2d_with_index_wrapper(
x,
kernel_size=[],
strides=[],
paddings=[],
global_pooling=False,
adaptive=False,
ceil_mode=False,
):
dilations = [1, 1]
return paddle._C_ops.max_pool2d_with_index(
x,
kernel_size,
strides,
paddings,
dilations,
global_pooling,
adaptive,
ceil_mode,
)
class TestCase4(TestMaxPoolWithIndex_Op):
def init_test_case(self):
self.op_type = "max_pool2d_with_index"
self.python_api = max_pool2d_with_index_wrapper
self.pool_forward_naive = max_pool2D_forward_naive
self.shape = [2, 3, 7, 7]
self.ksize = [3, 3]
self.strides = [1, 1]
self.paddings = [1, 1]
def init_global(self):
self.global_pool = True
class TestCase5(TestCase4):
def init_global(self):
self.global_pool = False
class TestCase6(TestMaxPoolWithIndex_Op):
def init_test_case(self):
self.op_type = "max_pool2d_with_index"
self.python_api = max_pool2d_with_index_wrapper
self.pool_forward_naive = max_pool2D_forward_naive
self.shape = [2, 3, 7, 7]
self.ksize = [3, 3]
self.strides = [2, 2]
self.paddings = [0, 0]
def init_global(self):
self.global_pool = True
class TestCase7(TestCase6):
def init_global(self):
self.global_pool = False
class TestCastAdaptive2d(TestCase6):
def init_adaptive(self):
self.adaptive = True
# ----------------max_pool2d_with_index_fp16----------------
def create_test_fp16_class(parent):
@unittest.skipIf(
not (core.is_compiled_with_cuda() or is_custom_device()),
"core is not compiled with CUDA",
)
class TestMaxPool2dFP16(parent):
def init_dtype(self):
self.dtype = np.float16
def test_check_output(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
if core.is_float16_supported(place):
self.check_output_with_place(place)
def test_check_grad(self):
place = get_device_place()
if core.is_float16_supported(place):
self.check_grad_with_place(place, {'X'}, ['Out'])
cls_name = "{}_{}".format(parent.__name__, "FP16OP")
TestMaxPool2dFP16.__name__ = cls_name
globals()[cls_name] = TestMaxPool2dFP16
create_test_fp16_class(TestCase4)
create_test_fp16_class(TestCase5)
create_test_fp16_class(TestCase6)
create_test_fp16_class(TestCase7)
create_test_fp16_class(TestCastAdaptive2d)
# ----------------max_pool2d_with_index_bf16----------------
def create_test_bf16_class(parent):
@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 do not support bfloat16",
)
class TestMaxPool2dBF16(parent):
def init_dtype(self):
self.dtype = np.uint16
def get_numeric_grad(self, place, check_name):
scope = core.Scope()
self._check_grad_helper()
op = create_op(
scope, self.op_type, self.inputs, self.outputs, self.attrs
)
return get_numeric_gradient(
place, scope, op, self.inputs_fp32, check_name, ['Out']
)
def test_check_output(self):
place = get_device_place()
if core.is_bfloat16_supported(place):
self.check_output_with_place(place)
def test_check_grad(self):
place = get_device_place()
numeric_grads = self.get_numeric_grad(place, 'X')
if core.is_bfloat16_supported(place):
self.check_grad_with_place(
place, {'X'}, ['Out'], user_defined_grads=[numeric_grads]
)
cls_name = "{}_{}".format(parent.__name__, "BF16OP")
TestMaxPool2dBF16.__name__ = cls_name
globals()[cls_name] = TestMaxPool2dBF16
create_test_bf16_class(TestCase4)
create_test_bf16_class(TestCase5)
create_test_bf16_class(TestCase6)
create_test_bf16_class(TestCase7)
create_test_bf16_class(TestCastAdaptive2d)
def skip_unit_test():
if is_custom_device():
return False
return (
not core.is_compiled_with_cuda()
or not core.is_compiled_with_cudnn_frontend()
or paddle.device.cuda.get_device_capability()[0] < 8
)
@unittest.skipIf(
skip_unit_test(),
"Only support Ampere or later devices; "
"Paddle should be built with WITH_CUDNN_FRONTEND=ON.",
)
class TestMaxPool2dV2Op(OpTest):
def setUp(self):
self.init_layout()
self.init_test_case()
self.init_global()
self.init_adaptive()
self.init_dtype()
if self.is_bfloat16_op():
input = np.random.random(self.shape).astype(np.float32)
input = convert_uint16_to_float(
convert_float_to_uint16(np.round(input * 100.0, 2))
)
else:
input = np.random.random(self.shape).astype(self.dtype)
input = np.round(input * 100.0, 2)
output, _ = self.pool_forward_naive(
input,
self.ksize,
self.strides,
self.paddings,
self.global_pool,
self.adaptive,
)
if self.is_bfloat16_op():
output = output.astype(np.float32)
else:
output = output.astype(self.dtype)
self.attrs = {
'strides': self.strides,
'paddings': self.paddings,
'kernel_size': self.ksize,
'data_format': self.data_format,
'global_pooling': self.global_pool,
'adaptive': self.adaptive,
}
if self.data_format == 'NHWC':
input = input.transpose((0, 2, 3, 1))
output = output.transpose((0, 2, 3, 1))
saved_idx = np.zeros(shape=output.shape, dtype=np.int32)
if self.is_bfloat16_op():
self.inputs = {
'x': convert_float_to_uint16(
input, data_format=self.data_format
)
}
self.outputs = {
'out': convert_float_to_uint16(
output, data_format=self.data_format
),
'saved_idx': saved_idx,
}
self.inputs_fp32 = {'x': input}
else:
self.inputs = {'x': input}
self.outputs = {'out': output, 'saved_idx': saved_idx}
def init_layout(self):
self.data_format = "NHWC"
def init_dtype(self):
self.dtype = np.float32
def test_check_output(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
self.check_output_with_place(
place, no_check_set=['saved_idx'], check_dygraph=False
)
def test_check_grad(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
self.check_grad_with_place(
place,
{'x'},
['out'],
max_relative_error=0.05,
check_dygraph=False,
)
def init_test_case(self):
self.op_type = "max_pool2d_v2"
self.pool_forward_naive = max_pool2D_forward_naive
self.shape = [2, 3, 7, 7]
self.ksize = [3, 3]
self.strides = [1, 1]
self.paddings = [1, 1]
def init_global(self):
self.global_pool = True
def init_adaptive(self):
self.adaptive = False
class TestCase8(TestMaxPool2dV2Op):
def init_global(self):
self.global_pool = False
def test_check_grad(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
self.check_grad_with_place(
place,
{'x'},
['out'],
max_relative_error=0.5,
check_dygraph=False,
)
class TestCase9(TestMaxPool2dV2Op):
def init_test_case(self):
self.op_type = "max_pool2d_v2"
self.pool_forward_naive = max_pool2D_forward_naive
self.shape = [2, 3, 7, 7]
self.ksize = [3, 3]
self.strides = [2, 2]
self.paddings = [0, 0]
def init_global(self):
self.global_pool = True
class TestCase10(TestCase9):
def init_global(self):
self.global_pool = False
def create_test_fp16_class(parent):
class TestMaxPool2dV2FP16(parent):
def init_dtype(self):
self.dtype = np.float16
def test_check_output(self):
if core.is_compiled_with_cuda() or is_custom_device():
place = get_device_place()
if core.is_float16_supported(place):
self.check_output_with_place(
place, no_check_set=['saved_idx'], check_dygraph=False
)
def test_check_grad(self):
place = get_device_place()
if core.is_float16_supported(place):
self.check_grad_with_place(
place, {'x'}, ['out'], check_dygraph=False
)
cls_name = "{}_{}".format(parent.__name__, "FP16OP")
TestMaxPool2dV2FP16.__name__ = cls_name
globals()[cls_name] = TestMaxPool2dV2FP16
create_test_fp16_class(TestMaxPool2dV2Op)
create_test_fp16_class(TestCase8)
create_test_fp16_class(TestCase9)
create_test_fp16_class(TestCase10)
def create_test_bf16_class(parent):
@unittest.skipIf(
skip_unit_test() or not core.is_bfloat16_supported(get_device_place()),
"core is not compiled with CUDA and do not support bfloat16",
)
class TestMaxPool2dV2BF16(parent):
def init_dtype(self):
self.dtype = np.uint16
def get_numeric_grad(self, place, check_name):
scope = core.Scope()
self._check_grad_helper()
op = create_op(
scope, self.op_type, self.inputs, self.outputs, self.attrs
)
return get_numeric_gradient(
place,
scope,
op,
self.inputs_fp32,
check_name,
['out'],
delta=0.005,
)
def test_check_output(self):
place = get_device_place()
if core.is_bfloat16_supported(place):
self.check_output_with_place(
place, no_check_set=['saved_idx'], check_dygraph=False
)
def test_check_grad(self):
place = get_device_place()
numeric_grads = self.get_numeric_grad(place, 'x')
if core.is_bfloat16_supported(place):
self.check_grad_with_place(
place,
{'x'},
['out'],
user_defined_grads=[numeric_grads],
check_dygraph=False,
)
cls_name = "{}_{}".format(parent.__name__, "BF16OP")
TestMaxPool2dV2BF16.__name__ = cls_name
globals()[cls_name] = TestMaxPool2dV2BF16
create_test_bf16_class(TestMaxPool2dV2Op)
create_test_bf16_class(TestCase8)
create_test_bf16_class(TestCase9)
create_test_bf16_class(TestCase10)
if __name__ == '__main__':
unittest.main()