1 subroutine utsoqu ( somare, a1, a2, a3, a4,
2 > sa1a2, sa2a3, sa3a4, sa4a1 )
3 c ______________________________________________________________________
7 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
9 c Version originale enregistree le 18 juin 1996 sous le numero 96036
10 c aupres des huissiers de justice Simart et Lavoir a Clamart
11 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
12 c aupres des huissiers de justice
13 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
15 c HOMARD est une marque deposee d'Electricite de France
21 c ______________________________________________________________________
23 c UTilitaire - SOmmets d'un QUadrangle
25 c Remarque : cela suppose que les aretes a1, a2, a3, a4 soient donnees
26 c dans l'ordre standard
27 c Remarque : programme semblable a utoraq
28 c ______________________________________________________________________
30 c . nom . e/s . taille . description .
31 c .____________________________________________________________________.
32 c . somare . e .2*nbaret. numeros des extremites d'arete .
33 c .a1,..,a4. e . 1 . les numeros d'aretes du quadrangle .
34 c . sa1a2 . s . 1 . sommet commun aux aretes a1 et a2 .
35 c . sa2a3 . s . 1 . sommet commun aux aretes a2 et a3 .
36 c . sa3a4 . s . 1 . sommet commun aux aretes a3 et a4 .
37 c . sa4a1 . s . 1 . sommet commun aux aretes a4 et a1 .
38 c ______________________________________________________________________
41 c 0. declarations et dimensionnement
44 c 0.1. ==> generalites
54 integer a1, a2, a3, a4
55 integer sa1a2, sa2a3, sa3a4, sa4a1
57 c 0.4. ==> variables locales
60 c 1. recherche des sommets
73 if ( somare(1,a1).eq.somare(1,a2) .or.
74 > somare(1,a1).eq.somare(2,a2) ) then
75 c le 1er noeud de l'arete 1 est un sommet de a2 ;
76 c donc le 2nd noeud de l'arete 1 est un sommet de a4
82 c le 1er noeud de l'arete 1 n'est pas un sommet de a2 ;
83 c donc c'est qu'il est un des sommets de a4
84 c donc le 2nd noeud de l'arete 1 est un sommet de a2
90 if ( somare(1,a3).eq.somare(1,a4) .or.
91 > somare(1,a3).eq.somare(2,a4) ) then
92 c le 1er noeud de l'arete 3 est un sommet de a4 ;
93 c donc le 2nd noeud de l'arete 3 est un sommet de a2
98 c le 1er noeud de l'arete 3 n'est pas un sommet de a4 ;
99 c donc c'est qu'il est un des sommets de a2
100 c donc le 2nd noeud de l'arete 3 est un sommet de a4