Files
paddlepaddle--paddle/test/legacy_test/test_fill_diagonal_tensor_op.py
2026-07-13 12:40:42 +08:00

211 lines
6.4 KiB
Python

# Copyright (c) 2019 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.base import core
def fill_diagonal_ndarray(x, value, offset=0, dim1=0, dim2=1):
"""Fill value into the diagonal of x that offset is ${offset} and the coordinate system is (dim1, dim2)."""
strides = x.strides
shape = x.shape
if dim1 > dim2:
dim1, dim2 = dim2, dim1
assert 0 <= dim1 < dim2 <= 2
assert len(x.shape) == 3
dim_sum = dim1 + dim2
dim3 = len(x.shape) - dim_sum
if offset >= 0:
diagdim = min(shape[dim1], shape[dim2] - offset)
diagonal = np.lib.stride_tricks.as_strided(
x[:, offset:] if dim_sum == 1 else x[:, :, offset:],
shape=(shape[dim3], diagdim),
strides=(strides[dim3], strides[dim1] + strides[dim2]),
)
else:
diagdim = min(shape[dim2], shape[dim1] + offset)
diagonal = np.lib.stride_tricks.as_strided(
x[-offset:, :] if dim_sum in [1, 2] else x[:, -offset:],
shape=(shape[dim3], diagdim),
strides=(strides[dim3], strides[dim1] + strides[dim2]),
)
diagonal[...] = value
return x
def fill_gt(x, y, offset, dim1, dim2):
if dim1 > dim2:
dim1, dim2 = dim2, dim1
offset = -offset
xshape = x.shape
yshape = y.shape
if len(xshape) != 3:
perm_list = []
unperm_list = [0] * len(xshape)
idx = 0
for i in range(len(xshape)):
if i != dim1 and i != dim2:
perm_list.append(i)
unperm_list[i] = idx
idx += 1
perm_list += [dim1, dim2]
unperm_list[dim1] = idx
unperm_list[dim2] = idx + 1
x = np.transpose(x, perm_list)
y = y.reshape(-1, yshape[-1])
nxshape = x.shape
x = x.reshape((-1, xshape[dim1], xshape[dim2]))
out = fill_diagonal_ndarray(x, y, offset, 1, 2)
if len(xshape) != 3:
out = out.reshape(nxshape)
out = np.transpose(out, unperm_list)
return out
class TensorFillDiagTensor_Test(OpTest):
def setUp(self):
self.op_type = "fill_diagonal_tensor"
self.python_api = paddle.tensor.manipulation.fill_diagonal_tensor
self.init_kernel_type()
x = np.random.random((10, 10)).astype(self.dtype)
y = np.random.random((10,)).astype(self.dtype)
dim1 = 0
dim2 = 1
offset = 0
out = fill_gt(x, y, offset, dim1, dim2)
self.inputs = {"X": x, "Y": y}
self.outputs = {'Out': out}
self.attrs = {"offset": offset, "dim1": dim1, "dim2": dim2}
def init_kernel_type(self):
self.dtype = np.float64
def test_check_output(self):
self.check_output(check_pir=True)
def test_check_grad(self):
self.check_grad(['X'], 'Out', check_pir=True)
class TensorFillDiagTensor_Test2(TensorFillDiagTensor_Test):
def setUp(self):
self.op_type = "fill_diagonal_tensor"
self.python_api = paddle.tensor.manipulation.fill_diagonal_tensor
self.init_kernel_type()
x = np.random.random((2, 20, 25)).astype(self.dtype)
y = np.random.random((2, 20)).astype(self.dtype)
dim1 = 2
dim2 = 1
offset = -3
out = fill_gt(x, y, offset, dim1, dim2)
self.inputs = {"X": x, "Y": y}
self.outputs = {'Out': out}
self.attrs = {"offset": offset, "dim1": dim1, "dim2": dim2}
def init_kernel_type(self):
self.dtype = np.float32
class TensorFillDiagTensor_Test3(TensorFillDiagTensor_Test):
def setUp(self):
self.op_type = "fill_diagonal_tensor"
self.python_api = paddle.tensor.manipulation.fill_diagonal_tensor
self.init_kernel_type()
x = np.random.random((2, 20, 20, 3)).astype(self.dtype)
y = np.random.random((2, 3, 18)).astype(self.dtype)
dim1 = 1
dim2 = 2
offset = 2
out = fill_gt(x, y, offset, dim1, dim2)
self.inputs = {"X": x, "Y": y}
self.outputs = {'Out': out}
self.attrs = {"offset": offset, "dim1": dim1, "dim2": dim2}
def init_kernel_type(self):
self.dtype = np.float16
class TensorFillDiagTensorFP16OP(TensorFillDiagTensor_Test):
def init_kernel_type(self):
self.dtype = np.float16
@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 TensorFillDiagTensorBF16(OpTest):
def setUp(self):
self.op_type = "fill_diagonal_tensor"
self.python_api = paddle.tensor.manipulation.fill_diagonal_tensor
self.init_kernel_type()
self.init_config()
self.init_input_output()
def init_kernel_type(self):
self.dtype = np.uint16
def init_config(self):
self.x = np.random.random((10, 10)).astype(np.float32)
self.y = np.random.random((10,)).astype(np.float32)
self.dim1 = 0
self.dim2 = 1
self.offset = 0
def init_input_output(self):
out = fill_gt(self.x, self.y, self.offset, self.dim1, self.dim2)
self.inputs = {
"X": convert_float_to_uint16(self.x),
"Y": convert_float_to_uint16(self.y),
}
self.outputs = {'Out': convert_float_to_uint16(out)}
self.attrs = {
"offset": self.offset,
"dim1": self.dim1,
"dim2": self.dim2,
}
def test_check_output(self):
place = get_device_place()
self.check_output_with_place(place, check_pir=True)
def test_check_grad(self):
place = get_device_place()
self.check_grad_with_place(place, ['X'], 'Out', check_pir=True)
if __name__ == '__main__':
paddle.enable_static()
unittest.main()