Sources Esope | Cast3M

prin35
C PRIN35    SOURCE    CHAT      05/01/13    02:27:21     5004
c**************************************************************************
      SUBROUTINE PRIN35(B,BPR,QEN,Q6,Q6T)
c**************************************************************************
      IMPLICIT INTEGER(I-N)
      real*8 b(4),bpr(3),qen(3,3),q6(6,6),q6t(6,6)
      integer i,j,k,l,nrot
      qen(1,1) = b(1)
      qen(2,2) = b(2)
      qen(3,3) = b(3)
      qen(2,1) = b(4)
      qen(1,2) = b(4)
      qen(3,1) = 0.d0
      qen(3,2) = 0.d0
      qen(1,3) = 0.d0
      qen(2,3) = 0.d0
      call eig35(qen,bpr,nrot)
      do i=1,3
        do j=1,3
          k=1+mod(i,3)
          l=1+mod(j,3)
          q6(i,j)     = qen(i,j)**2
          q6t(i,j)    = qen(j,i)**2
          q6(i+3,j)   = qen(i,j)*qen(k,j)
          q6t(i+3,j)  = qen(j,i)*qen(j,k)
          q6(j,i+3)   = qen(j,i)*qen(j,k)*2.D0
          q6t(j,i+3)  = qen(i,j)*qen(k,j)
          q6(j+3,i+3) = qen(j,i)*qen(l,k)+qen(j,k)*qen(l,i)
          q6t(j+3,i+3)=(qen(i,j)*qen(k,l)+qen(k,j)*qen(i,l))/2d0
        end do
      end do
      return
      end