[PYTHON/FORTRAN(f2py)] IMSL:DQDAGS <-> scipty.integrate.dblquad()

[squallmrs]

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]

[kango]

bonjour,

quel est le type de func dans:

Code :

scipy.integrate.dblquad(func , A, B, 0e00,1e-08, RESULT, ERREST)

il semblerait que ce soit un float et non une fonction.