1! { dg-do compile }
2! { dg-options "-O2 -ffast-math -funroll-loops -ftree-vectorize -g" }
3  SUBROUTINE block_15_1_1_1(kbd,kbc,kad,kac,pbd,pbc,pad,pac,prim,scale)
4    INTEGER, PARAMETER :: dp=8
5    REAL(KIND=dp) :: kbd(1*1), kbc(1*1), kad(15*1), kac(15*1), pbd(1*1), &
6      pbc(1*1), pad(15*1), pac(15*1), prim(15*1*1*1), scale
7    INTEGER                                  :: ma, mb, mc, md, p_index
8      DO md = 1,1
9        DO mc = 1,1
10          DO mb = 1,1
11            DO ma = 1,15
12              p_index=p_index+1
13              tmp = scale*prim(p_index)
14              ks_bd = ks_bd + tmp* pac((mc-1)*15+ma)
15            END DO
16            kbd((md-1)*1+mb) = kbd((md-1)*1+mb) - ks_bd
17          END DO
18        END DO
19      END DO
20  END SUBROUTINE block_15_1_1_1
21