scipy积分器顺序参数,手动控制时间步长

时间:2019-03-17 01:48:54

标签: python scipy python-3.5 qutip

我正在使用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,unknown_istate_msg))) UserWarning:zvode:检测到非法输入。 (请参阅打印的消息。) 追溯(最近一次通话):   文件“ /home/justin/Methodology_Rotation1/Quantum_Neuron_Sim/Redfield_Simulation/redfield_test.py”在第107行中     Cor1_z =关联_2op_1t(HN.full_liquid_hamiltonian,HN.coupled_state,tau,[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行,相关_2op_1t     求解器=求解器,参数=参数,选项=选项)[0]   _correlation_2t中的文件“ /usr/local/lib/python3.5/dist-packages/qutip/correlation.py”,第1054行     args = args,options = options)   _correlation_me_2t中的文件“ /usr/local/lib/python3.5/dist-packages/qutip/correlation.py”,行1105     [b_op],args = _args,options = options   在mesolve中的文件“ /usr/local/lib/python3.5/dist-packages/qutip/mesolve.py”,第321行     进度条)   _mesolve_list_str_td中的文件“ /usr/local/lib/python3.5/dist-packages/qutip/mesolve.py”,第731行     return _generic_ode_solve(r,rho0,tlist,e_ops,opt,progress_bar)   _generic_ode_solve中的文件“ /usr/local/lib/python3.5/dist-packages/qutip/mesolve.py”,行1013     引发Exception(“ ODE集成错误:尝试增加” 异常:ODE集成错误:尝试通过增加Options类中的nsteps参数来增加允许的子步骤数。

我知道这是scipy集成商试图计算比可接受的时间步长更多的时间步长的问题,所以我想知道以下几点:

  1. 在scipy积分器中,order参数做什么?我问是因为我想知道将其更改为默认值之外的其他值是否可以帮助我解决问题。

  2. 是否可以在scipy积分器中手动替代自动时间步长计算机制?我想做的是设置积分器,以便与其根据要尝试集成的函数的属性来增加时间步长,还可以使用我自己确定的恒定时间步长。 / p>

0 个答案:

没有答案