From fd470ee623f3a9585a5be3c6558fd2e99cf9af98 Mon Sep 17 00:00:00 2001 From: Joshua Lampert Date: Fri, 14 Jun 2024 16:00:08 +0200 Subject: [PATCH] remove l_hat2 since it this case is impossible --- ext/SummationByPartsOperatorsOptimForwardDiffExt.jl | 10 +++------- 1 file changed, 3 insertions(+), 7 deletions(-) diff --git a/ext/SummationByPartsOperatorsOptimForwardDiffExt.jl b/ext/SummationByPartsOperatorsOptimForwardDiffExt.jl index c20b7ddb..920ad83d 100644 --- a/ext/SummationByPartsOperatorsOptimForwardDiffExt.jl +++ b/ext/SummationByPartsOperatorsOptimForwardDiffExt.jl @@ -176,13 +176,9 @@ end D_plus_one_half_trunc = trunc(D_plus_one_half) if D_plus_one_half == D_plus_one_half_trunc int_D_plus_one_half = trunc(Int, D_plus_one_half_trunc) - l_hat1 = N + int_D_plus_one_half - 1 - l_hat2 = N - int_D_plus_one_half - if 1 <= l_hat1 <= i - 1 - daij_dsigmak[i, j, k] -= V[l_hat1, j] - end - if 1 <= l_hat2 <= i - 1 - daij_dsigmak[i, j, k] -= V[l_hat2, j] + l_hat = N - int_D_plus_one_half + if 1 <= l_hat <= i - 1 + daij_dsigmak[i, j, k] -= V[l_hat, j] end end end