Télécharger crech2.eso

Retour à la liste

Numérotation des lignes :

  1. C CRECH2 SOURCE CHAT 09/10/12 21:15:01 6526
  2. SUBROUTINE CRECH2(KCHPOI,mvectx,mvecri,icle)
  3. *
  4. * cree uun cpoint à partir d'un vecteur issu d'une resolution
  5. * a partir des info stockees dans mvecri (qui lui meme est rempli par
  6. * crech3
  7. *
  8. * si ICLE=1 on fait un vecteur primal ; si =2 un dual
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8(A-H,O-Z)
  11. -INC CCOPTIO
  12. -INC SMCHPOI
  13. -INC SMVECTD
  14. SEGACT,mvecri
  15. NSOUPO= MELZON(/1)
  16. NAT=1
  17. SEGINI,MCHPOI
  18. IFOPOI=IFOMOD
  19. MTYPOI=' '
  20. MOCHDE=' '
  21. C
  22. C **** ON VA FABRIQUER LES SEGMENTS MSOUPO POUR LES REMPLIR IL FAUT
  23. C **** CONNAITRE LES INCONNUES DU SOUS CHAMPS ET L'OBJET GEOMETRIQUE
  24. C **** SUPPORT
  25. C
  26. DO 100 I=1,NSOUPO
  27. C
  28. C **** ON CHERCHE D'ABORD LA LISTE DES INCONNUES A PARTIR DE LA
  29. C **** VALEUR DE IBINN ET ON REMPLIT NOCOMP
  30. C
  31. nc= NUMNOM(i)
  32. n=numnoe(i)
  33. SEGINI,MSOUPO
  34. IPCHP(I)=MSOUPO
  35. IB=0
  36. DO 14 J=1,NC
  37. IF(ICLE.EQ.1)NOCOMP(J)= NAMINC(I,j)
  38. IF(ICLE.EQ.2)NOCOMP(J)= NAMDUA(I,j)
  39. NOHARM(J)=NUMHAR(I,J)
  40. 14 CONTINUE
  41. C
  42. C **** ON CHERCHE COMBIEN DE NOEUD DANS L'OBJET MELEME,ON LE CREE
  43. C
  44.  
  45.  
  46. SEGINI,MPOVAL
  47. IGEOC=melzon(I)
  48. ipoval=mpoval
  49. 100 CONTINUE
  50. * remplissage des vpocha
  51. mvectd=mvectx
  52. segact mvectd
  53. do ia=1,vectbb(/1)
  54. izo= NUMZON(ia)
  55. msoupo=ipchp(izo)
  56. mpoval=ipoval
  57. i=NUNOLO(ia)
  58. j=NUINLO(ia)
  59. vpocha(i,j)=vectbb(ia)
  60. enddo
  61. * desactivation du chpoint
  62. do ia= 1,ipchp(/1)
  63. msoupo=ipchp(ia)
  64. mpoval=ipoval
  65. segdes mpoval
  66. segdes msoupo
  67. enddo
  68.  
  69.  
  70. SEGDES,mvecri
  71. KCHPOI=MCHPOI
  72. SEGDES,MCHPOI
  73. RETURN
  74. END
  75.  
  76.  
  77.  

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