diff --git a/src/delaporte.f90 b/src/delaporte.f90 index f199f0f..416fcc7 100644 --- 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