1! { dg-do compile }
2! { dg-options "-O2 -ftree-vectorize" }
3      subroutine gen3delem(nel,ial,ifix,xta,xnoref,dd,jact,nelshell)
4      real*8 xnoref(3),xta(3,100),xn1(3,100)
5      if(nel.gt.0) then
6        do j=1,nel
7        enddo
8        do
9        enddo
10      endif
11      do
12        if(ifix.eq.0) then
13          do j=nelshell,nel
14            if(ial(j).eq.0) then
15            endif
16          enddo
17        endif
18        do j=nelshell,nel
19        enddo
20        do j=1,3
21          xnoref(j)=xnoref(j)/dd
22        enddo
23        xn1(2,jact)=xnoref(3)*xta(1,jact)-xnoref(1)*xta(3,jact)
24        xn1(3,jact)=xnoref(1)*xta(2,jact)-xnoref(2)*xta(1,jact)
25        call foo(xn1)
26      enddo
27      end
28