2 if ( etanp1.ge.31 .and. etanp1.le.34 ) then
3 c ................. .................
7 c ................. ===> . . . .
11 c ................. .................
13 c quadrangle predecoupe en 4 quad et decoupe en trois triangles
16 f1hp = -filqua(quhnp1)
31 if (etanp1 .eq. 31) then
39 cgn write(6,*) 'etanp1', etanp1
40 cgn write(6,*) 'f1hn+2=',f1hn+2,
41 cgn > 'mod(anhetr(f1hn+2),10)',mod(anhetr(f1hn+2),10)
42 if ( mod(anhequ(f1hn),10).eq.0 ) then
45 elseif ( mod(anhequ(f1hn),10).eq.etanp1 ) then
55 cgn write(ulsort,*) '_7h A codret', codret
56 cgn write (ulsort,texte(langue,4)) 'n ', trhn
57 cgn write (ulsort,texte(langue,5)) 'n ', etan
58 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
59 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
61 cgn write(6,*) 'etanp1', etanp1
62 cgn write(6,*) 'mod(anhetr(f1hn+3),10)',mod(anhetr(f1hn+3),10)
63 if ( mod(anhequ(f1hn+1),10).eq.0 ) then
66 elseif ( mod(anhequ(f1hn+1),10).eq.etanp1 ) then
76 cgn write (ulsort,*) '_5h B codret', codret
77 cgn write (ulsort,texte(langue,4)) 'n ', trhn
78 cgn write (ulsort,texte(langue,5)) 'n ', etan
79 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
80 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
83 if ( g1.eq.0 .and. d1.eq.0 ) then
85 do 2611 , nrofon = 1 , nbfonc
87 c dans le triangle NT1
89 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf4cn)
91 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
92 > vafoen(nrofon,q4,prf2cn))
94 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf3cn)
96 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
97 > vafoen(nrofon,q6,prf4cn))
99 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
100 > vafoen(nrofon,q7,prf3cn))
102 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
103 > vafoen(nrofon,q4,prf4cn))
105 c dans le triangle NT2
107 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf3cn)
109 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q2,prf4cn)
111 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf2cn)
113 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
114 > vafoen(nrofon,q7,prf3cn))
116 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf2cn)
118 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
119 > vafoen(nrofon,q4,prf3cn))
121 c dans le triangle NT3
123 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf4cn)
125 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf1cn)
127 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf1cn)
129 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
130 > vafoen(nrofon,q2,prf4cn))
132 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf1cn)
134 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
135 > vafoen(nrofon,q6,prf4cn))
139 elseif ( g1.ne.0 .and. d1.eq.0 ) then
140 do 26111 , nrofon = 1 , nbfonc
142 c dans le triangle NT1
144 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf4cn)
146 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
147 > vafoen(nrofon,q4,prf2cn))
149 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf3cn)
151 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf4cn)
152 > +vatren(nrofon,q6,prfg1n))
154 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
155 > vafoen(nrofon,q7,prf3cn))
157 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)
158 > +vafoen(nrofon,q4,prf4cn))
160 c dans le triangle NT2
162 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf3cn)
164 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q2,prf4cn)
166 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf2cn)
168 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
169 > vafoen(nrofon,q7,prf3cn))
171 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf2cn)
173 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
174 > vafoen(nrofon,q4,prf3cn))
177 c dans le triangle NT3
179 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf4cn)
181 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
183 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfg2n)
185 vatrtt(nrofon,q4t,f3cp) = unstr*(vatren(nrofon,q1,prfg1n)
186 > +vatren(nrofon,q1,prfg3n)
187 > +vafoen(nrofon,q2,prf4cn))
189 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q2,prfg1n)
190 > +vatren(nrofon,q2,prfg2n)
191 > +vatren(nrofon,q3,prfg3n))
193 vatrtt(nrofon,q6t,f3cp) = unsde*(vatren(nrofon,q6,prfg1n)
194 > +vafoen(nrofon,q6,prf4cn))
198 elseif ( g1.eq.0 .and. d1.ne.0 ) then
199 do 26112 , nrofon = 1 , nbfonc
201 c dans le triangle NT1
203 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf4cn)
205 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf1cn)
206 > +vatren(nrofon,q2,prfd3n))
208 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf3cn)
210 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)
211 > +vatren(nrofon,q6,prf4cn))
213 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf3cn)
214 > +vatren(nrofon,q6,prfd1n))
216 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)
217 > +vafoen(nrofon,q4,prf4cn))
219 c dans le triangle NT2
221 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf3cn)
223 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q3,prfg2n)
225 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
227 vatrtt(nrofon,q4t,f2cp) = unsde*(vatren(nrofon,q6,prfd1n)
228 > +vafoen(nrofon,q7,prf3cn))
230 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfd1n)
231 > +vatren(nrofon,q5,prfd2n)
232 > +vatren(nrofon,q3,prfd3n))
234 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf3cn)
235 > +vatren(nrofon,q3,prfd1n)
236 > +vatren(nrofon,q1,prfd2n))
238 c dans le triangle NT3
240 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf4cn)
242 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf1cn)
244 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf1cn)
246 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
247 > vafoen(nrofon,q2,prf4cn))
249 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf1cn)
251 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
252 > vafoen(nrofon,q6,prf4cn))
256 elseif ( g1.ne.0 .and. d1.ne.0 ) then
257 do 26113 , nrofon = 1 , nbfonc
259 c dans le triangle NT1
261 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf4cn)
263 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)
264 > +vatren(nrofon,q2,prfd3n))
266 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf3cn)
268 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf4cn)
269 > +vatren(nrofon,q6,prfg1n))
271 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf3cn)
272 > +vatren(nrofon,q6,prfd1n))
274 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)
275 > +vafoen(nrofon,q4,prf4cn))
277 c dans le triangle NT2
279 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf3cn)
281 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
283 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
285 vatrtt(nrofon,q4t,f2cp) = unsde*(vatren(nrofon,q6,prfd1n)
286 > +vafoen(nrofon,q7,prf3cn))
288 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfd1n)
289 > +vatren(nrofon,q2,prfd2n)
290 > +vatren(nrofon,q3,prfd3n))
292 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf3cn)
293 > +vatren(nrofon,q3,prfd1n)
294 > +vatren(nrofon,q1,prfd2n))
296 c dans le triangle NT3
298 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf4cn)
300 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
302 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfg2n)
304 vatrtt(nrofon,q4t,f3cp) = unstr*(vatren(nrofon,q1,prfg1n)
305 > +vatren(nrofon,q1,prfg3n)
306 > +vafoen(nrofon,q2,prf4cn))
308 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q2,prfg1n)
309 > +vatren(nrofon,q2,prfg2n)
310 > +vatren(nrofon,q3,prfg3n))
312 vatrtt(nrofon,q6t,f3cp) = unsde*(vatren(nrofon,q6,prfg1n)
313 > +vafoen(nrofon,q6,prf4cn))
318 c quadrangle predecoupe en 4 quad et decoupe en trois triangles
321 elseif (etanp1 .eq. 32) then
324 prf1cn = prfcan(f1cn)
325 f4cn = nqueca(f1hn+3)
326 prf4cn = prfcan(f4cn)
329 cgn write(6,*) 'etanp1', etanp1
330 cgn write(6,*) 'f1hn+2=',f1hn+2,
331 cgn > 'mod(anhetr(f1hn+1),10)',mod(anhetr(f1hn+1),10)
332 if ( mod(anhequ(f1hn+1),10).eq.0 ) then
333 f2cn = nqueca(f1hn+1)
334 prf2cn = prfcan(f2cn)
335 elseif ( mod(anhequ(f1hn+1),10).eq.etanp1 ) then
345 cgn write(ulsort,*) '_7h A codret', codret
346 cgn write (ulsort,texte(langue,4)) 'n ', trhn
347 cgn write (ulsort,texte(langue,5)) 'n ', etan
348 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
349 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
351 cgn write(6,*) 'etanp1', etanp1
352 cgn write(6,*) 'mod(anhetr(f1hn+2),10)',mod(anhetr(f1hn+2),10)
353 if ( mod(anhequ(f1hn+2),10).eq.0 ) then
354 f3cn = nqueca(f1hn+2)
355 prf3cn = prfcan(f3cn)
356 elseif ( mod(anhequ(f1hn+2),10).eq.etanp1 ) then
366 cgn write (ulsort,*) '_5h B codret', codret
367 cgn write (ulsort,texte(langue,4)) 'n ', trhn
368 cgn write (ulsort,texte(langue,5)) 'n ', etan
369 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
370 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
374 if ( g1.eq.0 .and. d1.eq.0 ) then
376 do 2612 , nrofon = 1 , nbfonc
377 c dans le triangle NT1
379 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf1cn)
381 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
382 > vafoen(nrofon,q4,prf3cn))
384 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf4cn)
386 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
387 > vafoen(nrofon,q7,prf2cn))
389 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
390 > vafoen(nrofon,q7,prf4cn))
392 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
393 > vafoen(nrofon,q2,prf4cn))
395 c dans le triangle NT2
397 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf4cn)
399 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf3cn)
401 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf3cn)
403 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
404 > vafoen(nrofon,q7,prf4cn))
406 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf3cn)
408 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
409 > vafoen(nrofon,q4,prf4cn))
411 c dans le triangle NT3
413 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf1cn)
415 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf2cn)
417 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf2cn)
419 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
420 > vafoen(nrofon,q4,prf2cn))
422 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf2cn)
424 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
425 > vafoen(nrofon,q7,prf2cn))
429 elseif ( g1.ne.0 .and. d1.eq.0 ) then
430 do 26121 , nrofon = 1 , nbfonc
432 c dans le triangle NT1
434 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf1cn)
436 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
437 > vafoen(nrofon,q4,prf3cn))
439 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf4cn)
441 vatrtt(nrofon,q4t,f1cp) = unsde*(vatren(nrofon,q6,prfg1n)+
442 > vafoen(nrofon,q6,prf1cn))
444 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
445 > vafoen(nrofon,q7,prf4cn))
447 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
448 > vafoen(nrofon,q2,prf4cn))
450 c dans le triangle NT2
452 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf4cn)
454 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf3cn)
456 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf3cn)
458 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
459 > vafoen(nrofon,q7,prf4cn))
461 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf3cn)
463 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
464 > vafoen(nrofon,q4,prf4cn))
466 c dans le triangle NT3
468 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf1cn)
470 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
472 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfg2n)
474 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf1cn)
475 > +vatren(nrofon,q1,prfg3n)
476 > +vatren(nrofon,q1,prfg1n))
478 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q3,prfg3n)
479 > +vatren(nrofon,q2,prfg2n)
480 > +vatren(nrofon,q1,prfg1n))
482 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
483 > vafoen(nrofon,q6,prfg1n))
487 elseif ( g1.eq.0 .and. d1.ne.0 ) then
488 do 26122 , nrofon = 1 , nbfonc
490 c dans le triangle NT1
492 vatrtt(nrofon,q1t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
493 > vatren(nrofon,q2,prfd3n))
495 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
496 > vatren(nrofon,q2,prfd3n))
498 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf4cn)
500 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q7,prf2cn)+
501 > vafoen(nrofon,q6,prf1cn))
503 vatrtt(nrofon,q5t,f1cp) = unsde*(vatren(nrofon,q6,prfd1n)+
504 > vafoen(nrofon,q7,prf4cn))
506 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
507 > vafoen(nrofon,q2,prf4cn))
509 c dans le triangle NT2
511 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf4cn)
513 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
515 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
517 vatrtt(nrofon,q4t,f2cp) = unsde*(vatren(nrofon,q6,prfd1n)+
518 > vafoen(nrofon,q7,prf4cn))
520 vatrtt(nrofon,q5t,f2cp) = unstr*(vafoen(nrofon,q4,prf4cn)
521 > +vatren(nrofon,q3,prfd1n)
522 > +vatren(nrofon,q1,prfd2n))
524 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q3,prfd3n)
525 > +vatren(nrofon,q2,prfd2n)
526 > +vatren(nrofon,q2,prfd1n))
528 c dans le triangle NT3
530 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf1cn)
532 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf2cn)
534 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf2cn)
536 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
537 > vafoen(nrofon,q4,prf2cn))
539 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf2cn)
541 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
542 > vafoen(nrofon,q7,prf2cn))
546 elseif ( g1.ne.0 .and. d1.ne.0 ) then
547 do 26123 , nrofon = 1 , nbfonc
549 c dans le triangle NT1
551 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf1cn)
553 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
554 > vatren(nrofon,q2,prfd3n))
556 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf4cn)
558 vatrtt(nrofon,q4t,f1cp) = unsde*(vatren(nrofon,q6,prfg1n)+
559 > vafoen(nrofon,q6,prf1cn))
561 vatrtt(nrofon,q5t,f1cp) = unsde*(vatren(nrofon,q6,prfd1n)+
562 > vafoen(nrofon,q7,prf4cn))
564 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
565 > vafoen(nrofon,q2,prf4cn))
567 c dans le triangle NT2
569 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf4cn)
571 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
573 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
575 vatrtt(nrofon,q4t,f2cp) = unsde*(vatren(nrofon,q6,prfd1n)+
576 > vafoen(nrofon,q7,prf4cn))
578 vatrtt(nrofon,q5t,f2cp) = unstr*(vafoen(nrofon,q3,prfd3n)
579 > +vatren(nrofon,q2,prfd2n)
580 > +vatren(nrofon,q2,prfd1n))
582 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf4cn)
583 > +vatren(nrofon,q3,prfd1n)
584 > +vatren(nrofon,q1,prfd2n))
585 c dans le triangle NT3
587 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf1cn)
589 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
591 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfg2n)
593 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf1cn)
594 > +vatren(nrofon,q1,prfg3n)
595 > +vatren(nrofon,q1,prfg1n))
597 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q3,prfg3n)
598 > +vatren(nrofon,q2,prfg2n)
599 > +vatren(nrofon,q1,prfg1n))
601 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
602 > vafoen(nrofon,q6,prfg1n))
607 c quadrangle predecoupe en 4 quad et decoupe en trois triangles
610 elseif (etanp1 .eq. 33) then
613 prf1cn = prfcan(f1cn)
614 f2cn = nqueca(f1hn+1)
615 prf2cn = prfcan(f2cn)
618 cgn write(6,*) 'etanp1', etanp1
619 cgn write(6,*) 'f1hn+2=',f1hn+2,
620 cgn > 'mod(anhetr(f1hn+2),10)',mod(anhetr(f1hn+2),10)
621 if ( mod(anhequ(f1hn+2),10).eq.0 ) then
622 f3cn = nqueca(f1hn+2)
623 prf3cn = prfcan(f3cn)
624 elseif ( mod(anhequ(f1hn+2),10).eq.etanp1 ) then
634 cgn write(ulsort,*) '_7h A codret', codret
635 cgn write (ulsort,texte(langue,4)) 'n ', trhn
636 cgn write (ulsort,texte(langue,5)) 'n ', etan
637 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
638 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
640 cgn write(6,*) 'etanp1', etanp1
641 cgn write(6,*) 'mod(anhetr(f1hn+3),10)',mod(anhetr(f1hn+3),10)
642 if ( mod(anhequ(f1hn+3),10).eq.0 ) then
643 f4cn = nqueca(f1hn+3)
644 prf4cn = prfcan(f4cn)
645 elseif ( mod(anhequ(f1hn+3),10).eq.etanp1 ) then
655 cgn write (ulsort,*) '_5h B codret', codret
656 cgn write (ulsort,texte(langue,4)) 'n ', trhn
657 cgn write (ulsort,texte(langue,5)) 'n ', etan
658 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
659 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
663 if ( g1.eq.0 .and. d1.eq.0 ) then
665 do 2613 , nrofon = 1 , nbfonc
666 c dans le triangle NT1
668 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf2cn)
670 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
671 > vafoen(nrofon,q4,prf4cn))
673 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf1cn)
675 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
676 > vafoen(nrofon,q7,prf3cn))
678 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
679 > vafoen(nrofon,q6,prf4cn))
681 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
682 > vafoen(nrofon,q4,prf2cn))
684 c dans le triangle NT2
686 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf1cn)
688 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf4cn)
690 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf4cn)
692 vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
694 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf4cn)
696 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
697 > vafoen(nrofon,q2,prf4cn))
699 c dans le triangle NT3
701 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf2cn)
703 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf3cn)
705 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf3cn)
707 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
708 > vafoen(nrofon,q4,prf3cn))
710 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf3cn)
712 vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
716 elseif ( g1.ne.0 .and. d1.eq.0 ) then
717 do 26131 , nrofon = 1 , nbfonc
719 c dans le triangle NT1
721 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf2cn)
723 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
724 > vafoen(nrofon,q4,prf4cn))
726 vatrtt(nrofon,q3t,f1cp) = vatren(nrofon,q3,prfg2n)
728 vatrtt(nrofon,q4t,f1cp) = unstr*(vafoen(nrofon,q2,prf2cn)+
729 > vatren(nrofon,q1,prfg1n)+
730 > vatren(nrofon,q1,prfg3n))
732 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
733 > vafoen(nrofon,q6,prf4cn))
735 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
736 > vatren(nrofon,q6,prfg1n))
738 c dans le triangle NT2
740 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf1cn)
742 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf4cn)
744 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf4cn)
746 vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
748 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf4cn)
750 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
751 > vafoen(nrofon,q2,prf4cn))
753 c dans le triangle NT3
755 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf2cn)
757 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
759 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q3,prfg2n)
761 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf2cn)+
762 > vatren(nrofon,q1,prfg1n)+
763 > vatren(nrofon,q1,prfg3n))
765 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q2,prfg1n)+
766 > vatren(nrofon,q2,prfg2n)+
767 > vatren(nrofon,q3,prfg3n))
769 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
770 > vatren(nrofon,q6,prfg1n))
774 elseif ( g1.eq.0 .and. d1.ne.0 ) then
775 do 26132 , nrofon = 1 , nbfonc
777 c dans le triangle NT1
779 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf2cn)
781 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
782 > vatren(nrofon,q2,prfd3n))
784 vatrtt(nrofon,q3t,f1cp) = vatren(nrofon,q3,prfg2n)
786 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
787 > vatren(nrofon,q7,prf3cn))
789 vatrtt(nrofon,q5t,f1cp) = unstr*(vatren(nrofon,q3,prfg3n)+
790 > vatren(nrofon,q2,prfg1n)+
791 > vatren(nrofon,q2,prfg2n))
793 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
794 > vatren(nrofon,q6,prfg1n))
796 c dans le triangle NT2
798 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf1cn)
800 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
802 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
804 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
805 > vatren(nrofon,q6,prfd1n))
807 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfd2n)+
808 > vatren(nrofon,q2,prfd1n)+
809 > vatren(nrofon,q3,prfd3n))
811 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf1cn)+
812 > vatren(nrofon,q3,prfd1n)+
813 > vatren(nrofon,q1,prfd2n))
815 c dans le triangle NT3
817 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf2cn)
819 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf3cn)
821 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf3cn)
823 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
824 > vafoen(nrofon,q4,prf3cn))
826 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf3cn)
828 vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
832 elseif ( g1.ne.0 .and. d1.ne.0 ) then
833 do 26133 , nrofon = 1 , nbfonc
834 c dans le triangle NT1
836 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf2cn)
838 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
839 > vatren(nrofon,q2,prfd3n))
841 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf1cn)
843 vatrtt(nrofon,q4t,f1cp) = unstr*(vafoen(nrofon,q2,prf2cn)+
844 > vatren(nrofon,q1,prfg1n)+
845 > vatren(nrofon,q1,prfg3n))
847 vatrtt(nrofon,q5t,f1cp) = unstr*(vatren(nrofon,q3,prfg3n)+
848 > vatren(nrofon,q2,prfg1n)+
849 > vatren(nrofon,q2,prfg2n))
851 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
852 > vatren(nrofon,q6,prfg1n))
854 c dans le triangle NT2
856 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf1cn)
858 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
860 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
862 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q7,prf1cn)+
863 > vatren(nrofon,q6,prfd1n))
865 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfd2n)+
866 > vatren(nrofon,q2,prfd1n)+
867 > vatren(nrofon,q3,prfd3n))
869 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf1cn)+
870 > vatren(nrofon,q3,prfd1n)+
871 > vatren(nrofon,q1,prfd2n))
873 c dans le triangle NT3
875 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf2cn)
877 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfg3n)
879 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q3,prfg2n)
881 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf2cn)+
882 > vatren(nrofon,q1,prfg1n)+
883 > vatren(nrofon,q1,prfg3n))
885 vatrtt(nrofon,q5t,f3cp) = unstr*(vatren(nrofon,q2,prfg1n)+
886 > vatren(nrofon,q2,prfg2n)+
887 > vatren(nrofon,q3,prfg3n))
889 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf2cn)+
890 > vatren(nrofon,q6,prfg1n))
895 c quadrangle predecoupe en 4 quad et decoupe en trois triangles
898 elseif (etanp1 .eq. 34) then
900 f2cn = nqueca(f1hn+1)
901 prf2cn = prfcan(f2cn)
902 f3cn = nqueca(f1hn+2)
903 prf3cn = prfcan(f3cn)
906 cgn write(6,*) 'etanp1', etanp1
907 cgn write(6,*) 'f1hn+2=',f1hn+2,
908 cgn > 'mod(anhetr(f1hn),10)',mod(anhetr(f1hn),10)
909 if ( mod(anhequ(f1hn),10).eq.0 ) then
911 prf1cn = prfcan(f1cn)
912 elseif ( mod(anhequ(f1hn),10).eq.etanp1 ) then
922 cgn write(ulsort,*) '_7h A codret', codret
923 cgn write (ulsort,texte(langue,4)) 'n ', trhn
924 cgn write (ulsort,texte(langue,5)) 'n ', etan
925 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
926 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
928 cgn write(6,*) 'etanp1', etanp1
929 cgn write(6,*) 'mod(anhetr(f1hn+3),10)',mod(anhetr(f1hn+3),10)
930 if ( mod(anhequ(f1hn+3),10).eq.0 ) then
931 f4cn = nqueca(f1hn+3)
932 prf4cn = prfcan(f4cn)
933 elseif ( mod(anhequ(f1hn+3),10).eq.etanp1 ) then
943 cgn write (ulsort,*) '_5h B codret', codret
944 cgn write (ulsort,texte(langue,4)) 'n ', trhn
945 cgn write (ulsort,texte(langue,5)) 'n ', etan
946 cgn write (ulsort,texte(langue,4)) 'n+1', trhnp1
947 cgn write (ulsort,texte(langue,5)) 'n+1', etanp1
951 if ( g1.eq.0 .and. d1.eq.0 ) then
953 do 2614 , nrofon = 1 , nbfonc
954 c dans le triangle NT1
956 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf3cn)
958 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
959 > vafoen(nrofon,q2,prf4cn))
961 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf2cn)
963 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
964 > vafoen(nrofon,q7,prf4cn))
966 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
967 > vafoen(nrofon,q7,prf2cn))
969 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
970 > vafoen(nrofon,q4,prf3cn))
972 c dans le triangle NT2
974 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf2cn)
976 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf1cn)
978 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf1cn)
980 vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
982 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf1cn)
984 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q2,prf1cn)+
985 > vafoen(nrofon,q4,prf2cn))
987 c dans le triangle NT3
989 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf3cn)
991 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf4cn)
993 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf4cn)
995 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
996 > vafoen(nrofon,q4,prf4cn))
998 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf4cn)
1000 vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
1004 elseif ( g1.ne.0 .and. d1.eq.0 ) then
1005 do 26141 , nrofon = 1 , nbfonc
1007 c dans le triangle NT1
1009 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf3cn)
1011 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q2,prf4cn)+
1012 > vatren(nrofon,q2,prfg3n))
1014 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf2cn)
1016 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
1017 > vafoen(nrofon,q7,prf4cn))
1019 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf2cn)+
1020 > vatren(nrofon,q6,prfg1n))
1022 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
1023 > vafoen(nrofon,q4,prf3cn))
1025 c dans le triangle NT2
1027 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf2cn)
1029 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfg3n)
1031 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfg2n)
1033 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q7,prf2cn)+
1034 > vatren(nrofon,q6,prfg1n))
1036 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfg1n)+
1037 > vatren(nrofon,q2,prfg2n)+
1038 > vatren(nrofon,q3,prfg3n))
1040 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf2cn)+
1041 > vatren(nrofon,q1,prfg2n)+
1042 > vatren(nrofon,q3,prfg1n))
1044 c dans le triangle NT3
1046 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf3cn)
1048 vatrtt(nrofon,q2t,f3cp) = vafoen(nrofon,q1,prf4cn)
1050 vatrtt(nrofon,q3t,f3cp) = vafoen(nrofon,q2,prf4cn)
1052 vatrtt(nrofon,q4t,f3cp) = unsde*(vafoen(nrofon,q2,prf3cn)+
1053 > vafoen(nrofon,q4,prf4cn))
1055 vatrtt(nrofon,q5t,f3cp) = vafoen(nrofon,q5,prf4cn)
1057 vatrtt(nrofon,q6t,f3cp) = vatrtt(nrofon,q4t,f1cp)
1061 elseif ( g1.eq.0 .and. d1.ne.0 ) then
1062 do 26142 , nrofon = 1 , nbfonc
1064 c dans le triangle NT1
1066 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf3cn)
1068 vatrtt(nrofon,q2t,f1cp) = unsde*(vafoen(nrofon,q4,prf1cn)+
1069 > vatren(nrofon,q3,prfd2n))
1071 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf2cn)
1073 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
1074 > vatren(nrofon,q6,prfd1n))
1076 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q6,prf1cn)+
1077 > vafoen(nrofon,q7,prf2cn))
1079 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
1080 > vafoen(nrofon,q4,prf3cn))
1082 c dans le triangle NT2
1084 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf2cn)
1086 vatrtt(nrofon,q2t,f2cp) = vafoen(nrofon,q4,prf1cn)
1088 vatrtt(nrofon,q3t,f2cp) = vafoen(nrofon,q1,prf1cn)
1090 vatrtt(nrofon,q4t,f2cp) = vatrtt(nrofon,q5t,f1cp)
1092 vatrtt(nrofon,q5t,f2cp) = vafoen(nrofon,q8,prf1cn)
1094 vatrtt(nrofon,q6t,f2cp) = unsde*(vafoen(nrofon,q1,prf1cn)+
1095 > vafoen(nrofon,q4,prf2cn))
1097 c dans le triangle NT3
1099 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf3cn)
1101 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfd3n)
1103 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfd2n)
1105 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf3cn)+
1106 > vatren(nrofon,q1,prfd1n)+
1107 > vatren(nrofon,q1,prfd3n))
1109 vatrtt(nrofon,q5t,f3cp) = unstr*(vafoen(nrofon,q2,prfd1n)+
1110 > vatren(nrofon,q2,prfd2n)+
1111 > vatren(nrofon,q3,prfd3n))
1113 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
1114 > vatren(nrofon,q6,prfd1n))
1118 elseif ( g1.ne.0 .and. d1.ne.0 ) then
1119 do 26143 , nrofon = 1 , nbfonc
1120 c dans le triangle NT1
1122 vatrtt(nrofon,q1t,f1cp) = vafoen(nrofon,q1,prf3cn)
1124 vatrtt(nrofon,q2t,f1cp) = unsde*(vatren(nrofon,q3,prfg2n)+
1125 > vatren(nrofon,q2,prfd3n))
1127 vatrtt(nrofon,q3t,f1cp) = vafoen(nrofon,q1,prf2cn)
1129 vatrtt(nrofon,q4t,f1cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
1130 > vatren(nrofon,q6,prfg1n))
1132 vatrtt(nrofon,q5t,f1cp) = unsde*(vafoen(nrofon,q7,prf2cn)+
1133 > vatren(nrofon,q6,prfd1n))
1135 vatrtt(nrofon,q6t,f1cp) = unsde*(vafoen(nrofon,q2,prf2cn)+
1136 > vafoen(nrofon,q4,prf3cn))
1138 c dans le triangle NT2
1140 vatrtt(nrofon,q1t,f2cp) = vafoen(nrofon,q1,prf2cn)
1142 vatrtt(nrofon,q2t,f2cp) = vatren(nrofon,q2,prfd3n)
1144 vatrtt(nrofon,q3t,f2cp) = vatren(nrofon,q3,prfd2n)
1146 vatrtt(nrofon,q4t,f2cp) = unsde*(vafoen(nrofon,q7,prf2cn)+
1147 > vatren(nrofon,q6,prfg1n))
1149 vatrtt(nrofon,q5t,f2cp) = unstr*(vatren(nrofon,q2,prfg1n)+
1150 > vatren(nrofon,q2,prfg2n)+
1151 > vatren(nrofon,q3,prfg3n))
1153 vatrtt(nrofon,q6t,f2cp) = unstr*(vafoen(nrofon,q4,prf2cn)+
1154 > vatren(nrofon,q1,prfg2n)+
1155 > vatren(nrofon,q3,prfg1n))
1157 c dans le triangle NT3
1159 vatrtt(nrofon,q1t,f3cp) = vafoen(nrofon,q1,prf3cn)
1161 vatrtt(nrofon,q2t,f3cp) = vatren(nrofon,q2,prfd3n)
1163 vatrtt(nrofon,q3t,f3cp) = vatren(nrofon,q3,prfd2n)
1165 vatrtt(nrofon,q4t,f3cp) = unstr*(vafoen(nrofon,q2,prf3cn)+
1166 > vatren(nrofon,q1,prfd1n)+
1167 > vatren(nrofon,q1,prfd3n))
1169 vatrtt(nrofon,q5t,f3cp) = unstr*(vafoen(nrofon,q2,prfd1n)+
1170 > vatren(nrofon,q2,prfd2n)+
1171 > vatren(nrofon,q3,prfd3n))
1173 vatrtt(nrofon,q6t,f3cp) = unsde*(vafoen(nrofon,q6,prf3cn)+
1174 > vatren(nrofon,q6,prfd1n))