+2015-11-08 Joost VandeVondele <vondele@gcc.gnu.org>
+
+ * gfortran.dg/PR67518.f90: move from here...
+ * gfortran.dg/graphite/PR67518.f90: to here.
+ * gfortran.dg/PR53852.f90: move from here...
+ * gfortran.dg/graphite/PR53852.f90: to here.
+
2015-11-08 Eric Botcazou <ebotcazou@adacore.com>
* gcc.dg/sso-1.c: New test.
+++ /dev/null
-! { dg-do compile }
-! { dg-options "-floop-nest-optimize -O2 -ffast-math" }
-! PR53852 : compile time / memory hog
-SUBROUTINE build_d_tensor_gks(d5f,v,d5)
- INTEGER, PARAMETER :: dp=8
- REAL(KIND=dp), DIMENSION(3, 3, 3, 3, 3), &
- INTENT(OUT) :: d5f
- REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: v
- REAL(KIND=dp), INTENT(IN) :: d5
- INTEGER :: k1, k2, k3, k4, k5
- REAL(KIND=dp) :: w
-
- d5f = 0.0_dp
- DO k1=1,3
- DO k2=1,3
- DO k3=1,3
- DO k4=1,3
- DO k5=1,3
- d5f(k5,k4,k3,k2,k1)=d5f(k5,k4,k3,k2,k1)+ &
- v(k1)*v(k2)*v(k3)*v(k4)*v(k5)*d5
- ENDDO
- w=v(k1)*v(k2)*v(k3)*d4
- d5f(k1,k2,k3,k4,k4)=d5f(k1,k2,k3,k4,k4)+w
- d5f(k1,k2,k4,k3,k4)=d5f(k1,k2,k4,k3,k4)+w
- d5f(k1,k4,k2,k3,k4)=d5f(k1,k4,k2,k3,k4)+w
- d5f(k4,k1,k2,k3,k4)=d5f(k4,k1,k2,k3,k4)+w
- d5f(k1,k2,k4,k4,k3)=d5f(k1,k2,k4,k4,k3)+w
- d5f(k1,k4,k2,k4,k3)=d5f(k1,k4,k2,k4,k3)+w
- d5f(k4,k1,k2,k4,k3)=d5f(k4,k1,k2,k4,k3)+w
- d5f(k1,k4,k4,k2,k3)=d5f(k1,k4,k4,k2,k3)+w
- d5f(k4,k1,k4,k2,k3)=d5f(k4,k1,k4,k2,k3)+w
- d5f(k4,k4,k1,k2,k3)=d5f(k4,k4,k1,k2,k3)+w
- ENDDO
- ENDDO
- ENDDO
- ENDDO
-END SUBROUTINE build_d_tensor_gks
+++ /dev/null
-! { dg-do compile }
-! { dg-options "-floop-nest-optimize -O2 -ffast-math" }
-! PR67518 : isl: position out of bounds
-MODULE ao_util
- INTEGER, PARAMETER :: dp=8
-CONTAINS
- FUNCTION exp_radius(l,alpha,threshold,prefactor,epsin) RESULT(radius)
- REAL(KIND=dp), INTENT(IN) :: alpha, threshold, prefactor
- REAL(KIND=dp), INTENT(IN), OPTIONAL :: epsin
- DO
- IF (iter.gt.maxiter) THEN
- CALL stop_program(routineN,moduleN,1,"exceeded")
- ENDIF
- ENDDO
- CALL stop_program(routineN,moduleN,1,"exceeded")
- END FUNCTION exp_radius
- FUNCTION exp_radius_very_extended(la_min,la_max,lb_min,lb_max,pab,o1,o2,ra,rb,rp,&
- zetp,eps,prefactor,cutoff,epsin) RESULT(radius)
- REAL(KIND=dp), DIMENSION(:, :), &
- OPTIONAL, POINTER :: pab
- REAL(KIND=dp), INTENT(IN) :: ra(3), rb(3), rp(3), zetp, &
- eps, prefactor, cutoff
- REAL(KIND=dp) :: bini, binj, coef(0:20), &
- epsin_local, polycoef(0:60), &
- rad_b, s1, s2
- IF (PRESENT(pab)) THEN
- ENDIF
- DO lxa=0,la_max
- DO lxb=0,lb_max
- coef(0:la_max+lb_max)=0.0_dp
- DO i=0,lxa
- DO j=0,lxb
- coef(lxa+lxb-i-j)=coef(lxa+lxb-i-j) + bini*binj*s1*s2
- ENDDO
- ENDDO
- DO i=0,lxa+lxb
- polycoef(i)=MAX(polycoef(i),coef(i))
- ENDDO
- ENDDO
- ENDDO
- DO i=0,la_max+lb_max
- radius=MAX(radius,exp_radius(i,zetp,eps,polycoef(i),epsin_local) )
- ENDDO
- END FUNCTION exp_radius_very_extended
-END MODULE ao_util
--- /dev/null
+! { dg-do compile }
+! { dg-options "-floop-nest-optimize -O2 -ffast-math" }
+! PR53852 : compile time / memory hog
+SUBROUTINE build_d_tensor_gks(d5f,v,d5)
+ INTEGER, PARAMETER :: dp=8
+ REAL(KIND=dp), DIMENSION(3, 3, 3, 3, 3), &
+ INTENT(OUT) :: d5f
+ REAL(KIND=dp), DIMENSION(3), INTENT(IN) :: v
+ REAL(KIND=dp), INTENT(IN) :: d5
+ INTEGER :: k1, k2, k3, k4, k5
+ REAL(KIND=dp) :: w
+
+ d5f = 0.0_dp
+ DO k1=1,3
+ DO k2=1,3
+ DO k3=1,3
+ DO k4=1,3
+ DO k5=1,3
+ d5f(k5,k4,k3,k2,k1)=d5f(k5,k4,k3,k2,k1)+ &
+ v(k1)*v(k2)*v(k3)*v(k4)*v(k5)*d5
+ ENDDO
+ w=v(k1)*v(k2)*v(k3)*d4
+ d5f(k1,k2,k3,k4,k4)=d5f(k1,k2,k3,k4,k4)+w
+ d5f(k1,k2,k4,k3,k4)=d5f(k1,k2,k4,k3,k4)+w
+ d5f(k1,k4,k2,k3,k4)=d5f(k1,k4,k2,k3,k4)+w
+ d5f(k4,k1,k2,k3,k4)=d5f(k4,k1,k2,k3,k4)+w
+ d5f(k1,k2,k4,k4,k3)=d5f(k1,k2,k4,k4,k3)+w
+ d5f(k1,k4,k2,k4,k3)=d5f(k1,k4,k2,k4,k3)+w
+ d5f(k4,k1,k2,k4,k3)=d5f(k4,k1,k2,k4,k3)+w
+ d5f(k1,k4,k4,k2,k3)=d5f(k1,k4,k4,k2,k3)+w
+ d5f(k4,k1,k4,k2,k3)=d5f(k4,k1,k4,k2,k3)+w
+ d5f(k4,k4,k1,k2,k3)=d5f(k4,k4,k1,k2,k3)+w
+ ENDDO
+ ENDDO
+ ENDDO
+ ENDDO
+END SUBROUTINE build_d_tensor_gks
--- /dev/null
+! { dg-do compile }
+! { dg-options "-floop-nest-optimize -O2 -ffast-math" }
+! PR67518 : isl: position out of bounds
+MODULE ao_util
+ INTEGER, PARAMETER :: dp=8
+CONTAINS
+ FUNCTION exp_radius(l,alpha,threshold,prefactor,epsin) RESULT(radius)
+ REAL(KIND=dp), INTENT(IN) :: alpha, threshold, prefactor
+ REAL(KIND=dp), INTENT(IN), OPTIONAL :: epsin
+ DO
+ IF (iter.gt.maxiter) THEN
+ CALL stop_program(routineN,moduleN,1,"exceeded")
+ ENDIF
+ ENDDO
+ CALL stop_program(routineN,moduleN,1,"exceeded")
+ END FUNCTION exp_radius
+ FUNCTION exp_radius_very_extended(la_min,la_max,lb_min,lb_max,pab,o1,o2,ra,rb,rp,&
+ zetp,eps,prefactor,cutoff,epsin) RESULT(radius)
+ REAL(KIND=dp), DIMENSION(:, :), &
+ OPTIONAL, POINTER :: pab
+ REAL(KIND=dp), INTENT(IN) :: ra(3), rb(3), rp(3), zetp, &
+ eps, prefactor, cutoff
+ REAL(KIND=dp) :: bini, binj, coef(0:20), &
+ epsin_local, polycoef(0:60), &
+ rad_b, s1, s2
+ IF (PRESENT(pab)) THEN
+ ENDIF
+ DO lxa=0,la_max
+ DO lxb=0,lb_max
+ coef(0:la_max+lb_max)=0.0_dp
+ DO i=0,lxa
+ DO j=0,lxb
+ coef(lxa+lxb-i-j)=coef(lxa+lxb-i-j) + bini*binj*s1*s2
+ ENDDO
+ ENDDO
+ DO i=0,lxa+lxb
+ polycoef(i)=MAX(polycoef(i),coef(i))
+ ENDDO
+ ENDDO
+ ENDDO
+ DO i=0,la_max+lb_max
+ radius=MAX(radius,exp_radius(i,zetp,eps,polycoef(i),epsin_local) )
+ ENDDO
+ END FUNCTION exp_radius_very_extended
+END MODULE ao_util