- Dec 30, 2023
-
-
Florian Atteneder authored
* update ref test data srhd_simple_wave_avmda * fix up dispatch of register_variables; fix up smoothing on initial data * mda: use inplace operations * cons2prim: disable simd; disable atmosphering warning * update blastwave1/2 av mda examples * adjust tau when energy density exceeds lower bound * don't recompute max speed for timestep * determine initial viscosity using all dynamic variables as indicators * add blastwave1 example * add blastwave2 examples * fixup blastwave2 test; don't smooth FVElement mesh initial data * commit conservative_fixing experiments (not active) * cons2prim: temporarily error in case speed is too large * disable conservative_fixing again * change fixedprms in setup * fix domain assertion for blastwave1 test * add initial data smoother based on Bernsteining * do not recompute viscosity on initial data * modify examples * implement fd derivative again * implement covariant regularization * add prms to Project type * add example parfiles * update ref test data srhd_simple_wave_fv the difference is quite pronounced * use maximum of max_v instead of v for time step estimation * update ref test data srhd_simple_wave_avmda * compute cellmax_v properly; compute maxspeed in equation callback * implement sod shock tube initial data
-
- Dec 24, 2023
-
-
Florian Atteneder authored
* insert line number nodes for better stack traces
-
- Dec 19, 2023
-
-
Florian Atteneder authored
* add ref test srhd_simple_wave_avmda * undo some code deletion * special case FV method on entropy production computation * init EP for all setups * add entropy tci computation that can work with precomputed entropy production * add entropy viscosity computation that can work with precomputed entropy production * compute entropy production in equation callback * add ref test srhd_simple_wave_fv reset ref data for srhd_simple_wave_fv
-
Florian Atteneder authored
* update ref test increase resolution ref test increase tol for ref test and lower resolution * add ref test srhd_simple_wave_fv * implement FV version for SRHD eqs * simplify timestep computation * initial data setup functions need to take a mesh
-
- Dec 18, 2023
-
-
Florian Atteneder authored
* fixup speed computation burgers eq * fix upload logic for diffoutput uploading * update ref test data * migrate all implementations to use new broadcasting interface * evict broadcast_volume_2!
-
Florian Atteneder authored
* remove cons2prim unit test * add srhd_simple_wave.toml * update ref test burger_sine_aventropy * undo experiments in ev.jl * undo changes to advection_sine ref test * fix up entropy production computation * fix entropy variables computation according to formulae from EFL and EHL papers * WIP * make AV regularization run * disable conservative fixing as it causes wrong results in srhd_smooth test; update ref test data * evict old broadcasting interface * define entropy variables for polytrope * compute fluxes after interpolating bdry data * optimize mda computation * fix ramping function for mda * enable conservative_fixing every to make ref tests pass -- however, ref test restults don't look correct update srhd_smooth ref test data * move cons2prim in rhs down * use only D as HRSC indicator * fix simple wave initial data and disable conservative fixing as it kills the simple wave propagation * make bdry conditions work for shock tube test * document print_idxs for @with_signature * add mda_n_low, mda_n_high parameters that allow to control the truncation interval for the MDA heuristic * fix logic in combining mda indicators from multiple vars * enable TCIs for SRHD * SRHD: allow to compute MDA AV from multiple vars and combine them with a max function * SRHD: make srhd_smooth test run with AV (though not perfect) * SRHD: add ldg_nflux * update ref test srhd_smooth * add more srhd parfiles * split conservative fixing from cons2prim; fixup cons2prim so that symmetric shock tubes are actually symmetric * update 0th order smoother * compute viscosity for initial data * use mix of conserved and primitives for flux computation * cleanup equations * cleanup cons2prim * implement av entropy for 1d eq -- wip * utilize kastaun cons2prim in llf * remove cons2prim_insane and friends * fix type instability in cons2prim caused by infamous closure bug fixup cons2prim * add Kastaun's cons2prim with conservative fixing * rename variables * remove comment * convert llf to use new broadcasting interface * convert cons2prim to new broadcasting * convert entropy_variables, but it is todo * remove unused * convert flux to new broadcast interface * convert speed to new broadcast interface * clean up module file
-
- Dec 14, 2023
-
-
Florian Atteneder authored
* update euler_isentropic_flow ref test to use interpolate_aligned_ts * reorganize testconfig.toml * update advection_sine_glgl ref test to use interpolate_aligned_ts * update advection2d_bump_threshold ref test to use aligned_ts * remove unused toml files * fix testparfiles to not list global options * update heateq_sine ref test to use aligned_ts * update burgers2d_gaussian_entropyav ref test to use aligned_ts * update burgers2d_gaussian ref test to use aligned_ts * update srhd_smooth ref test to use aligned_ts * update euler_isentropic_flow ref test to use aligned_ts * update advection2d_sine ref test to use aligned_ts * update fv_euler_isentropic_flow ref test to use aligned_ts update output.h5
-
Florian Atteneder authored
* fixup success check * add interpolate.h5 output for heateq_gaussian ref test * use aligned_ts for heateq_gaussian ref test * update ref test advection_sine to include interpolate.h5 output * upgrade IntegrationTests to also compare interpolate.h5 outputs (if available)
-
Florian Atteneder authored
* remove whitespaces * use mktempdir instead of custom mockdir gymnastics * trigger a TimeAlignedCB also when t+dt == next_t
-
Florian Atteneder authored
* InterpolationCallback can only interpolate fields that are registered as dynamic, static or rhs * add check for registered variable in get_variable * add a few @inline decorations * add isregistered method to cache and variablesdict types * throw proper error message when unregistered variable is requested * add test for InterpolationCallback * fixup callback interface docs * enable interpolation options of Output parameter section * add parameters for interpolation output * fixup InterpolatioCallback * Invocation of CallbackSet should report the call status I think initially the return value of all callbacks was meant to indicate whether that callback was actually called. I think what we need now is that callbacks return their status, e.g. whether they failed or not. * remove superfluous Base.show overloads * add InterpolationCallback use refactored meshinterpolator * introduce AbstractSaveCallback
-
Florian Atteneder authored
* update tests * evict old plot functionality evict remainders
-
Florian Atteneder authored
* update dg1dTests manifest * refactor MeshInterpolator to accept arrays for interpolation * use ranges to construct the cells' extends, because they ensure that the endpoints are preserved
-
Florian Atteneder authored
* add sod_shock_tube examples * add examples mono and navier-stokes regularization of euler_isentropic_flow problem * remove dead code * update ref test burgers_sine_aventropy * undo deletion * remove dead code * beautify log details when unequal number of output times is encountered * refactor printing of failed test details * im in a rush * fix ref test log info * remove unused var * revert callbacks * add invalid assertion in glgl integration method * remove bernsteining from AV regularization * add 0th order smoother * remove lower threshold on entropy av * wip -- tinker with applying Bernsteining to the LDG variables * also forward project to rhs * remove unnecessary computations from callback * put lower threshold on entropy AV * inline bernsteining workhorse * use inplace ops for entropy TCI * add ref test fv_euler_isentropic_flow * fixup hrsc initialization for EulerEq * fix args for Project constructor * update deps; put compat bound on PrettyTables.jl * add ref test euler_isentropic_flow_regulate-ns -- disabled for now * add dispatch for av regularization in rhs!; wip implementation of new regularizations * add a namedtuple for parameters to Project struct * fix call to broadcast in callbacks.jl * add av_regularization parameter * dispatch TODO on Symbol
-
- Dec 11, 2023
-
-
Florian Atteneder authored
* IntegrationTests: print Inf norm if mismatch is found * update refdata burgers_sine_avmda * update refdata burgers_sine_avdb * update refdata burgers_sine_aventropy * update refdata burgers_sine_avtci * fixup constructor call * more examples * update output times; fixup gitignore * utilize Sandbox to parse expressions from parfiles * fixup ScalarEq bump initial data * add examples folder with parfiles that are not ref tested * specialize timestep for AV scheme * implement 1st order fd approx for derivative of regularization term * enable callbacks * add prms field to ScalarEq's project struct also add av_derivative_scheme parameter to ScalarEq project to be used subsequentely * simplify dispatch of rhs!
-
- Dec 06, 2023
-
-
Florian Atteneder authored
* update some ref tests - euler_isentropic_flow - euler_sod_shock_tube_bernstein needed because previously we employed a timestep that was too large, because it only used the local cell width but not the smallest grid spacing * add ref test for FV euler sod shock tube * update ScalarEq so that FV version can use non-periodic bdry conditions * implement FV method for Euler equation - also refactor timestep to work with DG and FV elements - also refactored fv_update_step! methods to now accept buffers for bdry values
-
Florian Atteneder authored
* remove more Riemann solver stuff * remove non-existent export * evict riemannsolver.jl and broadcast_[volume,faces,bdry]! * evict Riemann solver related stuff from SRHD * evict Riemann solver related stuff from EulerEq * evict Riemann solver related stuff from ScalarEq * convert all in EulerEq remove unrelated file * convert all in ScalarEq * re-enable srhd_smooth test update srhd_smooth reftest * evict old broadcasting methods
-
- Dec 05, 2023
-
-
Florian Atteneder authored
Add a finite volume method for ScalarEq See merge request !75
-
Florian Atteneder authored
* add reftest fv_advection_sine * migrate all of ScalarEq to use broadcast_*_2! methods * add isperiodic method for meshes * fix grid spacing computation for SpectralElement * fixup parametric types to subtype Mesh1d/2d * make TODO error more comprehensive * fixup cfl usage; add abstractions for min grid spacing * enable add mesh kwarg to Evolution setup in order to enable FV evolution * add timestep! and rhs! for FVElement * do away with the State(...) syntax * add docs to dg_rhs methods * add FV stepper * limit dg_rhs.jl methods to Meshes with SpectralElement * add mesh kwarg to Evolution so that we can use a spearate stepper based on mesh type * add a scheme parameter to Mesh * add FVElement type
-
Florian Atteneder authored
Refactor CFL parameter usage See merge request !74
-
Florian Atteneder authored
* fix tests * remove an unreachable check can't be reached, because a error would be thrown earlier when multiplying a number with nothing * remove cfl field from Evolution struct * move cfl parameter from Mesh to Evolution section also update all reftest parfiles * remove CFL field from Mesh * apply CFL factor in evolve loop and not inside projects
-
Florian Atteneder authored
call the right interpolator for HRSC/smoother.jl See merge request !73
-
Florian Atteneder authored
* call the right interpolator
-
- Dec 04, 2023
-
-
Florian Atteneder authored
Allow dispatch on second argument See merge request !72
-
Florian Atteneder authored
* enable mesh argument in broadcasting methods * remove debug statement * fix docs * add unit test * disable euler_sod_shock_tube_entropyav ref test * fixup tests * enable mesh as second arg for legacy mode
-
Florian Atteneder authored
make main() accept numbers which correspond to testparfiles() list See merge request !70
-
Florian Atteneder authored
* make main() accept numbers which correspond to testparfiles() list
-
- Nov 30, 2023
-
-
Florian Atteneder authored
bump version See merge request !69
-
Florian Atteneder authored
* bump version
-
Florian Atteneder authored
Cleanup See merge request !68
-
Florian Atteneder authored
* rm deps BenchmarkTools, GLob, InteractiveUitls, LaTeXStrings, REPL; update dg1dTests' dependencies * remove Plots, PlotThemes, Interpolations
-
Florian Atteneder authored
improve a few things in the `make_movie.jl` script See merge request !63
-
Florian Atteneder authored
* improve a few things - avoid axis jittering by using the Inside() alignmode for axes - allow to plot cell data as constant functions - update manifest for v1.9
-
Florian Atteneder authored
Improve `HeatEq` See merge request !65
-
Florian Atteneder authored
* add heateq_sine reftest * rename reftest from heateq to heateq_gaussian * rename numflux methods * change time step heuristic to align with what nodal DG book suggests (p. 250) * add sine initial data for heateq
-
- Nov 29, 2023
-
-
Florian Atteneder authored
Stepper for time evolution See merge request !67
-
Florian Atteneder authored
* fix evolve! by wrapping stepper result also in callback previously we copied the stepper's result into the initial data buffer but now we just swap the vectors we now push a wrapper for the dynamic state vector (like we also use it for the rhs_fn, timestep_fn) to the fullstep callback so that we don't have to do the wrapping in projects * teach pushfirst to CallbackSet * avoid adding cfl twice * disable the timestep-retaking tests as they are now deprecated * grab cfl parameter from mesh parameters * fix while condition to be aware of symbol return values * fix up evolve test * fix time stepping with Evolution type - make a copy of initial data to initialize up1, because we need to allocate memory anyways and the stepper will also swap it before stepping - apply cfl to timestep * fix construction of Evolution type - fix sanity checks for callbacks - add sanity check for cfl - introduce TimeStepAlgorithm super type and move all algorithm types into types.jl - add docstring to Evolution type - update tests * update unit tests * introduce a state full type Evolution so that we can evolve step by step
-
- Nov 27, 2023
-
-
Florian Atteneder authored
Fix CI See merge request !66
-
Florian Atteneder authored
* fix undef var * fix vector signature * try to capture and log full backtrace in ref tests * revert some changes * ignore v2.2.3 in PrettyTables in IntegrationTests * what is the real error we see? * wrap parameter printing into try-finally, because PrettyTables.jl causes the CI to crash here * also upload errors.log * log the captured error to a file * also resolve before running ref tests
-
- Nov 16, 2023
-
-
Florian Atteneder authored
convert error to warning when unknown parameters are found See merge request !64
-
Florian Atteneder authored
* convert error to warning when unknown parameters are found
-