Numérotation des lignes :

C DEFELA    SOURCE    AM        09/12/08    21:15:10     6582      SUBROUTINE DEFELA(DEP,EPSILO,NSTRS,SIGF,IFOUR,KERRE)C_______________________________________________________________________CC CALCUL DE LA DEFORMATION ELASTIQUE EN FIN DE PASC_______________________________________________________________________C                                                                            C      IMPLICIT INTEGER(I-N)      IMPLICIT REAL*8(A-H,O-Z)      DIMENSION DEP(NSTRS,NSTRS)      DIMENSION EPSILO(NSTRS)      DIMENSION SIGF(NSTRS)      DIMENSION TEMP1(3,3),TEMP2(10,10)CC MATRICE DEP A INVERSERC      NDI =10      CALL ZERO(TEMP2,NDI,NSTRS)CC      IF (IFOUR.EQ.-2) THEN        CALL ZERO(TEMP1,3,3)        DO 10 I=1,3          DO 10 J=1,3            TEMP1(I,J)=DEP(I,J)   10   CONTINUEC        CALL INVMA1(TEMP1,3,3,ISING)        IF (ISING.EQ.1) THEN          WRITE(*,*) 'Matrice singuliere dans defela'          KERRE = 99          RETURN        ENDIFC        TEMP2(1,1)=TEMP1(1,1)        TEMP2(1,2)=TEMP1(1,2)        TEMP2(2,1)=TEMP1(2,1)        TEMP2(2,2)=TEMP1(2,2)        TEMP2(3,3)=TEMP1(3,3)      ELSE        DO 30 I=1,NSTRS          DO 30 J=1,NSTRS            TEMP2(I,J)=DEP(I,J)   30 CONTINUE        CALL INVMA1(TEMP2,NDI,NSTRS,ISING)        IF (ISING.EQ.1) THEN          WRITE(*,*) 'Matrice singuliere dans defela'          KERRE = 99          RETURN        ENDIF      ENDIFCC CALCUL DU VECTEUR DEFORMATION ELASTIQUE EN FIN DE PASC      DO 50 I=1,NSTRS        EPSILO(I)=0.D0        DO 50 J=1,NSTRS          EPSILO(I)=EPSILO(I)+TEMP2(I,J)*SIGF(J)   50 CONTINUEC      END

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