Discussion :

[Merlin_II]

Bonjour,

j'essaye difficilement de migrer de Fortran a Python, mais je sens que ma logique ne suit pas.

Je fais une boucle dans laquelle je fit 3 parametres, et je souhaite rentrer ces 3 parametres dans un tableau, voici un extrait du code:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
...
rvs = []
...
    for k in range(1,4):
        ri,vi = init(x,y)
        si=1.
        guess=[ri,vi,si]
        pars, cov = curve_fit(f=fonc, xdata=x, ydata=y, p0=guess, bounds=(binf, bsup))
        rvs[k-1].append(pars[0],pars[1],pars[2])     <=======================================================
        stdevs = np.sqrt(np.diag(cov))
...

L'erreur est au niveau de la fleche, avec le message:

Code :

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

1
2
3
 
rvs[k-1].append(pars[0],pars[1],pars[2])
IndexError: list index out of range

pourtant le premier indice de mon tableau est bien zero (j'ai cru comprendre que ce message d'erreur correspondait a des indices de tableau qui n'existaient pas).
J'ai essaye aussi de creer 3 tableaux sur le meme model, avec une "affectation" simple (rt=[] et rt[k-1].append(pars[0]) et meme au tout depart avec rt[k-1]=pars[0], mais ca ne fonctionne pas...
Si vous pouvez me mettre sur la voie ce serait super!

[Sve@r]

Bonjour

pourtant le premier indice de mon tableau est bien zero (j'ai cru comprendre que ce message d'erreur correspondait a des indices de tableau qui n'existaient pas).

Effectivement le message "list out of range" indique que tu dépasses la limite de ton tableau (en plus ou en moins car en python tu as droit aux index négatifs qui commencent alors à -1 pour le dernier élément). Et il est probable que ce message vienne de rvs[k-1] car rvs est un tableau vide (donc ne contient aucun élément pouvant être récupéré avec son indice)
Tu as commis 2 erreurs

• tu n'as pas étudié suffisement append() qui s'applique sur la liste et non pas sur un de ces éléments (sauf si cet élément est lui-même une liste) => rvs.append(truc) et non rvs[i].append(truc)
• append() ne prend qu'un paramètre en compte. Si tu veux rajouter 3 trucs, il te faut faire 3 append() à suivre. Ou alors utiliser extend() qui lui accepte un itérable en paramètre pour étendre la liste sur laquelle il est appliqué =>rvs.extend((pars[0], pars[1], pars[2])) ou rvs.extend(pars[0:3]) ou rvs.extend(pars) si "pars" n'a que 3 éléments ou encore plus simplement rvs=list(pars) si rvs est créé from scratch

Accessoirement je ne pige pas pourquoi tu fais varier k de 1 à 3 pour taper dans k-1 alors que tu aurais pu faire varier k de 0 à 2 directement...

meme au tout depart avec rt[k-1]=pars[0], mais ca ne fonctionne pas...

Interdit. Tu ne peux modifier par affectation directe qu'un élément déjà existant (donc ayant déjà un indice). L'agrandissement d'un tableau ne peut se faire qu'en utilisant des méthodes précises du tableau (append, extend, insert)

Et en Python, les tableaux (sous-entendu "tableaux classiques des autres langages") se nomment "listes". Car il y a d'autres outils permettant de représenter plus ou moins les tableaux: les listes (ok les tableaux classiques donc), les tuples (listes non modifiables), les ensembles (tableaux ne contenant aucun élément dupliqué) et les dictionnaires (tableaux associatifs).

[Merlin_II]

Bonjour Sve@r,
c'est effectivement ca, merci beaucoup, cela fonctionne maintenant!
Je n'ai pas utiliser le extend (mais je note son existence!), j'ai donc mis rvs.append(pars) puisqu'effectivement pars n'a que 3 elements.
Le passage Fortran (77!) a Pyhton reste difficile pour moi, et cette aide m'est tres precieuse pour ne pas me decourager et finalement rester a Fortran
Merci beaucoup, je clos la discussion!

[Merlin_II]

Bonjour Wiztricks,
desole, la reponse est apparue entre-temps. Effectivement ca fait plus concis, mais a mon niveau, j'ai davantage besoin de clarte que de concision ^^' Mais ca ira mieux dans quelques temps avec ces aspects!
Merci beaucoup!!

[Sve@r]

Je n'ai pas utiliser le extend (mais je note son existence!), j'ai donc mis rvs.append(pars) puisqu'effectivement pars n'a que 3 elements.

Non, erreur !!!
rvs.append(X) ajoute "X" au tableau "rvs". Si "X" est lui-même un tableau, alors rvs deviendra un tableau 2D (tableau à un élément où cet élément est lui-même un tableau)
Dans ton cas, c'est (enfin il me semble) c'est rvs.extend(pars) qu'il faut écrire, ce qui se traduira alors par une extension 1D prise depuis le contenu de pars. Et la solution de wiztricks (qui m'était sortie de la tête tellement il y a de façons différentes de faire les choses en Python) fonctionnera aussi.

[Merlin_II]

Bonjour Sve@r,
oui, ca fonctionnait avec append mais effectivement, ca me faisait des tableaux de tableaux, donc pas facile de retrouver mes petits ensuite. Avec extend, j'ai exactement ce qu'il me faut, rvs est bien un vecteur 1D
Merci beaucoup!!

[wiztricks]

Salut,

Si rvs est une liste, rvs[k-1] en serait le (k-1)^ième élément.
Or si cette liste est vide, elle n'a aucun élément... et peut probable qu'il soit une liste.

Après si le but est d'ajouter à cette liste les 3 premiers éléments de pars à chaque itération rvs += pars[:3] sera plus concis.

- W