! { dg-options "-O2 -floop-interchange" }

      subroutine linel(icmdl,stre,anisox)
      real*8 stre(6),tkl(3,3),ekl(3,3),anisox(3,3,3,3)
            do m1=1,3
               do m2=1,m1
                  do m3=1,3
                     do m4=1,3
                        tkl(m1,m2)=tkl(m1,m2)+
     &                       anisox(m1,m2,m3,m4)*ekl(m3,m4)
                     enddo
                  enddo
               enddo
            enddo
            stre(1)=tkl(1,1)
      end