Compare commits

...

32 Commits

Author SHA1 Message Date
de3da1f6a1 [obsidian] vault backup: 2025-08-01 13:05:00[
Some checks failed
Build / build (push) Failing after 14m32s
2025-08-01 13:05:00 +09:00
3d7c9c9048 [obsidian] vault backup: 2025-07-31 19:25:07[
Some checks failed
Build / build (push) Failing after 14m32s
2025-07-31 19:25:07 +09:00
b32ef1d348 [obsidian] vault backup: 2025-07-31 18:24:57[
Some checks failed
Build / build (push) Failing after 14m44s
2025-07-31 18:24:57 +09:00
bb453aba6b [obsidian] vault backup: 2025-07-31 14:24:30[
All checks were successful
Build / build (push) Successful in 6m45s
2025-07-31 14:24:30 +09:00
bab3afd790 [obsidian] vault backup: 2025-07-31 12:24:15[
All checks were successful
Build / build (push) Successful in 9m35s
2025-07-31 12:24:15 +09:00
17df20a502 [obsidian] vault backup: 2025-07-26 10:59:10
Some checks failed
Build / build (push) Failing after 10m35s
2025-07-26 10:59:10 +09:00
51fa440e02 [obsidian] vault backup: 2025-07-25 21:06:47[
All checks were successful
Build / build (push) Successful in 11m19s
2025-07-25 21:06:47 +09:00
245acb784e [obsidian] vault backup: 2025-07-24 22:59:15[
All checks were successful
Build / build (push) Successful in 11m40s
2025-07-24 22:59:16 +09:00
baa447e054 [obsidian] vault backup: 2025-07-24 18:23:12[
All checks were successful
Build / build (push) Successful in 10m8s
2025-07-24 18:23:12 +09:00
b595629042 [obsidian] vault backup: 2025-07-22 23:44:33[
All checks were successful
Build / build (push) Successful in 8m41s
2025-07-22 23:44:33 +09:00
08a1cc5b30 [obsidian] vault backup: 2025-07-22 22:44:10[
Some checks failed
Build / build (push) Failing after 10m25s
2025-07-22 22:44:10 +09:00
81fbd200a8 [obsidian] vault backup: 2025-07-22 18:09:04[
Some checks failed
Build / build (push) Failing after 10m43s
2025-07-22 18:09:04 +09:00
6ef107b199 [obsidian] vault backup: 2025-07-22 17:09:00[
All checks were successful
Build / build (push) Successful in 9m15s
2025-07-22 17:09:00 +09:00
b3646b2c71 [obsidian] vault backup: 2025-07-21 21:43:09[
All checks were successful
Build / build (push) Successful in 9m45s
2025-07-21 21:43:09 +09:00
0eba0d32f6 Merge remote-tracking branch 'home/v4' into v4
All checks were successful
Build / build (push) Successful in 8m3s
2025-07-20 20:41:07 +09:00
e4d9573ffb [obsidian] vault backup: 2025-07-20 12:34:46[
Some checks failed
Build / build (push) Failing after 14m46s
2025-07-20 12:34:46 +09:00
101d6203a3 [obsidian] vault backup: 2025-07-20 11:34:26[
Some checks failed
Build / build (push) Failing after 10m7s
2025-07-20 11:34:27 +09:00
a5e33b30b3 [obsidian] vault backup: 2025-07-20 10:34:06[
All checks were successful
Build / build (push) Successful in 17m35s
2025-07-20 10:34:06 +09:00
43da672850 [obsidian] vault backup: 2025-07-20 02:31:41[
Some checks failed
Build / build (push) Failing after 12m24s
2025-07-20 02:31:42 +09:00
3d4c6fcf5c [obsidian] vault backup: 2025-07-18 18:14:10[ 2025-07-18 18:14:10 +09:00
c81997e527 [obsidian] vault backup: 2025-07-15 18:11:12[
All checks were successful
Build / build (push) Successful in 7m19s
2025-07-15 18:11:12 +09:00
cbe09d0852 [obsidian] vault backup: 2025-07-15 17:14:39[
All checks were successful
Build / build (push) Successful in 9m10s
2025-07-15 17:14:39 +09:00
160db1948f [obsidian] vault backup: 2025-07-15 17:11:03[
Some checks are pending
Build / build (push) Has started running
2025-07-15 17:11:03 +09:00
011edf7bed [obsidian] vault backup: 2025-07-11 20:10:23[
All checks were successful
Build / build (push) Successful in 8m35s
2025-07-11 20:10:24 +09:00
aa71637fe6 [obsidian] vault backup: 2025-07-11 13:08:46[
All checks were successful
Build / build (push) Successful in 7m55s
2025-07-11 13:08:46 +09:00
3dbbb6da30 [obsidian] vault backup: 2025-07-11 00:45:49[
All checks were successful
Build / build (push) Successful in 7m55s
2025-07-11 00:45:50 +09:00
bcfb639da7 [obsidian] vault backup: 2025-07-09 19:27:38[
All checks were successful
Build / build (push) Successful in 11m2s
2025-07-09 19:27:38 +09:00
23b72c51ab [obsidian] vault backup: 2025-07-08 17:15:15[
All checks were successful
Build / build (push) Successful in 8m27s
2025-07-08 17:15:15 +09:00
f0d85878da [obsidian] vault backup: 2025-07-08 16:12:02[ 2025-07-08 16:12:02 +09:00
9ceccfeec5 [obsidian] vault backup: 2025-07-05 15:48:36[
All checks were successful
Build / build (push) Successful in 8m7s
2025-07-05 15:48:36 +09:00
9ceda0f1a1 Merge remote-tracking branch 'home/v4' into v4
All checks were successful
Build / build (push) Successful in 9m9s
2025-07-04 13:17:31 +09:00
a47cd93b67 [obsidian] vault backup: 2025-07-04 06:05:53[ 2025-07-04 06:05:53 +09:00
21 changed files with 429 additions and 5 deletions

View File

@@ -0,0 +1,23 @@
#programming
主にゲームエンジンなどで採用されるプログラミングのパラダイム。
最近のUnityでも内部的に採用されている。
[[Rust]]だと[[Bevy]]が有名。
[[オブジェクト指向]]と比べると、データのメモリ分布が、オブジェクトごとに並ぶのではなく個別のメンバー変数ごとに並ぶことになり、CPUのメモリキャッシュに乗りやすいなどの利点がある。
[\[Rust\] ECSアーキテクチャ \[bevy\_ecs\] | DevelopersIO](https://dev.classmethod.jp/articles/ecs-rust-bevy/)
[Intro to ECS - Unofficial Bevy Cheat Book](https://bevy-cheatbook.github.io/programming/ecs-intro.html)
## データのモデリング方法としてなにがうれしいのか
オブジェクト指向と比べた時の利点が基本的にさっきのような、パフォーマンス面での利点が強調されることが多い。ただ、何かパフォーマンスのためにせっかくプログラミング言語が用意してくれたデータのモデリング技法を犠牲にしているような気がしてならず、あんまり旨味がよくわからなかった。
ただ、Bevyの設計をいろいろ読んでいると、モデリングとしての核心は、「集合の中から特定の要素を持つもの部分集合を指定して挙動を個別に操作する依存性がない操作同士はどういう順番で実行してもいい」ということかなと思った。

7
content/Gluon.md Normal file
View File

@@ -0,0 +1,7 @@
#programming-language
[[Rust]]で書かれた、[[Lua]]のような埋め込みを想定しつつ静的型付けの言語。
モジュールの読み込みの文法とか複雑な型情報をRustと相互でやり取りする時のやりかたとか、参考になる
[Gluon](https://gluon-lang.org/)

9
content/Hackett.md Normal file
View File

@@ -0,0 +1,9 @@
#programming-language
[[Racket]]上で実装された、静的型付けなうえで、型安全なマクロシステムを搭載したプログラミング言語
[1 The Hackett Guide](https://lexi-lambda.github.io/hackett/guide.html)
["Hackett: a metaprogrammable Haskell" by Alexis King - YouTube](https://www.youtube.com/watch?v=5QQdI3P7MdY)

View File

@@ -4,3 +4,9 @@ date: 2024-10-04 13:31
#person
[[Not Art&Tech - On the role of Media Theory at Universities of Applied Art, Technology and Art and Technology.]]
---
[Turing Complete User](https://contemporary-home-computing.org/turing-complete-user/)
[[橋本麦]]さんによる試訳 [『チューリング完全ユーザー』試訳 - baku](https://scrapbox.io/glisp/%E3%80%8E%E3%83%81%E3%83%A5%E3%83%BC%E3%83%AA%E3%83%B3%E3%82%B0%E5%AE%8C%E5%85%A8%E3%83%A6%E3%83%BC%E3%82%B6%E3%83%BC%E3%80%8F%E8%A9%A6%E8%A8%B3)

View File

@@ -0,0 +1,13 @@
#compiler-design
[Red-Green Trees: an Overview. Six months ago I dug into Roslyns… | by Bayastan | Jun, 2025 | Medium](https://medium.com/@krendelia2021/red-green-trees-an-overview-17bae2d84e8c)
[[抽象構文木]]というか[[具象構文木]]の実装方法の一つ。
構文の途中で差し込まれたコメントなどのトリビアなどを保持して、完全なソースコードに戻せるけど実行効率などをよくしたままにできるやり方。
部分的に構文木を書き換えたりするとソースコードのロケーション情報がずれたりするので、文字幅とオフセットを別々に持たせるというやり方をしている
ながいけどこれがわかりやすい
[Ruby Parser開発日誌 (19) - 最高の構文木の設計 2024年版 - かねこにっき](https://yui-knk.hatenablog.com/entry/2024/08/23/113543)

9
content/Rune.md Normal file
View File

@@ -0,0 +1,9 @@
#programming-language
[The Rune Programming Language](https://rune-rs.github.io/)
[[Rust]]で書かれた埋め込みプログラミング言語。
Rustの文法ほぼそのままに動的型付けの言語をつくるとこうなる、みたいな感じ
[[Moonbit]]とか[[Rhai]]も近い感じRhaiは直接影響を受けたとのこと

View File

@@ -30,7 +30,12 @@ https://github.com/mimium-org/mimium-rs
- [[mimiumのMIRコンパイル過程を真面目に考える]]
- [[lambda-mmm(実用版)]]
- [[mimiumグローバル環境評価について]]
- [[多段階計算を命令型VMインストラクションで表現したい]]
## マクロ
- [[mimiumと多段階計算]]
### ランタイム周りについて
@@ -46,6 +51,13 @@ https://github.com/mimium-org/mimium-rs
ふと思ったけど、[[SuperCollider]]や[[PureData]]と比べると、これらの言語は組み込みに使おうと思うとLinuxが動く環境を想定することになる[[Heavy]]はそれを全く別の処理系作ることで対応してたけど)。シーケンサとかスケジューラーがあるような、[[Faust]]だと難しいタイプのプログラムを[[Arduino]]とかに持っていくには向いているのではないか([[Extempore]]だって仕組み的に言えばそうかもしれないけど)
- できれば教育用途とかに持ち込めるのが一番いい
- ブートストラップできるといい
- いろいろな人がmimiumで拡張やライブラリを書けるようになってからが本当の本番
- 言語自体の拡張機能をその言語上でたくさん作れるとよい
- しかし、極端に自由度が高いとそれはそれで参入障壁が高い
[[mimiumでのシーケンサ]]
[[mimiumでのライブコーディングエンジン]]

View File

@@ -7,4 +7,32 @@ date: 2025-07-03 18:34
二重インクルードの問題とかもあるので、先に分割コンパイルの仕組みを整える方がいい?
型の環境の中に変数解決が必要になってくる
内部的には幽霊、Opaque、バリアントで分かれててもいいかもな
```rust
enum Variant{
Phantom(Symbol),
Opaque(Symbol,TypeNodeId),
Union(Box<Self>,Box<Self>)
}
enum Type{
...
Variant(Variant)
}
```
Foldでネストするくらいなら`Vec<Box<Self>>`に収められるかな?
パーサーとASTとしてはmatch節も加えないと意味ない
で、型推論に加えなきゃいけないコードは、
- 型コンストラクタを関数として使えるように型→値環境に追加し、App式の評価で使えるように
- match式の型推論
matchのコード生成がジャンプテーブルとか出てくるとめんどくさいんだよなあ

View File

@@ -19,3 +19,6 @@ ChucKではエフェクトのテールが更新時にぶちぎれる問題があ
各トラックごとのエフェクトのライブ切り替えとかも実現しようと思えばできるかな
...これ、結局[[SuperCollider]]のJITLibと同じことかもな
[jitlib\_basic\_concepts\_01 \| SuperCollider 3.14.0-dev Help](https://doc.sccode.org/Tutorials/JITLib/jitlib_basic_concepts_01.html)

View File

@@ -0,0 +1,71 @@
#memo #mimium
- [[多段階計算を命令型VMインストラクションで表現したい]]
- 一時期考えていたが、あんまり筋が良くないのでやめた
[Implement multi-stage computation intepreter by tomoyanonymous · Pull Request #136 · mimium-org/mimium-rs · GitHub](https://github.com/mimium-org/mimium-rs/pull/136)
シンプルに構文木レベルでのインタプリタを別途作ることで実現できそう。
## 課題
- 組み込みの関数の型情報がマクロ用レベル0とVM用レベル1で分かれているので、両方にあることを別々に宣言しなければならない
- 逆に、ベクターのappend・removeや文字列操作など、メモリを確保する操作はレベル0でのみ利用できる組み込み関数として制限すると、便利かもしれない
- 逆に、`self``delay`はレベル1でのみ利用できる関数とする
- レベル01両方で利用できるpersistentなモジュールは、上記2種類のいずれも使わず、かつescapeとBracketも利用しないものであればよい、ということになる
- もともとディレイの最大サイズはリテラルで指定しなければいけないという問題があったが、これを`make_delay(size:float)-> <(float,float)->float>`レベル0で最大時間を指定すると、入力と遅延時間をとる関数のコードを返すという関数にできるかも
- ただ、これやるとディレイが絡む関数はすべてレベル0定義になるのかな
- なんかそれよりは、[[数値プリミティブ型に常に範囲をつける]]とかのほうが筋がいいかも
```rust
fn fbdelay(max_time){
`|input,time,fb| {
(input + self*fb, time) |> make_delay!(max_time)
}
}
```
- こういうfbdelayを2個以上つかうときのコードはどうなるのかな
- というか、fbdelayで生成されたコードをMIRに持ってくときにどうなるのか
- Value型を単なるCodeじゃなくてDelayとして特別な値にリダクションすればいいのかな
[[otopoiesis]]でパラメーターを生成するのにも使える?
```rust
fn synth(freq,gate){
osc(freq)*gate
}
fn synth_module(freq:()->float,gate:()->float){
`||{ synth(freq(),gate()) }
}
```
これのモジュールを評価すると、freq,gateがUIに現れるという感じでできるのかなそして、ここでUIの範囲制限をするためにも数値型が範囲を持っていた方がいいということになりそう
うーん、サンクを手動で使わず表現できるような何かが欲しいなあ
```rust
fn wrap_module(param:Param, synth:(Param)->float ){
`{ | | synth(param |> invoke) }
}
```
Param型はinvokeでuiからの値を取れる、サンクをアンラップするようなメソッドを持つ型クラスに属している、という感じで、ジェネリクスが実装出来たらいけそうね
## マクロを提供するプラグイン
プラグインも、各関数が返す型と同時にどのレベルでその関数が使えるかをコンパイラに渡す必要がある
### ファイル読み込みはどうなる?
[[mimiumのファイルIO]]の話。マクロの展開より先に型の評価が行われるので、例えばファイルを読み込んでからチャンネル数が分かる、みたいな場合型情報を読み込む段階でファイルを読みこまなければならない、がその評価のためにはマクロを実行して文字列の値を受け取らなければならない、、、ということで、結局[[依存型]]がないと多段階計算でも対応できない。
結局、読み込んで返す値として次元数の値を持つ配列を返すしかないのかも

View File

@@ -147,7 +147,7 @@ myugen({ default_v <- freq = 200.0 })
---
## 部分型と型クラス
## 部分型と[[型クラス]]
[[構造的部分型]]を採用するつもり。

View File

@@ -123,6 +123,26 @@ fn reverse(origin:Region)->Region{
- 再生前prepareToPlay)
- 信号再生時process
## 多段階計算と組み合わせる
[[mimiumと多段階計算]]で、それなりに多段階計算の実装が間に合ってきた。
FadeinOutのようなリージョン→リージョンの関数はステージ0の計算と考えることができる。
また、Generator系も、基本的には周波数や音量といったUIパラメーターは、ステージ0での評価時にUIを生成して値を受け取るチャンネルを作り、ステージ1再生中にそのUIからの値を受け取るという方式で捉えられる
```rust
fn audiofx(param1=100,param2=200){
`|input|{... }
}
fn gen_component(){
Param{..} |> sinosc
}
```
`Param`はジェネリックな
---
以下は昔に考えていたこと

View File

@@ -0,0 +1,8 @@
---
date: 2025-07-05 15:45
---
#research
会議が煮詰まってくると、机にお菓子を投げ込んでくれるデバイス
https://ipsj.ixsq.nii.ac.jp/records/183464

View File

@@ -0,0 +1,10 @@
#compiler-design
- ある程度コンパイラのコードベースが大きくなってくると、認知負荷がでかくなってくる
- かといって、コードをコンパクトに保とうとすると、一箇所の変更がしづらくなっていく
- かといって、頑張ってインターフェースを切って分離性を高めていくと、全体の把握はしづらくなる
- そりゃ全体を把握しなくいても継続的に開発できるようにするのが目的なのでそうなんだけど
- 例えば、シンタックスツリーに[[Red Green Syntax Tree]]をあとから使おうとすると、さすがに書き換えが大変
- だが、Language Serverを作ろうと思ったらいつかは必要
- これをインクリメンタルに開発するのは無理だよなあ
-

View File

@@ -0,0 +1,41 @@
---
date: 2025-07-22 17:00
---
#programming-language
一般的なコードフォーマッターの実装について。
[How to write a code formatter](https://yorickpeterse.com/articles/how-to-write-a-code-formatter/)
1行当たり最大文字数がこのくらい、としたときに、どういう戦略で折り返すかは結構難しい問題。
文字列→トークン→ASTという順番で変換されるので、この逆順が良いのかとも思ったが、
AST→フォーマッタ用の専用の木構造みたいな中間表現を一度挟んだほうが賢いのかもしれない
[[Rust]]の[pretty](https://docs.rs/pretty/latest/pretty/)クレートがコンビネーターとして定義してあるやつ
[[Tree-sitter]]を使った汎用フォーマッター[Topiary](https://topiary.tweag.io/)とかいうのもある
---
Exprの途中に差し込まれたトリビア主にコメントをどうやって抽出するか
ExprNodeIdに対するSecondary MapがSpanに対して作れているのだから、Trailng Triviaとしてコメントを保持するのは一応できるか
パーサコンビネータでどうにか処理できるもんか?
```rust
fn parse_expr_top<Output>()->impl Parser<Token,Output,Error>{
not(comment()).padded_by(comment().repeated())
}
```
[Parser in chumsky - Rust](https://docs.rs/chumsky/latest/chumsky/trait.Parser.html#method.map_with)
`map_with`使えばいけるかしら
Stateにトリビアを書き込んでおけばいいのか

View File

@@ -0,0 +1,14 @@
---
date: 2025-07-09 19:14
---
#programming-language
言語のライブラリなどをリモートリポジトリから取得、管理する仕組み。
バージョンの依存性解決なども機能のうちに含む。
[[Suwa Takashi]]さんの記事がパッケージマネージャとは何かをよく解説してくれている
[パッケージマネージャを自作するときに考えること - gfnweb](https://gfngfn.github.io/ja/posts/2023-02-15-on-creating-package-managers/)

15
content/型クラス.md Normal file
View File

@@ -0,0 +1,15 @@
#programming-language
[[Rust]]におけるTraitとかに近いもの
あるメソッド群を持つジェネリックな型の分類
パラメトリックなジェネリクスに対して、対象が広すぎるものを、アドホック多相的に制限する
["Hackett: a metaprogrammable Haskell" by Alexis King - YouTube](https://youtu.be/5QQdI3P7MdY)
この動画の説明わかりやすかった12:10
式をもとに型を生成するのが[[型推論]]、型情報をもとに式を生成するのがジェネリクスというループ
[[Hackett]]

View File

@@ -32,9 +32,18 @@ http://www.cs.tsukuba.ac.jp/~kam/lecture/gairon2-2012/gairon2.pdf
λ○□:両方を統合
[[SATySFi]]:組版言語だが、マクロのシステムとして多段階計算が導入されている。
[多段階計算の型システムの基礎 - gfnweb](https://gfngfn.github.io/ja/posts/2022-05-12-slides-type-system-matsuri-2020/)([[Suwa Takashi]])
[SATySFi の多段階計算入門](https://sankantsu.hatenablog.com/entry/2022/08/19/215024)
## 用語
**Cross-Stage Persistence**基本的にステージ0/1でのコードはそれぞれステージ0/1の中でしか使えない。何らかの方法で両方のステージで跨いで使える仕組みを作れると便利
**Run**プリミティブ:
## [[依存型]]との組み合わせ
[A Dependently Typed Multi-Stage Calculus](https://arxiv.org/abs/1908.02035)

View File

@@ -18,9 +18,9 @@ VMは最初、マクロ評価時に自分の現在の評価ステージを0と
VMは、命令記録モードと実行モードの2種類で、現在の評価ステージが0のときに実行モードになる
mainプログラムの先頭はマクロ展開のため必ずincl_levelからスタート命令記録モードから開始
mainプログラムの先頭はマクロ展開のため必ずincr_levelからスタート命令記録モードから開始
命令記録モードでは、decl_levelがでてくるまでバッファに実行した命令をコピーし続ける
命令記録モードでは、decr_levelがでてくるまでバッファに実行した命令をコピーし続ける
decl_levelでレベル0になったら...うーん
```ocaml
@@ -62,3 +62,7 @@ endescape
ret 1
```
VMの構造にメタプログラミングを埋め込むよりも、MIR生成段階でなにかしらのサンクを呼び出すということにして一時的に空にしておく、みたいなやり方のほうが素直かもな

View File

@@ -0,0 +1,122 @@
---
date: 2025-07-24 17:54
---
#memo #mathematics
例えば、[[双方向プログラミング]]的なもので、数値をスライダーで調整できるようになっていた場合、実際のところは5/7とかわかりやすい有理数だったりとか、整数に近い値を表したいのにめちゃくちゃ長い小数点の数値が残ったりする。
これを、適当なグリッドにスナップするUIの一つとして、許容できる誤差範囲と整数の複雑さを指定して分数に変換するのが良いのではないか。
ChatGPTに聞いた。[[連分数展開]]というのを使うと良いらしい。[[Rust]]のコードを書いてもらった。
```rust
/// 任意の実数 x を近似する収束分数・半収束分数を求める
/// x: 近似したい実数
/// eps: 絶対誤差許容値
/// max_den: 分母の上限
/// max_num: 分子の上限
fn rational_approx(
x: f64,
eps: f64,
max_den: u64,
max_num: u64,
) -> Vec<(u64, u64)> {
// 連分数展開の a_k リストを構築
let mut a: Vec<u64> = Vec::new();
let mut r = x;
for _ in 0..64 {
let ak = r.floor() as u64;
a.push(ak);
let frac = r - (ak as f64);
if frac.abs() < f64::EPSILON {
break;
}
r = 1.0 / frac;
}
let mut candidates: Vec<(u64, u64)> = Vec::with_capacity(a.len());
let (mut p_nm2, mut q_nm2) = (0u128, 1u128);
let (mut p_nm1, mut q_nm1) = (1u128, 0u128);
for &ak in &a {
// 収束分数
let p = ak as u128 * p_nm1 + p_nm2;
let q = ak as u128 * q_nm1 + q_nm2;
if q > max_den as u128 { break; }
let num = p as u64;
let den = q as u64;
if num <= max_num {
let approx = (num as f64) / (den as f64);
if (x - approx).abs() <= eps {
candidates.push((num, den));
}
}
// 半収束分数
for d in 1..=(ak / 2) {
let ap = ak - d;
let p2 = ap as u128 * p_nm1 + p_nm2;
let q2 = ap as u128 * q_nm1 + q_nm2;
if q2 <= max_den as u128 && p2 <= max_num as u128 {
let num2 = p2 as u64;
let den2 = q2 as u64;
let approx2 = (num2 as f64) / (den2 as f64);
if (x - approx2).abs() <= eps {
candidates.push((num2, den2));
}
}
}
p_nm2 = p_nm1;
q_nm2 = q_nm1;
p_nm1 = p;
q_nm1 = q;
}
candidates.sort_by(|&(p1, q1), &(p2, q2)| {
let e1 = (x - (p1 as f64)/(q1 as f64)).abs();
let e2 = (x - (p2 as f64)/(q2 as f64)).abs();
e1.partial_cmp(&e2).unwrap()
});
candidates.dedup();
candidates
}
/// x を percent% 精度で近似するラッパー
/// x: 近似したい実数
/// percent: 相対誤差許容値(%
/// max_num: 分母、分子の上限
fn rational_approx_pct(
x: f64,
percent: f64,
max_num: u64,
) -> Vec<(u64, u64)> {
let eps = (percent / 100.0) * x.abs();
rational_approx(x, eps, max_num, max_num)
}
fn main() {
let x = 2.7182818284;
let percent = 1.0; // 1%
let max_num = 500;
let results = rational_approx_pct(x, percent, max_num);
println!("近似候補percent={}%:", percent);
for (p, q) in results {
let approx = p as f64 / q as f64;
let err = (x - approx).abs();
println!("{}/{} = {:.8}, 誤差 {} ({:.4}%差)",
p, q, approx, err, err / x.abs() * 100.0);
}
}
```
```
近似候補percent=1%:
193/71 = 2.71830986, 誤差 0.00002803075492963103 (0.0010%差)
106/39 = 2.71794872, 誤差 0.00033311045128181505 (0.0123%差)
87/32 = 2.71875000, 誤差 0.00046817160000012237 (0.0172%差)
68/25 = 2.72000000, 誤差 0.0017181716000003178 (0.0632%差)
49/18 = 2.72222222, 誤差 0.003940393822222443 (0.1450%差)
19/7 = 2.71428571, 誤差 0.003996114114285465 (0.1470%差)
```
悪くなさそうここから5個ぐらいまでを提示するとして、