FR 2954876 A1 20110701 FR 0906391 A 20091230 PROCEDE ET SYSTEME POUR DETERMINER ET OPTIMISER LA DISPONIBILITE D'UN RESEAU DE COMMUNICATION L'objet de la présente invention concerne un procédé et un système pour déterminer la probabilité de pouvoir communiquer (disponibilité) entre des noeuds, dans un réseau de télécommunications dont l'architecture initiale est connue. Les résultats de la disponibilité sont, par exemple, mémorisés dans un fichier relié à un dispositif permettant d'adapter l'architecture du réseau, c'est-à-dire le nombre et la qualité des noeuds et les liens associés, au besoin d'un utilisateur ou à une interface Homme machine en fonction des exigences du fonctionnement du réseau requis préalablement par un utilisateur. Définitions utilisées dans la description de l'invention La sûreté de fonctionnement d'un réseau de communication s'appuie généralement sur les notions suivantes : • Disponibilité C'est la probabilité qu'un composant ou un système soit en état de marche à un instant donné t. • Fiabilité : C'est la probabilité qu'un composant ou un système fonctionne sur l'intervalle de temps [0,t]. • Maintenabilité : C'est l'aptitude d'une entité à être maintenue ou remise en état de fonctionnement (la prise en considération de l'aptitude du système de soutien à maintenir ou à remettre en état l'entité est du domaine du soutien logistique intégré). • Sécurité : C'est l'aptitude d'un produit à respecter, pendant toutes les phases de vie, un niveau acceptable de risques d'accident susceptible de causer une agression du personnel ou une dégradation majeure du produit ou de son environnement.30 Dans le domaine des réseaux de télécommunications, il s'avère indispensable de connaître au mieux la disponibilité ou l'état de fonctionnement de composants participant à la transmission des données ou des informations. Ainsi, il est souhaitable de connaître l'état de fonctionnement ou disponibilité de noeuds et des liens qui vont participer au relayage des informations. Différentes méthodes sont décrites dans l'art antérieur. La plupart de ces méthodes font appel à des outils complexes qui requièrent de disposer de puissance de calcul importante. La figure 1 représente un exemple de réseau constitué de plusieurs noeuds Ni reliés entre eux par des liens de communication li, pouvant être un média physique ou encore tout autre média connu et utilisé dans le domaine des communications sans fil. Les noeuds sont définis par : un identifiant unique, représenté par un nombre entre 0 et 999 pour les noeuds. un indicateur de bon fonctionnement, qui permet de définir si l'élément est en état de bon fonctionnement ou si au contraire celui-ci est en panne. une valeur de disponibilité au sens général, c'est à dire une valeur entre 0 et 1 qui peut représenter une disponibilité intrinsèque, une disponibilité opérationnelle, une probabilité de bon fonctionnement sur un temps donné, etc. un type (abonné ou transit). Certaines méthodes ne traitent que sur les noeuds abonnés (le calcul analytique par exemple). le réseau étant constitué d'une couche logique et d'une couche physique il y a donc un indicateur définissant l'appartenance à telle ou telle couche. Les liaisons entre noeuds sont définies par exemple avec: un identifiant unique dont le nombre minimal est, par exemple, 1000. une disponibilité. un poids qui représente en réalité le coût technique d'une liaison. un indicateur qui définit le type de la liaison (liaison physique ou logique). II est connu de l'art antérieur de déterminer de manière analytique la disponibilité de noeuds et de liaisons par un calcul analytique. Il s'agit de calculer analytiquement la probabilité d'avoir au moins k routes disponibles entre deux points parmi n routes. Les probabilités pour les routes d'être disponibles n'étant pas indépendantes les unes des autres, il convient d'appliquer les formules de probabilités correspondantes. La probabilité pour une route d'être disponible est égale au produit des disponibilités des éléments qui la composent. En revanche, pour 2 routes A et B qui relient deux points Rd et Ra, la probabilité d'avoir au moins une (k=1) route disponible entre Rd et Ra se note : P(AuB) = P(A) + P(B) ù P(AnB) soit la somme des probabilités de fonctionnement des routes A et B, à laquelle on retranche la probabilité que les deux routes soient en bon 15 fonctionnement simultanément. La généralisation est donnée par la formule de Poincaré (au moins 1 parmi n) : m m j-1 m j 1 1-1 A(t)=E p(Pj)EE p(PiPj+EEE p(PkPiPj)+...+(ù1)mp(P1P2...Pm) j=l j=2 1=1 j=3 i=2 k=l 20 où p(Pj) représente la probabilité pour qu'à l'instant t les éléments x de l'ensemble de routes Pj, de probabilité qx(t) soient en marche, donc p(Pj) = fl qx (t) .Formule de Poincaré. xej On remarque que lorsqu'on calcule la probabilité de l'intersection des événements Ai (« la route i est disponible » ou encore « tous les éléments 25 qui composent la route i sont disponibles ») et Aj (« la route j est disponible » ou encore « tous les éléments qui composent la route j sont disponibles »), il faut prendre en compte le fait que Ai et Aj ne sont pas forcément disjoints (« les routes i et j ne sont pas forcément totalement disjointes » ou encore « il peut exister des éléments qui appartiennent en même temps à i et à j »). On calcule donc la probabilité de l'intersection des événements Ai et Aj (« tous les éléments de la route i et tous les éléments de la route j sont disponibles » sans compter deux fois les éléments qui sont dans i et dans j) et non pas le simple produit des probabilités des routes. Pour rendre le calcul possible en temps de calcul, il est connu d'utiliser le filtrage de routes. Ainsi, la sélection des routes est un critère déterminant dans les résultats obtenus par le calcul analytique. Pour chacun de ces critères, il est demandé à l'utilisateur le nombre de routes qu'il souhaite calculer : • Sélection des routes les plus disjointes auxquelles s'ajoutent les plus courtes au sens de l'algorithme de routage réseau (RIP) • Sélection des routes les plus disjointes • Sélection des routes les plus disjointes auxquelles s'ajoutent les moins lourdes au sens de l'algorithme de routage réseau (OSPF) • Sélection de routes basée sur une couche logique • Sélection d'une route forcée. Le calcul analytique utilisé selon l'art antérieur présente notamment comme inconvénient d'être coûteux en temps de calcul. Calcul de disponibilité par simulation de type Monte-Carlo pour un 20 calcul de bout en bout II est aussi connu de l'art antérieur qu'effectuer un calcul de disponibilité entre deux noeuds du réseau par simulation de type Monte-Carlo est une méthode faisant appel à l'ensemble du réseau. Pour ce faire, il faut dans un premier temps générer une série de 25 tirs ou de simulations. Lors de chaque tir ou simulation, on génère des nombres aléatoires compris entre 0 et 1, et une loi de probabilité (loi uniforme) est utilisée. Ainsi pour chaque pas de simulation, on affecte à chaque noeud et chaque liaison un nombre généré aléatoirement. La disponibilité des noeuds Ni et des liaisons li est comparée avec les nombres 30 aléatoirement générés, si ceux-ci sont inférieurs à la disponibilité on considère que les liaisons ou les noeuds fonctionnent correctement. La disponibilité pour un calcul sur une paire du réseau est donc donnée par le rapport entre le nombre de succès et le nombre total de simulation. Pour vérifier la disponibilité d'un chemin entre 2 noeuds, on exécute un algorithme de parcours en profondeur qui va être explicité plus loin dans la description. L'invention repose sur une nouvelle approche de calcul de disponibilité depuis une extrémité du réseau vers plusieurs autres extrémités, qui prend donc en compte un noeud de départ et un ou plusieurs noeuds destinataires du réseau de communications, ainsi que les liens existants entre lesdits noeuds. L'invention concerne un procédé pour déterminer et/ou optimiser la disponibilité d'un réseau de communications caractérisé en ce qu'il exécute au moins les étapes suivantes : Entrée : soit un Graphe g G(N, L, p, q), Nombre de simulations nbSimulations, tableau de paires, avec N le nombre de noeuds, L le nombre de liaisons, p la probabilité des liaisons p (p,, p2, P3, ... ,PN), q le probabilité des noeuds q = (q1, q2, q3, ... qL), et deux fonctions La fonction RandUniform0l () qui permet de générer des nombres aléatoires grâce au générateur « Marsenne Twister ». La fonction DFS_IsRouteExiste() qui parcourt le graphe en profondeur pour vérifier l'existence d'un chemin entre 2 noeuds du graphe. Sortie : Une approximation de la disponibilité D Pour i allant de 1 à nbSimulations //Initialisation 0 Pour j allant de 1 à N On génère un nombre aléatoire avec RandUniform0l () pour chaque Noeuds du réseau et on le compare avec sa probabilité pi Fin pour Pour k allant de 1 à L On génère un nombre aléatoire avec RandUniform01()pour chaque liaisons du réseau et on le compare avec sa probabilité qk, valeur initiale du noeud Fin pour Faire On cherche ensuite l'existence d'un chemin communiquant entre les 2 noeuds que l'on veut étudier avec DFS_IsRouteExisteO. Tant que (il reste des paires à traiter) Si (toutes les paires de noeuds ont pu communiquer) alors X *-- X + 1 Fin si Fin Pour Retourner la valeur de disponibilité D = X/nbSimulations, Mémoriser ladite valeur de disponibilité dans un fichier. D'autres caractéristiques et avantages de la présente invention apparaîtront mieux à la lecture de la description qui suit d'un exemple de réalisation du procédé et du système selon l'invention donnée à titre illustratif et nullement limitatif annexé des figures qui représentent • La figure 1, un exemple de réseaux comprenant plusieurs noeuds connectés, 1T,... • La figure 2, Une interface machine développée pour la mise en oeuvre du procédé, • La figure 3, une interface possible intégrant une fonctionnalité de convergence, • La figure 4, un graphe explicatif dans le cas d'un calcul de disponibilité globale, et • La figure 5, une interface Homme machine pour le calcul de la disponibilité globale, et • La figure 6, un diagramme des étapes exécutées par la mise en oeuvre du logiciel selon l'invention piloté par un utilisateur, Le procédé pour déterminer la disponibilité d'une extrémité ou d'un noeud du réseau vers plusieurs autres noeuds va notamment utiliser un nombre de simulations donné. Pour chaque simulation à chaque noeud Ni et à chaque lien lij, un nombre généré aléatoirement va être affecté. Hypothèse Pour le calcul de disponibilité selon l'invention « UnVersPlusieurs » on ne sélectionne pas de route au préalable, mais le procédé s'appuie sur l'ensemble du réseau. Le procédé selon l'invention va s'intéresser à un calcul de disponibilité des noeuds Ni et/ou des liens lij entre un premier noeud du réseau et plusieurs noeuds du réseau (contexte de point à multipoints). II va ensuite comparer la disponibilité (probabilité) des noeuds et des liaisons calculée en mettant en oeuvre les étapes spécifiques du procédé selon l'invention avec la valeur des nombres générés aléatoirement. Si le nombre généré aléatoirement est inférieur à la disponibilité de la liaison ou du noeud alors on la marque comme OK, donc en état de bon fonctionnement (permettant le passage de l'information), dans le cas contraire on le considère comme défaillant, NOK (ne laissant pas passer l'information). Le procédé vérifie ensuite en fonction des liaisons et noeuds ouverts (en état de marche) s'il existe un chemin menant du point d'origine vers tous les points d'arrivée définis par l'utilisateur, et ce sans exception. Pour vérifier l'existence d'un chemin entre deux noeuds, le procédé utilise, par exemple, l'algorithme de parcours en profondeur d'un graphe connu de l'Homme du métier, DFS (Depth First Search). Cette méthode repose sur le marquage des noeuds voisins d'un noeud donné et ainsi se propager dans le graphe. Dans l'exploration, l'algorithme cherche à aller très vite "profondément" dans le graphe, en s'éloignant du sommet s de départ. La recherche sélectionne à chaque étape un sommet voisin du sommet marqué à l'étape précédente. S'il existe un chemin entre le point de départ et tous les points d'arrivée le procédé considère cela comme un succès. Dans le cas contraire il s'agit d'un échec. La disponibilité est alors donnée par le rapport entre le nombre de succès sur le nombre total de simulations. On obtient ainsi la disponibilité pour un noeud vers plusieurs autres noeuds. Le calcul de disponibilité selon l'invention repose notamment sur la génération de nombres aléatoires. II est possible d'utiliser comme générateur de nombres aléatoires le générateur (Pseudo Random Number Generator) qui se nomme le 'Marsenne Twister' connu de l'Homme du métier. Ce générateur aléatoire a été développé en 1996/1997, avec une mise à jour en 2002, par M. Makoto Matsumoto et Takuji Nishimura (Université d'Hiroshima). Description de l'étape de calcul de disponibilité selon l'invention d'un premier noeud vers plusieurs autres noeuds du réseau Fonctions utiles : o La fonction RandUniform01O qui permet de générer des nombres aléatoires grâce au générateur « Marsenne Twister ». o La fonction DFS_IsRouteExiste() qui parcourt le graphe en profondeur pour vérifier l'existence d'un chemin entre 2 noeuds du graphe. Soit un graphe g = G(N, L, p, q ) N : nombre de noeuds L : nombre de liaisons p : probabilité des liaisons p = (pi, p2: p3, ... ,PN) q : probabilité des noeuds q = (ch, q2, q3, ... , qL)30 Entrée Un Graphe g = G(N, L, p, q), Nombre de simulations nbSimulations, tableau de paires. Sortie : Une approximation de la disponibilité D Pour i allant de 1 à nbSimulations //Initialisation X.ù 0 Pour j allant de 1 à N On génère un nombre aléatoire avec RandUniform01() pour chaque Noeuds du réseau et on le compare avec sa probabilité pi Fin pour Pour k allant de 1 à L On génère un nombre aléatoire avec RandUniform01()pour chaque liaisons du réseau et on le compare avec sa probabilité qk, valeur initiale du noeud Fin pour Faire On cherche ensuite l'existence d'un chemin communiquant entre les 2 noeuds que l'on veut étudier avec DFS_IsRouteExiste(). Tant que (il reste des paires à traiter) Si (toutes les paires de noeuds ont pu communiquer) alors X 4-- X+1 Fin si Fin Pour Retourner la valeur de disponibilité D = X/nbSimulations La valeur trouvée de la disponibilité D pour un noeud vers plusieurs autres noeuds du réseau est, par exemple, mémorisée dans un fichier qui est en liaison avec une interface Homme Machine et qui permet l'affichage des données trouvées sur une console ou bien dans un fichier en relation avec un processeur qui va permettre d'interagir et modifier un ou plusieurs éléments (noeud, liaison) à partir du réseau initial sur lequel les étapes ont été exécutées dans le cas où la disponibilité qui a été trouvée ne convient pas au vu des exigences données par l'utilisateur du réseau. Une interface homme machine développée pour la mise en oeuvre du procédé selon l'invention est représentée à la figure 2. Elle comporte notamment des boutons permettant le lancement de 10 la méthode d'une extrémité vers plusieurs autres en calcul rapide et avec étude de la convergence telle que détaillée ci-après. Exemple de l'application du procédé en relation avec la figure 1 Avec les différentes valeurs des noeuds et des liaisons. Noeuds Laisons C T OK 0,9 L 2 0K 0,8 C T OK 0,9 L 5 0K 0,8 C T OK 0,9 L 6 0K 0,8 C T OK 0,9 L 16 0K 0,8 C T OK 0,9 L 3 0K 0,8 C T OK 0,9 L 4 0K 0,8 C T OK 0,9 L 7 0K 0,8 C T OK 0,9 L 8 0K 0,8 C T OK 0,9 L 9 0K 0,8 C T OK 0,9 L 10 0K 0,8 C T OK 0,9 L 11 0K 0,8 C T OK 0,9 L 12 0K 0,8 C T OK 0,9 L 13 0K 0,8 C T OK 0,9 L 14 0K 0,8 C T OK 0,9 L 15 0K 0,8 T OK 0,9 L 16 0K 0,8 C T OK 0,9 L 10 0K 0,8 C T OK 0,9 L 13 0K 0,8 C T OK 0,9 L 34 0K 0,8 C T OK 0,9 L 33 0K 0,8 C T 0 K0 , 9 L 7 0K 0,8 C T OK 0,9 L 9 0K 0,8 C T 0 K0 , 9 L 10 0K 0,8 C T OK 0,9 L 7 0K 0,8 C T OK 0,9 L 14 0K 0,8 C T OK 0,9 L 11 0K 0,8 C T OK 0,9 L 12 0K 0,8 C T OK 0,9 L 9 0K 0,8 C T 0 K0 , 9 L 3 0K 0,8 C T 0 K0 , 9 L 13 0K 0,8 C T OK 0,9 L 15 0K 0,8 C T OK 0,9 L 11 0K 0,8 C T 0 K0 , 9 L 12 0K 0,8 C T OK 0,9 L 14 0K 0,8 C A OK 1 L 5 0K 0,8 C A OK 1 L 16 0K 0,8 L 6 0K 0,8 L 15 0K 0,8 L 16 0K 0,8 L 1 0K 1 L 17 0K 1 Selon une variante de mise en oeuvre de l'invention, le procédé peut comporter une étape supplémentaire qui permet la convergence dans les résultats de calcul de disponibilité décrit précédemment. La figure 3 schématise une interface possible intégrant cette fonctionnalité. Au niveau de cette interface on trouve deux fonctionnalités permettant le calcul de un vers plusieurs avec ou sans convergence, ainsi que la possibilité de faire le calcul de disponibilité de bout en bout avec étude de la convergence. La convergence est définie par des critères d'arrêts qui peuvent être les suivants : le nombre simulations et la précision qu'un utilisateur souhaite. Selon un exemple de mise en oeuvre, le procédé comporte une 10 étape d'étude de la convergence des résultats obtenus afin d'affiner les résultats. Pour calculer la disponibilité d'un réseau de communication avec étude de la convergence, il est possible de fixer, par exemple, deux critères d'arrêt de la méthode. Le premier critère peut être choisi en fonction d'un 15 nombre d'itérations afin d'éviter de se trouver dans une configuration de boucle infinie, un deuxième critère peut prendre en compte la précision du résultat. Celui-ci permet de donner un résultat valide à plusieurs chiffres après la virgule (le nombre de chiffres après la virgule dépendant de la précision qu-' . : ° • = i : s le champ prévu à cet effet 20 dans l'interface- figure 5). Etant donné que les champs renseignés par un utilisateur au niveau de l'interface sont des critères d'arrêt dans le code, les valeurs d'initialisation de l'algorithme sont fixées en dur au préalable à sa mise en oeuvre. On commence donc le calcul, par exemple, avec 1000 simulations, 25 la valeur 1000 est donc modifiable qu'à partir.du code source. Pour lancer le calcul l'utilisateur va entrer un nombre de simulations supérieur à 1000 et une précision comprise en ] 0 ; 0.1 ]. Pour s'assurer de la précision, le calcul de disponibilité de "Un vers plusieurs" va être exécuté 2 fois, pour un même nombre de simulations. 30 Lors de la première itération le nombre de simulations et donc de 1000. Une fois les 2 calculs effectués on compare les résultats. Si la soustraction de ces 2 nombres en valeur absolue est inférieure à la précision choisie alors le résultat est valide, on arrête donc le calcul. Si la précision n'est par contre pas atteinte, alors on multiplie le nombre de simulations par 10 et on recommence les calculs pour effectuer à nouveau la 5 comparaison. On s'arrête lorsque l'un des critères d'arrêt est atteint (les 2 peuvent être atteints en même temps). Un exemple chiffré de la mise en oeuvre du procédé selon 10 l'invention avec convergence est donné ci-après. En utilisant la méthode de calcul de CMCS Un Vers Plusieurs() avec étude de la convergence selon le procédé selon l'invention les résultats suivants sont obtenus : Critère d'arrêt : nombre de simulations 1 000 000 000 (109), précision 15 0,00001 ^ En tenant compte des noeuds de départ et d'arrivée. Paire Nb simulations Disponibilité Nb rejeux critère d'arrêt atteint Temps de calcul 1 X 10000000 0.132278200000000 5 Précision Oh 2mn 5s 1 X 100000000 0.132253660000000 6 Précision Oh 20mn 10s 1 X 10000000 0.132282700000000 5 Précision Oh 2mn 5s 1 X 10000000 0.132261600000000 5 Précision Oh 2mn 6s Une autre variante de mise en oeuvre du procédé selon l'invention 20 concerne le cas de calcul de disponibilité entre tous les points ou noeuds du réseau afin de pouvoir exprimer une disponibilité globale du réseau. A partir de cette valeur, il sera alors possible de vérifier si cette disponibilité correspond bien aux exigences requises par un utilisateur. La disponibilité globale est définie de la manière suivante: 25 La disponibilité globale est définie comme étant le fait qu'un noeud de départ s puisse communiquer avec tous les autres noeuds du réseau. s communique avec t1 .... tä _i ( pour un réseau comportant n noeuds). o Un noeud est capable de communiquer avec tous les autres. pour tout j tel que j # i, le couple (n;, n;) est connecté il est possible d'exprimer ceci de manière plus général en disant que la disponibilité globale est la probabilité que tout noeud ni, ni communique avec tous les autres noeuds. o tous les noeuds sont capables de communiquer entre eux. Pour tout i , j, tel que i # j le couple (n;, ni) communique Ces 2 notions sont équivalentes en effet : Pour tout i , j tel que i # j, pour tout j tel que j # i, le couple le couple (ni,nj) communique (n;, ni) communique => pour tout i, j, i # j le couple (ni,ni) communique donc pour i fixé, pour tout j # i on a bien (ni,ni) connecté. La figure 6 représente les différentes étapes pour la mise en oeuvre du procédé selon l'invention. 20 Les étapes mises en oeuvre dans les différentes variantes de réalisation exposées précédemment sont par exemple implémentées en utilisant le langage C++ ou tout autre langage de programmation qui permet une manipulation des différents objets ; noeuds, liaisons, réseau. 25 L'interface Homme machine peut être une interface avec fenêtres de saisie sous Windows, des boutons d'action, des menus déroulant, une utilisation de dispositifs de pointage. Le procédé selon l'invention offre un outil permettant de déterminer 30 différents types de disponibilité dans un réseau de communications, facile d'utilisation, présentant une interface graphique. Il répond notamment aux besoins suivants : • Participer aux efforts de conception des réseaux fixes ou faiblement mobiles, grâce aux grandeurs de sûreté de fonctionnement annoncées, • Justifier le choix d'une architecture réseau face à des exigences de sûreté de fonctionnement requis par les utilisateurs, REVENDICATIONS1 û Procédé pour déterminer et/ou optimiser la disponibilité d'un réseau de communications, c'est-à-dire la probabilité pour des noeuds de communiquer entre eux dans ledit réseau de communication, ledit procédé étant exécuté sur un processeur caractérisé en ce qu'il exécute au moins les étapes suivantes : Définir en Entrée : un Graphe g = G(N, L, p, q), Nombre de simulations nbSimulations, tableau de paires, avec N le nombre de noeuds, L le nombre de liaisons, p la probabilité des liaisons p = (pi, p2, P3, ... ,PN), q le probabilité des noeuds q = (qI, q2, q3, ... , qL), et utiliser une fonction RandUniform0l () qui permet de générer des nombres aléatoires grâce au générateur « Marsenne Twister ». utiliser une fonction DFS_IsRouteExiste() qui parcourt le graphe en profondeur pour vérifier l'existence d'un chemin entre 2 noeuds du graphe. Déterminer en Sortie : Une approximation de la valeur de disponibilité D en exécutant les étapes suivantes : Pour i allant de 1 à nbSimulations //Initialisation X 0 Pour j allant de 1 à N On génère un nombre aléatoire avec RandUniform0l() pour chaque Noeud du réseau et on le compare avec sa probabilité pi Fin pour Pour k allant de 1 à L On génère un nombre aléatoire avec RandUniformOlOpour chaque liaisons du réseau et on le compare avec sa probabilité qk, valeur initiale du noeud Fin pour Faire 30On cherche ensuite l'existence d'un chemin communiquant entre les 2 noeuds que l'on veut étudier avec DFS 1sRouteExiste(). Tant que (il reste des paires à traiter) Si (toutes les paires de noeuds ont pu communiquer) alors X •ù X+1 Fin si Fin Pour Retourner la valeur de disponibilité D = X/nbSimulations, Mémoriser ladite valeur de disponibilité D pour un noeud vers un autre noeud dans un fichier, afin d'interagir et modifier un ou plusieurs éléments (noeud, liaison) de l'architecture du réseau ou valider ladite architecture. 2 ù Procédé selon la revendication 1 caractérisé en ce qu'il comporte une étape de détermination de la disponibilité globale du réseau, définie comme étant le fait qu'un noeud source S peut communiquer avec tous les autres noeuds du réseau, ladite étape comportant une étape de calcul de disponibilité pour que toutes les paires de noeud du réseau soient connectées entre elles, en exécutant les étapes de la revendication 1 en partant d'un noeud source vers tous les autres noeuds. 3 ù Procédé selon l'une des revendications 1 ou 2 caractérisé en ce qu'il comporte une étape supplémentaire d'étude de convergence utilisant deux critères d'arrêt, le premier étant fixé en fonction d'un nombre d'itérations pour ne pas se trouver dans une configuration de boucle infinie, le deuxième critère prenant en compte la précision du résultat, ces critères étant paramétrés dans le code du logiciel selon l'invention, le procédé comportant au moins les étapes suivantes Soit X la valeur du nombre de simulations dans le code et Y l'intervalle de la précision requise,exécuter au moins 2 fois le calcul de disponibilité de Un vers plusieurs pour un même nombre de simulations selon les étapes de l'une des revendications 1 ou 2, une fois les 2 calculs effectués comparer les résultats, et si la soustraction de ces 2 nombres en valeur absolue est inférieure à la précision choisie alors arrêter le calcul. ii la précision n'est par contre pas atteinte, alors multiplier le nombre de simulations par un facteur a et recommencer les calculs pour effectuer à nouveau la comparaison, arrêter l'exécution des étapes lorsque l'un des critères d'arrêt est atteint, et interagir et modifier un ou plusieurs éléments (noeud, liaison) de l'architecture du réseau ou valider ladite architecture.15