Discussion :

[membreComplexe12]

Bonjour tous,
j'ai un problème que je trouve très bizarre, j'espère que vous pourrez m'expliquer qu'es ce qui ne va pas dans mon code car là je ne comprends plus du tout

Ci dessous vous allez voir dans ma fonction que j'ai un array qui s'appelle "Loadings", il est initialisé au début de la simulation et ensuite je n'y touche plus (je m'en sers juste comme entrée pour faire des petits calculs)

A mon grand étonnement lorsque je vais la vérification de son contenu (où je met le commentaire #verif1 en réalité il y a une fonction qui vérifie le contenu de Loading) je me rend compte qu 'il est modifié dans cette fonction

Et le truc que comprends encore moins c'est que si j'efface la ligne AllResultsHere[f][tp1][:]=np.copy(CurrentResultTp1) alors la matrice "Loadings" reste intacte comme je le souhaite

J'ai entendu dire qu'avec numpy on ne peux pas faire A=B (je ne sais pas pourquoi ) et qu'il faut faire A=np.copy(B) et c'est bien ce que j'ai fais donc je suis vraiment étonné d'avoir des soucis....

merci pour l'aide que vous pourrez me donner

bonne soirée

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
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
    def allsimulations(self,Param,systemO,inputO,plotO,finalPlot):
        Loadings=np.asarray(inputO.allExpLoadMatrix)
        Results=np.asarray(inputO.allExpRsltMatrix)
        LoadingsOpt=np.asarray(inputO.allOptloadMatrix)
        ResultsOpt=np.asarray(inputO.allOptResultMatrix)
        AllResultsHere=np.asarray(inputO.allExpLoadMatrix)#np.copy(Loadings)
        Times=np.asarray(inputO.allTimeVector)
        ResultsIni=np.asarray(inputO.allInitialResults)
        ### 
        ### 
        #verif1
        ### 
        ### 
        tempCol=[]
        for i in range(0,len(inputO.iniAuxVariables),1):
            tempCol.append(inputO.iniAuxVariables[i])
        ListOfAux=[]
        for i in range(0,len(Loadings),1):
            tempLine=[]
            for j in range(0,len(Loadings[i]),1):
                tempLine.append(tempCol)
            ListOfAux.append(tempLine)
        AuxVar=np.asarray(ListOfAux)
        ### 
        ### 
        #verif2
        ### 
        ### 
        for f in range(0,len(Loadings),1):
            for tp1 in range(0,len(Times[f]),1):
                if tp1==0:
                    CurrentResultTp1=np.copy(ResultsIni[f][0][:])
                    CurrentAuxVarTp1=np.copy(AuxVar[f][0][:])
                else:
                    timeT=Times[f][tp1-1]
                    timeTp1=Times[f][tp1]
                    dTime=timeTp1-Times[f][tp1-1]
                    LoadT=np.copy(Loadings[f][tp1-1][:])
                    LoadTp1=np.copy(Loadings[f][tp1][:])
                    DLoad=np.copy(LoadTp1-Loadings[f][tp1-1][:])
                    DLoadDt=DLoad/dTime
                    ResultT=np.copy(CurrentResultTp1)
                    ResultTp1=np.copy(CurrentResultTp1)
                    AuxVarT=np.copy(CurrentAuxVarTp1)
                    AuxVarTp1=np.copy(CurrentAuxVarTp1)
                    sizeR=len(CurrentResultTp1)
                    sizeA=len(CurrentAuxVarTp1)
                    CurrentResultTp1=runC.simulation(self,Param,sizeR,sizeA,timeT,timeTp1,dTime,LoadT,LoadTp1,DLoad,DLoadDt,ResultT,ResultTp1,AuxVarT,AuxVarTp1)
                ### 
                ### 
                #verif3
                ### 
                ### 
                AllResultsHere[f][tp1][:]=np.copy(CurrentResultTp1)
                ### 
                ### 
                #verif4
                ### 
                ### 
        ##############
        return AllResultsHere

[wiztricks]

Salut,

Pourquoi écrire:

Code :

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

        Loadings=np.asarray(inputO.allExpLoadMatrix)

alors que vous dites:

J'ai entendu dire qu'avec numpy on ne peux pas faire A=B (je ne sais pas pourquoi ) et qu'il faut faire A=np.copy(B) et c'est bien ce que j'ai fais donc je suis vraiment étonné d'avoir des soucis....

La documentation de numpy.asarray donne des exemples qui devraient vous alerter.

- W

[membreComplexe12]

merci beaucoup de ta réponse mais désolé je n'ai pas compris.

J'ai une liste qui se trouve dans une classe input et je veux copier le contenu dans une nouvelle variable "Loading" en faisant une conversion de "list" (biblio de base) vers "array" (numpy)

--> pour ceci np.asarray(a) est donc pas nécessaire ?
--> il faut donc faire np.copy(np.asarray(a)) ?

mais je ne comprends pas trop pourquoi on est obligé de faire ceci ?

c'est surement lié à cette phrase " If a is a subclass of ndarray, a base class ndarray is returned." mais je ne l'a comprends pas trop... pourrais tu l'éclaircir stp pour un débutant comme moi ?

au fait c'est quoi un "ndarray" ce n'est pas comme un "array" ?

[wiztricks]

c'est surement lié à cette phrase "f a is a subclass of ndarray, a base class ndarray is returned." mais je ne l'a comprends pas trop...
pourrais tu l'éclaircir stp pour un débutant comme moi ?

Code :

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

1
2
class A: pass
class B(A): pass

B est une sous classe de A.
La phrase dit que si on passe un B en argument, il retournera un A.
C'est de la gymnastique POO de base.

Mais le plus important est dans la phrase précédente: No copy is performed if the input is already an ndarray.
De plus, regardez les exemples, ils sont assez "parlants".

- W

[marco056]

Pour expliquer de façon simplifiée (d'après ce que j'ai compris), A= B font que A et B sont à la même adresse dans l'ordinateur. Modifier l'un affecte l'autre.
Pour avoir 2 adressages différents, il faut faire une copie.

[membreComplexe12]

merci pour votre aide

j'ai donc remplacé toutes les commandes :

Code :

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

A=np.asarray(B)

par des :

Code :

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

A=np.copy(np.asarray(B))

et le problème persiste....

[wiztricks]

Salut,

C'est que cela ne corrige qu'une partie du problème ou que la cause du problème est autre.
Sans fournir un exemple de code qui permette de reproduire votre problème, pas facile de savoir.

- W

[membreComplexe12]

merci mais à mon avis le problème vient bien de ce bout de code car :

- la variable Loading n'est défini et utilisée que dans cette fonction nulle part ailleur donc à priori il n'y a que les lignes que l'on voit ci dessus qui peuvent avoir un impact sur Loading ? comment une fonction n'ayant rien avoir avec celle là pourrait modifier le contenu d'une variable locale de cette fonction ?

- si je met cette ligne en commentaire "AllResultsHere[f][tp1][:]=np.copy(CurrentResultTp1)" mon problème est résolu

[membreComplexe12]

merci beaucoup de prendre le temps de m'aider.

Sans fournir un exemple de code qui permette de reproduire votre problème, pas facile de savoir.

Suite à cette remarque j'ai réduis mon programme afin d'avoir le cas le plus simple possible qui puisse te permettre de reproduire ce que j'observe.

En pièce jointe de ce message j'ai mis mes trois fichiers pythons, dans le prochain message je mettrais les fichiers d'entrées.

Si tu lances ce bout de code tu vas d'apercevoir que ma matrice "Loading" est plein de 0 alors que je ne l'a rempli pas de 0....
par contre si je commente la ligne dont je t'ai parlé tout à l'heure (dans le fichier runClass.py) alors ma matrice "Loading" garde bien un contenu correct...

si tu peux jeter un coup d'oeil à mon code et me dire d'où est le problème ça serait super, j'ai essayé de nombres choses et là je ne comprends plus rien

merci beaucoup