src/tool/AP_Conversion/pcsqu2_1.h

   1 c           .................         .................
   2 c           .               .         .      . .      .
   3 c           .               .         .     .   .     .
   4 c           .               .         .    .     .    .
   5 c           .               .  ===>   .   .       .   .
   6 c           .               .         .  .         .  .
   7 c           .               .         . .           . .
   8 c           .               .         ..             ..
   9 c           .................         .................
  10 c
  11 c
  12             f1hp = -filqua(quhnp1)
  13             f1cp = ntrsca(f1hp)
  14             f2cp = ntrsca(f1hp+1)
  15             f3cp = ntrsca(f1hp+2)
  16             prftrp(f1cp) = 1
  17             prftrp(f2cp) = 1
  18             prftrp(f3cp) = 1
  19 c
  20             q1t = 1
  21             q2t = 2
  22             q3t = 3
  23             q4t = 4
  24             q5t = 5
  25             q6t = 6
  26
  27 c         Pour un decoupage par l'arete numero 1 :
  28
  29             if ( etanp1.eq.31 ) then
  30 c
  31             do 2221 , nrofon = 1 , nbfonc
  32 c
  33 c           Pour le triangle NT1 :
  34 c
  35               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q4,prqucn)
  36 c
  37               vatrtt(nrofon,q2t,f1cp) = vafoen(nrofon,q5,prqucn)
  38 c
  39               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q3,prqucn)
  40 c
  41               vatrtt(nrofon,q4t,f1cp) =
  42      >       -trssz*
  43      >              (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  44      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  45      >        trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
  46      >        unsqu*vafoen(nrofon,q6,prqucn)+
  47      >        trsqu*vafoen(nrofon,q8,prqucn)
  48 c
  49               vatrtt(nrofon,q5t,f1cp) =
  50      >       -trssz*
  51      >              (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
  52      >               vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
  53      >        trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
  54      >        unsqu*vafoen(nrofon,q8,prqucn)+
  55      >        trsqu*vafoen(nrofon,q6,prqucn)
  56 c
  57               vatrtt(nrofon,q6t,f1cp) = vafoen(nrofon,q7,prqucn)
  58 c
  59 c           Pour le triangle NT2 :
  60 c
  61               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q3,prqucn)
  62 c
  63               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q5,prqucn)
  64 c
  65               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q2,prqucn)
  66 c
  67               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
  68 c
  69               vatrtt(nrofon,q5t,f2cp) =
  70      >        trshu*vafoen(nrofon,q2,prqucn)-
  71      >        unshu*vafoen(nrofon,q1,prqucn)+
  72      >        trsqu*vafoen(nrofon,q5,prqucn)
  73 c
  74               vatrtt(nrofon,q6t,f2cp) = vafoen(nrofon,q6,prqucn)
  75 c
  76 c           Pour le triangle NT3 :
  77 c
  78               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q4,prqucn)
  79 c
  80               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prqucn)
  81 c
  82               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q5,prqucn)
  83 c
  84               vatrtt(nrofon,q4t,f3cp) = vafoen(nrofon,q8,prqucn)
  85 c
  86               vatrtt(nrofon,q5t,f3cp) =
  87      >        trshu*vafoen(nrofon,q1,prqucn)-
  88      >        unshu*vafoen(nrofon,q2,prqucn)+
  89      >        trsqu*vafoen(nrofon,q5,prqucn)
  90 c
  91               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
  92 c
  93  2221       continue
  94
  95 c         Pour un decoupage par l'arete numero 2 :
  96
  97           elseif ( etanp1.eq.32 ) then
  98 c
  99             do 2222 , nrofon = 1 , nbfonc
 100 c
 101 c           Pour le triangle NT1 :
 102 c
 103               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prqucn)
 104 c
 105               vatrtt(nrofon,q2t,f1cp) = vafoen(nrofon,q6,prqucn)
 106 c
 107               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q4,prqucn)
 108 c
 109               vatrtt(nrofon,q4t,f1cp) =
 110      >       -trssz*
 111      >            (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 112      >             vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 113      >        trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
 114      >        unsqu*vafoen(nrofon,q7,prqucn)+
 115      >        trsqu*vafoen(nrofon,q5,prqucn)
 116 c
 117               vatrtt(nrofon,q5t,f1cp) =
 118      >       -trssz*
 119      >            (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 120      >             vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 121      >        trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
 122      >        unsqu*vafoen(nrofon,q5,prqucn)+
 123      >        trsqu*vafoen(nrofon,q7,prqucn)
 124 c
 125               vatrtt(nrofon,q6t,f1cp) = vafoen(nrofon,q8,prqucn)
 126 c
 127 c           Pour le triangle NT2 :
 128 c
 129               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q4,prqucn)
 130 c
 131               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q6,prqucn)
 132 c
 133               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q3,prqucn)
 134 c
 135               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
 136 c
 137               vatrtt(nrofon,q5t,f2cp) =
 138      >        trshu*vafoen(nrofon,q3,prqucn)-
 139      >        unshu*vafoen(nrofon,q2,prqucn)+
 140      >        trsqu*vafoen(nrofon,q6,prqucn)
 141 c
 142               vatrtt(nrofon,q6t,f2cp) = vafoen(nrofon,q7,prqucn)
 143 c
 144 c           Pour le triangle NT3 :
 145 c
 146               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prqucn)
 147 c
 148               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q2,prqucn)
 149 c
 150               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q6,prqucn)
 151 c
 152               vatrtt(nrofon,q4t,f3cp) = vafoen(nrofon,q5,prqucn)
 153 c
 154               vatrtt(nrofon,q5t,f3cp) =
 155      >        trshu*vafoen(nrofon,q2,prqucn)-
 156      >        unshu*vafoen(nrofon,q3,prqucn)+
 157      >        trsqu*vafoen(nrofon,q6,prqucn)
 158 c
 159               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 160 c
 161  2222       continue
 162
 163 c         Pour un decoupage par l'arete numero 3 :
 164
 165           elseif ( etanp1.eq.33 ) then
 166 c
 167             do 2223 , nrofon = 1 , nbfonc
 168 c
 169 c           Pour le triangle NT1 :
 170 c
 171               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q2,prqucn)
 172 c
 173               vatrtt(nrofon,q2t,f1cp) = vafoen(nrofon,q7,prqucn)
 174 c
 175               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prqucn)
 176 c
 177               vatrtt(nrofon,q4t,f1cp) =
 178      >       -trssz*
 179      >           (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 180      >            vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 181      >      trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
 182      >        unsqu*vafoen(nrofon,q8,prqucn)+
 183      >        trsqu*vafoen(nrofon,q6,prqucn)
 184 c
 185               vatrtt(nrofon,q5t,f1cp) =
 186      >       -trssz*
 187      >            (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 188      >             vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 189      >      trshu*(vafoen(nrofon,q5,prqucn)+vafoen(nrofon,q7,prqucn))+
 190      >        unsqu*vafoen(nrofon,q6,prqucn)+
 191      >        trsqu*vafoen(nrofon,q8,prqucn)
 192 c
 193               vatrtt(nrofon,q6t,f1cp) = vafoen(nrofon,q5,prqucn)
 194 c
 195 c           Pour le triangle NT2 :
 196 c
 197               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prqucn)
 198 c
 199               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q7,prqucn)
 200 c
 201               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q4,prqucn)
 202 c
 203               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
 204 c
 205               vatrtt(nrofon,q5t,f2cp) =
 206      >        trshu*vafoen(nrofon,q4,prqucn)-
 207      >        unshu*vafoen(nrofon,q3,prqucn)+
 208      >        trsqu*vafoen(nrofon,q7,prqucn)
 209 c
 210               vatrtt(nrofon,q6t,f2cp) = vafoen(nrofon,q8,prqucn)
 211 c
 212 c           Pour le triangle NT3 :
 213 c
 214               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q2,prqucn)
 215 c
 216               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q3,prqucn)
 217 c
 218               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q7,prqucn)
 219 c
 220               vatrtt(nrofon,q4t,f3cp) = vafoen(nrofon,q6,prqucn)
 221 c
 222               vatrtt(nrofon,q5t,f3cp) =
 223      >        trshu*vafoen(nrofon,q3,prqucn)-
 224      >        unshu*vafoen(nrofon,q4,prqucn)+
 225      >        trsqu*vafoen(nrofon,q7,prqucn)
 226 c
 227               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 228 c
 229  2223       continue
 230
 231 c         Pour un decoupage par l'arete numero 4 :
 232
 233           elseif ( etanp1.eq.34 ) then
 234 c
 235             do 2224 , nrofon = 1 , nbfonc
 236 c
 237 c           Pour le triangle NT1 :
 238 c
 239               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q3,prqucn)
 240 c
 241               vatrtt(nrofon,q2t,f1cp) = vafoen(nrofon,q8,prqucn)
 242 c
 243               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q2,prqucn)
 244 c
 245               vatrtt(nrofon,q4t,f1cp) =
 246      >       -trssz*
 247      >          (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 248      >           vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 249      >      trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
 250      >      unsqu*vafoen(nrofon,q5,prqucn)+
 251      >      trsqu*vafoen(nrofon,q7,prqucn)
 252 c
 253               vatrtt(nrofon,q5t,f1cp) =
 254      >       -trssz*
 255      >          (vafoen(nrofon,q1,prqucn)+vafoen(nrofon,q2,prqucn)+
 256      >           vafoen(nrofon,q3,prqucn)+vafoen(nrofon,q4,prqucn))+
 257      >      trshu*(vafoen(nrofon,q6,prqucn)+vafoen(nrofon,q8,prqucn))+
 258      >      unsqu*vafoen(nrofon,q7,prqucn)+
 259      >      trsqu*vafoen(nrofon,q5,prqucn)
 260 c
 261               vatrtt(nrofon,q6t,f1cp) = vafoen(nrofon,q6,prqucn)
 262 c
 263 c           Pour le triangle NT2 :
 264 c
 265               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q2,prqucn)
 266 c
 267               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q8,prqucn)
 268 c
 269               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prqucn)
 270 c
 271               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
 272 c
 273               vatrtt(nrofon,q5t,f2cp) =
 274      >        trshu*vafoen(nrofon,q1,prqucn)-
 275      >        unshu*vafoen(nrofon,q4,prqucn)+
 276      >        trsqu*vafoen(nrofon,q8,prqucn)
 277 c
 278               vatrtt(nrofon,q6t,f2cp) = vafoen(nrofon,q5,prqucn)
 279 c
 280 c           Pour le triangle NT3 :
 281 c
 282               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q3,prqucn)
 283 c
 284               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q4,prqucn)
 285 c
 286               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q8,prqucn)
 287 c
 288               vatrtt(nrofon,q4t,f3cp) = vafoen(nrofon,q7,prqucn)
 289 c
 290               vatrtt(nrofon,q5t,f3cp) =
 291      >        trshu*vafoen(nrofon,q4,prqucn)-
 292      >        unshu*vafoen(nrofon,q1,prqucn)+
 293      >        trsqu*vafoen(nrofon,q8,prqucn)
 294 c
 295               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 296 c
 297  2224       continue
 298 c
 299            endif