284 lines
8.3 KiB
Python
284 lines
8.3 KiB
Python
# 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.base import core
|
|
|
|
np.random.seed(10)
|
|
|
|
|
|
def logit(x, eps):
|
|
if eps:
|
|
x_min = np.minimum(x, 1.0 - eps)
|
|
x_max = np.maximum(x_min, eps)
|
|
return np.log(x_max / (1.0 - x_max))
|
|
else:
|
|
return np.where(
|
|
(x < 0.0) | (x > 1.0),
|
|
np.array(np.nan, dtype=x.dtype),
|
|
np.log(x / (1.0 - x)),
|
|
)
|
|
|
|
|
|
def logit_grad(x, eps=1e-8):
|
|
if eps:
|
|
tmp_x = np.select(
|
|
[x < eps, x > (1.0 - eps)], [x * 0.0, x * 0.0], default=-1.0
|
|
)
|
|
x_1 = 1.0 - x
|
|
_x = np.select([tmp_x == -1.0], [np.reciprocal(x * x_1)], default=0.0)
|
|
else:
|
|
tmp_x = np.select(
|
|
[x < 0.0, x > 1.0],
|
|
[np.array(np.nan, dtype=x.dtype), np.array(np.nan, dtype=x.dtype)],
|
|
default=-1.0,
|
|
)
|
|
x_1 = 1.0 - x
|
|
_x = np.select(
|
|
[tmp_x == -1.0],
|
|
[np.reciprocal(x * x_1)],
|
|
default=np.array(np.nan, dtype=x.dtype),
|
|
)
|
|
|
|
if _x.size == 0:
|
|
dout = np.full_like(x, fill_value=0.0)
|
|
else:
|
|
dout = np.full_like(x, fill_value=1.0 / _x.size)
|
|
dx = dout * _x
|
|
return dx
|
|
|
|
|
|
class TestLogitOp(OpTest):
|
|
def setUp(self):
|
|
self.op_type = 'logit'
|
|
self.python_api = paddle.logit
|
|
self.set_attrs()
|
|
x = np.random.uniform(-1.0, 1.0, self.shape).astype(self.dtype)
|
|
out = logit(x, self.eps)
|
|
self.x_grad = logit_grad(x, self.eps)
|
|
self.inputs = {'X': x}
|
|
self.outputs = {'Out': out}
|
|
self.attrs = {'eps': self.eps}
|
|
|
|
def set_attrs(self):
|
|
self.dtype = np.float64
|
|
self.shape = [120]
|
|
self.eps = 1e-8
|
|
|
|
def test_check_output(self):
|
|
self.check_output(check_pir=True, check_symbol_infer=False)
|
|
|
|
def test_check_grad(self):
|
|
self.check_grad(
|
|
['X'],
|
|
['Out'],
|
|
user_defined_grads=[self.x_grad],
|
|
check_pir=True,
|
|
)
|
|
|
|
|
|
class TestLogitOpFp32(TestLogitOp):
|
|
def set_attrs(self):
|
|
self.dtype = np.float32
|
|
self.shape = [120]
|
|
self.eps = 1e-8
|
|
|
|
def test_check_output(self):
|
|
self.check_output(check_pir=True)
|
|
|
|
|
|
class TestLogitOpFp16(TestLogitOp):
|
|
def set_attrs(self):
|
|
self.dtype = np.float16
|
|
self.shape = [120]
|
|
self.eps = 1e-8
|
|
|
|
def test_check_output(self):
|
|
self.check_output(check_pir=True)
|
|
|
|
|
|
@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 TestLogitOpBf16(OpTest):
|
|
def setUp(self):
|
|
self.op_type = 'logit'
|
|
self.python_api = paddle.logit
|
|
self.set_attrs()
|
|
x = np.random.uniform(-0.5, 0.5, self.shape).astype(np.float32)
|
|
out = logit(x, self.eps)
|
|
self.x_grad = logit_grad(x, self.eps)
|
|
self.inputs = {'X': convert_float_to_uint16(x)}
|
|
self.outputs = {'Out': convert_float_to_uint16(out)}
|
|
self.attrs = {'eps': self.eps}
|
|
|
|
def set_attrs(self):
|
|
self.dtype = np.uint16
|
|
self.shape = [120]
|
|
self.eps = 1e-8
|
|
|
|
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, check_pir=True, check_symbol_infer=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'],
|
|
user_defined_grads=[self.x_grad],
|
|
check_pir=True,
|
|
)
|
|
|
|
|
|
class TestLogitShape(TestLogitOp):
|
|
def set_attrs(self):
|
|
self.dtype = np.float64
|
|
self.shape = [2, 60]
|
|
self.eps = 1e-8
|
|
|
|
|
|
class TestLogitEps(TestLogitOp):
|
|
def set_attrs(self):
|
|
self.dtype = np.float32
|
|
self.shape = [120]
|
|
self.eps = 1e-8
|
|
|
|
|
|
class TestLogit_ZeroSize(TestLogitOp):
|
|
def set_attrs(self):
|
|
self.dtype = np.float64
|
|
self.shape = [2, 0]
|
|
self.eps = 1e-8
|
|
|
|
|
|
class TestLogitAPI(unittest.TestCase):
|
|
def init_data(self):
|
|
self.x_shape = [120]
|
|
self.x_dtype = "float32"
|
|
|
|
def setUp(self):
|
|
self.init_data()
|
|
self.x = np.random.uniform(-1.0, 1.0, self.x_shape).astype(self.x_dtype)
|
|
self.place = get_device_place()
|
|
|
|
def check_api(self, eps=1e-8):
|
|
ref_out = logit(self.x, eps)
|
|
# test static api
|
|
with paddle.static.program_guard(paddle.static.Program()):
|
|
x = paddle.static.data(
|
|
name='x', shape=self.x_shape, dtype=self.x_dtype
|
|
)
|
|
y = paddle.logit(x, eps)
|
|
exe = paddle.static.Executor(self.place)
|
|
out = exe.run(feed={'x': self.x}, fetch_list=[y])
|
|
np.testing.assert_allclose(out[0], ref_out, rtol=1e-05)
|
|
# test dygrapg api
|
|
paddle.disable_static()
|
|
x = paddle.to_tensor(self.x, dtype=self.x_dtype)
|
|
y = paddle.logit(x, eps)
|
|
np.testing.assert_allclose(y.numpy(), ref_out, rtol=1e-05)
|
|
paddle.enable_static()
|
|
|
|
def check_api_grad(self, eps=1e-8):
|
|
ref_grad = logit_grad(self.x, eps)
|
|
numpy_tensor = np.ones(self.x_shape).astype(self.x_dtype)
|
|
# test dygrapg api
|
|
paddle.disable_static()
|
|
paddle_outgrad = paddle.to_tensor(numpy_tensor / numpy_tensor.size)
|
|
x = paddle.to_tensor(self.x, dtype=self.x_dtype)
|
|
x.stop_gradient = False
|
|
y = paddle.logit(x, eps)
|
|
x_grad = paddle.grad([y], [x], [paddle_outgrad])
|
|
np.testing.assert_allclose(x_grad[0].numpy(), ref_grad, rtol=1e-05)
|
|
paddle.enable_static()
|
|
|
|
def test_check_api(self):
|
|
paddle.enable_static()
|
|
for eps in [1e-6, 0.0]:
|
|
self.check_api(eps)
|
|
self.check_api_grad(eps)
|
|
|
|
def test_errors(self):
|
|
paddle.enable_static()
|
|
with paddle.static.program_guard(paddle.static.Program()):
|
|
x = paddle.static.data(name='X1', shape=[100], dtype='bool')
|
|
self.assertRaises(TypeError, paddle.logit, x)
|
|
|
|
x = paddle.static.data(name='X2', shape=[100], dtype='float32')
|
|
self.assertRaises(TypeError, paddle.logit, x, dtype='int32')
|
|
|
|
|
|
class TestLogitAPI_NAN_Val(unittest.TestCase):
|
|
def setUp(self):
|
|
self.init_input_output()
|
|
self.place = [paddle.CPUPlace()]
|
|
if core.is_compiled_with_cuda() or is_custom_device():
|
|
self.place.append(get_device_place())
|
|
|
|
def init_input_output(self):
|
|
self.x = [-0.1, 1.1, 2]
|
|
self.expect_out = [np.nan, np.nan, np.nan]
|
|
self.expect_x_grad = [np.nan, np.nan, np.nan]
|
|
|
|
def test_nan_val(self):
|
|
def _test_nan_val_with_place(place):
|
|
with paddle.base.dygraph.guard():
|
|
x = paddle.to_tensor(self.x, stop_gradient=False, place=place)
|
|
y = paddle.logit(x)
|
|
loss = y.sum()
|
|
loss.backward()
|
|
np.testing.assert_allclose(
|
|
y.numpy(), self.expect_out, rtol=1e-05
|
|
)
|
|
np.testing.assert_allclose(
|
|
x.grad.numpy(), self.expect_x_grad, rtol=1e-05
|
|
)
|
|
|
|
for place in self.place:
|
|
_test_nan_val_with_place(place)
|
|
|
|
|
|
class TestLogitAPICase1(unittest.TestCase):
|
|
def init_data(self):
|
|
self.x_shape = [120]
|
|
self.x_dtype = "float64"
|
|
|
|
|
|
class TestLogitAPICase2(unittest.TestCase):
|
|
def init_data(self):
|
|
self.x_shape = [120]
|
|
self.x_dtype = "float16"
|
|
|
|
|
|
if __name__ == "__main__":
|
|
unittest.main()
|