Publicité
+ Répondre à la discussion
Affichage des résultats 1 à 2 sur 2
  1. #1
    Invité régulier
    Inscrit en
    novembre 2008
    Messages
    38
    Détails du profil
    Informations forums :
    Inscription : novembre 2008
    Messages : 38
    Points : 8
    Points
    8

    Par défaut [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]

  2. #2
    Membre émérite
    Homme Profil pro
    Inscrit en
    décembre 2007
    Messages
    754
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 35
    Localisation : France

    Informations professionnelles :
    Secteur : Aéronautique - Marine - Espace - Armement

    Informations forums :
    Inscription : décembre 2007
    Messages : 754
    Points : 887
    Points
    887

    Par défaut

    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.

Liens sociaux

Règles de messages

  • Vous ne pouvez pas créer de nouvelles discussions
  • Vous ne pouvez pas envoyer des réponses
  • Vous ne pouvez pas envoyer des pièces jointes
  • Vous ne pouvez pas modifier vos messages
  •