Skip to content

Commit 933e069

Browse files
DanielDoehringDanielDoehring
committed
set_zero! instead of reset_du 
1 parent d0708d9 commit 933e069

File tree

16 files changed

+30
-31
lines changed

16 files changed

+30
-31
lines changed

‎src/solvers/dgmulti/dg.jl‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -154,7 +154,7 @@ function allocate_nested_array(uEltype, nvars, array_dimensions, dg)
154154
nvars))
155155
end
156156

157-
function reset_du!(du, dg::DGMulti, other_args...)
157+
function set_zero!(du, dg::DGMulti, other_args...)
158158
@threaded for i in eachindex(du)
159159
du[i] = zero(eltype(du))
160160
end
@@ -665,7 +665,7 @@ end
665665
function rhs!(du, u, t, mesh, equations,
666666
boundary_conditions::BC, source_terms::Source,
667667
dg::DGMulti, cache) where {BC, Source}
668-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
668+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
669669

670670
@trixi_timeit timer() "volume integral" begin
671671
calc_volume_integral!(du, u, mesh,

‎src/solvers/dgmulti/dg_parabolic.jl‎

Lines changed: 4 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -166,7 +166,7 @@ function calc_gradient!(gradients, u::StructArray, t, mesh::DGMultiMesh,
166166
boundary_conditions, dg::DGMulti, parabolic_scheme,
167167
cache, cache_parabolic)
168168
for dim in eachindex(gradients)
169-
reset_du!(gradients[dim], dg)
169+
set_zero!(gradients[dim], dg)
170170
end
171171

172172
calc_gradient_volume_integral!(gradients, u, mesh, equations, dg, cache,
@@ -303,7 +303,7 @@ function calc_viscous_fluxes!(flux_viscous, u, gradients, mesh::DGMultiMesh,
303303
equations::AbstractEquationsParabolic,
304304
dg::DGMulti, cache, cache_parabolic)
305305
for dim in eachdim(mesh)
306-
reset_du!(flux_viscous[dim], dg)
306+
set_zero!(flux_viscous[dim], dg)
307307
end
308308

309309
(; local_u_values_threaded) = cache_parabolic
@@ -430,7 +430,7 @@ function calc_divergence!(du, u::StructArray, t, flux_viscous, mesh::DGMultiMesh
430430
equations::AbstractEquationsParabolic,
431431
boundary_conditions, dg::DGMulti, parabolic_scheme, cache,
432432
cache_parabolic)
433-
reset_du!(du, dg)
433+
set_zero!(du, dg)
434434

435435
calc_divergence_volume_integral!(du, u, flux_viscous, mesh, equations, dg, cache,
436436
cache_parabolic)
@@ -474,7 +474,7 @@ function rhs_parabolic!(du, u, t, mesh::DGMultiMesh,
474474
equations_parabolic::AbstractEquationsParabolic,
475475
boundary_conditions, source_terms,
476476
dg::DGMulti, parabolic_scheme, cache, cache_parabolic)
477-
reset_du!(du, dg)
477+
set_zero!(du, dg)
478478

479479
@trixi_timeit timer() "transform variables" begin
480480
(; u_transformed, gradients, flux_viscous) = cache_parabolic

‎src/solvers/dgmulti/flux_differencing.jl‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -618,7 +618,7 @@ end
618618
# Also called by DGMultiFluxDiff{<:GaussSBP} solvers.
619619
function rhs!(du, u, t, mesh, equations, boundary_conditions::BC,
620620
source_terms::Source, dg::DGMultiFluxDiff, cache) where {Source, BC}
621-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
621+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
622622

623623
# this function evaluates the solution at volume and face quadrature points (which was previously
624624
# done in `prolong2interfaces` and `calc_volume_integral`)
@@ -664,7 +664,7 @@ end
664664
function rhs!(du, u, t, mesh, equations,
665665
boundary_conditions::BC, source_terms::Source,
666666
dg::DGMultiFluxDiffSBP, cache) where {BC, Source}
667-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
667+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
668668

669669
@trixi_timeit timer() "volume integral" calc_volume_integral!(du, u, mesh,
670670
have_nonconservative_terms(equations),

‎src/solvers/dgmulti/flux_differencing_gauss_sbp.jl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -585,7 +585,7 @@ end
585585
function rhs!(du, u, t, mesh, equations, boundary_conditions::BC,
586586
source_terms::Source, dg::DGMultiFluxDiff{<:GaussSBP},
587587
cache) where {Source, BC}
588-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
588+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
589589

590590
# this function evaluates the solution at volume and face quadrature points (which was previously
591591
# done in `prolong2interfaces` and `calc_volume_integral`)

‎src/solvers/dgsem_p4est/dg_2d_parabolic.jl‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -89,7 +89,7 @@ function rhs_parabolic!(du, u, t, mesh::Union{P4estMesh{2}, P4estMesh{3}},
8989
# TODO: parabolic; reconsider current data structure reuse strategy
9090

9191
# Reset du
92-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
92+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
9393

9494
# Calculate volume integral
9595
@trixi_timeit timer() "volume integral" begin
@@ -158,8 +158,8 @@ function calc_gradient!(gradients, u_transformed, t,
158158

159159
# Reset du
160160
@trixi_timeit timer() "reset gradients" begin
161-
reset_du!(gradients_x, dg, cache)
162-
reset_du!(gradients_y, dg, cache)
161+
set_zero!(gradients_x, dg, cache)
162+
set_zero!(gradients_y, dg, cache)
163163
end
164164

165165
# Calculate volume integral

‎src/solvers/dgsem_p4est/dg_3d_parabolic.jl‎

Lines changed: 3 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -35,9 +35,9 @@ function calc_gradient!(gradients, u_transformed, t,
3535

3636
# Reset du
3737
@trixi_timeit timer() "reset gradients" begin
38-
reset_du!(gradients_x, dg, cache)
39-
reset_du!(gradients_y, dg, cache)
40-
reset_du!(gradients_z, dg, cache)
38+
set_zero!(gradients_x, dg, cache)
39+
set_zero!(gradients_y, dg, cache)
40+
set_zero!(gradients_z, dg, cache)
4141
end
4242

4343
# Calculate volume integral

‎src/solvers/dgsem_p4est/dg_3d_parallel.jl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -29,7 +29,7 @@ function rhs!(du, u, t,
2929
end
3030

3131
# Reset du
32-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
32+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
3333

3434
# Calculate volume integral
3535
@trixi_timeit timer() "volume integral" begin

‎src/solvers/dgsem_structured/dg.jl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -40,7 +40,7 @@ function rhs!(du, u, t,
4040
boundary_conditions, source_terms::Source,
4141
dg::DG, cache) where {Source}
4242
# Reset du
43-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
43+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
4444

4545
# Calculate volume integral
4646
@trixi_timeit timer() "volume integral" begin

‎src/solvers/dgsem_tree/dg_1d_parabolic.jl‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -54,7 +54,7 @@ function rhs_parabolic!(du, u, t, mesh::TreeMesh{1},
5454
# TODO: parabolic; reconsider current data structure reuse strategy
5555

5656
# Reset du
57-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
57+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
5858

5959
# Calculate volume integral
6060
@trixi_timeit timer() "volume integral" begin
@@ -508,7 +508,7 @@ function calc_gradient!(gradients, u_transformed, t, mesh::TreeMesh{1},
508508

509509
# Reset du
510510
@trixi_timeit timer() "reset gradients" begin
511-
reset_du!(gradients, dg, cache)
511+
set_zero!(gradients, dg, cache)
512512
end
513513

514514
# Calculate volume integral

‎src/solvers/dgsem_tree/dg_2d.jl‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -111,7 +111,7 @@ function rhs!(du, u, t,
111111
boundary_conditions, source_terms::Source,
112112
dg::DG, cache) where {Source}
113113
# Reset du
114-
@trixi_timeit timer() "reset ∂u/∂t" reset_du!(du, dg, cache)
114+
@trixi_timeit timer() "reset ∂u/∂t" set_zero!(du, dg, cache)
115115

116116
# Calculate volume integral
117117
@trixi_timeit timer() "volume integral" begin

0 commit comments

Comments
 (0)