Télécharger la plaquette Cast3M


fr en

Télécharger invmat.eso

Retour à la liste

Numérotation des lignes : 

invmat

  1. C INVMAT    SOURCE    PV        20/09/26    21:17:28     10724          
  2.       SUBROUTINE INVMAT(MAT,MATRIK,IRET)
  3.       IMPLICIT INTEGER(I-N)
  4.       IMPLICIT REAL*8 (A-H,O-Z)
  5. c     **************************************
  6. C     * ROUTINE d inversion de MATRIK      *
  7. C     * Diagonale                          *
  8. C     * MAT : MATRIK entree                *
  9. C     * MATRIK : MATRIK :sortie            *
  10. C     * MATRIK est la matrice inverse de   *
  11. C     * MAT sous forme matrice elementaire *
  12. C     * est sous forme morse               *
  13. C     **************************************
  14.  
  16.       POINTEUR MAT.MATRIK
  17.  
  18.       SEGACT MAT
  19.       IF (MAT.IRIGEL(7,1).NE.5) THEN
  20.         WRITE(6,*) 'Impossible d inverser une'
  21.         WRITE(6,*) 'matrice non diagonale'
  22.         IRET=1
  23.         RETURN
  24.       END IF
  25.       SEGDES MAT
  26.  
  27.       CALL MCOPY(MAT,MATRIK)
  28.       SEGACT MATRIK*MOD
  29.  
  30.       IMATRI=IRIGEL(4,1)
  31.       SEGACT IMATRI
  32.  
  33.       NBSOUS=LIZAFM(/1)
  34.       NBME=LIZAFM(/2)
  35.  
  36.       DO I=1,NBSOUS
  37.         DO J=1,NBME
  38.            IZAFM=LIZAFM(I,J)
  39.            SEGACT IZAFM*MOD
  40.            NBEL=AM(/1)
  41.            DO K=1,NBEL
  42.              AM(K,1,1)=1/AM(K,1,1)
  43.            END DO
  44.            SEGDES IZAFM
  45.         END DO
  46.       END DO
  47.  
  48.       SEGDES IMATRI
  49.       SEGDES MATRIK
  50.  
  51.       END
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  

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