[PYTHON/FORTRAN(f2py)] IMSL:DQDAGS <-> scipty.integrate.dblquad()
Bonjour tout le monde,
Je continue mon projet ou je dois remplacer un vieu code fortran, par du python. Et je suis bloqué a cette erreur... que je ne comprends pas
Voici le code fortran :
Code:
1
2
|
CALL DQDAGS(FONC,A,B,0.D00,1.D-08,RESULT,ERREST) |
Voici le code python qui va remplacer le code fortran , avec son erreur :
Code:
1
2
3
4
5
|
RESULT = 0e00
ERREST = 0e00
scipy.integrate.dblquad(func , A, B, 0e00,1e-08, RESULT, ERREST) |
Et voici le resultat :
File "TRANSF.py", line 191, in ?
tr.CC()
File "TRANSF.py", line 130, in CC
TRANSF_SCHWARZ_CHRISTOFFEL.EVALUE(Transf)
File "/ecf/xa99222/HodoGraph/V1.0_Python/transf/TRANSF_SCHWARZ_CHRISTOFFEL.py", line 56, in EVALUE
Transf.zctran.lint_[i-1]= math.fabs(INTEGR(Transf.f, BORN1,BORN2))
File "/ecf/xa99222/HodoGraph/V1.0_Python/transf/TRANSF_SCHWARZ_CHRISTOFFEL.py", line 10, in INTEGR
integrate.dblquad(func , A, B, 0e00,1e-08, RESULT, ERREST)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadpack.py", line 351, in dblquad
return quad(_infunc,a,b,(func,gfun,hfun,args),epsabs=epsabs,epsrel=epsrel)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadpack.py", line 187, in quad
retval = _quad(func,a,b,args,full_output,epsabs,epsrel,limit,points)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadpack.py", line 251, in _quad
return _quadpack._qagse(func,a,b,args,full_output,epsabs,epsrel,limit)
File "/usr/lib64/python2.4/site-packages/scipy/integrate/quadpack.py", line 311, in _infunc
a = gfun(x)
TypeError: 'float' object is not callable
[4] + Fait ge /ecf/xa99222/HODO/V0.0/src/transf/integr.F &
[/code]
