#2931: Claude AI edits to TwoPuncturesX
Reporter: Steven R. Brandt
Status: new
Milestone:
Version:
Type: bug
Priority: major
Component:
Comment (by Steven R. Brandt):
I say we leave it as is and wait for the clankers to become smart enough to figure it out.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2931/claude-ai-edits-t…
#2929: GRHayLET/IllinoisGRMHD convert_IllinoisGRMHD_to_HydroBase schedule compatibility with VolumeIntegrals_*
Reporter: Maxwell Rizzo
Status: new
Milestone: ET_2026_05
Version: ET_2025_05
Type: bug
Priority: major
Component: EinsteinToolkit thorn
Comment (by Leonardo Rosa Werneck):
[PR #15](https://github.com/GRHayL/GRHayLET/pull/15) was merged. Please close this ticket if in fact issue was resolved.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2929/grhaylet-illinois…
#2936: CarpetX tsv output does not append files
Reporter: Maxwell Rizzo
Status: new
Milestone:
Version:
Type: enhancement
Priority: minor
Component: CarpetX
CarpetX’s TSV output, excluding the norms/ files, does not append to a single file per group in the way Carpet’s ASCII output does. Instead, it creates one file per group per iteration: `<group>.it000000.<x/y/z>.tsv` for grid functions and arrays, and `<group>.it000000.tsv` for scalars.
Appending subsequent iterations to the same file would make analysis and post-processing easier, while also reducing the amount of output files.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2936/carpetx-tsv-outpu…
#2916: Include Cottonmouth
Reporter: Beyhan Karakaş
Status: open
Milestone: ET_2026_05
Version:
Type: enhancement
Priority: major
Component: EinsteinToolkit thorn
Comment (by Zach Etienne):
Issues found (sorry for the uneven formatting):
```
Issue 1:
Missing symmetrization in the conformal Ricci tensor inside the constraints/analysis block.
Where the bug is
In `fun_bssn_cons`, the code uses
+ Delta[uc] * Gammat[la, lb, lc]
but in the BSSN Ricci decomposition the corresponding term is the symmetrized object
Delta^c * Gammat_(ab)c
which means
(1/2) * Delta^c * Gammat_abc
+ (1/2) * Delta^c * Gammat_bac.
Minimal fix
Replace
+ Delta[uc] * Gammat[la, lb, lc]
by
+ Rational(1, 2) * Delta[uc] * Gammat[la, lb, lc]
+ Rational(1, 2) * Delta[uc] * Gammat[lb, la, lc]
Note that your own `fun_bssn_rhs` block already contains the correct symmetrized pair and therefore serves as an internal consistency check.
```
2. bssnok.py
```
* Issue 2
- Severity high
- Description of the error
- The momentum-constraint matter source carries an extra factor of `w**2`. With `g_{ij} = (1 / w**2) * gt_{ij}`, one has `gamma^{ij} = w**2 * gt^{ij}`, but the conformal BSSN momentum constraint with free upper index uses `gt^{ij} * S_j`, not `gamma^{ij} * S_j`. The current code therefore overweights the source term by a factor of `w**2`.
- Filename + line number(s)
- Pasted file (filename not provided):638-648
- Original code snippet
- fun_bssn_cons.add_eqn(
MomCons[ua],
+ gt[ua, uc] * gt[ub, ud] * (
D(At[lc, ld], lb)
- Gammat[uk, lc, lb] * At[lk, ld]
- Gammat[uk, ld, lb] * At[lc, lk]
)
+ 6 * At[ua, ub] * cdphi[lb]
- Rational(2, 3) * gt[ua, ub] * D(trK, lb)
# Matter
- 8 * pi * w**2 * gt[ua, ub] * S[lb]
)
- Fixed code snippet.
- fun_bssn_cons.add_eqn(
MomCons[ua],
+ gt[ua, uc] * gt[ub, ud] * (
D(At[lc, ld], lb)
- Gammat[uk, lc, lb] * At[lk, ld]
- Gammat[uk, ld, lb] * At[lc, lk]
)
+ 6 * At[ua, ub] * cdphi[lb]
- Rational(2, 3) * gt[ua, ub] * D(trK, lb)
# Matter
- 8 * pi * gt[ua, ub] * S[lb]
)
- Confidence
- 96%
* Issue 3
- Severity low
- Description of the error
- The comments defining `Delta` are mathematically malformed. They describe the contraction with the wrong index pattern, which can mislead a reader into thinking `Delta^i` is built from `tilde gamma^{i,j} tilde Gamma^a_{ab}` or `tilde gamma^{ij} tilde Gamma^i_{jk}`. The correct definition is `Delta^i = tilde gamma^{jk} tilde Gamma^i_{jk}`. The executable code below the comments is correct; only the comments are wrong.
- Filename + line number(s)
- Pasted file (filename not provided):273-275
- Original code snippet
- # \tilde{\gamma}^{i, j} \tilde{\Gamma}^a_{a b}
# When \tilde{\Gamma}^{i} when it appears and its derivative are not needed,
# the substitution \tilde{\Gamma}^{i} \rightarrow \tilde{gamma}^{ij} \tilde{\Gamma}^{i}_{jk} = \Delta^i
- Fixed code snippet.
- # \Delta^i = \tilde{\gamma}^{jk} \tilde{\Gamma}^i_{jk}
# When \tilde{\Gamma}^{i} appears without derivatives, we replace it by
# \Delta^i = \tilde{\gamma}^{jk} \tilde{\Gamma}^{i}_{jk}
- Confidence
- 99%
```
```
Issue 4:
linear_wave.py:
File sets an unphysical default amplitude for the wave:
amplitude = cottonmouth_linear_wave_id.add_param(
"amplitude",
default=1.0,
desc="Linear wave amplitude"
)
A << 1 for a valid linearized wave.
```
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2916/include-cottonmou…
#2928: New Simfactory Option
Reporter: Steven R. Brandt
Status: new
Milestone:
Version:
Type: enhancement
Priority: major
Component: SimFactory
Changes (by Roland Haas):
responsible: [] (was )
assignee: Steven R. Brandt (was )
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2928/new-simfactory-op…
#2935: Kuibit: Support CarpetX tsv output (1D GF/Array, Norm/Scalar)
Reporter: Maxwell Rizzo
Status: new
Milestone:
Version:
Type: proposal
Priority: minor
Component:
Changes (by Maxwell Rizzo):
See issue on Kuibit's repository: https://github.com/Sbozzolo/kuibit/issues/49
To summarize the file output differences:
Carpet text output was either per variable or per group `.asc` ASCII files.
CarpetX text output currently is only per group `.tsv` TSV files.
CarpetX only appends the output to the `norms/` files, all other output files (scalars, GFs, Arrays) get their own file per group per iteration (per axis where applicable).
Lastly, the headers have changed, the CarpetX headers are a single line, as opposed to the multi-line Carpet headers.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2935/kuibit-support-ca…
#2935: Kuibit: Support CarpetX tsv output (1D GF/Array, Norm/Scalar)
Reporter: Maxwell Rizzo
Status: new
Milestone:
Version:
Type: proposal
Priority: minor
Component:
See issue on Kuibit's repository: https://github.com/Sbozzolo/kuibit/issues/49
To summarize the file output differences: Carpet text output was either per variable or per group `.asc` ASCII files. CarpetX text output currently is only per group `.tsv` TSV files. CarpetX only appends the output to the `norms/` files, all other output files (scalars, GFs, Arrays) get their own file per group per iteration (per axis where applicable). Lastly, the headers have changed, the CarpetX headers are a single line, as opposed to the multi-line Carpet headers.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2935/kuibit-support-ca…
#2175: Test "Single, stable neutron star" example
Reporter: Roland Haas
Status: open
Milestone: ET_2026_05
Version: development version
Type: task
Priority: major
Component: EinsteinToolkit website
Comment (by Beyhan Karakaş):
For the testing in future please check the "tov_20260325.tar.gz" and upload only the required files for each run which are less than a MB in total.
--
Ticket URL: https://bitbucket.org/einsteintoolkit/tickets/issues/2175/test-single-stabl…