Télécharger f_s2pi.procedur

Retour à la liste

Numérotation des lignes :

  1. * F_S2PI PROCEDUR CHAT 93/01/11 21:18:40 780
  2. 'DEBPROC' F_S2PI CHARM*'MAILLAGE' SCIRC*'CHPOINT' NN*'ENTIER'
  3. OO*'POINT' ;
  4. *************************************************************************
  5. * magnetostatique 2d analyse harmonique procedure appelee par *
  6. * DDFOUR
  7. *************************************************************************
  8. VALIN = 'REDU' SCIRC ( CHARM 'POINT' 'INITIAL' ) ;
  9. VALFI = 'REDU' SCIRC ( CHARM 'POINT' 'FINAL' ) ;
  10. V1 = 'MAXI' VALIN ;V2 = 'MAXI' VALFI ;
  11. K = -1 ;
  12. 'SI' (( V1 < V2 ) 'ET' ( NUM 'NEG' 1 )) ;
  13. CHARM = 'INVE' CHARM ;K = 1 ;
  14. 'FINSI' ;
  15. 1P = CHARM 'POINT' 'INITIAL' ;2P = CHARM 'POINT' 'FINAL' ;
  16. IP = 0 ;
  17. 'REPETER' BOUC (NN - 1 );
  18. IP = IP + 1 ;
  19. 'SI' ( 'EGA' IP 1 ) ;
  20. SOL2 = (SCIRC 'PLUS' (0. 0.)) * K ;
  21. 'SINON' ;
  22. SOL2 = (SOL2 'PLUS' (0. 0.)) * K ;
  23. 'FINS' ;
  24. MAI2 = 'EXTR' SOL2 MAILLAGE ;
  25. PDOU = CHARM 'POINT' 'FINAL' ;
  26. VALDOU = 'REDU' SCIRC PDOU ;
  27.  
  28. MAI2 = 'DEPLACER' MAI2 'SYME' 'DROITE' OO PDOU ;
  29. 'SI' ( 'EGA' IP 1 ) ;
  30. CHARM2= ( 'INVE' (CHARM 'SYME' 'DROITE' OO PDOU )) ;
  31. 'SINON' ;
  32. CHARM2= ( 'INVE' (CHARM2 'SYME' 'DROITE' OO PDOU )) ;
  33. 'FINSI' ;
  34. 'ELIM' .0001 CHARM2 MAI2 ;
  35. 'CONFONDRE' ( CHARM 'POINT' 'FINAL') ( CHARM2 'POINT' 'INITIAL' ) ;
  36. CHARM = CHARM 'ET' CHARM2 ;
  37. SCIRC = SCIRC + SOL2 - VALDOU ;
  38. K = K * -1 ;
  39. 'FIN' BOUC ;
  40. *
  41. SUPO = CHARM ;
  42. * TITRE ' SOLUTION SUR 2PI ';
  43. * EVV1 = EVOL ROUG CHPO SCIRC 'T' CHARM ;
  44. * titre ' evo '( dime evv1 ) 'supp ' ( supo nbno ) ;
  45. * trac supo ;
  46. * dess evv1 ;
  47. 'FINPROC' SCIRC SUPO ;
  48.  
  49.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales