Problème dans une partie d'une simulation (Python)

[invitec0cbc192]

Bonsoir,

Je me demandais si quelqu'un pouvait m'aider avec une simulation en Python (une petite partie du code où je me pose des questions) concernant le modèle d'Ising.
Je voulais notamment obtenir la magnétisation M, et cette dernière est définie comme :


d'où j'ai écrit dans le code, cette partie pour avoir M (et pas de problème, ça a fonctionné) :

sigmaxop = []
sites = []
for i in range(L):
    # \sum_i sigma^z(i)
    sigmaxop.append((sigmaz).tolist())
    sites.append([i])

où L est la longueur du réseau (càd le nb de spins).

Mais maintenant j'aimerais obtenir les fluctuations moyennes càd :


Et pour ça, évidemment, il "suffit" que je réécrive le code ci-dessus j'imagine pour obtenir les fluctuations mais je me trompe quand j'essaie quelque chose car je n'obtiens pas ces fluctuations Ça fait un certain temps que j'essaie mais en vain, c'est pourquoi je me tourne vers le forum.

Est-ce que quelqu'un aurait une idée de comment je peux modifier cette partie de code pour obtenir les fluctuations ?

Merci d'avance.
-----

[invite896757ff]

A mon avis tu n'a pas à réécrire le code qui te sort le tableau
Il faut l'imbriquer dans une boucle ou tu refera ton calcul (ta boucle i) avec une variation en entrée et mettre en tableau.
Ensuite il faudra comparer les tableaux pour en extraire la variation des résultats.
A moins qu'une instruction savante te fasse ça en une ligne.

[invite6c250b59]

d'où j'ai écrit dans le code, cette partie pour avoir M (et pas de problème, ça a fonctionné) :

Perso je n'arrive pas à comprendre l'utilité de ce bout de code. Est-ce que tu pourrais ajouter le code qui suit, celui qui précède, et dans un autre code/code ce que ce que tu as essayé qui n'a pas marché et comment tu as deviné que ça n'a pas marché?