Missed one if

aadler · Jun 9, 2024 · bbe69e3 · bbe69e3
1 parent 028762d
commit bbe69e3
Showing 1 changed file with 2 additions and 4 deletions.
diff --git a/src/delaporte.f90 b/src/delaporte.f90
@@ -160,9 +160,7 @@ subroutine ddelap_f(x, nx, a, na, b, nb, l, nl, lg, threads, pmfv) &
         do i = 1, nx
             pmfv(i) = ddelap_f_s(x(i), a(imk(i, na)), b(imk(i, nb)), &
             l(imk(i, nl)))
-            if (lg == 1) then
-                pmfv(i) = log(pmfv(i))
-            end if
+            if (lg == 1) pmfv(i) = log(pmfv(i))
         end do
         !$omp end parallel do simd
 
@@ -199,7 +197,7 @@ pure elemental function pdelap_f_s(q, alpha, beta, lambda) result(cdf)
             do i = 1_INT64, k
                 cdf = cdf + ddelap_f_s(real(i, c_double), alpha, beta, lambda)
             end do
-        cdf = cFPe(cdf)                           ! Clear floating point errors
+            cdf = cFPe(cdf)                       ! Clear floating point errors
         end if
 
     end function pdelap_f_s