diff --git a/gcc/ChangeLog.graphite b/gcc/ChangeLog.graphite
index f0b08515a85..5a8c8fe97eb 100644
--- a/gcc/ChangeLog.graphite
+++ b/gcc/ChangeLog.graphite
@@ -1,3 +1,10 @@
+2010-02-07  Sebastian Pop  <sebastian.pop@amd.com>
+
+	* gfortran.dg/graphite/id-19.f: New.
+	* gfortran.dg/graphite/pr14741.f90: New.
+	* gfortran.dg/graphite/pr41924.f90: New.
+	* gfortran.dg/graphite/run-id-2.f90: New.
+
 2010-02-07  Sebastian Pop  <sebastian.pop@amd.com>
 
 	PR middle-end/42988
diff --git a/gcc/testsuite/gfortran.dg/graphite/id-19.f b/gcc/testsuite/gfortran.dg/graphite/id-19.f
new file mode 100644
index 00000000000..e05f764b294
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/graphite/id-19.f
@@ -0,0 +1,15 @@
+      SUBROUTINE ECCODR(FPQR)
+      DIMENSION FPQR(25,25,25)
+      INTEGER P,Q,R
+      DIMENSION REC(73)
+      DO 150 P=1,N4MAX,2
+         QM2=-ONE
+         DO 140 Q=1,N4MAX,2
+            DO 130 R=1,N4MAX,2
+               IF(P.GT.1) THEN
+                  FPQR(P,Q,R)= QM2*FPQR(P,Q-2,R)*REC(P+Q-2+R)
+               END IF
+  130       RM2= RM2+TWO
+  140    QM2= QM2+TWO
+  150 PM2= PM2+TWO
+      END
diff --git a/gcc/testsuite/gfortran.dg/graphite/pr14741.f90 b/gcc/testsuite/gfortran.dg/graphite/pr14741.f90
new file mode 100644
index 00000000000..3fe1d690c38
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/graphite/pr14741.f90
@@ -0,0 +1,29 @@
+! { dg-options "-O3 -ffast-math -floop-interchange -floop-block -fdump-tree-graphite-all" }
+
+  INTEGER, PARAMETER :: N=1024
+  REAL*8 :: A(N,N), B(N,N), C(N,N)
+  REAL*8 :: t1,t2
+  A=0.1D0
+  B=0.1D0
+  C=0.0D0
+  CALL cpu_time(t1)
+  CALL mult(A,B,C,N)
+  CALL cpu_time(t2)
+  write(6,*) t2-t1,C(1,1)
+END program
+
+SUBROUTINE mult(A,B,C,N)
+  REAL*8 :: A(N,N), B(N,N), C(N,N)
+  INTEGER :: I,J,K,N
+  DO J=1,N
+     DO I=1,N
+        DO K=1,N
+           C(I,J)=C(I,J)+A(I,K)*B(K,J)
+        ENDDO
+     ENDDO
+  ENDDO
+END SUBROUTINE mult
+
+! { dg-final { scan-tree-dump-times "number of SCoPs: 1" 1 "graphite" { xfail *-*-* } } }
+! { dg-final { scan-tree-dump-times "will be loop blocked" 1 "graphite" { xfail *-*-* } } }
+! { dg-final { cleanup-tree-dump "graphite" } }
diff --git a/gcc/testsuite/gfortran.dg/graphite/pr41924.f90 b/gcc/testsuite/gfortran.dg/graphite/pr41924.f90
new file mode 100644
index 00000000000..f8dc8078e5c
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/graphite/pr41924.f90
@@ -0,0 +1,15 @@
+! { dg-options "-O2 -fgraphite-identity " }
+
+MODULE MAIN1
+ REAL , ALLOCATABLE :: HRVALD(:)
+END MODULE MAIN1
+
+SUBROUTINE VOLCALC()
+ USE MAIN1
+ INTEGER :: ITYP
+ LOGICAL :: WETSCIM
+
+ DO ITYP = 1 , 100
+    IF ( WETSCIM ) HRVALD(ITYP) = 0.0
+ ENDDO
+END SUBROUTINE VOLCALC
diff --git a/gcc/testsuite/gfortran.dg/graphite/run-id-2.f90 b/gcc/testsuite/gfortran.dg/graphite/run-id-2.f90
new file mode 100644
index 00000000000..c4fa1d06141
--- /dev/null
+++ b/gcc/testsuite/gfortran.dg/graphite/run-id-2.f90
@@ -0,0 +1,66 @@
+  IMPLICIT NONE
+  INTEGER, PARAMETER :: dp=KIND(0.0D0)
+  REAL(KIND=dp)      :: res
+
+  res=exp_radius_very_extended(  0    ,      1   ,       0      ,    1, &
+                               (/0.0D0,0.0D0,0.0D0/),&
+                               (/1.0D0,0.0D0,0.0D0/),&
+                               (/1.0D0,0.0D0,0.0D0/),&
+                                 1.0D0,1.0D0,1.0D0,1.0D0)
+  if (res.ne.1.0d0) call abort()
+
+CONTAINS
+
+ FUNCTION exp_radius_very_extended(la_min,la_max,lb_min,lb_max,ra,rb,rp,&
+                          zetp,eps,prefactor,cutoff) RESULT(radius)
+
+    INTEGER, INTENT(IN)                      :: la_min, la_max, lb_min, lb_max
+    REAL(KIND=dp), INTENT(IN)                :: ra(3), rb(3), rp(3), zetp, &
+                                                eps, prefactor, cutoff
+    REAL(KIND=dp)                            :: radius
+
+    INTEGER                                  :: i, ico, j, jco, la(3), lb(3), &
+                                                lxa, lxb, lya, lyb, lza, lzb
+    REAL(KIND=dp)                            :: bini, binj, coef(0:20), &
+                                                epsin_local, polycoef(0:60), &
+                                                prefactor_local, rad_a, &
+                                                rad_b, s1, s2
+
+    epsin_local=1.0E-2_dp
+
+    prefactor_local=prefactor*MAX(1.0_dp,cutoff)
+    rad_a=SQRT(SUM((ra-rp)**2))
+    rad_b=SQRT(SUM((rb-rp)**2))
+
+    polycoef(0:la_max+lb_max)=0.0_dp
+    DO lxa=0,la_max
+    DO lxb=0,lb_max
+       coef(0:la_max+lb_max)=0.0_dp
+       bini=1.0_dp
+       s1=1.0_dp
+       DO i=0,lxa
+          binj=1.0_dp
+          s2=1.0_dp
+          DO j=0,lxb
+             coef(lxa+lxb-i-j)=coef(lxa+lxb-i-j) + bini*binj*s1*s2
+             binj=(binj*(lxb-j))/(j+1)
+             s2=s2*(rad_b)
+          ENDDO
+          bini=(bini*(lxa-i))/(i+1)
+          s1=s1*(rad_a)
+       ENDDO
+       DO i=0,lxa+lxb
+          polycoef(i)=MAX(polycoef(i),coef(i))
+       ENDDO
+    ENDDO
+    ENDDO
+
+    polycoef(0:la_max+lb_max)=polycoef(0:la_max+lb_max)*prefactor_local
+    radius=0.0_dp
+    DO i=0,la_max+lb_max
+          radius=MAX(radius,polycoef(i)**(i+1))
+    ENDDO
+
+  END FUNCTION exp_radius_very_extended
+
+END