Discussion :

[sebsainte]

Bonjour,

Je suis amené à calculer des séries mathématiques pour des applications
scientifiques. (ce sont des sommes de termes qui dépendent d'un paramètre)

l'algo serait le suivant :

------------------------------------------------

somme = 0

alpha =150
beta=12
...


pour n de 0 à 100
somme = somme + f(n, alpha, beta, ...)
fin pour

Resultat = somme

------------------------------------------------


avec f une fonction qui dépend d'un entier n (entier sur lequel on
effectue la sommation) mais aussi d'autres paramètes qui eux sont fixés
avant la somme (alpha, beta, ...) .


je cherche à créer une fonction qui permet de réaliser directement la somme

ex :

Resultat = MaSuperSomme(f (n,alpha, beta, .... ), ParametreDeSomme = n, n=0, n=100)

Je suis en galère :

- pour passer la fonction f en paramètre de MaSuperSomme
- pour déterminer sur lequel des arguments de f on effecue la somme.


Merci pour vos idées !

Seb

[see++]

Salut,

Voici une solution.
c'est certainement pas la meilleure, mais ca te permettra peut être de te débloquer...

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
def toto(n,alpha,beta) :
    return n+alpha+beta
 
def MaSuperSomme(f,args,parametre,debut,fin):
    somme = 0
    for i in range(debut,fin) :
        args[parametre]=i
        fct = f.__name__+"("
        for cle,valeur in args.items(): fct+=cle+"="+str(valeur)+","
        fct=fct[:-1] # pour virer la virgule
        fct += ")"
        exec "somme+="+fct
    return somme
 
 
print MaSuperSomme(toto,{'n':0,'alpha':3,'beta':4},parametre='n',debut=0,fin=100)

Si tu as des questions sur le code, je suis là

See++

[sebsainte]

Bonjour see++

La solution que tu proposes est vraiment intéressante. Elle ouvre
un certain nombre d'applicatons sur les manipulations génériques
de fonctions.

Merci pour l'effort d'avoir construit un exemple fonctionnel

Le script est très clair

Merci et à bientôt sur le forum !

++

Seb

[see++]

De rien, ca fait toujours plaisir de pouvoir aider quelqu'un

A+ sur le forum

See++

[Matthieu Brucher]

On devrait pouvoir s'arranger sans passer par le exec.
Si tu passes ta fonction en argument, tu peux directement l'appeler, sans pb. Le réel problème reste les arguments.
Pour cela, tu peux faire qqch du style :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
def MaSuperFonction(f, debut, fin, **parameters):
  somme = 0 #
  for i in range(debut, fin):
    somme += f(parameters)
  return somme

l'intérêt du **parameters est que c'est un dictionnaire contenant toutes les paramètres passés à la fonction qui ne font pas partie des attributs définis explicitements. Ensuite, il te reste à faire un :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

print MaSuperFonction(f = maFonction, debut = 0, fin = 10, alpha = 30, beta = 45, gamma = 56, ...)

[Guigui_]

On peut aussi se servir des fonctions sum et map:
(y'a peut-être moyen d'être plus générique)

exemple:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
alpha =150
beta=12
 
def somme(f, debut, fin, a, b):
    return sum(map(f, range(debut, fin+1), (fin-debut+1) * [a], (fin-debut+1) * [b]))
 
def f(n, a, b): return n*a+b
 
print somme(f, 0, 10, a = alpha, b = beta)

Edit: finalement avec les listes compréhensions, ce sera plus générique:

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
alpha =150
beta=12
 
def somme(f, debut, fin, **param):
    return sum([f(i, **param) for i in range(debut, fin+1)])
 
def f(n, a, b): return n*a+b
 
print somme(f, 0, 10, a = alpha, b = beta)

[SRT__]

Bonjour, (j'ai changé de paramètres de connexion... anciennement sebsainte)
Je me suis finalement rapproché de fonctions présentes dans la librairie scipy,
qui procède globalement comme les deux dernières réponses (notamment pour les fonctions d'intégration)

le code est très proche de ces deux dernières

Ca marche très bien !

Merci pour vos réponses