Sources Esope | Cast3M

matgho
C MATGHO    SOURCE    CHAT      05/01/13    01:34:02     5004
CCC
C Matriz de Hook ampliada con 1
CCC
      SUBROUTINE MATGENHOOK (AMAT,NDIMX,NDIMS)
      IMPLICIT INTEGER(I-N)
      integer ndimx,ndims,ifl,i,j
      real*8    Amat(ndimx,ndimx),xnu,you
      common /elasdata/ you,xnu
c Inversa de la matriz de hook
      ifl=1
      call mathok(Amat,ndimx,ndims,ifl)
c Parte de las variables internas
      do i=ndims+1,ndimx
         do j=1,ndimx
            Amat(i,j)=0.D0
            Amat(j,i)=0.D0
         enddo
         Amat(i,i)=1.D0
      enddo
      return
      end