Télécharger vch1.eso

Retour à la liste

Numérotation des lignes :

  1. C VCH1 SOURCE PV 16/11/17 22:01:40 9180
  2. SUBROUTINE VCH1 (MMATRX,MVECTX,ISOLU,KRIGI)
  3. C
  4. C **** A PARTIR D UN OBJET DE TYPE MATRICE TRANSFORME UN VECTEUR EN
  5. C **** CHPOIN DE TYPE PREMIER MEMBRE
  6. C **** KRIGI EST LE POINTEUR DE LA RIGIDITE.SERT A SAVOIRE SI LA GEOM
  7. C **** EXISTE DEJA.
  8. C
  9. IMPLICIT INTEGER(I-N)
  10. -INC SMMATRI
  11. -INC SMCHPOI
  12. -INC SMELEME
  13. -INC CCOPTIO
  14. -INC SMVECTD
  15. -INC SMRIGID
  16. -INC TMTRAV
  17. C
  18. IF(IIMPI.EQ.3) WRITE(IOIMP,3000) MMATRX,MVECTX
  19. 3000 FORMAT(' OPERATEUR VCH1 : POINTEUR DE LA MATRICE=',I5,
  20. 1 ' POINTEUR DU VECTEUR=',I5)
  21. C
  22. C **** ACTIVATION DES SEGMENTS
  23. C
  24. MRIGID=KRIGI
  25. SEGACT,MRIGID
  26. IMGEOQ=IMGEO1
  27.  
  28. MMATRI=MMATRX
  29. SEGACT,MMATRI
  30. MILIGN=IILIGN
  31. SEGACT,MILIGN
  32. INC=IPNO(/1)
  33. SEGDES,MILIGN
  34. MINCPO=IINCPO
  35. SEGACT,MINCPO
  36. MIMIK=IIMIK
  37. SEGACT,MIMIK
  38. MHARK=IHARK
  39. SEGACT,MHARK
  40. NNIN=IMIK(/2)
  41. midua=iidua
  42. segact midua
  43. MVECTD=MVECTX
  44. SEGACT MVECTD
  45. C
  46. C **** CREATION D'UN SEGMENT DE TYPE MCHPOI(VOIR SMCHPOI)
  47. C
  48.  
  49. if( ivecri.eq.0) then
  50. NND=VECTBB(/1)
  51. MELEME=IGEOMA
  52. SEGACT,MELEME
  53. NNNOE=NUM(/2)
  54. SEGINI,MTRAV
  55. DO 40 ITYU=1,NNNOE
  56. IGEO(ITYU)=NUM(1,ITYU)
  57. 40 CONTINUE
  58. SEGDES,MELEME
  59. DO 41 ITYU=1,NNIN
  60. NHAR(ITYU)=IHAR(ITYU)
  61. 41 INCO(ITYU)=IMIK(ITYU)
  62. DO 23 I=1,NNNOE
  63. DO 23 J=1,NNIN
  64. IK=INCPO(J,I)
  65. IF(IK.EQ.0) GO TO 23
  66. IBIN(J,I)=ik
  67. BB(J,I)=VECTBB(IK)
  68. 23 CONTINUE
  69. segact mrigid*mod
  70. ichoa=ichole
  71. if(ichole.eq.0) then
  72. ichole = mmatri
  73. endif
  74. CALL CRECH3(MTRAV,ISOLU,nnd,mrigid)
  75. ichole=ichoa
  76. SEGSUP,MTRAV
  77. else
  78. lvecri=ivecri
  79. CALL CRECH2(ISOLU,mvectx,lvecri,1)
  80. endif
  81.  
  82. SEGDES,MINCPO
  83. SEGDES,MIMIK
  84. segdes midua
  85. SEGDES,MMATRI
  86. SEGDES,MHARK
  87. * SEGDES MVECTD
  88. SEGDES,MRIGID
  89. IF(IIMPI.EQ.3) WRITE(IOIMP,3001) ISOLU
  90. 3001 FORMAT(' SUBROUTINE VCH1 : POINTEUR DE L'' OBJET CHPOINT=',I5)
  91. RETURN
  92. END
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  

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