Télécharger triang.eso

Retour à la liste

Numérotation des lignes :

  1. C TRIANG SOURCE PV 15/11/25 21:15:16 8707
  2. SUBROUTINE TRIANG(KRIGI,PREC,istab)
  3. C=======================================================================
  4. C ASSEMBLE LES PETITES MATRICES RIGIDITE ET LES MET SOUS LA FORME
  5. C Lt.D.L.
  6. C IL LE POINTEUR DE LA MATRICE RESULTANTE DANS ICHOLE
  7. C=======================================================================
  8. C
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8 (A-H,O-Z)
  11. -INC SMRIGID
  12. -INC PPARAM
  13. -INC CCOPTIO
  14. IF(IIMPI.EQ.1)THEN
  15. CALL GIBTEM(XKT)
  16. INTERR(1)=XKT
  17. CALL ERREUR(-259)
  18. ENDIF
  19. IF(IIMPI.EQ.1)WRITE(IOIMP,10)
  20. 10 FORMAT(' L''IMPRESSION PRECDENTE EST AVANT ASSEM1 ')
  21. CALL ASSEM1(KRIGI,MMATRX,INUINX,ITOPOX,IMINIX,IPOX,INCTRX,IITOPX)
  22. IF(IERR.NE.0) RETURN
  23. IF(IIMPI.EQ.1)THEN
  24. CALL GIBTEM(XKT)
  25. INTERR(1)=XKT
  26. CALL ERREUR(-259)
  27. ENDIF
  28. IF(IIMPI.EQ.1)WRITE(IOIMP,11)
  29. 11 FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT ASSEM2')
  30. CALL ASSEM2(KRIGI,ITOPOX,INUINX,IMINIX,MMATRX,IPOX,INCTRX,IITOPX)
  31. IF(IERR.NE.0) RETURN
  32. IF(IIMPI.EQ.1)THEN
  33. CALL GIBTEM(XKT)
  34. INTERR(1)=XKT
  35. CALL ERREUR(-259)
  36. ENDIF
  37. IF(IIMPI.EQ.1)WRITE(IOIMP,12)
  38. 12 FORMAT(' L''IMPRESSION PRECEDENTE EST AVANT CHOLE ')
  39. IF(IERR.NE.0) GO TO 5000
  40. CALL CHOLE(MMATRX,PREC,istab)
  41. IF(IERR.NE.0) RETURN
  42. IF(IIMPI.EQ.1)THEN
  43. CALL GIBTEM(XKT)
  44. INTERR(1)=XKT
  45. CALL ERREUR(-259)
  46. ENDIF
  47. IF(IERR.NE.0) GO TO 5000
  48. IF(IIMPI.EQ.1) WRITE(IOIMP,13)
  49. 13 FORMAT(' L''IMPRESSION PRECEDENTE EST APRES CHOLE')
  50. MRIGID=KRIGI
  51. SEGACT MRIGID*MOD
  52. ICHOLE=MMATRX
  53. SEGDES MRIGID
  54. 5000 CONTINUE
  55. RETURN
  56. END
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.  
  64.  
  65.  
  66.  

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