From: "VandeVondele Joost" <joost.vandevondele@mat.ethz.ch>
To: "fortran@gcc.gnu.org" <fortran@gcc.gnu.org>,
"gcc-patches@gcc.gnu.org" <gcc-patches@gcc.gnu.org>
Subject: [PATCH] PR67518 and PR53852 -- add testcase.
Date: Tue, 03 Nov 2015 18:22:00 -0000 [thread overview]
Message-ID: <908103EDB4893A42920B21D3568BFD9339EFE356@MBX23.d.ethz.ch> (raw)
[-- Attachment #1: Type: text/plain, Size: 375 bytes --]
Attached testcases for two previously fixed PRs (and thanks to Dominique who was quicker for PR67982).
2015-11-03 Joost VandeVondele <vondele@gnu.gcc.org>
PR middle-end/53852
PR middle-end/67518
* gfortran.dg/PR67518.f90: New test.
* gfortran.dg/PR53852.f90: New test.
OK for trunk after finished bootstrap and testing ?
Joost
[-- Attachment #2: patch.prs --]
[-- Type: application/octet-stream, Size: 3649 bytes --]
Index: gcc/testsuite/gfortran.dg/PR67518.f90
===================================================================
--- gcc/testsuite/gfortran.dg/PR67518.f90 (revision 0)
+++ gcc/testsuite/gfortran.dg/PR67518.f90 (revision 0)
@@ -0,0 +1,45 @@
+! { 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
Index: gcc/testsuite/gfortran.dg/PR53852.f90
===================================================================
--- gcc/testsuite/gfortran.dg/PR53852.f90 (revision 0)
+++ gcc/testsuite/gfortran.dg/PR53852.f90 (revision 0)
@@ -0,0 +1,37 @@
+! { 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
next reply other threads:[~2015-11-03 18:22 UTC|newest]
Thread overview: 8+ messages / expand[flat|nested] mbox.gz Atom feed top
2015-11-03 18:22 VandeVondele Joost [this message]
2015-11-06 7:04 ` VandeVondele Joost
2015-11-06 9:25 ` Paul Richard Thomas
2015-11-06 9:45 ` VandeVondele Joost
2015-11-08 17:02 ` Paul Richard Thomas
2015-11-08 17:18 ` Andre Vehreschild
2015-11-08 19:05 ` VandeVondele Joost
2015-11-08 19:14 ` VandeVondele Joost
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=908103EDB4893A42920B21D3568BFD9339EFE356@MBX23.d.ethz.ch \
--to=joost.vandevondele@mat.ethz.ch \
--cc=fortran@gcc.gnu.org \
--cc=gcc-patches@gcc.gnu.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).