# coding:utf8

"""
Ce programme est un exemple d'utilisation du module ising.

Copiez-collez ce programme pour l'adapter !

Vous pouvez renommer vos fichiers TP1_X-Y.py pour la réponse à la question X.Y.
(par exemple TP1_II-2.py pour la question II.2. du TP1.).

Pour lancer ce script, tapez `python3 TP1_example.py` dans un terminal.
"""

# Import des modules
import numpy as np
import matplotlib.pyplot as plt
import ising_NOM_PRENOM as ising



# Définit les variables d'entrée
size = 64               # Taille du réseau
Ncorr = 1000            # Nombre d'itération de l'algorithme de Metropolis entre chaque mesure
N_iter = 3000            # Nombre de fois que l'on appelle la fonction move_config
T = 1.0                 # Température
Nplots = 30             # Nombre de configurations à afficher


''' QUESTION II. 2''' 
# Initialise la configuration sur un réseau 64 x 64
config = ising.initialize_config(size)

# Boucle de 0 à N_iter inclus (N_iter + 1 exclu)
for it in range(0, N_iter + 1):
    # Met à jour la configuration courante
    config = ising.move_config(config, T, Ncorr)
    # Affiche la configuration Nplots fois par simulation
    if it % (N_iter//Nplots) == 0: # l'opérateur // permet de faire une division entière, l'opérateur % est le modulo
        # Construit le titre
        title = "T = " + str(T) + ", Ncorr = " + str(Ncorr) + ", Pas de temps = " + str(it)
        # Affiche la configuration sur la figure 0, puis fait une pause de 0.1 secondes
        ising.plot_config(config, title=title, fignumber=0, pause=.1)


# Attends une entrée utilisateur pour terminer le programme.
# En Python2, remplacez input par raw_input.
# En console interactive, commentez la ligne.
input('Appuyez sur "Entrée" pour terminer le programme.')


''' QUESTION II. 3'''