! { dg-do compile }
|
! { dg-do compile }
|
! { dg-options "-O3 -ffast-math -fdump-tree-reassoc1" }
|
! { dg-options "-O3 -ffast-math -fdump-tree-reassoc1" }
|
! { dg-options "-O3 -ffast-math -fdump-tree-reassoc1 --param max-completely-peel-times=16" { target spu-*-* } }
|
! { dg-options "-O3 -ffast-math -fdump-tree-reassoc1 --param max-completely-peel-times=16" { target spu-*-* } }
|
subroutine anisonl(w,vo,anisox,s,ii1,jj1,weight)
|
subroutine anisonl(w,vo,anisox,s,ii1,jj1,weight)
|
integer ii1,jj1,i1,iii1,j1,jjj1,k1,l1,m1,n1
|
integer ii1,jj1,i1,iii1,j1,jjj1,k1,l1,m1,n1
|
real*8 w(3,3),vo(3,3),anisox(3,3,3,3),s(60,60),weight
|
real*8 w(3,3),vo(3,3),anisox(3,3,3,3),s(60,60),weight
|
!
|
!
|
! This routine replaces the following lines in e_c3d.f for
|
! This routine replaces the following lines in e_c3d.f for
|
! an anisotropic material
|
! an anisotropic material
|
!
|
!
|
do i1=1,3
|
do i1=1,3
|
iii1=ii1+i1-1
|
iii1=ii1+i1-1
|
do j1=1,3
|
do j1=1,3
|
jjj1=jj1+j1-1
|
jjj1=jj1+j1-1
|
do k1=1,3
|
do k1=1,3
|
do l1=1,3
|
do l1=1,3
|
s(iii1,jjj1)=s(iii1,jjj1)
|
s(iii1,jjj1)=s(iii1,jjj1)
|
& +anisox(i1,k1,j1,l1)*w(k1,l1)*weight
|
& +anisox(i1,k1,j1,l1)*w(k1,l1)*weight
|
do m1=1,3
|
do m1=1,3
|
s(iii1,jjj1)=s(iii1,jjj1)
|
s(iii1,jjj1)=s(iii1,jjj1)
|
& +anisox(i1,k1,m1,l1)*w(k1,l1)
|
& +anisox(i1,k1,m1,l1)*w(k1,l1)
|
& *vo(j1,m1)*weight
|
& *vo(j1,m1)*weight
|
& +anisox(m1,k1,j1,l1)*w(k1,l1)
|
& +anisox(m1,k1,j1,l1)*w(k1,l1)
|
& *vo(i1,m1)*weight
|
& *vo(i1,m1)*weight
|
do n1=1,3
|
do n1=1,3
|
s(iii1,jjj1)=s(iii1,jjj1)
|
s(iii1,jjj1)=s(iii1,jjj1)
|
& +anisox(m1,k1,n1,l1)
|
& +anisox(m1,k1,n1,l1)
|
& *w(k1,l1)*vo(i1,m1)*vo(j1,n1)*weight
|
& *w(k1,l1)*vo(i1,m1)*vo(j1,n1)*weight
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
enddo
|
|
|
return
|
return
|
end
|
end
|
|
|
! There should be 22 multiplications left after un-distributing
|
! There should be 22 multiplications left after un-distributing
|
! weigth, w(k1,l1), vo(i1,m1) and vo(j1,m1) on the innermost two
|
! weigth, w(k1,l1), vo(i1,m1) and vo(j1,m1) on the innermost two
|
! unrolled loops.
|
! unrolled loops.
|
|
|
! { dg-final { scan-tree-dump-times "\[0-9\] \\\* " 22 "reassoc1" } }
|
! { dg-final { scan-tree-dump-times "\[0-9\] \\\* " 22 "reassoc1" } }
|
! { dg-final { cleanup-tree-dump "reassoc1" } }
|
! { dg-final { cleanup-tree-dump "reassoc1" } }
|
|
|