параметр порядка scipy интегратора, ручное управление временными шагами - PullRequest
Новая
       20

параметр порядка scipy интегратора, ручное управление временными шагами

0 голосов
/ 17 марта 2019

Я использую Qutip для вычисления квантового спектра.Моя ОС - Linux Mint 18, моя версия scipy - 1.2.1, а версия Qutip - 4.4.0.Когда я запускаю вычисления, я получаю следующий вывод: 

Warning (from warnings module):

Файл "/usr/local/lib/python3.5/dist-packages/scipy/integrate/_ode.py", строка 1009 self.messages.get (istate, surprise_istate_msg))) UserWarning: zvode: Обнаружен неправильный ввод..[J_c, J_c.dag ()], I_z, I_z, solver = 'me', args = HN.args, options = opts) Файл "/usr/local/lib/python3.5/dist-packages/qutip/correlation.py ", строка 139, в correlation_2op_1t solver = solver, args = args, options = options) [0] Файл" /usr/local/lib/python3.5/dist-packages/qutip/correlation.py ", строка 1054, в _correlation_2t args = args, options = options) Файл "/usr/local/lib/python3.5/dist-packages/qutip/correlation.py", строка 1105, в _correlation_me_2t [b_op], args = _args, options =Файл настроек "/usr/local/lib/python3.5/dist-packages/qutip/mesolve.py", строка 321, в mesolve progress_bar) Файл "/usr/local/lib/python3.5/dist-packages/qutip/mesolve.py ", строка 731, в _mesolve_list_str_td возвращает _generic_ode_solve (r, rho0, tlist, e_ops, opt, progress_bar) файл" / usr / local / lib / python3.5 / dist-packages / qutip / mesolve.py ", строка 1013, в _generic_ode_solve повысить Исключение (" Ошибка интеграции ODE: попытаться увеличить "Исключение: Ошибка интеграции ODE: попытаться увеличить допустимое числоподшагов путем увеличения параметра nsteps в классе Options.

Я понимаю, что это проблема интегратора scipy, пытающегося вычислить большее количество временных шагов, чем допустимо, поэтому мне интересно следующее:

  1. Что делает параметр порядка в интеграторе scipy?Я спрашиваю, потому что мне интересно, может ли изменение его значения на значение, отличное от значения по умолчанию, помочь мне решить мою проблему.

  2. Есть ли способ переопределить механизм автоматического расчета шага по времени вручную вскучный интегратор?Я хотел бы настроить интегратор таким образом, чтобы вместо увеличения количества временных шагов, основанных на свойствах функции, которую он пытается интегрировать, он использовал постоянное количество временных шагов, которые я определяю сам.

...