fixing bug in src/Evolution/odeintsgUnified.f

scarrazza · Feb 1, 2017 · 5c84d93 · 5c84d93
1 parent 0a3fd94
commit 5c84d93
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/Evolution/odeintsgUnified.f b/src/Evolution/odeintsgUnified.f
@@ -419,7 +419,7 @@ subroutine derivssgUnifiedS1(t,M,dMdt)
                   do beta=alpha,nin(igrid)
                      dMdt(i,j,alpha,beta) = 0d0
                      do l=1,5
-                        do delta=0,nin(igrid)-alpha
+                        do delta=alpha,beta
                            dMdt(i,j,alpha,beta) = dMdt(i,j,alpha,beta)
      1                          + integ2(alpha,delta,i,l)
      2                          * M(l,j,delta,beta)
@@ -888,7 +888,7 @@ subroutine derivssgUnifiedS2(t,M,dMdt)
                   do beta=alpha,nin(igrid)
                      dMdt(i,j,alpha,beta) = 0d0
                      do l=1,2
-                        do delta=0,nin(igrid)-alpha
+                        do delta=alpha,beta
                            dMdt(i,j,alpha,beta) = dMdt(i,j,alpha,beta)
      1                          + integ2(alpha,delta,i,l)
      2                          * M(l,j,delta,beta)