Hacktoberfest 2026：メンテナが10月に向けて印を付けた、オープンで初心者向けの issue。 Hacktoberfest の issue を見る

meta(sciml): expose prepared operators to the SciML integrator stack

オープン
  • 何か問題がありますか？
#81 コメント 1 件 リアクション 0 件 担当者 0 名 GitHub で見る

まだ誰も着手していません。

評価

難易度
5/5
見積もり時間
1週間以上
初心者へのやさしさ
25/100
issue の種類
機能追加
明瞭さ
おおむね明確
活発さ
活発
技術スタック
julia
領域
backend-api-design

調査の方向性

Treat this as an umbrella rather than a first patch. Read README.md, CLAUDE.md, and DESIGN.md §7, then inspect sub-issues #86, #65, #73, and #87. It is done when the ordered extension work and its constraints have been resolved without adding SciML packages to [deps].

索引モデルが issue の本文から書いたものです。

説明

enhancement

Umbrella for handing MFO operators to OrdinaryDiffEq / LinearSolve. Nothing SciML exists today: no weak dep, no ext/…SciMLExt.jl, and every example hand-rolls its integrator. README.md no longer advertises OrdinaryDiffEq interop that is not there (#87 retracted it; README.md:150 now says there is no integration yet).

Status (2026-09-30): #86 and #87 are done. What remains is #65 first, then #73.

Sub-issues, in order

  1. #86 done (PR #89) — prepared-operator trait hygiene — PreparedOperator has no supertype and PreparedComposed / PreparedAdjoint declare no islinear / isconstant, so a wrapper forwarding traits gets the false default by accident. An invariant hazard independent of SciML; ship first.
  2. #65 — next — the jac_prototype adapter (DESIGN §7). Enough that ODEFunction(f; jac_prototype = sciml(prepare(L))) works and KrylovJL_GMRES calls the zero-alloc flat mul!. The downstream forcing function is the microvascular model's KenCarp47 + sparse-ILU detour.
  3. #73 — after #65 — SplitODEFunction for the linear half plus a per-model work-precision search over IMEX / exponential integrators (a benchmark artifact, not a package feature).
  4. #87 done (PR #90) — docs — fix the README claim, and document field-level apply! as the explicit-stepping idiom (the prepared mul! leg pays 15–50% in flat copies that only Krylov needs).

Related, not parented: #63 (stacked multi-species operator) is what the reaction half of a split problem wants to look like; #62 (Robin) is what makes MultigridPreconditioner applicable to the downstream wall-coupled system.

Addendum questions (from the comment below)

  • Q1 (what OrdinaryDiffEq requires of a state u, and the minimal interface that makes Field a legal state) now lives in #73's checklist.
  • Q3 (ghosts inline in a padded array vs. a flat interior with ghosts beside it in face buffers) is a DESIGN question, not a SciML one, and now lives as a checkbox in #79.
  • Q2 (whether a Field-valued Krylov workspace beats the two staging copies per matvec, CPU and GPU separately) is still an open benchmark question here.

Constraints

  • Extension only. SciMLOperators / OrdinaryDiffEq are weak deps, nothing enters [deps] (CLAUDE.md, DESIGN.md §7).
  • Do not route around the mul!(y, ::AbstractOperator, x::AbstractVector) guard; a convenience constructor must still go through prepare.
  • The affine boundary_rhs term needs a declared home in a split problem.

—
🤖 Claude tidied this one and moved the addendum questions to their proper homes.
claude-opus-5-5[1m] (high) · claude-code 2.1.285 · unreviewed

主要言語
Julia
スター
3
フォーク
1
平均マージ
1日 28分
マージ済み PR（30日）
6

環境構築

このプロジェクトには開発コンテナ、Dockerfile、コントリビューションガイドがありません。まず README を読み、一般的な手順ははじめてのコントリビューションガイドを参照してください。

はじめの一歩

  1. issue を最後まで読み、次にプロジェクトのコントリビューションガイドを読みます。
  2. 着手することを issue にコメントします — 二人が同じ作業をするのを防げます。
  3. リポジトリをフォークし、ブランチを切って変更します。
  4. issue 番号を参照したプルリクエストを送ります。

RallypointOne/MatrixFreeOperators.jl のほかの issue

RallypointOne/MatrixFreeOperators.jl の issue をすべて見る

似ている issue

Julia の issue をもっと見る

新しい issue をメールで受け取る

初心者向けの GitHub issue を短くまとめたダイジェスト。