From 8070dc52710e8426d55ee40339c719e8f6502179 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 23 Jul 2021 11:18:54 -0400 Subject: [PATCH] Rearranging loops for efficiency --- src/obl/obl.f90 | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/src/obl/obl.f90 b/src/obl/obl.f90 index 8792f2399..f027908f9 100644 --- a/src/obl/obl.f90 +++ b/src/obl/obl.f90 @@ -57,8 +57,8 @@ module subroutine obl_acc_pl(self, system) cb%aobl(i) = -sum(pl%Gmass(1:npl) * pl%aobl(i, 1:npl)) / cb%Gmass end do - do i = 1, NDIM - pl%ah(i, 1:npl) = pl%ah(i, 1:npl) + pl%aobl(i, 1:npl) - cb%aobl(i) + do i = 1, npl + pl%ah(:, i) = pl%ah(:, i) + pl%aobl(:, i) - cb%aobl(:) end do end associate @@ -89,8 +89,8 @@ module subroutine obl_acc_tp(self, system) aoblcb = cb%aoblend end if - do i = 1, NDIM - tp%ah(i, 1:ntp) = tp%ah(i, 1:ntp) + tp%aobl(i, 1:ntp) - aoblcb(i) + do i = 1, ntp + tp%ah(:, i) = tp%ah(:, i) + tp%aobl(:, i) - aoblcb(:) end do end associate