[Bug target/84499] New: [8 Regression] ICE on cp2k in rs6000_emit_le_vsx_store starting with r256656

["jakub at gcc dot gnu.org" <gcc-bugzilla at gcc dot gnu dot org>]

https://gcc.gnu.org/bugzilla/show_bug.cgi?id=84499

            Bug ID: 84499
           Summary: [8 Regression] ICE on cp2k in rs6000_emit_le_vsx_store
                    starting with r256656
           Product: gcc
           Version: 8.0
            Status: UNCONFIRMED
          Severity: normal
          Priority: P3
         Component: target
          Assignee: unassigned at gcc dot gnu.org
          Reporter: jakub at gcc dot gnu.org
  Target Milestone: ---

Following testcase ICEs starting with r256656 with:
-O1 -fstack-protector-strong -mcpu=power8
on powerpc64le-linux.

MODULE rh1546743
  PRIVATE
  PUBLIC :: foo
  INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(14, 200)
  INTEGER, PARAMETER :: int_8 = SELECTED_INT_KIND(10)
  COMPLEX(KIND=dp), PARAMETER :: gaussi = (0.0_dp, 1.0_dp)
  REAL(KIND=dp), PARAMETER :: twopi = 6.28318530717958647692528677_dp
  INTEGER, DIMENSION(:, :, :), ALLOCATABLE :: coset
  INTEGER, DIMENSION(:), ALLOCATABLE :: ncoset
CONTAINS
  SUBROUTINE foo(S_G, hmat, h_inv, vol, G_bounds_1, R_bounds_2, order)
    REAL(KIND=dp), DIMENSION(3, 3), INTENT(IN) :: hmat, h_inv
    REAL(KIND=dp), INTENT(IN) :: vol
    REAL(KIND=dp), DIMENSION(3, 3), INTENT(IN) :: G_bounds_1, R_bounds_2
    INTEGER :: l_max, sum_method, sum_order
    REAL(KIND=dp) :: alpha, beta, G_rad, gamma, R_rad
    REAL(KIND=dp), ALLOCATABLE, DIMENSION(:, :, :) :: S_G_tmp
    REAL(KIND=dp), DIMENSION(3) :: G_bound, R1, R2, R_bound
    SELECT CASE (sum_method)
    CASE (1)
      SELECT CASE (sum_order)
      CASE (1)
        CALL bar(S_G_tmp, l_max, m_max, n_max, R1, R2, alpha, beta, &
          gamma, hmat, h_inv, vol, R_bound, G_bound, R_rad, G_rad)
        CALL bar(S_G_tmp, n_max, m_max, l_max, -R2, -R1, gamma, beta, &
          alpha, hmat, h_inv, vol, R_bound, G_bound, R_rad, G_rad)
      END SELECT
    END SELECT
  END SUBROUTINE foo
  PURE SUBROUTINE bar(S_G, l_max, m_max, n_max, R1, R2, alpha, beta, gamma, &
    hmat, h_inv, vol, R_c, G_c, R_rad, G_rad)
    REAL(KIND=dp), DIMENSION(:, :, :), INTENT(OUT) :: S_G
    INTEGER, INTENT(IN) :: l_max, m_max, n_max
    REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: R1, R2
    REAL(KIND=dp), INTENT(IN) :: alpha, beta, gamma
    REAL(KIND=dp), DIMENSION(3, 3), INTENT(IN) :: hmat, h_inv
    REAL(KIND=dp), INTENT(IN) :: vol
    REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: R_c, G_c
    REAL(KIND=dp), INTENT(IN) :: R_rad, G_rad
    COMPLEX(KIND=dp), DIMENSION(ncoset(m_max), ncoset(n_max)) :: S_R
    COMPLEX(KIND=dp), DIMENSION(ncoset(l_max), ncoset(m_max), ncoset(n_max)) ::
S_G_c
    REAL(KIND=dp) :: G_rad_sq, R_rad_sq
    REAL(KIND=dp), DIMENSION(3) :: G, G_y, G_z
    REAL(KIND=dp), DIMENSION(3, 3) :: ht
    INTEGER :: gx, gy, gz, k
    DO gx = -FLOOR(G_c(1)), FLOOR(G_c(1))
      DO gy = -FLOOR(G_c(2)), FLOOR(G_c(2))
        G_y = ht(:, 2)*gy
        DO gz = -FLOOR(G_c(3)), FLOOR(G_c(3))
          DO k = 1, 3
            DO l = 1, l_max
            ENDDO
          ENDDO
          CALL baz(S_R, m_max, n_max, G, R2, beta, gamma, hmat, h_inv, vol,
R_c, R_rad_sq)
          DO lco = 1, ncoset(l_max)
          ENDDO
          DO nco = 1, ncoset(n_max)
            DO mco = 1, ncoset(m_max)
              DO lco = 1, ncoset(l_max)
                S_G_c(lco, mco, nco) =
exp_tot*i_pow(l+m+n)*(-1.0_dp)**m*S_R(mco, nco)
              ENDDO
            ENDDO
          ENDDO
        ENDDO
      ENDDO
    ENDDO
  END SUBROUTINE bar
  PURE SUBROUTINE baz(S_R, l_max, m_max, G, R, alpha, beta, hmat, h_inv, vol,
R_c, R_rad_sq)
    COMPLEX(KIND=dp), DIMENSION(:, :), INTENT(OUT) :: S_R
    INTEGER, INTENT(IN) :: l_max, m_max
    REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: G, R
    REAL(KIND=dp), INTENT(IN) :: alpha, beta
    REAL(KIND=dp), DIMENSION(3, 3), INTENT(IN) :: hmat, h_inv
    REAL(KIND=dp), INTENT(IN) :: vol
    REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: R_c
    REAL(KIND=dp), INTENT(IN) :: R_rad_sq
    INTEGER :: k, t, tz, sx, tx, ty
    COMPLEX(KIND=dp), DIMENSION(ncoset(l_max+m_max)) :: S_R_t
    REAL(KIND=dp), DIMENSION(-1:l_max+m_max+1, -1:l_max, -1:m_max) :: E1, E2,
E3
    REAL(KIND=dp), DIMENSION(3) :: R_l, R_r, Rp, Rx, Ry, Rz, s_shift
    REAL(KIND=dp), DIMENSION(3, 0:l_max+m_max) :: R_pow_t
    INTERFACE
      PURE SUBROUTINE boo(l, m, a, b, R1, R2, H_or_C_product, E)
         INTEGER, PARAMETER :: dp = SELECTED_REAL_KIND(14, 200)
         INTENT(IN) :: l, m
         REAL(KIND=dp), INTENT(IN) :: a, b, R1, R2
         INTEGER, INTENT(IN) :: H_or_C_product
         REAL(KIND=dp), DIMENSION(-1:l+m+1, -1:m), INTENT(OUT) :: E
      END SUBROUTINE
    END INTERFACE
    s_shift = MATMUL(h_inv, R)
    R_l = -R_c+s_shift
    DO sx = CEILING(R_l(1)), FLOOR(R_r(1))
      Rx = hmat(:, 1)*sx
      DO sy = CEILING(R_l(2)), FLOOR(R_r(2))
        Ry = hmat(:, 2)*sy
        DO sz = CEILING(R_l(3)), FLOOR(R_r(3))
          Rz = hmat(:, 3)*sz
          exp_tot = EXP(-c1*Rp_sq)*EXP(-gaussi*c2*DOT_PRODUCT(Rp, G))
          DO k = 1, 3
            R_pow_t(k, 0) = 1.0_dp
            DO t = 1, l_max+m_max
            ENDDO
          ENDDO
          DO tco = 1, ncoset(l_max+m_max)
            S_R_t(tco) = S_R_t(tco)+(-1.0_dp)**t*i_pow(t)*exp_tot
          ENDDO
        ENDDO
      ENDDO
    ENDDO
    CALL boo(l_max, m_max, alpha, beta, G(3), 0.0_dp, 1, E3)
    DO mco = 1, ncoset(m_max)
      DO lco = 1, ncoset(l_max)
        DO tx = 0, lx+mx
        DO ty = 0, ly+my
        DO tz = 0, lz+mz
          S_R(lco, mco) = S_R(lco, mco)+E1(tx, lx, mx)*E2(ty, ly, my)*E3(tz,
lz, mz)*S_R_t(tco)
        ENDDO
        ENDDO
        ENDDO
      ENDDO
    ENDDO
  END SUBROUTINE baz
  PURE ELEMENTAL FUNCTION i_pow(i)
    INTEGER, INTENT(IN) :: i
  END FUNCTION
END MODULE rh1546743