1! { dg-do compile { target powerpc*-*-* } }
2! { dg-additional-options "-O3 -mcpu=power7 -ffast-math -mveclibabi=mass" }
3
4      integer index(18),i,j,k,l,ipiv(18),info,ichange,neq,lda,ldb,
5     &  nrhs,iplas
6      real*8 ep0(6),al10(18),al20(18),dg0(18),ep(6),al1(18),
7     &  al2(18),dg(18),ddg(18),xm(6,18),h(18,18),ck(18),cn(18),
8     &  c(18),d(18),phi(18),delta(18),r0(18),q(18),b(18),cphi(18),
9     &  q1(18),q2(18),stri(6),htri(18),sg(18),r(42),xmc(6,18),aux(18),
10     &  t(42),gl(18,18),gr(18,18),ee(6),c1111,c1122,c1212,dd,
11     &  skl(3,3),xmtran(3,3),ddsdde(6,6),xx(6,18)
12      do
13         do i=1,18
14            htri(i)=dabs(sg(i))-r0(i)-ck(i)*(dg(i)/dtime)**(1.d0/cn(i))
15            do j=1,18
16            enddo
17         enddo
18         do
19                  if(i.ne.j) then
20                  gr(index(i),1)=htri(i)
21               endif
22            call dgesv(neq,nrhs,gl,lda,ipiv,gr,ldb,info)
23         enddo
24      enddo
25      end
26
27