src/tool/AP_Conversion/pcstr2_1.h

   1 c                   .                         .
   2 c                  . .                       ...
   3 c                 .   .                     . . .
   4 c                .     .                   .  .  .
   5 c               .       .      ===>       .   .   .
   6 c              .         .               .    .    .
   7 c             .           .             .     .     .
   8 c            .             .           .      .      .
   9 c           .................         .................
  10 c
  11 c         Pour un decoupage selon l'arete numero 1 :
  12 c
  13           elseif ( etanp1.eq.1 ) then
  14 #ifdef _DEBUG_HOMARD_
  15       write (ulsort,texte(langue,1)) 'Entree', 'pcstr2_1 etanp1=1'
  16 #endif
  17 c
  18             f1hp = filtri(trhnp1)
  19             f1cp = ntrsca(f1hp)
  20             f2cp = ntrsca(f1hp+1)
  21             prfcap(f1cp) = 1
  22             prfcap(f2cp) = 1
  23 c
  24             do 2221 , nrofon = 1 , nbfonc
  25 c
  26 c           Pour le fils aine :
  27 c
  28               vafott(nrofon,q1,f1cp) = vafoen(nrofon,q1,prtrcn)
  29 c
  30               vafott(nrofon,q2,f1cp) = vafoen(nrofon,q5,prtrcn)
  31 c
  32               vafott(nrofon,q3,f1cp) = vafoen(nrofon,q3,prtrcn)
  33 c
  34               vafott(nrofon,q4,f1cp) =
  35      >       -unshu*(vafoen(nrofon,q2,prtrcn)+vafoen(nrofon,q3,prtrcn))+
  36      >        unsde*(vafoen(nrofon,q4,prtrcn)+vafoen(nrofon,q6,prtrcn))+
  37      >         unsqu*vafoen(nrofon,q5,prtrcn)
  38 c
  39               vafott(nrofon,q5,f1cp) =
  40      >       -unshu*vafoen(nrofon,q2,prtrcn)+
  41      >        trshu*vafoen(nrofon,q3,prtrcn)+
  42      >        trsqu*vafoen(nrofon,q5,prtrcn)
  43 c
  44               vafott(nrofon,q6,f1cp) = vafoen(nrofon,q6,prtrcn)
  45 c
  46 c           Pour le triangle fils NF+1 :
  47 c
  48               vafott(nrofon,q1,f2cp) = vafoen(nrofon,q1,prtrcn)
  49 c
  50               vafott(nrofon,q2,f2cp) = vafoen(nrofon,q2,prtrcn)
  51 c
  52               vafott(nrofon,q3,f2cp) = vafoen(nrofon,q5,prtrcn)
  53 c
  54               vafott(nrofon,q4,f2cp) = vafoen(nrofon,q4,prtrcn)
  55 c
  56               vafott(nrofon,q5,f2cp) =
  57      >        trshu*vafoen(nrofon,q2,prtrcn)-
  58      >        unshu*vafoen(nrofon,q3,prtrcn)+
  59      >        trsqu*vafoen(nrofon,q5,prtrcn)
  60 c
  61               vafott(nrofon,q6,f2cp) = vafott(nrofon,q4,f1cp)
  62 c
  63  2221       continue
  64 c
  65 c         Pour un decoupage selon l'arete numero 2 :
  66 c
  67           elseif ( etanp1.eq.2 ) then
  68 #ifdef _DEBUG_HOMARD_
  69       write (ulsort,texte(langue,1)) 'Entree', 'pcstr2_1 etanp1=2'
  70 #endif
  71 c
  72             f1hp = filtri(trhnp1)
  73             f1cp = ntrsca(f1hp)
  74             f2cp = ntrsca(f1hp+1)
  75             prfcap(f1cp) = 1
  76             prfcap(f2cp) = 1
  77 c
  78             do 2222 , nrofon = 1 , nbfonc
  79 c
  80 c           Pour le fils aine :
  81 c
  82               vafott(nrofon,q1,f1cp) = vafoen(nrofon,q6,prtrcn)
  83 c
  84               vafott(nrofon,q2,f1cp) = vafoen(nrofon,q2,prtrcn)
  85 c
  86               vafott(nrofon,q3,f1cp) = vafoen(nrofon,q3,prtrcn)
  87 c
  88               vafott(nrofon,q4,f1cp) =
  89      >       -unshu*(vafoen(nrofon,q1,prtrcn)+vafoen(nrofon,q3,prtrcn))+
  90      >        unsde*(vafoen(nrofon,q4,prtrcn)+vafoen(nrofon,q5,prtrcn))+
  91      >         unsqu*vafoen(nrofon,q6,prtrcn)
  92 c
  93               vafott(nrofon,q5,f1cp) = vafoen(nrofon,q5,prtrcn)
  94 c
  95               vafott(nrofon,q6,f1cp) =
  96      >       -unshu*vafoen(nrofon,q1,prtrcn)+
  97      >        trshu*vafoen(nrofon,q3,prtrcn)+
  98      >        trsqu*vafoen(nrofon,q6,prtrcn)
  99 c
 100 c           Pour le triangle fils NF+1 :
 101 c
 102               vafott(nrofon,q1,f2cp) = vafoen(nrofon,q1,prtrcn)
 103 c
 104               vafott(nrofon,q2,f2cp) = vafoen(nrofon,q2,prtrcn)
 105 c
 106               vafott(nrofon,q3,f2cp) = vafoen(nrofon,q6,prtrcn)
 107 c
 108               vafott(nrofon,q4,f2cp) = vafoen(nrofon,q4,prtrcn)
 109 c
 110               vafott(nrofon,q5,f2cp) = vafott(nrofon,q4,f1cp)
 111 c
 112               vafott(nrofon,q6,f2cp) =
 113      >        trshu*vafoen(nrofon,q1,prtrcn)-
 114      >        unshu*vafoen(nrofon,q3,prtrcn)+
 115      >        trsqu*vafoen(nrofon,q6,prtrcn)
 116 c
 117  2222       continue
 118 c
 119 c         Pour un decoupage selon l'arete numero 3 :
 120 c
 121           elseif ( etanp1.eq.3 ) then
 122 #ifdef _DEBUG_HOMARD_
 123       write (ulsort,texte(langue,1)) 'Entree', 'pcstr2_1 etanp1=3'
 124 #endif
 125 c
 126             f1hp = filtri(trhnp1)
 127             f1cp = ntrsca(f1hp)
 128             f2cp = ntrsca(f1hp+1)
 129             prfcap(f1cp) = 1
 130             prfcap(f2cp) = 1
 131 c
 132             do 2223 , nrofon = 1 , nbfonc
 133 c
 134 c           Pour le fils aine :
 135 c
 136               vafott(nrofon,q1,f1cp) = vafoen(nrofon,q4,prtrcn)
 137 c
 138               vafott(nrofon,q2,f1cp) = vafoen(nrofon,q2,prtrcn)
 139 c
 140               vafott(nrofon,q3,f1cp) = vafoen(nrofon,q3,prtrcn)
 141 c
 142               vafott(nrofon,q4,f1cp) =
 143      >       -unshu*vafoen(nrofon,q1,prtrcn)+
 144      >        trshu*vafoen(nrofon,q2,prtrcn)+
 145      >        trsqu*vafoen(nrofon,q4,prtrcn)
 146 c
 147               vafott(nrofon,q5,f1cp) = vafoen(nrofon,q5,prtrcn)
 148 c
 149               vafott(nrofon,q6,f1cp) =
 150      >       -unshu*(vafoen(nrofon,q1,prtrcn)+vafoen(nrofon,q2,prtrcn))+
 151      >        unsde*(vafoen(nrofon,q5,prtrcn)+vafoen(nrofon,q6,prtrcn))+
 152      >         unsqu*vafoen(nrofon,q4,prtrcn)
 153 c
 154 c           Pour le triangle fils NF+1 :
 155 c
 156               vafott(nrofon,q1,f2cp) = vafoen(nrofon,q1,prtrcn)
 157 c
 158               vafott(nrofon,q2,f2cp) = vafoen(nrofon,q4,prtrcn)
 159 c
 160               vafott(nrofon,q3,f2cp) = vafoen(nrofon,q3,prtrcn)
 161 c
 162               vafott(nrofon,q4,f2cp) =
 163      >        trshu*vafoen(nrofon,q1,prtrcn)-
 164      >        unshu*vafoen(nrofon,q2,prtrcn)+
 165      >        trsqu*vafoen(nrofon,q4,prtrcn)
 166 c
 167               vafott(nrofon,q5,f2cp) = vafott(nrofon,q6,f1cp)
 168 c
 169               vafott(nrofon,q6,f2cp) = vafoen(nrofon,q6,prtrcn)
 170 c
 171  2223       continue