Discussion :

[jcmic]

Bonjour,

Je rencontre un petit problème :
A l'intérieur d'une fonction, je voudrais remplir un tableaux de N flottants à l'aide d'un certain algorithme.
N est une valeur entiére passée en paramétre d'entré de la fonction.
Le tableau lui meme est renvoyé par la fonction.

Soit y le tableau en question. J'ai été voir dans différents tutoriels et il m'a semblé que l'on déclarait ce tableau par l'instruction y[].
Voici alors le début du code posant problème (à l'intérieur de la fonction):[/I]

y=[]
h=(borne-x0)/N
y[0]=y0

for i in range (1,N-1) :
y[i]=y[i-1]+h*phi(y[i-1]) etc...


Lors de l'execution, j'obtient le message suivant :
"... y[0]=y0
IndexError: list assignment index out of range"

qui semble m'indiquer que je déborde la dimension du tableau dés l'instruction y[0]=y0.
Quelle doit-etre la bonne syntaxe pour définir un tableau de taille N et le remplir de 0 à N-1 svp ?
Par avance merci

jcmic

[OliSharkPython]

Bonjour,

Le message est explicite.
Quand vous faites y[0] = y0, cela suppose que votre liste y est de dimension 1 (les indices commençant à 0). Or votre liste est de dimension 0.

Code :

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

y = []; print len(y)

pour vous en convaincre !.

Pour résoudre votre problème, soit vous dimensionnez votre liste (tableau) directement à la dimension N dès le début (car la valeur est connue), soit vous ajouter à votre liste (tableau) les valeurs au fur et à que vous les calculez.

Pour la première solution :

Code :

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

1
2
3
4
5
 
h=(borne-x0)/N
y=[ y0 for i in range(N) ]
for i in xrange(1,N) :
 y[i]=y[i-1]+h*phi(y[i-1])

Pour la seconde solution :

Code :

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

1
2
3
4
5
 
h=(borne-x0)/N
y = [ y0 ]
for i in xrange(1,N) :
 y.append( y[-1] + h*phi(y[-1]) )

Il y a sûrement des améliorations à ce que je vous propose (en une ligne ou presque ?) et/ou d'autres solutions.

Cordialement,

[jcmic]

Ok, merci pour votre réponse claire.
Ce qui me manquait, relativement à la premiére solution, c'était effectivement la syntaxe pour dimensionner le tableau avant de le gérer.
Merci encore

Cordialement

jcmic

[mont29]

Deux petites remarques*:

*D’une part, il beaucoup plus rapide de définir tout de suite le tableau à sa taille finale, si celle-ci est connue, plutôt que de rajouter les éléments un à un (on gagne déjà la moitié du temps avec la méthode proposée ci-dessus, avec py 2.6)*!

*D’autre part, il est beaucoup plus rapide d’utiliser y = [y0]*N plutôt qu’une liste comprehension (y = [y0 for i in range(N)]), ici*(plus de dix fois plus, toujours en py 2.6)*!

Notez que py3, bien que globalement environ 25% plus rapide, conserve les mêmes ordres de grandeur…

[wiztricks]

Salut,

Python vient avec une bibliothèque array qui est probablement plus efficace pour gérer les types de base que les listes.

Doug Hellman a fait il y a longtemps un PyMOTW là dessus. Il y montre comment çà s'utilise.

Nous avons aussi les "array" de numpy mais c'est autre chose.

- W

[jcmic]

Merci, effectivement cette solution parait intéressante
Cordialement