Sources Esope | Cast3M

C MAMU21    SOURCE    PV        20/01/14    21:15:07     10498          
      subroutine mamu21(long,val11     ,val1,lpl,pt)
      implicit real*8 (a-h,o-z)
      real*8 val11(*),val1(*)
      real*8 pt(6,6)
*          id12=lpl1
*          id13=lpl1+lpl1+1
           id2=lpl
*          id3=lpl+lpl+1
           pt1=0.d0
           pt4=0.d0
           do 110 i=1,long
              xval11=val11(i)
*             xval12=val11(i+id12)
*              xval13=val11(i+id13)
              xval1=val1(i)
              xval2=val1(i+id2)
*             xval3=val1(i+id3)
              pt1=pt1+xval1*xval11
*             pt(2)=pt(2)+xval1*xval12
*             pt(3)=pt(3)+xval1*xval13
              pt4=pt4+xval2*xval11
*             pt(5)=pt(5)+xval2*xval12
*              pt(6)=pt(6)+xval2*xval13
*             pt(7)=pt(7)+xval3*xval11
*             pt(8)=pt(8)+xval3*xval12
*              pt(9)=pt(9)+xval3*xval13
 110       continue
           pt(1,1)=pt(1,1)+pt1
           pt(1,2)=pt(1,2)+pt4
 
 
 
           end