src/tool/Utilitaire/utprmi.F

   1       subroutine utprmi ( v1, v2, v3, res )
   2 c ______________________________________________________________________
   3 c
   4 c                             H O M A R D
   5 c
   6 c Outil de Maillage Adaptatif par Raffinement et Deraffinement d'EDF R&D
   7 c
   8 c Version originale enregistree le 18 juin 1996 sous le numero 96036
   9 c aupres des huissiers de justice Simart et Lavoir a Clamart
  10 c Version 11.2 enregistree le 13 fevrier 2015 sous le numero 2015/014
  11 c aupres des huissiers de justice
  12 c Lavoir, Silinski & Cherqui-Abrahmi a Clamart
  13 c
  14 c    HOMARD est une marque deposee d'Electricite de France
  15 c
  16 c Copyright EDF 1996
  17 c Copyright EDF 1998
  18 c Copyright EDF 2002
  19 c Copyright EDF 2020
  20 c ______________________________________________________________________
  21 c
  22 c   UTilitaire : PRoduit MIxte
  23 c   --           --      --
  24 c     effectue res = V1.(V2^V3)
  25 c ______________________________________________________________________
  26 c .        .     .        .                                            .
  27 c .  nom   . e/s . taille .           description                      .
  28 c .____________________________________________________________________.
  29 c .  V1    . e   .   3    . premier vecteur                            .
  30 c .  V2    . e   .   3    . deuxieme vecteur                           .
  31 c .  V3    . e   .   3    . troisieme vecteur                          .
  32 c . res    .  s  .   1    . Le resultat                                .
  33 c .____________________________________________________________________.
  34 c
  35 c====
  36 c 0. declarations et dimensionnement
  37 c====
  38 c
  39 c 0.1. ==> generalites
  40 c
  41       implicit none
  42       save
  43 c
  44 c 0.2. ==> communs
  45 c
  46 c
  47 c 0.3. ==> arguments
  48 c
  49       double precision v1(3), v2(3), v3(3), res
  50 c
  51 c 0.4. ==> variables locales
  52 c
  53       double precision vn(3)
  54 c
  55       vn(1) = v2(2)*v3(3) - v2(3)*v3(2)
  56       vn(2) = v2(3)*v3(1) - v2(1)*v3(3)
  57       vn(3) = v2(1)*v3(2) - v2(2)*v3(1)
  58 c
  59       res = v1(1)*vn(1) + v1(2)*vn(2) + v1(3)*vn(3)
  60 c
  61       end