1 subroutine utqpyr ( lapyra, qualit, qualij, volume,
2 > coonoe, somare, aretri,
3 > facpyr, cofapy, arepyr )
4 c ______________________________________________________________________
8 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
10 c Version originale enregistree le 18 juin 1996 sous le numero 96036
11 c aupres des huissiers de justice Simart et Lavoir a Clamart
12 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
13 c aupres des huissiers de justice
14 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
16 c HOMARD est une marque deposee d'Electricite de France
22 c ______________________________________________________________________
24 c UTilitaire : Qualite d'une PYRamide
26 c ______________________________________________________________________
28 c . Jacobien normalise
29 c ______________________________________________________________________
31 c . nom . e/s . taille . description .
32 c .____________________________________________________________________.
33 c . lapyra . e . 1 . numero de la pyramide a examiner .
34 c . qualit . s . 1 . qualite non definie .
35 c . qualij . s . 1 . qualite par le jacobien normalise .
36 c . volume . s . 1 . volume .
37 c . coonoe . e . nbnoto . coordonnees des noeuds .
39 c . somare . e .2*nbarto. numeros des extremites d'arete .
40 c . aretri . e .nbtrto*3. numeros des 3 aretes des triangles .
41 c . facpyr . e .nbpycf*5. numeros des 5 faces des pyramides .
42 c . cofapy . e .nbpycf*5. codes des faces des pyramides .
43 c . arepyr . e .nbpyca*8. numeros des 8 aretes des pyramides .
44 c ______________________________________________________________________
47 c 0. declarations et dimensionnement
50 c 0.1. ==> generalites
64 double precision qualit, qualij, volume
65 double precision coonoe(nbnoto,3)
68 integer somare(2,nbarto)
69 integer aretri(nbtrto,3)
70 integer facpyr(nbpycf,5), cofapy(nbpycf,5), arepyr(nbpyca,8)
72 c 0.4. ==> variables locales
76 integer listar(8), listso(5)
80 c 0.5. ==> initialisations
81 c ______________________________________________________________________
99 c x--------------------------------------------------------x
101 c La face f5 est le quadrangle.
102 c La face fi, i<5, est le triangle s'appuyant sur l'arete ai.
105 c 1. les aretes et les sommets
108 call utaspy ( lapyra,
109 > nbtrto, nbpycf, nbpyca,
111 > facpyr, cofapy, arepyr,
118 call utvopy ( coonoe, listso, volume )
127 c 4. qualite par le jacobien normalise
129 c 4.1. ==> Liens sommet/aretes
171 c 4.2. ==> fonction generique
174 daux = sqrt(2.d0)/2.d0
175 call utqjno ( iaux, aresom, daux,
176 > listar, listso, somare, coonoe,
178 cgn write(1,*) '==> qualij : ', qualij