Évacuation d'une salle

Évacuation d'une salle - Python - Programmation

Marsh Posté le 19-05-2015 à 13:41:39    

Bonjour,  
 
J'essaye de programmer, dans le cadre d'un projet en informatique, une fonction qui va donner le temps d'évacuation d'une salle remplie de personnes et d'obstacles (comme des tables par exemple) en nombre d'itérations. Elle prend en arguments un tableau (composé de 0 (les vides), de 1 (les personnes) et de 2 (les obstacles)) modélisant la salle et la position (iporte,jporte) de la porte. Pour chaque 1 je calcule la distance à la porte des 4 cases l'entourant grâce à la formule : distance = sqrt((i-iporte)²+(j-jporte)²). La case vide de distance minimale reçoit un 3 (pour que le 1 ne soit pas déplacé plusieurs fois en une boucle) et l'ancienne un 0. Une boucle rechange ensuite tous les 3 en 1. Si la case contient un obstacle (2) ou est vide (0), elle n'est pas modifiée. A chaque entrée de boucle en remet à 0 la porte. Lorsque la somme des cases du tableau est égale à la somme des obstacles (il n'y a plus aucun 1) la fonction s'arrête et renvoie le nombre d'itérations nécessaires à l'évacuation de la salle. J'ai essayé de coder cette fonction en Python mais elle ne s'arrête pas et je ne vois pas où est le problème. Voilà le code de la fonction, est-ce que vous pourriez m'aider à voir ce qui ne vas pas ? Merci d'avance !  
 

Code :
  1. import numpy as np
  2. import numpy.random as rd
  4. def evacuation (M,iporte,jporte):
  5.     Temps_evacuation=0
  6.     (L,l)=np.shape(M)
  7.     nombre_obstacles=0
  8.     for i in range (L):
  9.         for j in range (l):
  10.             if M[i,j]==2:
  11.                 nombre_obstacles+=1
  12.     while np.sum(M)>2*nombre_obstacles:
  13.         print M       
  14.         Temps_evacuation+=1
  15.         M[iporte,jporte]=0 #On remet la porte à 0
  16.         for i in range (L):
  17.             for j in range (l):
  18.                 if M[i,j]==1:
  19.                     distance=np.ones(4)*1000 #On génère un vecteur rempli d'un grand nombre pour que le nombre reste grand si on ne peut pas calculer la distance
  20.                     if i>=0 and i+1<=L and j>=0 and j+1<=l:
  21.                         if i>0 and M[i-1,j]==0:
  22.                             distance[0]=np.sqrt((i-1-iporte)**2+(j-jporte)**2) #case du haut
  23.                         if i<L-1 and M[i+1,j]==0:
  24.                             distance[1]=np.sqrt((i+1-iporte)**2+(j-jporte)**2) #case du bas
  25.                         if j>0 and M[i,j-1]==0:
  26.                             distance[2]=np.sqrt((i-iporte)**2+(j-1-jporte)**2) # case de gauche
  27.                         if j<l-1 and M[i,j+1]==0:
  28.                             distance[3]=np.sqrt((i-iporte)**2+(j+1-jporte)**2) # case de droite
  29.                             distancemin = min (distance)
  30.                             indice_case_valeur_min=[]
  31.                             for k in range (4): #On cherche l'indice associé à la valeur minimale
  32.                                 if distance[k]==distancemin:
  33.                                     indice_case_valeur_min.append(k)
  34.                             if len (indice_case_valeur_min)==2: #Si deux distances sont les mêmes on tire au hasard
  35.                                 indice=rd.randint (0,2)
  36.                                 if indice==0:
  37.                                     indice_case_valeur_min=indice_case_valeur_min[0]
  38.                                 if indice==1:
  39.                                     indice_case_valeur_min=indice_case_valeur_min[1]
  40.                             else :
  41.                                 indice_case_valeur_min=indice_case_valeur_min[0]
  42.                             if indice_case_valeur_min==0 and i>0 and M[i-1,j]!=2:
  43.                                 M[i-1,j]=3
  44.                                 M[i,j]=0
  45.                             if indice_case_valeur_min==1 and i<L-1 and M[i+1,j]!=2:
  46.                                 M[i+1,j]=3
  47.                                 M[i,j]=0
  48.                             if indice_case_valeur_min==2 and j>0 and M[i,j-1]!=2:
  49.                                 M[i,j-1]=3
  50.                                 M[i,j]=0
  51.                             if indice_case_valeur_min==3 and j<l-1 and M[i,j+1]!=2:
  52.                                 M[i,j+1]=3
  53.                                 M[i,j]=0
  54.                         if M[i,j]==2: #Les 2 sont inchangés
  55.                             M[i,j]=2
  56.                         if M[i,j]==0: #Les 0 sont inchangés
  57.                             M[i,j]=0
  58.         for i in range (L): #On rechange les 3 en 1
  59.             for j in range (l):
  60.                 if M[i,j]==3:
  61.                     M[i,j]=1
  62.     return Temps_evacuation


Message édité par sherlock1995 le 19-05-2015 à 13:45:49
Reply

Marsh Posté le 19-05-2015 à 13:41:39   

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed