RuntimeWarning: overflow encountered in exp

[zcc2xj]

C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\uvc_model.py:96: RuntimeWarning: overflow encountered in exp
e_k = np.exp(-gamma_k[i] * (ep_eq - ep_eq_init))
C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\uvc_model.py:99: RuntimeWarning: invalid value encountered in scalar multiply
kin_modulus += c_k[i] * e_k - flow_dir * gamma_k[i] * e_k * alpha_components[i]
Traceback (most recent call last):
File "D:\Python\Mat2Excel\UVC_model.py", line 225, in
UVC_Ti17_R033R06()
File "D:\Python\Mat2Excel\UVC_model.py", line 144, in UVC_Ti17_R033R06
x_sol = rpl.uvc_param_opt(x_0, data_files, x_log, fxn_log, find_initial_point=True, filter_data=False,
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\uvc_parameter_identification.py", line 64, in uvc_param_opt
[opt_results, dumper] = scipy_factory(x.reshape(-1), filtered_data, x_log_file, function_log_file,
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\scipy_constr_opt_factory.py", line 55, in scipy_factory
scipy_sol = opt.minimize(fun, np.array(x_0), method='trust-constr', jac=jac, hess=hess, bounds=bounds,
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_minimize.py", line 753, in minimize
res = _minimize_trustregion_constr(fun, x0, args, jac, hess, hessp,
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_trustregion_constr\minimize_trustregion_constr.py", line 536, in _minimize_trustregion_constr
_, result = tr_interior_point(
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_trustregion_constr\tr_interior_point.py", line 336, in tr_interior_point
z, state = equality_constrained_sqp(
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_trustregion_constr\equality_constrained_sqp.py", line 161, in equality_constrained_sqp
f_next, b_next = fun_and_constr(x_next)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_trustregion_constr\tr_interior_point.py", line 95, in function_and_constraints
f = self.fun(x)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_differentiable_functions.py", line 326, in fun
self._update_fun()
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_differentiable_functions.py", line 295, in _update_fun
fx = self._wrapped_fun(self.x)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\scipy\optimize_differentiable_functions.py", line 21, in wrapped
fx = fun(np.copy(x), *args)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\mat_model_error_nda.py", line 56, in value
error_ensemble = error_ensemble + self.error_fun(x1, cleanedTest)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\uvc_model.py", line 172, in error_single_test_uvc
model_output = uvc_return_mapping(x_sol, test_clean)
File "C:\ProgramData\anaconda3\envs\Mat2Excel\lib\site-packages\RESSPyLab\uvc_model.py", line 136, in uvc_return_mapping
raise RuntimeError('Return mapping did not converge in ' + str(maximum_iterations) + ' iterations.')
RuntimeError: Return mapping did not converge in 1000 iterations.

CODE:
data_files = ['./RESSPyLab_Data/reduced_data_e2_388_r0_33v2.csv']

x_0 = np.array([200000., 355., 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1, 0.1])

x_log = './RESSPyLab_output/reduced_data_e2_388_r0_33v2_x_log_upd.txt'
fxn_log = './RESSPyLab_output/reduced_data_e2_388_r0_33v2_fxn_log_upd.txt'

its = [30, 30, 40]

x_sol = rpl.uvc_param_opt(x_0, data_files, x_log, fxn_log, find_initial_point=True, filter_data=False,
                          step_iterations=its)

data = rpl.load_data_set(data_files)
rpl.uvc_data_plotter(x_sol[0], data, output_dir='./RESSPyLab_output/', file_name='uvc_reduced_data_e2_388_r0_33v2_plots',
                     plot_label='Fitted-UVC')

########

reduced_data_e2_388_r0_33v2.csv data file link to download

链接: https://pan.baidu.com/s/1E6K-C0sKjhlcoY-E-5_AVA?pwd=fe4e