Skip to content

Dust size evolution updates#13

Open
calebchoban wants to merge 111 commits into
pfhopkins:gizmo-cppfrom
calebchoban:dust_evo_updates
Open

Dust size evolution updates#13
calebchoban wants to merge 111 commits into
pfhopkins:gizmo-cppfrom
calebchoban:dust_evo_updates

Conversation

@calebchoban

Copy link
Copy Markdown
Contributor

Updates to grain size evolution model from the previous merge.

  • Simplified Config options, removing unused models and outputs.
  • Speed ups to coagulation/shattering

mikegrudic and others added 30 commits July 22, 2025 11:13
fix: zero-initialize all bytes in  at startup to mitigate known issue…
remove jaco files from build system until the codegen is fully automated
…e dust temperature values in setups with COOLING but without RT_INFRARED
jaco build system temporary fix
pfhopkins and others added 29 commits March 27, 2026 08:40
…timestep criterion to not bother with stars on large timesteps since not needed
Adds ~30 new test suites and updates test harness infrastructure.

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
…ead of max_time_bin_active

Instead of always bumping woken particles to the highest active timebin,
use the waker's own timebin minus a fixed offset (ceil(log2(WAKEUP)) bins)
so woken particles get a timestep proportional to dt_waker / WAKEUP.

Encodes waker's TimeBin+1 in P[j].wakeup for hydro wakeups (>0), and
uses -1 for non-hydro wakeups (sinks, merge/split, AGS, SIDM) which
retain the old max_time_bin_active behavior.

Co-Authored-By: Claude Opus 4.6 (1M context) <noreply@anthropic.com>
…truction and a more accurate Riemann solver. This is not yet fully debugged, but preliminary tests suggest it is working well. I will continue to test and debug this in the coming days, but I wanted to share the code in case anyone else wants to try it out or provide feedback.
pulling in just non-invasive testing infrastructure developments from development branch
This reverts commit 7f8e841.
PR to make a couple minor tweaks i accidentally coded through here, and one test, compiler-flag-locked, for 2nd-order reconstruction of RHD fluxes, which I am testing. can be deleted as needed and will not be on by default
… still single-cell-and-neighbor based, with more robust mach estimation, using OUTPUT_SHOCK_MACH_NUMBER flag. note currently doesnt work with EOS_GENERAL because of assumption of hydrodynamic shock jump conditions. need to free this to get an approximate estimator, even, because ideal use cases do heavily involve EOS_GENERAL cases like CRs
targeted wakeup timebin: bump woken particles to dt_waker/WAKEUP instead of max_time_bin_active
…tions and added precomputing polynomial coefficients.
…lified and old Elemental model has been removed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants