Using an implicit-explicit (IMEX) Runge-Kutta solver based on Ariadne with Trixi.jl
using Trixi
using Theseus
using CairoMakieNotes: You must disable both Polyester and LoopVectorization for Enzyme to be able to differentiate Trixi.jl.
LocalPreferences.jl
[Trixi]
loop_vectorization = false
backend = "static"@assert Trixi._PREFERENCE_THREADING !== :polyester
@assert !Trixi._PREFERENCE_LOOPVECTORIZATIONFirst call to load callbacks
trixi_include(@__MODULE__, joinpath(examples_dir(), "p4est_2d_dgsem", "elixir_navierstokes_vortex_street.jl"), sol = nothing, analysis_interval = 10);
stepsize_callback = StepsizeCallback(cfl = 1.0)
callbacks = CallbackSet(summary_callback, analysis_callback, alive_callback, stepsize_callback)
tspan = (0.0, 0.5)
trixi_include(@__MODULE__, joinpath(examples_dir(), "p4est_2d_dgsem", "elixir_navierstokes_vortex_street.jl"), sol = nothing, callbacks = callbacks, tspan = tspan);
###############################################################################[ Info: You just called `trixi_include`. Julia may now compile the code, please be patient.
[ Info: You just called `trixi_include`. Julia may now compile the code, please be patient.run the simulation
sol = solve(
ode,
Theseus.ARS443(); # ARS111, ARS222, ARS443
dt = 0.001, # solve needs some value here but it will be overwritten by the stepsize_callback
ode_default_options()..., callback = callbacks,
krylov_algo = :gmres,
);
████████╗██████╗ ██╗██╗ ██╗██╗
╚══██╔══╝██╔══██╗██║╚██╗██╔╝██║
██║ ██████╔╝██║ ╚███╔╝ ██║
██║ ██╔══██╗██║ ██╔██╗ ██║
██║ ██║ ██║██║██╔╝ ██╗██║
╚═╝ ╚═╝ ╚═╝╚═╝╚═╝ ╚═╝╚═╝
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ SemidiscretizationHyperbolicParabolic │
│ ═════════════════════════════════════ │
│ #spatial dimensions: ………………………… 2 │
│ mesh: ………………………………………………………………… P4estMesh{2, 2, Float64} │
│ hyperbolic equations: ……………………… CompressibleEulerEquations2D │
│ parabolic equations: ………………………… CompressibleNavierStokesDiffusion2D │
│ initial condition: ……………………………… initial_condition │
│ source terms: …………………………………………… nothing │
│ source terms parabolic: ………………… nothing │
│ solver: …………………………………………………………… DG │
│ parabolic solver: ………………………………… ViscousFormulationBassiRebay1 │
│ total #DOFs per field: …………………… 36128 │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ P4estMesh{2, 2, Float64} │
│ ════════════════════════ │
│ #trees: …………………………………………………………… 2258 │
│ current #cells: ……………………………………… 2258 │
│ polydeg: ………………………………………………………… 3 │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ CompressibleEulerEquations2D │
│ ════════════════════════════ │
│ #variables: ………………………………………………… 4 │
│ │ variable 1: …………………………………………… rho │
│ │ variable 2: …………………………………………… rho_v1 │
│ │ variable 3: …………………………………………… rho_v2 │
│ │ variable 4: …………………………………………… rho_e │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ DG{Float64} │
│ ═══════════ │
│ basis: ……………………………………………………………… LobattoLegendreBasis{Float64}(polydeg=3) │
│ mortar: …………………………………………………………… LobattoLegendreMortarL2{Float64}(polydeg=3) │
│ surface integral: ………………………………… SurfaceIntegralWeakForm │
│ │ surface flux: ……………………………………… FluxHLL(min_max_speed_davis) │
│ volume integral: …………………………………… VolumeIntegralWeakForm │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ AnalysisCallback │
│ ════════════════ │
│ interval: ……………………………………………………… 10 │
│ analyzer: ……………………………………………………… LobattoLegendreAnalyzer{Float64}(polydeg=6) │
│ │ error 1: …………………………………………………… l2_error │
│ │ error 2: …………………………………………………… linf_error │
│ │ integral 1: …………………………………………… entropy_timederivative │
│ save analysis to file: …………………… no │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ AliveCallback │
│ ═════════════ │
│ interval: ……………………………………………………… 1 │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ StepsizeCallback │
│ ════════════════ │
│ CFL Advective: ………………………………………… Returns{Float64}(1.0) │
│ CFL Diffusive: ………………………………………… Returns{Float64}(0.0) │
│ Interval: ……………………………………………………… 1 │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ Time integration │
│ ════════════════ │
│ Start time: ………………………………………………… 0.0 │
│ Final time: ………………………………………………… 0.5 │
│ time integrator: …………………………………… ARS443 │
│ adaptive: ……………………………………………………… false │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ Environment information │
│ ═══════════════════════ │
│ #threads: ……………………………………………………… 1 │
│ threading backend: ……………………………… static │
│ LoopVectorization: ……………………………… disabled │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 0 run time: 1.34920291e+00 s
Δt: 1.00000000e-03 └── GC time: 0.00000000e+00 s (0.000%)
sim. time: 0.00000000e+00 (0.000%) time/DOF/rhs!: NaN s
PID: Inf s
#DOFs per field: 36128 alloc'd memory: 709.642 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 8.96518598e-17 8.96518598e-17 0.00000000e+00 1.27756504e-15
Linf error: 4.44089210e-16 4.44089210e-16 0.00000000e+00 2.66453526e-15
∑∂S/∂U ⋅ Uₜ : -1.20647355e-03
────────────────────────────────────────────────────────────────────────────────────────────────────
┌ Warning: TODO forward zero-set of memorycopy used memset rather than runtime type
│ Caused by:
│ Stacktrace:
│ [1] copy
│ @ ./array.jl:350
│ [2] unaliascopy
│ @ ./abstractarray.jl:1516
│ [3] unalias
│ @ ./abstractarray.jl:1500
│ [4] broadcast_unalias
│ @ ./broadcast.jl:946
│ [5] preprocess
│ @ ./broadcast.jl:953
│ [6] preprocess_args
│ @ ./broadcast.jl:956
│ [7] preprocess_args
│ @ ./broadcast.jl:955
│ [8] preprocess
│ @ ./broadcast.jl:952
│ [9] preprocess_args
│ @ ./broadcast.jl:956
│ [10] preprocess_args
│ @ ./broadcast.jl:955
│ [11] preprocess
│ @ ./broadcast.jl:952
│ [12] preprocess_args
│ @ ./broadcast.jl:955
│ [13] preprocess
│ @ ./broadcast.jl:952
│ [14] override_bc_copyto!
│ @ ~/.julia/packages/Enzyme/C5rSn/src/compiler/interpreter.jl:818
│ [15] copyto!
│ @ ./broadcast.jl:925
│ [16] materialize!
│ @ ./broadcast.jl:883
│ [17] materialize!
│ @ ./broadcast.jl:880
│ [18] RKIMEX
│ @ ~/work/Ariadne.jl/Ariadne.jl/libs/Theseus/src/imex/imex.jl:33
└ @ Enzyme.Compiler ~/.julia/packages/Enzyme/C5rSn/src/rules/llvmrules.jl:812
┌ Warning: TODO forward zero-set of memorycopy used memset rather than runtime type
│ Caused by:
│ Stacktrace:
│ [1] copy
│ @ ./array.jl:350
│ [2] unaliascopy
│ @ ./abstractarray.jl:1516
│ [3] unalias
│ @ ./abstractarray.jl:1500
│ [4] broadcast_unalias
│ @ ./broadcast.jl:946
│ [5] preprocess
│ @ ./broadcast.jl:953
│ [6] preprocess_args
│ @ ./broadcast.jl:956
│ [7] preprocess_args
│ @ ./broadcast.jl:955
│ [8] preprocess
│ @ ./broadcast.jl:952
│ [9] preprocess_args
│ @ ./broadcast.jl:956
│ [10] preprocess_args
│ @ ./broadcast.jl:955
│ [11] preprocess
│ @ ./broadcast.jl:952
│ [12] override_bc_copyto!
│ @ ~/.julia/packages/Enzyme/C5rSn/src/compiler/interpreter.jl:818
│ [13] copyto!
│ @ ./broadcast.jl:925
│ [14] materialize!
│ @ ./broadcast.jl:883
│ [15] materialize!
│ @ ./broadcast.jl:880
│ [16] RKIMEX
│ @ ~/work/Ariadne.jl/Ariadne.jl/libs/Theseus/src/imex/imex.jl:33
└ @ Enzyme.Compiler ~/.julia/packages/Enzyme/C5rSn/src/rules/llvmrules.jl:812
#timesteps: 1 │ Δt: 4.6182e-03 │ sim. time: 4.6182e-03 (0.924%) │ run time: 1.0459e+02 s
#timesteps: 2 │ Δt: 4.4948e-03 │ sim. time: 9.1130e-03 (1.823%) │ run time: 1.1628e+02 s
#timesteps: 3 │ Δt: 4.4714e-03 │ sim. time: 1.3584e-02 (2.717%) │ run time: 1.2798e+02 s
#timesteps: 4 │ Δt: 4.4766e-03 │ sim. time: 1.8061e-02 (3.612%) │ run time: 1.3924e+02 s
#timesteps: 5 │ Δt: 4.4931e-03 │ sim. time: 2.2554e-02 (4.511%) │ run time: 1.5168e+02 s
#timesteps: 6 │ Δt: 4.5254e-03 │ sim. time: 2.7079e-02 (5.416%) │ run time: 1.6238e+02 s
#timesteps: 7 │ Δt: 4.5643e-03 │ sim. time: 3.1644e-02 (6.329%) │ run time: 1.7465e+02 s
#timesteps: 8 │ Δt: 4.5740e-03 │ sim. time: 3.6218e-02 (7.244%) │ run time: 1.8525e+02 s
#timesteps: 9 │ Δt: 4.5552e-03 │ sim. time: 4.0773e-02 (8.155%) │ run time: 1.9701e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 10 run time: 2.10130893e+02 s
Δt: 4.55501201e-03 └── GC time: 1.10088290e+01 s (5.239%)
sim. time: 4.53278729e-02 (9.066%) time/DOF/rhs!: 5.92387965e-06 s
PID: 1.15477748e-04 s
#DOFs per field: 36128 alloc'd memory: 921.417 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 4.64923431e-03 8.85283724e-03 4.72300190e-03 4.80992841e-02
Linf error: 7.16158036e-01 1.42833015e+00 6.28720604e-01 7.47949255e+00
∑∂S/∂U ⋅ Uₜ : -1.98606192e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 11 │ Δt: 4.5783e-03 │ sim. time: 4.9906e-02 (9.981%) │ run time: 2.1912e+02 s
#timesteps: 12 │ Δt: 4.6280e-03 │ sim. time: 5.4534e-02 (10.907%) │ run time: 2.3004e+02 s
#timesteps: 13 │ Δt: 4.6537e-03 │ sim. time: 5.9188e-02 (11.838%) │ run time: 2.3911e+02 s
#timesteps: 14 │ Δt: 4.6430e-03 │ sim. time: 6.3831e-02 (12.766%) │ run time: 2.4926e+02 s
#timesteps: 15 │ Δt: 4.6370e-03 │ sim. time: 6.8468e-02 (13.694%) │ run time: 2.5827e+02 s
#timesteps: 16 │ Δt: 4.5785e-03 │ sim. time: 7.3046e-02 (14.609%) │ run time: 2.6841e+02 s
#timesteps: 17 │ Δt: 4.5293e-03 │ sim. time: 7.7576e-02 (15.515%) │ run time: 2.7740e+02 s
#timesteps: 18 │ Δt: 4.4909e-03 │ sim. time: 8.2067e-02 (16.413%) │ run time: 2.8738e+02 s
#timesteps: 19 │ Δt: 4.4629e-03 │ sim. time: 8.6529e-02 (17.306%) │ run time: 2.9655e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 20 run time: 3.07040925e+02 s
Δt: 4.44400238e-03 └── GC time: 1.74888873e+01 s (5.696%)
sim. time: 9.09734380e-02 (18.195%) time/DOF/rhs!: 5.76131320e-06 s
PID: 5.35616910e-05 s
#DOFs per field: 36128 alloc'd memory: 1106.723 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 6.53365544e-03 1.26659983e-02 7.73995358e-03 6.75443615e-02
Linf error: 5.43630658e-01 1.03003896e+00 7.51581870e-01 5.56166488e+00
∑∂S/∂U ⋅ Uₜ : -1.20583267e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 21 │ Δt: 4.4324e-03 │ sim. time: 9.5406e-02 (19.081%) │ run time: 3.1584e+02 s
#timesteps: 22 │ Δt: 4.4259e-03 │ sim. time: 9.9832e-02 (19.966%) │ run time: 3.2483e+02 s
#timesteps: 23 │ Δt: 4.4223e-03 │ sim. time: 1.0425e-01 (20.851%) │ run time: 3.3497e+02 s
#timesteps: 24 │ Δt: 4.4198e-03 │ sim. time: 1.0867e-01 (21.735%) │ run time: 3.4398e+02 s
#timesteps: 25 │ Δt: 4.4168e-03 │ sim. time: 1.1309e-01 (22.618%) │ run time: 3.5393e+02 s
#timesteps: 26 │ Δt: 4.4126e-03 │ sim. time: 1.1750e-01 (23.501%) │ run time: 3.6302e+02 s
#timesteps: 27 │ Δt: 4.4068e-03 │ sim. time: 1.2191e-01 (24.382%) │ run time: 3.7199e+02 s
#timesteps: 28 │ Δt: 4.3996e-03 │ sim. time: 1.2631e-01 (25.262%) │ run time: 3.8212e+02 s
#timesteps: 29 │ Δt: 4.3915e-03 │ sim. time: 1.3070e-01 (26.140%) │ run time: 3.9111e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 30 run time: 4.02817700e+02 s
Δt: 4.38295054e-03 └── GC time: 2.47921974e+01 s (6.155%)
sim. time: 1.35083929e-01 (27.017%) time/DOF/rhs!: 5.83340333e-06 s
PID: 5.29345170e-05 s
#DOFs per field: 36128 alloc'd memory: 1102.546 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 7.67756568e-03 1.52760645e-02 9.94581650e-03 7.93270443e-02
Linf error: 5.02625662e-01 1.03367280e+00 8.28764923e-01 5.09487231e+00
∑∂S/∂U ⋅ Uₜ : -1.30655214e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 31 │ Δt: 4.3747e-03 │ sim. time: 1.3946e-01 (27.892%) │ run time: 4.1043e+02 s
#timesteps: 32 │ Δt: 4.3672e-03 │ sim. time: 1.4383e-01 (28.765%) │ run time: 4.2071e+02 s
#timesteps: 33 │ Δt: 4.3608e-03 │ sim. time: 1.4819e-01 (29.637%) │ run time: 4.2976e+02 s
#timesteps: 34 │ Δt: 4.3558e-03 │ sim. time: 1.5254e-01 (30.508%) │ run time: 4.3973e+02 s
#timesteps: 35 │ Δt: 4.3521e-03 │ sim. time: 1.5689e-01 (31.379%) │ run time: 4.4890e+02 s
#timesteps: 36 │ Δt: 4.3499e-03 │ sim. time: 1.6124e-01 (32.249%) │ run time: 4.5784e+02 s
#timesteps: 37 │ Δt: 4.3487e-03 │ sim. time: 1.6559e-01 (33.119%) │ run time: 4.6795e+02 s
#timesteps: 38 │ Δt: 4.3484e-03 │ sim. time: 1.6994e-01 (33.988%) │ run time: 4.7696e+02 s
#timesteps: 39 │ Δt: 4.3485e-03 │ sim. time: 1.7429e-01 (34.858%) │ run time: 4.8714e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 40 run time: 4.97247125e+02 s
Δt: 4.34863116e-03 └── GC time: 3.12804181e+01 s (6.291%)
sim. time: 1.78638705e-01 (35.728%) time/DOF/rhs!: 5.77952714e-06 s
PID: 5.21852837e-05 s
#DOFs per field: 36128 alloc'd memory: 1188.398 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 8.47801438e-03 1.71648658e-02 1.18338209e-02 8.74177339e-02
Linf error: 4.82193165e-01 1.10797823e+00 7.28017033e-01 4.84304661e+00
∑∂S/∂U ⋅ Uₜ : -1.19358050e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 41 │ Δt: 4.3484e-03 │ sim. time: 1.8299e-01 (36.597%) │ run time: 5.0600e+02 s
#timesteps: 42 │ Δt: 4.3474e-03 │ sim. time: 1.8733e-01 (37.467%) │ run time: 5.1499e+02 s
#timesteps: 43 │ Δt: 4.3456e-03 │ sim. time: 1.9168e-01 (38.336%) │ run time: 5.2497e+02 s
#timesteps: 44 │ Δt: 4.3428e-03 │ sim. time: 1.9602e-01 (39.205%) │ run time: 5.3409e+02 s
#timesteps: 45 │ Δt: 4.3392e-03 │ sim. time: 2.0036e-01 (40.072%) │ run time: 5.4305e+02 s
#timesteps: 46 │ Δt: 4.3349e-03 │ sim. time: 2.0470e-01 (40.939%) │ run time: 5.5320e+02 s
#timesteps: 47 │ Δt: 4.3302e-03 │ sim. time: 2.0903e-01 (41.805%) │ run time: 5.6214e+02 s
#timesteps: 48 │ Δt: 4.3252e-03 │ sim. time: 2.1335e-01 (42.670%) │ run time: 5.7236e+02 s
#timesteps: 49 │ Δt: 4.3203e-03 │ sim. time: 2.1767e-01 (43.535%) │ run time: 5.8135e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 50 run time: 5.92880402e+02 s
Δt: 4.31568887e-03 └── GC time: 3.85350969e+01 s (6.500%)
sim. time: 2.21988498e-01 (44.398%) time/DOF/rhs!: 5.81264533e-06 s
PID: 5.28518172e-05 s
#DOFs per field: 36128 alloc'd memory: 1084.039 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 9.06400123e-03 1.86724634e-02 1.35000328e-02 9.31248152e-02
Linf error: 5.03004190e-01 1.10957733e+00 7.17579661e-01 5.07149968e+00
∑∂S/∂U ⋅ Uₜ : -1.22996211e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 51 │ Δt: 4.3115e-03 │ sim. time: 2.2630e-01 (45.260%) │ run time: 6.0043e+02 s
#timesteps: 52 │ Δt: 4.3079e-03 │ sim. time: 2.3061e-01 (46.122%) │ run time: 6.1041e+02 s
#timesteps: 53 │ Δt: 4.3049e-03 │ sim. time: 2.3491e-01 (46.983%) │ run time: 6.1934e+02 s
#timesteps: 54 │ Δt: 4.3027e-03 │ sim. time: 2.3922e-01 (47.843%) │ run time: 6.2918e+02 s
#timesteps: 55 │ Δt: 4.3010e-03 │ sim. time: 2.4352e-01 (48.703%) │ run time: 6.3779e+02 s
#timesteps: 56 │ Δt: 4.3000e-03 │ sim. time: 2.4782e-01 (49.563%) │ run time: 6.4586e+02 s
#timesteps: 57 │ Δt: 4.2994e-03 │ sim. time: 2.5212e-01 (50.423%) │ run time: 6.5447e+02 s
#timesteps: 58 │ Δt: 4.2992e-03 │ sim. time: 2.5642e-01 (51.283%) │ run time: 6.6179e+02 s
#timesteps: 59 │ Δt: 4.2992e-03 │ sim. time: 2.6071e-01 (52.143%) │ run time: 6.6909e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 60 run time: 6.77929466e+02 s
Δt: 4.29925400e-03 └── GC time: 4.31844753e+01 s (6.370%)
sim. time: 2.65013456e-01 (53.003%) time/DOF/rhs!: 5.63418945e-06 s
PID: 4.69930703e-05 s
#DOFs per field: 36128 alloc'd memory: 1252.737 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 9.53070797e-03 1.98856287e-02 1.49628554e-02 9.73582091e-02
Linf error: 4.32755742e-01 1.04262767e+00 7.27908089e-01 4.26683140e+00
∑∂S/∂U ⋅ Uₜ : -1.11519114e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 61 │ Δt: 4.2993e-03 │ sim. time: 2.6931e-01 (53.863%) │ run time: 6.8512e+02 s
#timesteps: 62 │ Δt: 4.2993e-03 │ sim. time: 2.7361e-01 (54.722%) │ run time: 6.9295e+02 s
#timesteps: 63 │ Δt: 4.2990e-03 │ sim. time: 2.7791e-01 (55.582%) │ run time: 7.0148e+02 s
#timesteps: 64 │ Δt: 4.2985e-03 │ sim. time: 2.8221e-01 (56.442%) │ run time: 7.0933e+02 s
#timesteps: 65 │ Δt: 4.2979e-03 │ sim. time: 2.8651e-01 (57.301%) │ run time: 7.1675e+02 s
#timesteps: 66 │ Δt: 4.2970e-03 │ sim. time: 2.9080e-01 (58.161%) │ run time: 7.2499e+02 s
#timesteps: 67 │ Δt: 4.2959e-03 │ sim. time: 2.9510e-01 (59.020%) │ run time: 7.3229e+02 s
#timesteps: 68 │ Δt: 4.2948e-03 │ sim. time: 2.9940e-01 (59.879%) │ run time: 7.4059e+02 s
#timesteps: 69 │ Δt: 4.2936e-03 │ sim. time: 3.0369e-01 (60.738%) │ run time: 7.4795e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 70 run time: 7.57785698e+02 s
Δt: 4.29240354e-03 └── GC time: 4.94311952e+01 s (6.523%)
sim. time: 3.07981065e-01 (61.596%) time/DOF/rhs!: 5.76608950e-06 s
PID: 4.41184285e-05 s
#DOFs per field: 36128 alloc'd memory: 1190.559 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 9.90399293e-03 2.08711478e-02 1.62492564e-02 1.00514287e-01
Linf error: 3.98488709e-01 1.04898027e+00 7.35060870e-01 3.89853935e+00
∑∂S/∂U ⋅ Uₜ : -1.17807394e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 71 │ Δt: 4.2913e-03 │ sim. time: 3.1227e-01 (62.454%) │ run time: 7.6422e+02 s
#timesteps: 72 │ Δt: 4.2902e-03 │ sim. time: 3.1656e-01 (63.313%) │ run time: 7.7189e+02 s
#timesteps: 73 │ Δt: 4.2892e-03 │ sim. time: 3.2085e-01 (64.170%) │ run time: 7.8069e+02 s
#timesteps: 74 │ Δt: 4.2882e-03 │ sim. time: 3.2514e-01 (65.028%) │ run time: 7.8844e+02 s
#timesteps: 75 │ Δt: 4.2873e-03 │ sim. time: 3.2943e-01 (65.885%) │ run time: 7.9675e+02 s
#timesteps: 76 │ Δt: 4.2864e-03 │ sim. time: 3.3371e-01 (66.743%) │ run time: 8.0421e+02 s
#timesteps: 77 │ Δt: 4.2854e-03 │ sim. time: 3.3800e-01 (67.600%) │ run time: 8.1244e+02 s
#timesteps: 78 │ Δt: 4.2844e-03 │ sim. time: 3.4228e-01 (68.457%) │ run time: 8.1991e+02 s
#timesteps: 79 │ Δt: 4.2833e-03 │ sim. time: 3.4657e-01 (69.313%) │ run time: 8.2723e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 80 run time: 8.36117462e+02 s
Δt: 4.28121331e-03 └── GC time: 5.41306511e+01 s (6.474%)
sim. time: 3.50847983e-01 (70.170%) time/DOF/rhs!: 5.65210129e-06 s
PID: 4.32712003e-05 s
#DOFs per field: 36128 alloc'd memory: 1266.272 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 1.02147848e-02 2.16637967e-02 1.73659533e-02 1.02893942e-01
Linf error: 3.71854479e-01 1.05297438e+00 7.20233874e-01 3.56770023e+00
∑∂S/∂U ⋅ Uₜ : -1.14063852e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 81 │ Δt: 4.2790e-03 │ sim. time: 3.5513e-01 (71.025%) │ run time: 8.4311e+02 s
#timesteps: 82 │ Δt: 4.2769e-03 │ sim. time: 3.5940e-01 (71.881%) │ run time: 8.5051e+02 s
#timesteps: 83 │ Δt: 4.2748e-03 │ sim. time: 3.6368e-01 (72.736%) │ run time: 8.5784e+02 s
#timesteps: 84 │ Δt: 4.2729e-03 │ sim. time: 3.6795e-01 (73.590%) │ run time: 8.6618e+02 s
#timesteps: 85 │ Δt: 4.2709e-03 │ sim. time: 3.7222e-01 (74.445%) │ run time: 8.7353e+02 s
#timesteps: 86 │ Δt: 4.2691e-03 │ sim. time: 3.7649e-01 (75.298%) │ run time: 8.8195e+02 s
#timesteps: 87 │ Δt: 4.2674e-03 │ sim. time: 3.8076e-01 (76.152%) │ run time: 8.8933e+02 s
#timesteps: 88 │ Δt: 4.2657e-03 │ sim. time: 3.8502e-01 (77.005%) │ run time: 8.9656e+02 s
#timesteps: 89 │ Δt: 4.2641e-03 │ sim. time: 3.8929e-01 (77.858%) │ run time: 9.0484e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 90 run time: 9.13737020e+02 s
Δt: 4.26264722e-03 └── GC time: 5.96644629e+01 s (6.530%)
sim. time: 3.93551528e-01 (78.710%) time/DOF/rhs!: 5.65335619e-06 s
PID: 4.28800074e-05 s
#DOFs per field: 36128 alloc'd memory: 1065.174 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 1.04832924e-02 2.23447434e-02 1.83272747e-02 1.04748173e-01
Linf error: 4.04043565e-01 1.06669114e+00 7.11255592e-01 3.91773525e+00
∑∂S/∂U ⋅ Uₜ : -1.14330464e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 91 │ Δt: 4.2612e-03 │ sim. time: 3.9781e-01 (79.563%) │ run time: 9.1969e+02 s
#timesteps: 92 │ Δt: 4.2599e-03 │ sim. time: 4.0207e-01 (80.415%) │ run time: 9.2710e+02 s
#timesteps: 93 │ Δt: 4.2587e-03 │ sim. time: 4.0633e-01 (81.266%) │ run time: 9.3536e+02 s
#timesteps: 94 │ Δt: 4.2575e-03 │ sim. time: 4.1059e-01 (82.118%) │ run time: 9.4271e+02 s
#timesteps: 95 │ Δt: 4.2564e-03 │ sim. time: 4.1485e-01 (82.969%) │ run time: 9.5097e+02 s
#timesteps: 96 │ Δt: 4.2554e-03 │ sim. time: 4.1910e-01 (83.820%) │ run time: 9.5837e+02 s
#timesteps: 97 │ Δt: 4.2544e-03 │ sim. time: 4.2336e-01 (84.671%) │ run time: 9.6660e+02 s
#timesteps: 98 │ Δt: 4.2535e-03 │ sim. time: 4.2761e-01 (85.522%) │ run time: 9.7401e+02 s
#timesteps: 99 │ Δt: 4.2526e-03 │ sim. time: 4.3186e-01 (86.372%) │ run time: 9.8127e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 100 run time: 9.91016554e+02 s
Δt: 4.25181873e-03 └── GC time: 6.50785007e+01 s (6.567%)
sim. time: 4.36113158e-01 (87.223%) time/DOF/rhs!: 5.63664668e-06 s
PID: 4.26913810e-05 s
#DOFs per field: 36128 alloc'd memory: 890.974 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 1.07266368e-02 2.29091724e-02 1.91465963e-02 1.06238167e-01
Linf error: 3.67962100e-01 1.07731365e+00 7.12625840e-01 3.57594740e+00
∑∂S/∂U ⋅ Uₜ : -1.13540624e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 101 │ Δt: 4.2510e-03 │ sim. time: 4.4036e-01 (88.073%) │ run time: 9.9710e+02 s
#timesteps: 102 │ Δt: 4.2503e-03 │ sim. time: 4.4461e-01 (88.923%) │ run time: 1.0044e+03 s
#timesteps: 103 │ Δt: 4.2496e-03 │ sim. time: 4.4886e-01 (89.773%) │ run time: 1.0117e+03 s
#timesteps: 104 │ Δt: 4.2490e-03 │ sim. time: 4.5311e-01 (90.623%) │ run time: 1.0190e+03 s
#timesteps: 105 │ Δt: 4.2484e-03 │ sim. time: 4.5736e-01 (91.472%) │ run time: 1.0272e+03 s
#timesteps: 106 │ Δt: 4.2478e-03 │ sim. time: 4.6161e-01 (92.322%) │ run time: 1.0346e+03 s
#timesteps: 107 │ Δt: 4.2473e-03 │ sim. time: 4.6586e-01 (93.171%) │ run time: 1.0420e+03 s
#timesteps: 108 │ Δt: 4.2468e-03 │ sim. time: 4.7010e-01 (94.021%) │ run time: 1.0493e+03 s
#timesteps: 109 │ Δt: 4.2463e-03 │ sim. time: 4.7435e-01 (94.870%) │ run time: 1.0576e+03 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 110 run time: 1.06648144e+03 s
Δt: 4.24591874e-03 └── GC time: 6.86620694e+01 s (6.438%)
sim. time: 4.78595613e-01 (95.719%) time/DOF/rhs!: 5.56852562e-06 s
PID: 4.16901822e-05 s
#DOFs per field: 36128 alloc'd memory: 1042.108 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 1.09501166e-02 2.34077070e-02 1.98398506e-02 1.07485665e-01
Linf error: 3.73765631e-01 1.08279798e+00 7.10650160e-01 3.31172185e+00
∑∂S/∂U ⋅ Uₜ : -1.15721767e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 111 │ Δt: 4.2456e-03 │ sim. time: 4.8284e-01 (96.568%) │ run time: 1.0725e+03 s
#timesteps: 112 │ Δt: 4.2453e-03 │ sim. time: 4.8709e-01 (97.417%) │ run time: 1.0808e+03 s
#timesteps: 113 │ Δt: 4.2450e-03 │ sim. time: 4.9133e-01 (98.266%) │ run time: 1.0881e+03 s
#timesteps: 114 │ Δt: 4.2448e-03 │ sim. time: 4.9558e-01 (99.115%) │ run time: 1.0964e+03 s
#timesteps: 115 │ Δt: 4.2447e-03 │ sim. time: 4.9982e-01 (99.964%) │ run time: 1.1038e+03 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 116 run time: 1.10933496e+03 s
Δt: 1.79026594e-04 └── GC time: 7.15266108e+01 s (6.448%)
sim. time: 5.00000000e-01 (100.000%) time/DOF/rhs!: 5.53709366e-06 s
PID: 3.93899533e-05 s
#DOFs per field: 36128 alloc'd memory: 1169.880 MiB
#elements: 2258
Variable: rho rho_v1 rho_v2 rho_e
L2 error: 1.10537553e-02 2.36419928e-02 2.01463055e-02 1.08021750e-01
Linf error: 3.79655915e-01 1.08633228e+00 7.08879871e-01 3.29414734e+00
∑∂S/∂U ⋅ Uₜ : -1.16165019e-04
────────────────────────────────────────────────────────────────────────────────────────────────────
────────────────────────────────────────────────────────────────────────────────────────────────────
Trixi.jl simulation finished. Final time: 0.5 Time steps: 116 (accepted), 116 (total)
────────────────────────────────────────────────────────────────────────────────────────────────────
───────────────────────────────────────────────────────────────────────────────────────
Trixi.jl Time Allocations
─────────────────────── ────────────────────────
Tot / % measured: 1109s / 89.0% 231GiB / 83.8%
Section ncalls time %tot avg alloc %tot avg
───────────────────────────────────────────────────────────────────────────────────────
parabolic rhs! 4.83k 984s 99.6% 204ms 193GiB 99.8% 40.9MiB
prolong2interfaces 4.83k 929s 94.0% 192ms 190GiB 98.3% 40.3MiB
calculate gradient 4.83k 20.1s 2.0% 4.17ms 3.33KiB 0.0% 0.71B
volume integral 4.83k 9.92s 1.0% 2.05ms 0.00B 0.0% 0.00B
surface integral 4.83k 4.77s 0.5% 987μs 0.00B 0.0% 0.00B
Jacobian 4.83k 1.65s 0.2% 341μs 0.00B 0.0% 0.00B
interface flux 4.83k 1.55s 0.2% 321μs 0.00B 0.0% 0.00B
prolong2interfaces 4.83k 1.49s 0.2% 308μs 0.00B 0.0% 0.00B
reset gradients 4.83k 662ms 0.1% 137μs 0.00B 0.0% 0.00B
boundary flux 4.83k 41.5ms 0.0% 8.58μs 0.00B 0.0% 0.00B
prolong2boundaries 4.83k 31.0ms 0.0% 6.41μs 0.00B 0.0% 0.00B
~calculate gradient~ 4.83k 17.9ms 0.0% 3.71μs 3.33KiB 0.0% 0.71B
prolong2mortars 4.83k 2.75ms 0.0% 568ns 0.00B 0.0% 0.00B
mortar flux 4.83k 2.68ms 0.0% 554ns 0.00B 0.0% 0.00B
prolong2boundaries 4.83k 14.1s 1.4% 2.93ms 2.84GiB 1.5% 616KiB
volume integral 4.83k 7.88s 0.8% 1.63ms 0.00B 0.0% 0.00B
calculate viscous fluxes 4.83k 7.16s 0.7% 1.48ms 0.00B 0.0% 0.00B
interface flux 4.83k 1.50s 0.2% 311μs 0.00B 0.0% 0.00B
transform variables 4.83k 1.43s 0.1% 296μs 0.00B 0.0% 0.00B
surface integral 4.83k 1.10s 0.1% 227μs 0.00B 0.0% 0.00B
Jacobian 4.83k 839ms 0.1% 174μs 0.00B 0.0% 0.00B
reset ∂u/∂t 4.83k 354ms 0.0% 73.2μs 0.00B 0.0% 0.00B
boundary flux 4.83k 62.2ms 0.0% 12.9μs 0.00B 0.0% 0.00B
~parabolic rhs!~ 4.83k 48.7ms 0.0% 10.1μs 5.38KiB 0.0% 1.14B
prolong2mortars 4.83k 16.5ms 0.0% 3.41μs 894KiB 0.0% 190B
mortar flux 4.83k 2.28ms 0.0% 472ns 0.00B 0.0% 0.00B
analyze solution 13 2.11s 0.2% 162ms 420MiB 0.2% 32.3MiB
rhs! 580 2.07s 0.2% 3.57ms 4.78KiB 0.0% 8.44B
volume integral 580 1.08s 0.1% 1.86ms 0.00B 0.0% 0.00B
interface flux 580 488ms 0.0% 842μs 0.00B 0.0% 0.00B
prolong2interfaces 580 179ms 0.0% 309μs 0.00B 0.0% 0.00B
Jacobian 580 141ms 0.0% 243μs 0.00B 0.0% 0.00B
surface integral 580 133ms 0.0% 229μs 0.00B 0.0% 0.00B
reset ∂u/∂t 580 30.6ms 0.0% 52.8μs 0.00B 0.0% 0.00B
boundary flux 580 13.4ms 0.0% 23.2μs 0.00B 0.0% 0.00B
prolong2boundaries 580 3.40ms 0.0% 5.86μs 0.00B 0.0% 0.00B
~rhs!~ 580 2.26ms 0.0% 3.89μs 4.78KiB 0.0% 8.44B
prolong2mortars 580 205μs 0.0% 353ns 0.00B 0.0% 0.00B
mortar flux 580 197μs 0.0% 339ns 0.00B 0.0% 0.00B
source terms 580 54.1μs 0.0% 93.3ns 0.00B 0.0% 0.00B
calculate dt 117 53.3ms 0.0% 456μs 16.5KiB 0.0% 144B
───────────────────────────────────────────────────────────────────────────────────────This page was generated using Literate.jl.