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: ……………………………………………… ParabolicFormulationBassiRebay1 │
│ 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_total │
└──────────────────────────────────────────────────────────────────────────────────────────────────┘
┌──────────────────────────────────────────────────────────────────────────────────────────────────┐
│ 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 Hyperbolic: …………………………………………………… Returns{Float64}(1.0) │
│ CFL Parabolic: ……………………………………………………… 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: 8.04609457e-01 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: 729.424 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 1 │ Δt: 4.6182e-03 │ sim. time: 4.6182e-03 (0.924%) │ run time: 8.0936e+01 s
#timesteps: 2 │ Δt: 4.4948e-03 │ sim. time: 9.1130e-03 (1.823%) │ run time: 9.2974e+01 s
#timesteps: 3 │ Δt: 4.4714e-03 │ sim. time: 1.3584e-02 (2.717%) │ run time: 1.0504e+02 s
#timesteps: 4 │ Δt: 4.4766e-03 │ sim. time: 1.8061e-02 (3.612%) │ run time: 1.1774e+02 s
#timesteps: 5 │ Δt: 4.4931e-03 │ sim. time: 2.2554e-02 (4.511%) │ run time: 1.3064e+02 s
#timesteps: 6 │ Δt: 4.5254e-03 │ sim. time: 2.7079e-02 (5.416%) │ run time: 1.4282e+02 s
#timesteps: 7 │ Δt: 4.5643e-03 │ sim. time: 3.1644e-02 (6.329%) │ run time: 1.5460e+02 s
#timesteps: 8 │ Δt: 4.5740e-03 │ sim. time: 3.6218e-02 (7.244%) │ run time: 1.6671e+02 s
#timesteps: 9 │ Δt: 4.5552e-03 │ sim. time: 4.0773e-02 (8.155%) │ run time: 1.7883e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 10 run time: 1.91871364e+02 s
Δt: 4.55501201e-03 └── GC time: 1.22060986e+01 s (6.362%)
sim. time: 4.53278729e-02 (9.066%) time/DOF/rhs!: 6.15115548e-06 s
PID: 1.05671261e-04 s
#DOFs per field: 36128 alloc'd memory: 697.007 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.0265e+02 s
#timesteps: 12 │ Δt: 4.6280e-03 │ sim. time: 5.4534e-02 (10.907%) │ run time: 2.1299e+02 s
#timesteps: 13 │ Δt: 4.6537e-03 │ sim. time: 5.9188e-02 (11.838%) │ run time: 2.2341e+02 s
#timesteps: 14 │ Δt: 4.6430e-03 │ sim. time: 6.3831e-02 (12.766%) │ run time: 2.3360e+02 s
#timesteps: 15 │ Δt: 4.6370e-03 │ sim. time: 6.8468e-02 (13.694%) │ run time: 2.4320e+02 s
#timesteps: 16 │ Δt: 4.5785e-03 │ sim. time: 7.3046e-02 (14.609%) │ run time: 2.5364e+02 s
#timesteps: 17 │ Δt: 4.5293e-03 │ sim. time: 7.7576e-02 (15.515%) │ run time: 2.6382e+02 s
#timesteps: 18 │ Δt: 4.4909e-03 │ sim. time: 8.2067e-02 (16.413%) │ run time: 2.7339e+02 s
#timesteps: 19 │ Δt: 4.4629e-03 │ sim. time: 8.6529e-02 (17.306%) │ run time: 2.8361e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 20 run time: 2.94093854e+02 s
Δt: 4.44400238e-03 └── GC time: 2.01853263e+01 s (6.864%)
sim. time: 9.09734380e-02 (18.195%) time/DOF/rhs!: 6.14122510e-06 s
PID: 5.65021167e-05 s
#DOFs per field: 36128 alloc'd memory: 868.600 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.0367e+02 s
#timesteps: 22 │ Δt: 4.4259e-03 │ sim. time: 9.9832e-02 (19.966%) │ run time: 3.1389e+02 s
#timesteps: 23 │ Δt: 4.4223e-03 │ sim. time: 1.0425e-01 (20.851%) │ run time: 3.2347e+02 s
#timesteps: 24 │ Δt: 4.4198e-03 │ sim. time: 1.0867e-01 (21.735%) │ run time: 3.3383e+02 s
#timesteps: 25 │ Δt: 4.4168e-03 │ sim. time: 1.1309e-01 (22.618%) │ run time: 3.4408e+02 s
#timesteps: 26 │ Δt: 4.4126e-03 │ sim. time: 1.1750e-01 (23.501%) │ run time: 3.5362e+02 s
#timesteps: 27 │ Δt: 4.4068e-03 │ sim. time: 1.2191e-01 (24.382%) │ run time: 3.6401e+02 s
#timesteps: 28 │ Δt: 4.3996e-03 │ sim. time: 1.2631e-01 (25.262%) │ run time: 3.7435e+02 s
#timesteps: 29 │ Δt: 4.3915e-03 │ sim. time: 1.3070e-01 (26.140%) │ run time: 3.8373e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 30 run time: 3.95094817e+02 s
Δt: 4.38295054e-03 └── GC time: 2.89063259e+01 s (7.316%)
sim. time: 1.35083929e-01 (27.017%) time/DOF/rhs!: 6.20392945e-06 s
PID: 5.58267930e-05 s
#DOFs per field: 36128 alloc'd memory: 733.126 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.0444e+02 s
#timesteps: 32 │ Δt: 4.3672e-03 │ sim. time: 1.4383e-01 (28.765%) │ run time: 4.1402e+02 s
#timesteps: 33 │ Δt: 4.3608e-03 │ sim. time: 1.4819e-01 (29.637%) │ run time: 4.2449e+02 s
#timesteps: 34 │ Δt: 4.3558e-03 │ sim. time: 1.5254e-01 (30.508%) │ run time: 4.3473e+02 s
#timesteps: 35 │ Δt: 4.3521e-03 │ sim. time: 1.5689e-01 (31.379%) │ run time: 4.4431e+02 s
#timesteps: 36 │ Δt: 4.3499e-03 │ sim. time: 1.6124e-01 (32.249%) │ run time: 4.5473e+02 s
#timesteps: 37 │ Δt: 4.3487e-03 │ sim. time: 1.6559e-01 (33.119%) │ run time: 4.6492e+02 s
#timesteps: 38 │ Δt: 4.3484e-03 │ sim. time: 1.6994e-01 (33.988%) │ run time: 4.7450e+02 s
#timesteps: 39 │ Δt: 4.3485e-03 │ sim. time: 1.7429e-01 (34.858%) │ run time: 4.8493e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 40 run time: 4.95740691e+02 s
Δt: 4.34863116e-03 └── GC time: 3.77068762e+01 s (7.606%)
sim. time: 1.78638705e-01 (35.728%) time/DOF/rhs!: 6.20338348e-06 s
PID: 5.56302322e-05 s
#DOFs per field: 36128 alloc'd memory: 1003.627 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.0445e+02 s
#timesteps: 42 │ Δt: 4.3474e-03 │ sim. time: 1.8733e-01 (37.467%) │ run time: 5.1501e+02 s
#timesteps: 43 │ Δt: 4.3456e-03 │ sim. time: 1.9168e-01 (38.336%) │ run time: 5.2529e+02 s
#timesteps: 44 │ Δt: 4.3428e-03 │ sim. time: 1.9602e-01 (39.205%) │ run time: 5.3469e+02 s
#timesteps: 45 │ Δt: 4.3392e-03 │ sim. time: 2.0036e-01 (40.072%) │ run time: 5.4519e+02 s
#timesteps: 46 │ Δt: 4.3349e-03 │ sim. time: 2.0470e-01 (40.939%) │ run time: 5.5542e+02 s
#timesteps: 47 │ Δt: 4.3302e-03 │ sim. time: 2.0903e-01 (41.805%) │ run time: 5.6485e+02 s
#timesteps: 48 │ Δt: 4.3252e-03 │ sim. time: 2.1335e-01 (42.670%) │ run time: 5.7527e+02 s
#timesteps: 49 │ Δt: 4.3203e-03 │ sim. time: 2.1767e-01 (43.535%) │ run time: 5.8548e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 50 run time: 5.95920228e+02 s
Δt: 4.31568887e-03 └── GC time: 4.57223138e+01 s (7.673%)
sim. time: 2.21988498e-01 (44.398%) time/DOF/rhs!: 6.14475292e-06 s
PID: 5.53628027e-05 s
#DOFs per field: 36128 alloc'd memory: 1062.959 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.0558e+02 s
#timesteps: 52 │ Δt: 4.3079e-03 │ sim. time: 2.3061e-01 (46.122%) │ run time: 6.1584e+02 s
#timesteps: 53 │ Δt: 4.3049e-03 │ sim. time: 2.3491e-01 (46.983%) │ run time: 6.2536e+02 s
#timesteps: 54 │ Δt: 4.3027e-03 │ sim. time: 2.3922e-01 (47.843%) │ run time: 6.3571e+02 s
#timesteps: 55 │ Δt: 4.3010e-03 │ sim. time: 2.4352e-01 (48.703%) │ run time: 6.4555e+02 s
#timesteps: 56 │ Δt: 4.3000e-03 │ sim. time: 2.4782e-01 (49.563%) │ run time: 6.5417e+02 s
#timesteps: 57 │ Δt: 4.2994e-03 │ sim. time: 2.5212e-01 (50.423%) │ run time: 6.6318e+02 s
#timesteps: 58 │ Δt: 4.2992e-03 │ sim. time: 2.5642e-01 (51.283%) │ run time: 6.7086e+02 s
#timesteps: 59 │ Δt: 4.2992e-03 │ sim. time: 2.6071e-01 (52.143%) │ run time: 6.7948e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 60 run time: 6.88846103e+02 s
Δt: 4.29925400e-03 └── GC time: 5.41882582e+01 s (7.867%)
sim. time: 2.65013456e-01 (53.003%) time/DOF/rhs!: 6.19640454e-06 s
PID: 5.12926993e-05 s
#DOFs per field: 36128 alloc'd memory: 759.227 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.9628e+02 s
#timesteps: 62 │ Δt: 4.2993e-03 │ sim. time: 2.7361e-01 (54.722%) │ run time: 7.0522e+02 s
#timesteps: 63 │ Δt: 4.2990e-03 │ sim. time: 2.7791e-01 (55.582%) │ run time: 7.1333e+02 s
#timesteps: 64 │ Δt: 4.2985e-03 │ sim. time: 2.8221e-01 (56.442%) │ run time: 7.2229e+02 s
#timesteps: 65 │ Δt: 4.2979e-03 │ sim. time: 2.8651e-01 (57.301%) │ run time: 7.2994e+02 s
#timesteps: 66 │ Δt: 4.2970e-03 │ sim. time: 2.9080e-01 (58.161%) │ run time: 7.3860e+02 s
#timesteps: 67 │ Δt: 4.2959e-03 │ sim. time: 2.9510e-01 (59.020%) │ run time: 7.4718e+02 s
#timesteps: 68 │ Δt: 4.2948e-03 │ sim. time: 2.9940e-01 (59.879%) │ run time: 7.5487e+02 s
#timesteps: 69 │ Δt: 4.2936e-03 │ sim. time: 3.0369e-01 (60.738%) │ run time: 7.6333e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 70 run time: 7.72052126e+02 s
Δt: 4.29240354e-03 └── GC time: 6.06552779e+01 s (7.856%)
sim. time: 3.07981065e-01 (61.596%) time/DOF/rhs!: 6.07908662e-06 s
PID: 4.59756914e-05 s
#DOFs per field: 36128 alloc'd memory: 823.216 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.8020e+02 s
#timesteps: 72 │ Δt: 4.2902e-03 │ sim. time: 3.1656e-01 (63.313%) │ run time: 7.8829e+02 s
#timesteps: 73 │ Δt: 4.2892e-03 │ sim. time: 3.2085e-01 (64.170%) │ run time: 7.9710e+02 s
#timesteps: 74 │ Δt: 4.2882e-03 │ sim. time: 3.2514e-01 (65.028%) │ run time: 8.0530e+02 s
#timesteps: 75 │ Δt: 4.2873e-03 │ sim. time: 3.2943e-01 (65.885%) │ run time: 8.1386e+02 s
#timesteps: 76 │ Δt: 4.2864e-03 │ sim. time: 3.3371e-01 (66.743%) │ run time: 8.2156e+02 s
#timesteps: 77 │ Δt: 4.2854e-03 │ sim. time: 3.3800e-01 (67.600%) │ run time: 8.3013e+02 s
#timesteps: 78 │ Δt: 4.2844e-03 │ sim. time: 3.4228e-01 (68.457%) │ run time: 8.3789e+02 s
#timesteps: 79 │ Δt: 4.2833e-03 │ sim. time: 3.4657e-01 (69.313%) │ run time: 8.4629e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 80 run time: 8.55135744e+02 s
Δt: 4.28121331e-03 └── GC time: 6.71936197e+01 s (7.858%)
sim. time: 3.50847983e-01 (70.170%) time/DOF/rhs!: 6.04870721e-06 s
PID: 4.59076538e-05 s
#DOFs per field: 36128 alloc'd memory: 736.432 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.6288e+02 s
#timesteps: 82 │ Δt: 4.2769e-03 │ sim. time: 3.5940e-01 (71.881%) │ run time: 8.7061e+02 s
#timesteps: 83 │ Δt: 4.2748e-03 │ sim. time: 3.6368e-01 (72.736%) │ run time: 8.7926e+02 s
#timesteps: 84 │ Δt: 4.2729e-03 │ sim. time: 3.6795e-01 (73.590%) │ run time: 8.8694e+02 s
#timesteps: 85 │ Δt: 4.2709e-03 │ sim. time: 3.7222e-01 (74.445%) │ run time: 8.9562e+02 s
#timesteps: 86 │ Δt: 4.2691e-03 │ sim. time: 3.7649e-01 (75.298%) │ run time: 9.0324e+02 s
#timesteps: 87 │ Δt: 4.2674e-03 │ sim. time: 3.8076e-01 (76.152%) │ run time: 9.1190e+02 s
#timesteps: 88 │ Δt: 4.2657e-03 │ sim. time: 3.8502e-01 (77.005%) │ run time: 9.2031e+02 s
#timesteps: 89 │ Δt: 4.2641e-03 │ sim. time: 3.8929e-01 (77.858%) │ run time: 9.2813e+02 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 90 run time: 9.36698752e+02 s
Δt: 4.26264722e-03 └── GC time: 7.37447571e+01 s (7.873%)
sim. time: 3.93551528e-01 (78.710%) time/DOF/rhs!: 6.05977187e-06 s
PID: 4.50661953e-05 s
#DOFs per field: 36128 alloc'd memory: 971.258 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.4455e+02 s
#timesteps: 92 │ Δt: 4.2599e-03 │ sim. time: 4.0207e-01 (80.415%) │ run time: 9.5295e+02 s
#timesteps: 93 │ Δt: 4.2587e-03 │ sim. time: 4.0633e-01 (81.266%) │ run time: 9.6087e+02 s
#timesteps: 94 │ Δt: 4.2575e-03 │ sim. time: 4.1059e-01 (82.118%) │ run time: 9.6946e+02 s
#timesteps: 95 │ Δt: 4.2564e-03 │ sim. time: 4.1485e-01 (82.969%) │ run time: 9.7717e+02 s
#timesteps: 96 │ Δt: 4.2554e-03 │ sim. time: 4.1910e-01 (83.820%) │ run time: 9.8560e+02 s
#timesteps: 97 │ Δt: 4.2544e-03 │ sim. time: 4.2336e-01 (84.671%) │ run time: 9.9346e+02 s
#timesteps: 98 │ Δt: 4.2535e-03 │ sim. time: 4.2761e-01 (85.522%) │ run time: 1.0021e+03 s
#timesteps: 99 │ Δt: 4.2526e-03 │ sim. time: 4.3186e-01 (86.372%) │ run time: 1.0098e+03 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 100 run time: 1.01944787e+03 s
Δt: 4.25181873e-03 └── GC time: 8.12799139e+01 s (7.973%)
sim. time: 4.36113158e-01 (87.223%) time/DOF/rhs!: 6.14393635e-06 s
PID: 4.57228502e-05 s
#DOFs per field: 36128 alloc'd memory: 763.462 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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: 1.0264e+03 s
#timesteps: 102 │ Δt: 4.2503e-03 │ sim. time: 4.4461e-01 (88.923%) │ run time: 1.0351e+03 s
#timesteps: 103 │ Δt: 4.2496e-03 │ sim. time: 4.4886e-01 (89.773%) │ run time: 1.0428e+03 s
#timesteps: 104 │ Δt: 4.2490e-03 │ sim. time: 4.5311e-01 (90.623%) │ run time: 1.0514e+03 s
#timesteps: 105 │ Δt: 4.2484e-03 │ sim. time: 4.5736e-01 (91.472%) │ run time: 1.0591e+03 s
#timesteps: 106 │ Δt: 4.2478e-03 │ sim. time: 4.6161e-01 (92.322%) │ run time: 1.0676e+03 s
#timesteps: 107 │ Δt: 4.2473e-03 │ sim. time: 4.6586e-01 (93.171%) │ run time: 1.0754e+03 s
#timesteps: 108 │ Δt: 4.2468e-03 │ sim. time: 4.7010e-01 (94.021%) │ run time: 1.0840e+03 s
#timesteps: 109 │ Δt: 4.2463e-03 │ sim. time: 4.7435e-01 (94.870%) │ run time: 1.0917e+03 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 110 run time: 1.10138473e+03 s
Δt: 4.24591874e-03 └── GC time: 8.81322545e+01 s (8.002%)
sim. time: 4.78595613e-01 (95.719%) time/DOF/rhs!: 6.02825955e-06 s
PID: 4.52726496e-05 s
#DOFs per field: 36128 alloc'd memory: 706.872 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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.1083e+03 s
#timesteps: 112 │ Δt: 4.2453e-03 │ sim. time: 4.8709e-01 (97.417%) │ run time: 1.1168e+03 s
#timesteps: 113 │ Δt: 4.2450e-03 │ sim. time: 4.9133e-01 (98.266%) │ run time: 1.1246e+03 s
#timesteps: 114 │ Δt: 4.2448e-03 │ sim. time: 4.9558e-01 (99.115%) │ run time: 1.1322e+03 s
#timesteps: 115 │ Δt: 4.2447e-03 │ sim. time: 4.9982e-01 (99.964%) │ run time: 1.1409e+03 s
────────────────────────────────────────────────────────────────────────────────────────────────────
Simulation running 'CompressibleEulerEquations2D' with DGSEM(polydeg=3)
────────────────────────────────────────────────────────────────────────────────────────────────────
#timesteps: 116 run time: 1.14610234e+03 s
Δt: 1.79026594e-04 └── GC time: 9.11840317e+01 s (7.956%)
sim. time: 5.00000000e-01 (100.000%) time/DOF/rhs!: 5.96456573e-06 s
PID: 4.11142909e-05 s
#DOFs per field: 36128 alloc'd memory: 832.532 MiB
#elements: 2258 device memory: 0.000 MiB
Variable: rho rho_v1 rho_v2 rho_e_total
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: 1146s / 92.3% 224GiB / 86.4%
Section ncalls time %tot avg alloc %tot avg
─────────────────────────────────────────────────────────────────────────────────────────
parabolic rhs! 4.83k 1053s 99.6% 218ms 193GiB 99.8% 40.9MiB
prolong2interfaces 4.83k 996s 94.1% 206ms 190GiB 98.3% 40.3MiB
calculate gradient 4.83k 21.5s 2.0% 4.44ms 3.33KiB 0.0% 0.71B
volume integral 4.83k 10.2s 1.0% 2.12ms 0.00B 0.0% 0.00B
surface integral 4.83k 5.91s 0.6% 1.22ms 0.00B 0.0% 0.00B
interface flux 4.83k 1.72s 0.2% 356μs 0.00B 0.0% 0.00B
Jacobian 4.83k 1.72s 0.2% 355μs 0.00B 0.0% 0.00B
prolong2interfaces 4.83k 1.32s 0.1% 273μs 0.00B 0.0% 0.00B
reset gradients 4.83k 446ms 0.0% 92.4μs 0.00B 0.0% 0.00B
boundary flux 4.83k 38.8ms 0.0% 8.03μs 0.00B 0.0% 0.00B
prolong2boundaries 4.83k 30.0ms 0.0% 6.22μs 0.00B 0.0% 0.00B
~calculate gradient~ 4.83k 24.4ms 0.0% 5.05μs 3.33KiB 0.0% 0.71B
prolong2mortars 4.83k 2.72ms 0.0% 563ns 0.00B 0.0% 0.00B
mortar flux 4.83k 2.45ms 0.0% 507ns 0.00B 0.0% 0.00B
prolong2boundaries 4.83k 14.9s 1.4% 3.09ms 2.84GiB 1.5% 616KiB
volume integral 4.83k 8.36s 0.8% 1.73ms 0.00B 0.0% 0.00B
calculate parabolic fluxes 4.83k 7.16s 0.7% 1.48ms 0.00B 0.0% 0.00B
interface flux 4.83k 1.65s 0.2% 343μs 0.00B 0.0% 0.00B
transform variables 4.83k 1.54s 0.1% 319μs 0.00B 0.0% 0.00B
surface integral 4.83k 1.49s 0.1% 309μs 0.00B 0.0% 0.00B
Jacobian 4.83k 866ms 0.1% 179μs 0.00B 0.0% 0.00B
reset ∂u/∂t 4.83k 226ms 0.0% 46.7μs 0.00B 0.0% 0.00B
~parabolic rhs!~ 4.83k 63.7ms 0.0% 13.2μs 5.67KiB 0.0% 1.20B
boundary flux 4.83k 50.0ms 0.0% 10.3μs 0.00B 0.0% 0.00B
prolong2mortars 4.83k 14.2ms 0.0% 2.94μs 894KiB 0.0% 190B
mortar flux 4.83k 2.99ms 0.0% 618ns 0.00B 0.0% 0.00B
source terms parabolic 4.83k 580μs 0.0% 120ns 0.00B 0.0% 0.00B
analyze solution 13 2.18s 0.2% 168ms 420MiB 0.2% 32.3MiB
rhs! 580 1.96s 0.2% 3.37ms 4.78KiB 0.0% 8.44B
volume integral 580 969ms 0.1% 1.67ms 0.00B 0.0% 0.00B
interface flux 580 431ms 0.0% 744μs 0.00B 0.0% 0.00B
surface integral 580 210ms 0.0% 362μs 0.00B 0.0% 0.00B
prolong2interfaces 580 167ms 0.0% 287μs 0.00B 0.0% 0.00B
Jacobian 580 141ms 0.0% 243μs 0.00B 0.0% 0.00B
reset ∂u/∂t 580 18.9ms 0.0% 32.5μs 0.00B 0.0% 0.00B
boundary flux 580 12.5ms 0.0% 21.6μs 0.00B 0.0% 0.00B
prolong2boundaries 580 3.17ms 0.0% 5.46μs 0.00B 0.0% 0.00B
~rhs!~ 580 3.09ms 0.0% 5.33μs 4.78KiB 0.0% 8.44B
prolong2mortars 580 308μs 0.0% 532ns 0.00B 0.0% 0.00B
mortar flux 580 253μs 0.0% 436ns 0.00B 0.0% 0.00B
source terms 580 44.8μs 0.0% 77.3ns 0.00B 0.0% 0.00B
calculate dt 117 35.9ms 0.0% 307μs 16.5KiB 0.0% 144B
─────────────────────────────────────────────────────────────────────────────────────────This page was generated using Literate.jl.