Sources Esope | Cast3M

x33x6
C X33X6     SOURCE    PV090527  23/01/27    21:16:10     11574          
c*************************************************************************************
c     passage matrice 33 symetrique vecteur 6 
      subroutine x33x6(x33,x6)
      real*8 x6(6),x33(3,3)
      integer i,k,l
      do i=1,6
       call indce0(i,k,l)
       if(i.le.3)then
        x6(i)=x33(k,l)
       else
        x6(i)=0.5*(x33(k,l)+x33(l,k))
       end if
      end do
      return
      end