项目文件夹

文件
2026-07-13 12:40:42 +08:00

782 行
25 KiB
Python

# Copyright (c) 2024 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
from copy import deepcopy
import numpy as np
from op_test import (
OpTest,
get_device,
get_device_place,
get_devices,
get_places,
is_custom_device,
)
import paddle
from paddle import base
from paddle.framework import core
RTOL = 1e-06
ATOL = 1e-06
def nadam_step(inputs, attributes, dtype='float32'):
param = inputs['param']
grad = inputs['grad']
lr = inputs['learning_rate']
# accumulators
momentum_decay_pow = inputs['momentum_decay_pow']
beta2_pow = inputs['beta2_pow']
mu_product = inputs['mu_product']
moment1 = inputs['moment1']
moment2 = inputs['moment2']
# attrs
epsilon = attributes['epsilon']
beta1 = attributes['beta1']
beta2 = attributes['beta2']
momentum_decay = attributes['momentum_decay']
momentum_decay_pow *= 0.96
beta2_pow *= beta2
mu_t = beta1 * (1.0 - 0.5 * (momentum_decay_pow**momentum_decay))
mu_t_1 = beta1 * (
1.0
- 0.5 * (momentum_decay_pow**momentum_decay) * (0.96**momentum_decay)
)
mu_product *= mu_t
mu_product_t_1 = mu_product * mu_t_1
moment1 = beta1 * moment1 + (1.0 - beta1) * grad
moment2 = beta2 * moment2 + (1.0 - beta2) * grad * grad
moment1_hat = mu_t_1 * moment1 / (1.0 - mu_product_t_1) + (
1.0 - mu_t
) * grad / (1.0 - mu_product)
moment2_hat = moment2 / (1.0 - beta2_pow)
param = param - lr * moment1_hat / (np.sqrt(moment2_hat) + epsilon)
# get accumulators
return (
param.astype(dtype),
momentum_decay_pow.astype(dtype),
beta2_pow.astype(dtype),
mu_product.astype(dtype),
moment1.astype(dtype),
moment2.astype(dtype),
)
def nadam_wrapper(
param,
grad,
lr,
momentum_decay_pow,
beta2_pow,
mu_product,
moment1,
moment2,
master_param=None,
beta1=0.9,
beta2=0.999,
epsilon=1e-8,
momentum_decay=0.004,
multi_precision=False,
):
_, _, _, _, _, _, _ = paddle._C_ops.nadam_(
param,
grad,
lr,
momentum_decay_pow,
beta2_pow,
mu_product,
moment1,
moment2,
master_param,
beta1,
beta2,
epsilon,
momentum_decay,
multi_precision,
)
class TestNAdamOp(OpTest):
def _init_param(self):
self.beta1 = 0.78
self.beta2 = 0.915
self.epsilon = 1e-8
self.momentum_decay = 0.004
def setUp(self):
'''Test NAdam Op with supplied attributes'''
np.random.seed(2024)
self.op_type = "nadam"
self.python_api = nadam_wrapper
self.python_out_sig = ['out']
param = np.random.uniform(-1, 1, (102, 105)).astype("float32")
grad = np.random.uniform(-1, 1, (102, 105)).astype("float32")
learning_rate = np.array(0.003).astype("float32")
self._init_param()
# accumulators
momentum_decay_pow = (np.ones((102, 105)) * (0.96**3)).astype("float32")
# use beta1 to fake mu_product
mu_product = (np.ones((102, 105)) * (self.beta1**3)).astype("float32")
beta2_pow = (np.ones((102, 105)) * (self.beta2**3)).astype("float32")
moment1 = np.random.uniform(-1, 1, (102, 105)).astype("float32")
# The second moment is positive
moment2 = np.random.random((102, 105)).astype("float32")
self.inputs = {
"param": param,
"grad": grad,
"momentum_decay_pow": momentum_decay_pow,
"beta2_pow": beta2_pow,
"mu_product": mu_product,
"moment1": moment1,
"moment2": moment2,
"learning_rate": learning_rate,
}
self.attrs = {
"epsilon": self.epsilon,
"beta1": self.beta1,
"beta2": self.beta2,
"momentum_decay": self.momentum_decay,
}
(
param_out,
momentum_decay_pow_out,
beta2_pow_out,
mu_product_out,
moment1_out,
moment2_out,
) = nadam_step(deepcopy(self.inputs), deepcopy(self.attrs))
self.outputs = {
"param_out": param_out,
"momentum_decay_pow_out": momentum_decay_pow_out,
"beta2_pow_out": beta2_pow_out,
"mu_product_out": mu_product_out,
"moment1_out": moment1_out,
"moment2_out": moment2_out,
}
def test_check_output(self):
self.check_output(check_pir=True, rtol=RTOL, atol=ATOL)
class TestNAdamOpWithDefault(TestNAdamOp):
def _init_param(self):
self.beta1 = 0.9
self.beta2 = 0.999
self.epsilon = 1.0e-8
self.momentum_decay = 0.004
@unittest.skipIf(
not (core.is_compiled_with_cuda() or is_custom_device()),
"core is not compiled with CUDA",
)
class TestNAdamOpGPU(TestNAdamOp):
def test_check_output(self):
self.check_output_with_place(
get_device_place(), check_pir=True, rtol=RTOL, atol=ATOL
)
class TestNAdamOpMultipleSteps(TestNAdamOp):
num_steps = 10
def test_check_output(self):
for _ in range(self.num_steps):
(
param_out,
momentum_decay_pow_out,
beta2_pow_out,
mu_product_out,
moment1_out,
moment2_out,
) = nadam_step(deepcopy(self.inputs), deepcopy(self.attrs))
self.outputs = {
"param_out": param_out,
"momentum_decay_pow_out": momentum_decay_pow_out,
"beta2_pow_out": beta2_pow_out,
"mu_product_out": mu_product_out,
"moment1_out": moment1_out,
"moment2_out": moment2_out,
}
# Verify output for this step
self.check_output()
# Output of this step becomes input for next step
self.inputs['param'] = param_out
self.inputs['momentum_decay_pow'] = momentum_decay_pow_out
self.inputs['beta2_pow'] = beta2_pow_out
self.inputs['mu_product'] = mu_product_out
self.inputs['moment1'] = moment1_out
self.inputs['moment2'] = moment2_out
# Randomize gradient for next step
self.inputs['grad'] = np.random.uniform(-1, 1, (102, 105)).astype(
"float32"
)
class TestNAdamAPI(unittest.TestCase):
def test_nadam_dygraph(self):
paddle.disable_static()
value = np.arange(26).reshape(2, 13).astype("float32")
a = paddle.to_tensor(value)
linear = paddle.nn.Linear(13, 5)
nadam = paddle.optimizer.NAdam(
learning_rate=0.01,
parameters=linear.parameters(),
weight_decay=0.01,
)
for _ in range(2):
out = linear(a)
out.backward()
nadam.step()
nadam.clear_gradients()
def test_nadam_apply_gradients(self):
paddle.disable_static()
value = np.arange(26).reshape(2, 13).astype("float32")
a = paddle.to_tensor(value)
linear = paddle.nn.Linear(13, 5)
nadam = paddle.optimizer.NAdam(
learning_rate=0.01,
parameters=linear.parameters(),
weight_decay=0.01,
)
for _ in range(2):
out = linear(a)
loss = paddle.mean(out)
param_grads = nadam.backward(loss)
nadam.apply_gradients(param_grads)
nadam.clear_gradients()
def test_nadam_static(self):
paddle.enable_static()
place = base.CPUPlace()
shape = [2, 3, 8, 8]
exe = base.Executor(place)
train_prog = base.Program()
startup = base.Program()
with (
base.program_guard(train_prog, startup),
base.unique_name.guard(),
):
data = paddle.static.data(name="data", shape=shape)
hidden = paddle.static.nn.fc(x=data, size=10)
loss = paddle.mean(hidden)
beta1 = 0.85
beta2 = 0.95
opt = paddle.optimizer.NAdam(
learning_rate=1e-5,
beta1=beta1,
beta2=beta2,
weight_decay=0.01,
epsilon=1e-8,
)
opt.minimize(loss)
exe.run(startup)
data_np = np.random.random(shape).astype('float32')
rets = exe.run(train_prog, feed={"data": data_np}, fetch_list=[loss])
assert rets[0] is not None
paddle.disable_static()
def test_pir_nadam(self):
with paddle.pir_utils.IrGuard():
place = base.CPUPlace()
shape = [2, 3, 8, 8]
exe = base.Executor(place)
train_prog = paddle.static.Program()
startup = paddle.static.Program()
with (
paddle.static.program_guard(train_prog, startup),
base.unique_name.guard(),
):
data = paddle.static.data(name="data", shape=shape)
hidden = paddle.static.nn.fc(x=data, size=10)
loss = paddle.mean(hidden)
beta1 = 0.85
beta2 = 0.95
opt = paddle.optimizer.NAdam(
learning_rate=1e-5,
beta1=beta1,
beta2=beta2,
weight_decay=0.01,
epsilon=1e-8,
)
opt.minimize(loss)
exe.run(startup)
data_np = np.random.random(shape).astype('float32')
rets = exe.run(
train_prog, feed={"data": data_np}, fetch_list=[loss]
)
assert rets[0] is not None
def test_nadam_invalid_input(self):
paddle.disable_static()
linear = paddle.nn.Linear(10, 10)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
learning_rate=-1, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, beta1=-1, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, beta2=-1, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, beta1=2.0, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, beta2=2.0, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, epsilon=-1, parameters=linear.parameters()
)
with self.assertRaises(ValueError):
_ = paddle.optimizer.NAdam(
0.1, momentum_decay=-1, parameters=linear.parameters()
)
class TestNAdamAPIWeightDecay(unittest.TestCase):
def test_weight_decay_int(self):
paddle.disable_static()
value = np.arange(26).reshape(2, 13).astype("float32")
a = paddle.to_tensor(value)
linear = paddle.nn.Linear(13, 5)
nadam = paddle.optimizer.NAdam(
learning_rate=0.01,
parameters=linear.parameters(),
weight_decay=1,
)
for _ in range(2):
out = linear(a)
out.backward()
nadam.step()
nadam.clear_gradients()
class TestNAdamAPIGroup(TestNAdamAPI):
def test_nadam_dygraph(self):
paddle.disable_static()
value = np.arange(26).reshape(2, 13).astype("float32")
a = paddle.to_tensor(value)
linear_1 = paddle.nn.Linear(13, 5)
linear_2 = paddle.nn.Linear(5, 3)
nadam = paddle.optimizer.NAdam(
learning_rate=0.01,
parameters=[
{'params': linear_1.parameters()},
{'params': linear_2.parameters(), 'weight_decay': 0.001},
],
weight_decay=0.01,
)
for _ in range(2):
out = linear_1(a)
out = linear_2(out)
out.backward()
nadam.step()
nadam.clear_gradients()
class TestNAdamMultiPrecision(unittest.TestCase):
def _test_nadam_dygraph_place_amp(self, place, use_amp=False):
paddle.disable_static()
paddle.seed(10)
paddle.set_device(place)
input = paddle.randn((5, 5))
model = paddle.nn.Linear(5, 5)
optimizer = paddle.optimizer.NAdam(
parameters=[
{
'params': model.parameters(),
'weight_decay': 0.001,
'beta1': 0.1,
'beta2': 0.99,
}
],
)
optimizer._multi_precision = use_amp
for _ in range(2):
if place == get_device() and use_amp:
model = paddle.amp.decorate(models=model, level='O2')
scaler = paddle.amp.GradScaler(init_loss_scaling=1024)
if place == get_device() and use_amp:
with paddle.amp.auto_cast(level='O2'):
output = model(input)
loss = paddle.mean(output)
scaled = scaler.scale(loss)
scaled.backward()
scaler.step(optimizer)
optimizer.clear_grad()
else:
output = model(input)
loss = paddle.mean(output)
loss.backward()
optimizer.step()
optimizer.clear_grad()
def test_main(self):
for place in get_devices():
use_amp_list = [True, False]
for use_amp in use_amp_list:
self._test_nadam_dygraph_place_amp(place, use_amp)
class TestNdamaxMultiPrecision2_0(unittest.TestCase):
def dygraph_nadam_mp(self, mp, use_amp):
paddle.disable_static()
paddle.seed(100)
paddle.set_device(get_device())
input = paddle.randn((2, 2))
model = paddle.nn.Linear(2, 2)
optimizer = paddle.optimizer.NAdam(0.1, parameters=model.parameters())
optimizer._multi_precision = mp
if use_amp:
model = paddle.amp.decorate(models=model, level='O2')
scaler = paddle.amp.GradScaler(init_loss_scaling=1024)
for idx in range(5):
if use_amp:
with paddle.amp.auto_cast(level='O2'):
output = model(input)
loss = paddle.mean(output)
scaled = scaler.scale(loss)
scaled.backward()
scaler.minimize(optimizer, scaled)
optimizer.clear_grad()
else:
output = model(input)
loss = paddle.mean(output)
loss.backward()
optimizer.step()
optimizer.clear_grad()
return output, model.parameters()
def static_nadam_mp(self, mp, use_amp):
paddle.enable_static()
paddle.seed(2024)
with paddle.pir_utils.OldIrGuard():
exe = paddle.static.Executor(get_device_place())
train_program = paddle.static.Program()
startup_program = paddle.static.Program()
optimizer = paddle.optimizer.NAdam(0.1)
optimizer._multi_precision = mp
if use_amp:
optimizer = paddle.static.amp.decorate(
optimizer,
init_loss_scaling=128.0,
use_dynamic_loss_scaling=True,
use_pure_fp16=True,
use_fp16_guard=False,
)
with paddle.static.program_guard(train_program, startup_program):
if use_amp:
data = paddle.static.data(
shape=[2, 2], name='X', dtype='float16'
)
else:
data = paddle.static.data(
shape=[2, 2], name='X', dtype='float32'
)
hidden = paddle.static.nn.fc(x=data, size=10)
loss = paddle.mean(hidden)
optimizer.minimize(loss)
exe.run(startup_program)
np.random.seed(2024)
if use_amp:
optimizer.amp_init(
place=get_device_place(),
scope=paddle.static.global_scope(),
)
x = np.random.random(size=(2, 2)).astype('float16')
else:
x = np.random.random(size=(2, 2)).astype('float32')
out = []
for idx in range(5):
(loss_data,) = exe.run(
train_program, feed={"X": x}, fetch_list=[loss.name]
)
out.append(loss_data)
return out
def pir_nadam_mp(self, mp, use_amp):
paddle.enable_static()
with paddle.pir_utils.IrGuard():
paddle.seed(2024)
exe = paddle.static.Executor(get_device_place())
train_program = paddle.static.Program()
startup_program = paddle.static.Program()
with paddle.static.program_guard(train_program, startup_program):
model = paddle.nn.Linear(2, 10)
optimizer = paddle.optimizer.NAdam(
0.1, parameters=model.parameters()
)
if use_amp:
data = paddle.static.data(
shape=[2, 2], name='X', dtype='float16'
)
model, optimizer = paddle.amp.decorate(
models=model,
optimizers=optimizer,
level='O2',
master_weight=mp,
)
scaler = paddle.amp.GradScaler(init_loss_scaling=128.0)
with paddle.amp.auto_cast(
level='O2', dtype="float16", use_promote=True
):
output = model(data)
loss = paddle.mean(output)
scaled = scaler.scale(loss)
scaler.minimize(optimizer, scaled)
else:
data = paddle.static.data(
shape=[2, 2], name='X', dtype='float32'
)
output = model(data)
loss = paddle.mean(output)
optimizer.minimize(loss)
exe.run(startup_program)
np.random.seed(2024)
if use_amp:
x = np.random.random(size=(2, 2)).astype('float16')
else:
x = np.random.random(size=(2, 2)).astype('float32')
out = []
for idx in range(5):
(loss_data,) = exe.run(
train_program, feed={"X": x}, fetch_list=[loss]
)
out.append(loss_data)
return out
def static_nadam_amp_o2_without_scaler(self):
paddle.enable_static()
paddle.seed(2024)
with paddle.pir_utils.IrGuard():
train_program = paddle.static.Program()
startup_program = paddle.static.Program()
with paddle.static.program_guard(train_program, startup_program):
exe = paddle.static.Executor(get_device_place())
linear = paddle.nn.Linear(2, 10)
optimizer = paddle.optimizer.NAdam(
0.1, parameters=linear.parameters()
)
linear, optimizer = paddle.amp.decorate(
optimizers=optimizer,
models=linear,
level='O2',
)
data = paddle.static.data(
shape=[2, 2], name='X', dtype='float32'
)
with paddle.amp.auto_cast(
level='O2', dtype='float16', use_promote=True
):
out = linear(data)
loss = paddle.mean(out)
optimizer.minimize(loss)
exe.run(startup_program)
np.random.seed(2024)
x = np.random.random(size=(2, 2)).astype('float32')
out = []
for idx in range(5):
(loss_data,) = exe.run(
train_program, feed={"X": x}, fetch_list=[loss]
)
out.append(loss_data)
return out
def test_main(self):
if not (paddle.is_compiled_with_cuda() or is_custom_device()):
return
"Test dygraph mode"
output1_dy, params1_dy = self.dygraph_nadam_mp(use_amp=True, mp=True)
output2_dy, params2_dy = self.dygraph_nadam_mp(use_amp=False, mp=False)
np.testing.assert_allclose(
output1_dy.astype('float32').numpy(),
output2_dy.astype('float32').numpy(),
rtol=1e-05,
atol=0.1,
)
for idx in range(len(params1_dy)):
np.testing.assert_allclose(
params1_dy[idx].astype('float32').numpy(),
params2_dy[idx].astype('float32').numpy(),
rtol=1e-05,
atol=0.1,
)
"Test static mode"
output1_st = self.static_nadam_mp(use_amp=True, mp=True)
output2_st = self.static_nadam_mp(use_amp=False, mp=False)
output3_st = self.static_nadam_amp_o2_without_scaler()
for idx in range(len(output1_st)):
np.testing.assert_allclose(
output1_st[idx].astype('float32'),
output2_st[idx].astype('float32'),
rtol=1e-05,
atol=0.1,
)
"Test pir mode"
output1_pir = self.pir_nadam_mp(use_amp=True, mp=True)
output2_pir = self.pir_nadam_mp(use_amp=False, mp=False)
for idx in range(len(output1_st)):
np.testing.assert_allclose(
output1_pir[idx].astype('float32'),
output2_pir[idx].astype('float32'),
rtol=1e-05,
atol=0.1,
)
for idx in range(len(output1_st)):
np.testing.assert_allclose(
output1_st[idx].astype('float32'),
output3_st[idx].astype('float32'),
rtol=1e-05,
atol=0.1,
)
class TestNAdamGroupWithLR(TestNAdamAPI):
def test_nadam(self):
paddle.disable_static()
value = np.arange(26).reshape(2, 13).astype("float32")
a = paddle.to_tensor(value)
linear_1 = paddle.nn.Linear(13, 5)
linear_2 = paddle.nn.Linear(5, 3)
nadam = paddle.optimizer.NAdam(
learning_rate=paddle.optimizer.lr.PiecewiseDecay(
boundaries=[3, 6], values=[0.1, 0.2, 0.3]
),
parameters=[
{
'params': linear_1.parameters(),
'learning_rate': 0.1,
},
{
'params': linear_2.parameters(),
'weight_decay': 0.001,
},
],
weight_decay=0.01,
)
for _ in range(2):
out = linear_1(a)
out = linear_2(out)
out.backward()
nadam.step()
nadam.clear_gradients()
def main_test_func(place, dtype):
paddle.enable_static()
main = base.Program()
startup = base.Program()
with (
base.program_guard(main, startup),
base.scope_guard(base.Scope()),
):
x = paddle.static.data(name='x', shape=[None, 13], dtype=dtype)
y = paddle.static.data(name='y', shape=[None, 1], dtype=dtype)
y_predict = paddle.static.nn.fc(x, size=1)
cost = paddle.nn.functional.square_error_cost(input=y_predict, label=y)
avg_cost = paddle.mean(cost)
nadam_optimizer = paddle.optimizer.NAdam(0.01)
nadam_optimizer.minimize(avg_cost)
fetch_list = [avg_cost]
train_reader = list(
zip(
np.random.rand(101, 13),
np.random.randint(12, size=(101, 1)),
)
)
feeder = base.DataFeeder(place=place, feed_list=[x, y])
exe = base.Executor(place)
exe.run(base.default_startup_program())
for data in train_reader:
exe.run(main, feed=feeder.feed([data]), fetch_list=fetch_list)
paddle.disable_static()
class NAdamFp32Test(unittest.TestCase):
def setUp(self):
self.dtype = 'float32'
def test_main(self):
for p in get_places():
main_test_func(p, self.dtype)
class NAdamFp64Test(NAdamFp32Test):
def setUp(self):
self.dtype = 'float64'
if __name__ == "__main__":
unittest.main()