src/tool/AP_Conversion/pcsqu2_6.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 c
  27 c    quadrangle predecoupe en 4 quad et decoupe en trois triangles
  28 c    a l'arete a1
  29
  30             if ( etanp1.eq.31 ) then
  31
  32             do 2511 , nrofon = 1 , nbfonc
  33 c
  34 c    dans le triangle NT1
  35 c
  36               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf4cn)
  37 c
  38               vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
  39      >                                        vafoen(nrofon,q4,prf2cn))
  40 c
  41               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf3cn)
  42 c
  43               vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
  44      >                                        vafoen(nrofon,q6,prf4cn))
  45 c
  46               vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
  47      >                                        vafoen(nrofon,q7,prf3cn))
  48 c
  49               vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
  50      >                                        vafoen(nrofon,q4,prf4cn))
  51 c
  52 c    dans le triangle NT2
  53 c
  54               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf3cn)
  55 c
  56               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf2cn)
  57 c
  58               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf2cn)
  59 c
  60               vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
  61      >                                        vafoen(nrofon,q7,prf3cn))
  62 c
  63               vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf2cn)
  64 c
  65               vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
  66      >                                        vafoen(nrofon,q4,prf3cn))
  67 c
  68 c    dans le triangle NT3
  69 c
  70               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf4cn)
  71 c
  72               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf1cn)
  73 c
  74               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf1cn)
  75 c
  76               vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
  77      >                                        vafoen(nrofon,q2,prf4cn))
  78 c
  79               vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf1cn)
  80 c
  81               vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
  82      >                                        vafoen(nrofon,q6,prf4cn))
  83 c
  84  2511       continue
  85 c
  86 c    quadrangle predecoupe en 4 quad et decoupe en trois triangles
  87 c    a l'arete a2
  88 c
  89             elseif ( etanp1.eq.32 ) then
  90 c
  91             do 2512 , nrofon = 1 , nbfonc
  92 c    dans le triangle NT1
  93 c
  94               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf1cn)
  95 c
  96               vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
  97      >                                        vafoen(nrofon,q4,prf3cn))
  98 c
  99               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf4cn)
 100 c
 101               vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
 102      >                                        vafoen(nrofon,q7,prf2cn))
 103 c
 104               vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
 105      >                                        vafoen(nrofon,q7,prf4cn))
 106 c
 107               vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
 108      >                                        vafoen(nrofon,q2,prf4cn))
 109 c
 110 c    dans le triangle NT2
 111 c
 112               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf4cn)
 113 c
 114               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf3cn)
 115 c
 116               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf3cn)
 117 c
 118               vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
 119      >                                        vafoen(nrofon,q7,prf4cn))
 120 c
 121               vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf3cn)
 122 c
 123               vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
 124      >                                        vafoen(nrofon,q4,prf4cn))
 125 c
 126 c    dans le triangle NT3
 127 c
 128               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf1cn)
 129 c
 130               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf2cn)
 131 c
 132               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf2cn)
 133 c
 134               vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
 135      >                                        vafoen(nrofon,q4,prf2cn))
 136 c
 137               vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf2cn)
 138 c
 139               vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
 140      >                                        vafoen(nrofon,q7,prf2cn))
 141 c
 142  2512       continue
 143 c
 144 c    quadrangle predecoupe en 4 quad et decoupe en trois triangles
 145 c    a l'arete a3
 146 c
 147             elseif ( etanp1.eq.33 ) then
 148 c
 149             do 2513 , nrofon = 1 , nbfonc
 150 c    dans le triangle NT1
 151 c
 152               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf2cn)
 153 c
 154               vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
 155      >                                        vafoen(nrofon,q4,prf4cn))
 156 c
 157               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf1cn)
 158 c
 159               vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
 160      >                                        vafoen(nrofon,q7,prf3cn))
 161 c
 162               vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
 163      >                                        vafoen(nrofon,q6,prf4cn))
 164 c
 165               vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
 166      >                                        vafoen(nrofon,q4,prf2cn))
 167 c
 168 c    dans le triangle NT2
 169 c
 170               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf1cn)
 171 c
 172               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf4cn)
 173 c
 174               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf4cn)
 175 c
 176               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
 177 c
 178               vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf4cn)
 179 c
 180               vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
 181      >                                        vafoen(nrofon,q2,prf4cn))
 182 c
 183 c    dans le triangle NT3
 184 c
 185               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf2cn)
 186 c
 187               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf3cn)
 188 c
 189               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf3cn)
 190 c
 191               vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
 192      >                                        vafoen(nrofon,q4,prf3cn))
 193 c
 194               vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf3cn)
 195 c
 196               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 197 c
 198  2513       continue
 199 c
 200 c    quadrangle predecoupe en 4 quad et decoupe en trois triangles
 201 c    a l'arete a4
 202 c
 203             elseif ( etanp1.eq.34 ) then
 204 c
 205             do 2514 , nrofon = 1 , nbfonc
 206 c    dans le triangle NT1
 207 c
 208               vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf3cn)
 209 c
 210               vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
 211      >                                         vafoen(nrofon,q2,prf4cn))
 212 c
 213               vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf2cn)
 214 c
 215               vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
 216      >                                        vafoen(nrofon,q7,prf4cn))
 217 c
 218               vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
 219      >                                        vafoen(nrofon,q7,prf2cn))
 220 c
 221               vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
 222      >                                        vafoen(nrofon,q4,prf3cn))
 223 c
 224 c    dans le triangle NT2
 225 c
 226               vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf2cn)
 227 c
 228               vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf1cn)
 229 c
 230               vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf1cn)
 231 c
 232               vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
 233 c
 234               vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf1cn)
 235 c
 236               vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
 237      >                                        vafoen(nrofon,q4,prf2cn))
 238 c
 239 c    dans le triangle NT3
 240 c
 241               vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf3cn)
 242 c
 243               vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf4cn)
 244 c
 245               vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf4cn)
 246 c
 247               vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
 248      >                                        vafoen(nrofon,q4,prf4cn))
 249 c
 250               vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf4cn)
 251 c
 252               vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
 253 c
 254  2514       continue
 255 c
 256             endif