RuntimeError: невозможно проанализировать аргументы при использовании scipy.optimize.bisect - PullRequest
Новая
       29

RuntimeError: невозможно проанализировать аргументы при использовании scipy.optimize.bisect

0 голосов
/ 05 июля 2018

У меня проблемы с кодом, использующим корневой искатель scipy.optimize.bisect 1 . При использовании этой процедуры я получаю следующую ошибку: 

Traceback (most recent call last):
File "project1.py", line 100, in <module>
zero1=bisect(rho_function1,x[0],crit,maxiter=1e6,full_output=True)
File "/home/irya/anaconda3/envs/hubble/lib/python3.6/site-packages/scipy        /optimize/zeros.py", line 287, in bisect
r = _zeros._bisect(f,a,b,xtol,rtol,maxiter,args,full_output,disp)
RuntimeError: Unable to parse arguments

Код, который я использую, следующий: 

import numpy as np
from scipy.optimize import bisect
gamma_default=-0.8
gamma_used=-2.0
def rho_function(x,P_0=5.0e3,B_0=3.0e-6,gamma=gamma_default):
    """function used to solve rho. Here, B_0[G] (default 3.0 uG), _0[K/cm^3] and 
     gamma are constants. The variable x is mean to be rho/rho_0"""
     kb=1.3806e-16 # boltzmann  constant in erg/k
     f= B_0**2/(8.0*np.pi*kb)*(x**(4./3)-1) + P_0*(x**gamma-1)
     return f,P_0,B_0,gamma
def rho_function1(x):
     P_0=5.0e3
     B_0=3.0e-6
     gamma=gamma_default
     """function used to solve rho. Here, B_0[G] (default 3.0 uG),  P_0[K/cm^3] and 
     gamma are constants. The variable x is mean to be rho/rho_0"""
     kb=1.3806e-16 # boltzmann  constant in erg/k
     f= B_0**2/(8.0*np.pi*kb)*(x**(4./3)-1) + P_0*(x**gamma-1)
     return f
def rho_prime_function(x,P_0=5.0e3,B_0=3.0e-6,gamma=gamma_default):
     """Derivative of the rho_function"""
     kb=1.3806e-16 # boltzmann  constant in erg/k
     f=B_0**2/(6.0*np.pi*kb)*x**(1./3) + P_0*gamma*x**(gamma-1)
     return f
def rho_2nd_prime_function(x,P_0=5.0e3,B_0=3.0e-6,gamma=gamma_default):
     kb=1.3806e-16 # boltzmann  constant in erg/k
     f=B_0**2/(18.0*np.pi*kb)*x**(-2./3) + P_0*gamma*(gamma-1)*x**   (gamma-2)
     return f
def magnetic_term(x,B_0=3.0e-6):
     """"Magnetic term of the rho function"""
     kb=1.3806e-16 # boltzmann  constant in erg/k
     f=B_0**2/(8.0*np.pi*kb)*(x**(4./3)-1) 
     return f
def TI_term(x,P_0=5.0e3,gamma=gamma_default):
     f=P_0*(x**gamma-1)
     return f

x=np.arange(0.8,2,0.01)

f,P_0_out,B_0_out,gamma_out=rho_function(x,gamma=gamma_used)
f_prime=rho_prime_function(x,gamma=gamma_used)
b_term=magnetic_term(x)
ti_term=TI_term(x,gamma=gamma_used)

kb=1.3806e-16 
crit=(-B_0_out**2/(6.*np.pi*kb*P_0_out*gamma_out))**(3./(3.*gamma_out-4))
print("crit =",crit)
print("Using interval: a=",x[0],", b=",crit)
#using bisect
zero1=bisect(rho_function1,x[0],crit,maxiter=1e6,full_output=True)

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

Приветствия

1 Ответ

0 голосов
/ 06 июля 2018

Параметр maxiter должен быть целым числом. 1e6 это поплавок. Используйте 

maxiter=int(1e6)

или 

maxiter=1000000
...