Skip to content

Conversation

@melilly
Copy link
Owner

@melilly melilly commented Jun 3, 2025

No description provided.

zhichen3 and others added 30 commits March 4, 2025 19:54
also remove the old version of the subch inputs
Temperature was 0, causing tempi to be infinity during the eos call.
I'm able to run wdmerger to the step # 3.
cellsize is not set at this point. So manually set it here.
…d of amrex (#3048)

This pr adds explicit evaluation of the thermal diffusion source term in castro, without needing to use amrex.

In particular, this seems to work well for 2D spherical geometry.
…t test (#3047)

this pr adds problem_bc_fill.H, which allows to use inflow, i.e. Dirichlet, boundary condition by using the analytic solution for different problems.
Even though this is intended to use for spherical 2d input file, but any input file theoretically can use this as long as its constant conductivity.
we had a ** due to line wrapping
* dataset processing can run in parallel

* added code to make interactive plots of each dataset that show the raw slice data, the averaged profile, and the positions for each metric
  
* metric thresholds are calculated using both the local maximum and global maximum at the same time (under <field>[<fac>%] and <field>[<fac>%]_gmax, respectively)
  
* the local maximum for each field is stored in the output file, in the max_<field> column
  
* if a field value never goes below the threshold for a metric, its position will be stored as NaN
  
* --[xyz]lim arguments now work, despite some yt bugs
  
* allow ignoring plotfiles after a certain time with --tmax
  
* default FRB resolution is set to match the finest level, and is automatically adjusted if custom limits are used
…ffusion test (#3054)

Updates problem_bc_fill to do either face-centered or cell centered depending on where
castro.diffuse_use_amrex_mlmg is used.

Gets 2nd order accuracy (with L-inf error) with max_level=0 for 2d spherica

To test theta dependence, the analytic solution of spherical-bessel function and legendre polynomial is used.
Currently only castro evaluation works.

Also added an 2d spherical input file for noncenter Gaussian, a convergence test for this still needs to be done in the future.
if we mistakenly target "main" and then switch to "development",
then check-ifdefs, codespell, and docs never run.  This changes
the action to look like the others so these should hopefully run.
some of the inputs were redudant
also renamed a few
updated the README.md
finally added a dt option to the profiles.py script
we fill this in Microphysics, but it was not being captured in job_info
* 1D inputs and new model

* update README

* update 1d lower boundary
r is the reduction parameter but we were also using it
in the reduction loop for the radius (for 2D spherical)
…3061)

this shows the mass fractions behind the detonation
we can't to .ptp() anymore, but instead should use np.ptp()
some loop counters were also uninitialized.
this also fixes a few compiler warnings.
this was accidentally deleted in #3068
also remove some old inputs files
zhichen3 and others added 30 commits October 24, 2025 15:43
This was only an issue if we were doing 4th order and we were initializing the data as cell-averages directly.
Bumps [actions/upload-artifact](https://github.com/actions/upload-artifact) from 4 to 5.
- [Release notes](https://github.com/actions/upload-artifact/releases)
- [Commits](actions/upload-artifact@v4...v5)

---
updated-dependencies:
- dependency-name: actions/upload-artifact
  dependency-version: '5'
  dependency-type: direct:production
  update-type: version-update:semver-major
...

Signed-off-by: dependabot[bot] <support@github.com>
Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com>
this adds Ryan's RNAAS to the section on shock burning
this now allows Ne20 and also the C/O ratio can be something other
than 1/2.  Each core X is specified separately now.
this uses a more robust solve for NSE
This reads the data from the location of peak enuc and uses pynucastro to visualize
the network flow
Bumps [actions/checkout](https://github.com/actions/checkout) from 5 to 6.
- [Release notes](https://github.com/actions/checkout/releases)
- [Changelog](https://github.com/actions/checkout/blob/main/CHANGELOG.md)
- [Commits](actions/checkout@v5...v6)

---
updated-dependencies:
- dependency-name: actions/checkout
  dependency-version: '6'
  dependency-type: direct:production
  update-type: version-update:semver-major
...

Signed-off-by: dependabot[bot] <support@github.com>
Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com>
this small function can take 10% of the runtime of a simulation
Bumps [actions/upload-artifact](https://github.com/actions/upload-artifact) from 5 to 6.
- [Release notes](https://github.com/actions/upload-artifact/releases)
- [Commits](actions/upload-artifact@v5...v6)

---
updated-dependencies:
- dependency-name: actions/upload-artifact
  dependency-version: '6'
  dependency-type: direct:production
  update-type: version-update:semver-major
...

Signed-off-by: dependabot[bot] <support@github.com>
Co-authored-by: dependabot[bot] <49699333+dependabot[bot]@users.noreply.github.com>
this gives problems on some machines that don't have a python3 command
(Castro.H -> timestep.H -> Rotation.H -> Castro.H)
this uses a better buffer size, adds the "ash" field, and cleans
up some sizing
this speeds up copy_cons_to_burn by about 10x.
this uses Array4.cellData() to reducing some indexing
and uses a constexpr_for
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.

8 participants