Skip to content

examples: correct the Heun claim and list rfm_fisher_simplex - #16

Merged
arclabs561 merged 2 commits into
mainfrom
examples/accuracy-20261103
Nov 4, 2026
Merged

arclabs561 merged 2 commits into
mainfrom
examples/accuracy-20261103

Conversation

@arclabs561

@arclabs561 arclabs561 commented Nov 4, 2026 •

Copy link
Copy Markdown
Owner

Two small corrections to the examples.

  • ode_comparison ended with "Heun's radius stays near 1.0 (symplectic-like behavior)". Its own output shows Heun's radius at 1.2106 after 10 steps and 1.0246 after 20, growing every step. Heun is not symplectic. The closing text now says both methods gain radius, about h^2/2 per step for Euler and h^4/8 for Heun.
  • rfm_fisher_simplex is registered under the riemannian feature and run by CI but was missing from the examples README. It is now listed. fynch_bregman_energy printed its inputs with ndarray's debug format (shape=[3], strides=[1], layout=CFcf ...); it now prints plain lists.

Checked: cargo fmt --check, cargo clippy --release -D warnings on both changed examples, and both examples run with exit 0.

Henry Wallace added 2 commits November 3, 2026 22:02
The closing text said Heun's radius "stays near 1.0 (symplectic-like
behavior)". The example's own output contradicts it: Heun's radius is 1.2106
at 10 steps and 1.0246 at 20, and it grows every step. Heun is not symplectic;
on dx/dt = -y, dy/dt = x its per-step radius factor is sqrt(1 + h^4/4) against
Euler's sqrt(1 + h^2), so it gains radius about h^2/4 times slower.

The text now says both methods gain radius and gives the two rates. Run with
`cargo run --release --example ode_comparison`; clippy -D warnings and fmt
pass.

Claude-Session: https://claude.ai/code/session_01UJZpXYJLkeY3RtMaiJnqSR
…ts as lists

rfm_fisher_simplex is registered with the `riemannian` feature and run by CI,
but the examples README did not list it. Add it to the feature table and the
run block.

fynch_bregman_energy printed its inputs with `{:?}` on an ndarray, which emits
`shape=[3], strides=[1], layout=CFcf (0xf), const ndim=1`. Print the values
as a plain list instead.

Claude-Session: https://claude.ai/code/session_01UJZpXYJLkeY3RtMaiJnqSR
@arclabs561
arclabs561 merged commit 89ae9ff into main Nov 4, 2026
2 checks passed
@arclabs561
arclabs561 deleted the examples/accuracy-20261103 branch November 4, 2026 07:15
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.

1 participant