src/tool/AP_Conversion/pcsqu2_3.h

   1 c           .................         .................
   2 c           .      . .      .         .               .
   3 c           .     .   .     .         .               .
   4 c           .    .     .    .         .               .
   5 c           .   .       .   .  ===>   .               .
   6 c           .  .         .  .         .               .
   7 c           . .           . .         .               .
   8 c           ..             ..         .               .
   9 c           .................         .................
  10 c
  11             qucnp1 = nqusca(quhnp1)
  12             prfcap(qucnp1) = 1
  13 c
  14 c quadrangle predecoupe en 3 triangles par l'arete A1
  15 c
  16             if ( etan.eq.31 ) then
  17 c
  18             do 2311 , nrofon = 1 , nbfonc
  19 c
  20               vafott(nrofon,q1,qucnp1) = vatren(nrofon,q2t,prf3cn)
  21               vafott(nrofon,q2,qucnp1) = vatren(nrofon,q3t,prf2cn)
  22               vafott(nrofon,q3,qucnp1) =
  23      >                    unsde * ( vatren(nrofon,q1t,prf2cn)
  24      >                            + vatren(nrofon,q3t,prf1cn) )
  25               vafott(nrofon,q4,qucnp1) =
  26      >                    unsde * ( vatren(nrofon,q1t,prf3cn)
  27      >                            + vatren(nrofon,q1t,prf1cn) )
  28               vafott(nrofon,q5,qucnp1) =
  29      >                    unstr * ( vatren(nrofon,q2t,prf1cn)
  30      >                            + vatren(nrofon,q2t,prf2cn)
  31      >                            + vatren(nrofon,q3t,prf3cn))
  32               vafott(nrofon,q6,qucnp1) = vatren(nrofon,q6t,prf2cn)
  33               vafott(nrofon,q7,qucnp1) = vatren(nrofon,q6t,prf1cn)
  34               vafott(nrofon,q8,qucnp1) = vatren(nrofon,q4t,prf3cn)
  35 c
  36  2311       continue
  37 c
  38 c quadrangle predecoupe en 3 triangles par l'arete A2
  39 c
  40             elseif ( etan.eq.32 ) then
  41 c
  42             do 2312 , nrofon = 1 , nbfonc
  43 c
  44               vafott(nrofon,q1,qucnp1) =
  45      >                    unsde * ( vatren(nrofon,q1t,prf3cn)
  46      >                            + vatren(nrofon,q1t,prf1cn) )
  47               vafott(nrofon,q2,qucnp1) = vatren(nrofon,q2t,prf3cn)
  48               vafott(nrofon,q3,qucnp1) = vatren(nrofon,q3t,prf2cn)
  49               vafott(nrofon,q4,qucnp1) =
  50      >                    unsde * ( vatren(nrofon,q1t,prf2cn)
  51      >                            + vatren(nrofon,q3t,prf1cn) )
  52               vafott(nrofon,q5,qucnp1) = vatren(nrofon,q4t,prf3cn)
  53               vafott(nrofon,q6,qucnp1) =
  54      >                    unstr * ( vatren(nrofon,q2t,prf1cn)
  55      >                            + vatren(nrofon,q2t,prf2cn)
  56      >                            + vatren(nrofon,q3t,prf3cn))
  57               vafott(nrofon,q7,qucnp1) = vatren(nrofon,q6t,prf2cn)
  58               vafott(nrofon,q8,qucnp1) = vatren(nrofon,q6t,prf1cn)
  59 c
  60  2312        continue
  61 c
  62 c quadrangle predecoupe en 3 triangles par l'arete A3
  63 c
  64             elseif ( etan.eq.33 ) then
  65 c
  66             do 2313 , nrofon = 1 , nbfonc
  67
  68               vafott(nrofon,q1,qucnp1) =
  69      >                    unsde * ( vatren(nrofon,q1t,prf2cn)
  70      >                            + vatren(nrofon,q3t,prf1cn) )
  71               vafott(nrofon,q2,qucnp1) =
  72      >                    unsde * ( vatren(nrofon,q1t,prf3cn)
  73      >                            + vatren(nrofon,q1t,prf1cn) )
  74               vafott(nrofon,q3,qucnp1) = vatren(nrofon,q2t,prf3cn)
  75               vafott(nrofon,q4,qucnp1) = vatren(nrofon,q3t,prf2cn)
  76               vafott(nrofon,q5,qucnp1) = vatren(nrofon,q6t,prf1cn)
  77               vafott(nrofon,q6,qucnp1) = vatren(nrofon,q4t,prf3cn)
  78               vafott(nrofon,q7,qucnp1) =
  79      >                    unstr * ( vatren(nrofon,q2t,prf1cn)
  80      >                            + vatren(nrofon,q2t,prf2cn)
  81      >                            + vatren(nrofon,q3t,prf3cn))
  82               vafott(nrofon,q8,qucnp1) = vatren(nrofon,q6t,prf2cn)
  83 c
  84  2313       continue
  85 c
  86 c quadrangle predecoupe en 3 triangles par l'arete A4
  87 c
  88             elseif ( etan.eq.34 ) then
  89 c
  90             do 2314 , nrofon = 1 , nbfonc
  91               vafott(nrofon,q1,qucnp1) = vatren(nrofon,q3t,prf2cn)
  92               vafott(nrofon,q2,qucnp1) =
  93      >                    unsde * ( vatren(nrofon,q1t,prf2cn)
  94      >                            + vatren(nrofon,q3t,prf1cn) )
  95               vafott(nrofon,q3,qucnp1) =
  96      >                    unsde * ( vatren(nrofon,q1t,prf3cn)
  97      >                            + vatren(nrofon,q1t,prf1cn) )
  98               vafott(nrofon,q4,qucnp1) = vatren(nrofon,q2t,prf3cn)
  99               vafott(nrofon,q5,qucnp1) = vatren(nrofon,q6t,prf2cn)
 100               vafott(nrofon,q6,qucnp1) = vatren(nrofon,q6t,prf1cn)
 101               vafott(nrofon,q7,qucnp1) = vatren(nrofon,q4t,prf3cn)
 102               vafott(nrofon,q8,qucnp1) =
 103      >                    unstr * ( vatren(nrofon,q2t,prf1cn)
 104      >                            + vatren(nrofon,q2t,prf2cn)
 105      >                            + vatren(nrofon,q3t,prf3cn))
 106 c
 107  2314       continue
 108 c
 109             endif