ML demo: a U-Net in SQL, for einfold to optimize convolutions against

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

メンテナーはふだん 3 日以内に返信

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

評価

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

調査の方向性

Start with benchmarks/ml/CONTRIBUTING.md and compare the structure of 01_attention.py and 02_nanogpt.py. Build a readable U-Net demo with a JAX reference, timed SQL steps across hosts in _hosts.py, and checks for forward results and input and weight gradients; completion means those comparisons agree to about 1e-15. einfold's shifted-key optimizer work is tracked separately, and results belong on #285.

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

説明

enhancement ml

🪭 Claude (einfold)

A request for a new demo in benchmarks/ml, following its CONTRIBUTING.md, which already lists "a small convolutional network over a grid, where a convolution is a join on shifted coordinates" among the demos that would fit. The proposal: a U-Net (Ronneberger, Fischer and Brox, U-Net: Convolutional Networks for Biomedical Image Segmentation, 2015), the encoder–decoder CNN with skip connections that geoscience uses for downscaling, nowcasting and segmentation.

einfold, XQL Systems' query optimizer that rewrites SQL plans into array algorithms, would use it the way it uses 01_attention.py: as the workload it optimizes against.

Why this model

  • It is the convolution shape einfold doesn't detect yet. A convolution is a fold over a join whose keys are shifted, i.y = o.y + w.dy, not equal columns. einfold's next design step for it is shifted keys in its fold IR, so that the dense kernels take convolutions directly. Strides and dilations become a multiplier on the offset. Padding needs no special case: border pixels simply have fewer taps, which is SQL's own result when the shifted rows don't exist. A demo gives that work a correctness reference (JAX) and a benchmark, as attention did for einfold's fused pipelines.
  • It exercises more than convolutions.
    • Convolution, then bias, then ReLU, is a fold pipeline, which einfold's M2 fuses.
    • Downsampling (a strided convolution, or a max-pool, which is MAX over shifted keys) and upsampling (a transposed convolution) are the same join with o.y = 2·i.y + dy and its inverse.
    • Skip connections concatenate channels from two levels: a UNION ALL, or a join on the pixel.
  • It fits xarray-sql's audience. A U-Net on an ERA5 field, such as downscaling a coarse temperature grid or predicting a field one step ahead, is a model geoscientists train on gridded data. Its input is an xarray Dataset registered as a table of (time, y, x, channel).
  • ddx gets new gradients to check. The gradient of a convolution with respect to its input is a transposed convolution, and with respect to its weights is a correlation of the input with the output gradient: the same shifted-key shape again.

Suggested shape

The pieces, as 02_nanogpt.py builds its model from functions named after the reference's:

  • conv(x, w, stride, dilation): SUM(x.v * w.v) over x JOIN w ON x.c = w.ci AND x.y = o.y * stride + w.dy * dilation AND …, grouped by the output's (n, y, x, co). One way to write o is a join against a coordinate table, or x.y - w.dy grouped directly;
  • relu, bias: derived factors;
  • down(x): a stride-2 convolution, or a 2 × 2 max-pool;
  • up(x): a stride-2 transposed convolution, or nearest-neighbour JOIN ON o.y / 2 = i.y;
  • skip(a, b): channels of two levels, concatenated;
  • unet(depth, width): two or three levels, 8–32 channels, 3 × 3 kernels, on a grid of about 64 × 64 to 256 × 256.

Checked against the same network in JAX: jax.lax.conv_general_dilated, or Flax's nn.Conv, with 64-bit floats, so they agree to about 1e-15. Check the forward pass, and ddx's gradients with respect to the input and to every weight table.

The interesting knobs for einfold's measurements:

  • grid size;
  • channel widths, which make the convolutions dense or near-dense;
  • depth, which is the number of folds in a pipeline;
  • stride, against pooling.

As in the other demos: the same SQL on every host in _hosts.py, each step timed, with its error against JAX.

Not asked here

Speed. The demo should be written to be read, as the others are. einfold's side, detecting and running shifted-key folds, will be tracked in einfold. Results go on #285.

主要言語
Python
スター
142
フォーク
19
平均マージ
3日 47分
マージ済み PR（30日）
8

環境構築

はじめの一歩

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

xqlsystems/xarray-sql のほかの issue

xqlsystems/xarray-sql の issue をすべて見る

似ている issue

Python の issue をもっと見る

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

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