黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2020年08月29日(土)
そりゃもー、数学者さん達と無関係に、誰かの権威に頼ることなく、「それは間違い」と断言できないようだと困る。 twitter.com/genkuroki/stat...
タグ：
posted at 00:04:01
#Julia言語
MLStyle.jl
Julia functional programming infrastructures and metaprogramming facilities
Topics:
pattern-matching
algebraic-data-types
syntax-sugar
scalability
statically-generated github.com/thautwarm/MLSt...
タグ： Julia言語
posted at 00:39:10
Implementing #rstats Notebook in #VSCode: github.com/Ikuyadeu/vscod... still at very early stage pic.twitter.com/QjcWNXerWu
posted at 01:23:50
#Julia言語 楕円曲線暗号で有名な Edwards curves と呼ばれる楕円曲線のプロット。このアニメーションは単位円を含んでいる。
Edwards curveによる楕円曲線の表示では加法公式は場合分けのない非常に綺麗な形をしているので、高速に計算できてかつ安全性も高い。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/3mjbiP8IdU
タグ： Julia言語
posted at 01:26:40
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
パナマの友人から「この曲が日本の作品だってのはわかってるんだけど、作品名がどうしてもわからないんだ！知っていたら教えてくれ。もう何年も探しているんだ」って連絡が来たから教えてあげたらすごく喜んでた
名もわからない作品を耳コピで楽譜を完成させたらしいけど、その熱意と演奏に感動したよ pic.twitter.com/kw8V8EZMLH
タグ：
posted at 08:10:46
#Julia言語 おめでとうございます。
このJuliaとJupyter notebookの連携に苦労している人は多いと思う。kernel deadが止まらなくなる失敗がある。
「み」の字さんのように、シンプルに pkg> add IJulia して、Juliaが~/.julia以下にインストールしたPython環境のjupyterを使うとそうなり難い。続く twitter.com/misaki70533995...
タグ： Julia言語
posted at 08:18:22
#Julia言語 すでにインストール済みのjupyterを使う場合には
julialang.github.io/IJulia.jl/stab...
にも書いてあるように、
julia> ENV["JUPYTER"] = "使いたいjupyterのフルパス"
としてから
pkg> add IJulia
pkg> build IJulia
とする必要があります。kernel deadが続く場合の原因がおそらくこれ。
タグ： Julia言語
posted at 08:22:10
#Julia言語 私による解説はここにある↓
nbviewer.jupyter.org/github/genkuro...
Julia v1.1.0 の Windows 8.1 へのインストール
ちょっと時代遅れになっている部分もありますが、色々書いてあるので、参考になることもあると思います。
タグ： Julia言語
posted at 08:24:18
#Julia言語
using Plots
x = y = range(-2, 2, length=51)
f(x, y) = x^2 - y^2
z = @. f(x', y)
wireframe(x, y, z)
で行けます。
comrehensionより、
z = @. f(x', y)
や
z = f.(x', y)
や
z = @. x'^2 - y^2
のbroadcastがシンプルに書けて便利です。
wireframe(x, y, f)
でも行ける。 twitter.com/misaki70533995... pic.twitter.com/UnQdbcz5dj
タグ： Julia言語
posted at 08:40:36
非公開
タグ：
posted at xx:xx:xx
#Julia言語
pyplot() バックエンドを使った場合のアニメーション。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/qvjpZAZETL
タグ： Julia言語
posted at 09:10:24
Anyone aware of work on a Julia equivalent?
Notebook directly in Vscode would be cool.
#julialang twitter.com/renkun_ken/sta...
タグ： julialang
posted at 09:21:47
田中秀臣さん”MMTには権力関係が中心を占めています。需要供給の市場メカニズムは基本的に不在と考えていいでしょう。実証というよりは規範的で、MMTが依拠する「現実」は解釈学的な境地“
市場メカニズム、期待が欠けるのは人間の生活を軽視した印象。どおりでMMフレンズになじめない^_^ twitter.com/shinchanchi/st...
タグ：
posted at 09:23:29
#Julia言語
plotlyjs() バックエンド製のsurface(x, y, f)のアニメーション pic.twitter.com/acsMOGFwnN
タグ： Julia言語
posted at 09:26:44
MMTは税や財政にビルトインされた「仕掛け」でインフレ率をうまくコントロールする考えらしい。
レイの入門本を見ても価格統制の話題に触れられていたし、MMフレンズは、実証もモデルも不充分なのに
「魔法使いの先生より上手く物価コントロールできるもん！」
とのたまう”我賢し“な魔法使いの弟子デシ
タグ：
posted at 09:32:27
田中秀臣さん“「緊縮財政が悪い」というのは一面正しいのだけれども、それは、別にMMTだけのオリジナルのロジックでなく、経済学者も非経済学者も正統も異端も全員が共通して論じていること“
まさに。既存の議論より有益な知見が見当たらないMM…
タグ：
posted at 09:38:57
@hidetomitanaka @KoheiMorinaga とても興味深く勉強になりました。期待や市場メカニズムと縁遠いMMフレンズとは友達になれそうもないですが、経済安定化のためにマクロ経済政策が重要であり、その先には人間の生活や人生があることを忘れてはならないと改めて思いました。 twitter.com/shinchanchi/st...
タグ：
posted at 09:45:37
参考にしたのはこの黒木さんの半年前の解説。ちなみに言及されている記事のkatsu1110氏は僕ではありません。
twitter.com/genkuroki/stat...
タグ：
posted at 12:41:58
ごまふあざらし(GomahuAzaras @MathSorcerer
Python の psutil 依存だけれどCPU利用状況をリアルタイムにモニタリングする機能を作れた． pic.twitter.com/Pb6huDJhYR
タグ：
posted at 13:26:36
「信頼区間はその区間にデータが含まれる確率じゃないですよ」と言った後に「それはベイズ確信区間」というところまでがセットのまずい説明。そこから「ベイズの方が直感的で分かりやすいからベイズが良いですよ」と続く場合が多い。数年やってみて話が違うってなりました。
タグ：
posted at 13:28:05
Dockerにrust, latex, julia, pyhthon環境整えて、theiaとjupyterlabでUI立てて、nginxでリバースプロクシした。これでコンテナに引き篭れる。
タグ：
posted at 13:34:46
Comparison of the Wasserstein, Hellinger, Kullback-Leibler and reverse KL on the space of Gaussian distributions. djalil.chafai.net/blog/2020/01/2... pic.twitter.com/BtekhkMZGK
タグ：
posted at 14:00:00
自宅のパソコンで数値計算していて、メモリ100%、ディスクアクセスも100%になった。さぁ、ここからが君の本領だ！とマウスを握ったが、マウスポインタすら出てこない。
タグ：
posted at 14:18:39
@genkuroki 掛け算順序や漢字のとめはらいなど小学校限定ルールでバツにされることに対して「本当は（大人になったら）どっちでも良い」を教える事が出来ても「どっちでも良いのに先生はなぜバツにするのか？」という子供の疑問に答えられず「先生は教え方の教科書通りにしか教えられないらしい」と説明しました。
タグ：
posted at 14:36:26
@genkuroki 「学校では先生に合わせてマルをもらって来なさい」と指導したことはありませんが「納得いかなくても堪えてやって」と言ったように記憶してます。今から振り返ってもどうすればよかったのか、よくわかりません。
ただ、算数から数学に変わってどんどん伸びていったので、悪くはなかったのかなと。
タグ：
posted at 14:49:38
非公開
タグ：
posted at xx:xx:xx
Really interesting blog post about Julia's code generation and compilation latency improvements:
julialang.org/blog/2020/08/i...
my random thoughts on this follows ...
タグ：
posted at 16:09:44
Roughly speaking, the current approach to improve the compiler latency problem is fixing type-unstabilities that invalidate already generated code caches.
Well written code is often type-stable by itself, but still we often need additional type annotations that help the compiler
タグ：
posted at 16:09:44
generate more type-stable cache.
Julia's flexible design allows us to balance the trade-off between human vs. machine needs, and here we (humans) are willing to help machine (the compiler) by additional, sensible annotations.
Having said that, the annotations usually cumbersome
タグ：
posted at 16:09:45
and kind of against the Julia's philosophy; they can spoil the generics of code when scattered around without sensible manner.
So what matters would be good tooling. We need such tools that help us make type annotations easily and correctly, ideally in an automated way.
タグ：
posted at 16:09:45
SnoopCompile.jl is one of the examples of such an attempt if I understand correctly. But I think there still remains huge possible solution space in this area. The recent work to expose APIs to interact with Julia's native compiler can be helpful for example.
タグ：
posted at 16:09:45
Juliaのcompiler latencyに関するブログポストを読んだので、雑な感想を書いていきます:
julialang.org/blog/2020/08/i... twitter.com/kdwkshh/status...
タグ：
posted at 16:20:51
簡単にまとめると 、型不安定なコードはすでに生成された(それに依存する)コードキャッシュを"invalidate"してしまうため、コンパイラのlatencyを低めるためにJuliaユーザが取りうるアプローチは、コードの型安定性を高めることになる。
よく書かれたコードは往々にして型安定だけど、それでもやはり
タグ：
posted at 16:20:51
コンパイラにより型安定なコードを生成させるための追加的な型注釈が必要になる場合がある。
Juliaはプログラマが人間と機械の間の(面倒さと効率の)トレードオフのバランスを取ることができるように設計されていて、ここで人間側は追加的な型注釈を施すことによって機械(コンパイラ)にヒントを与える
タグ：
posted at 16:20:52
とはいえ、型注釈は基本的に面倒だし適当に使われた場合はJuliaコードの持ちうる多相性を損なってしまうので、そうした型注釈を正確に簡単にできれば自動化されたやり方で行えるようにするtoolingが重要になる。SnoopCompile.jlはその例の1つだけど、このあたりにはまだまだ色々な選択肢がありそう。
タグ：
posted at 16:20:52
非公開
タグ：
posted at xx:xx:xx
@yukisan65 おお！そういうナマの声は参考になります。
コメントどうもありがとうございます。
親よりも当事者の子供の側が結局色々分かっているというようなこともあるのではないかと思いました。
タグ：
posted at 17:07:41
非公開
タグ：
posted at xx:xx:xx
@mat_der_D 嫌われそうですがPythonやりはしめました。Jupyter に感激したからです。2日でここまで来たので次に同報メール送信とかして、Pandusをやります。基本早いのが好きですからジュリアとか、Vとかをとても注目してます。
タグ：
posted at 17:31:21
ロジスティック方程式のベクトル場と4次ルンゲ＝クッタ法でx=0.1からスタートした場合の動き。Juliaでプロット。 pic.twitter.com/mHc08tCZ5X
タグ：
posted at 17:32:18
#Julia言語
Juliaは複数の意味で「過学習者殺し」
* 1-based index
* 試行錯誤時にjuliaの再起動を繰り返さない
どうしても0-indexを使いたい場合には
github.com/JuliaArrays/Of...
OffsetArrays.jl
Fortran-like arrays with arbitrary,
zero or negative starting indices.
を使ってみるとか。 twitter.com/chohcooh2/stat...
タグ： Julia言語
posted at 17:34:50
#Julia言語 Edwards曲線(楕円曲線のいち表示)の
x² + y² = 1 + m x²y²
は
(x, y) = (cd(u|m), sn(u|m)), 0≤u≤4K(m)
でパラメトライズされる.
添付動画では x² + y² = 1 + m x²y² を m=1+50～1-50についてプロット。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/P1PeCUDHjW
タグ： Julia言語
posted at 18:14:09
#数楽 Jacobiのsn函数がsinの一般化になっていることは有名。Edwards曲線
x² + y² = 1 + m x²y²
(x, y) = (cd(u|m), sn(u|m)), 0≤u≤4K(m)
の立場でのcosの一般化はJacobiのcd函数になる。
円周率は第一種楕円積分K(m)の2倍に一般化される。
github.com/nolta/Elliptic...
タグ： 数楽
posted at 18:17:32
@mat_der_D ありがとうございます。
とにかくオブジェクト指向が嫌いでしたのでやってませんでした。
でも開発環境？IDE？エディタ？わかりませんが Jupyter ですね、これは面白い。これならすぐできると思いましてやり始めました。今週からです
タグ：
posted at 18:19:31
例の Julia で書いたら Python より約25倍速くなったコード載せました
はてなブログに投稿しました #はてなブログ
Julia が Python より約25倍速いとあるコード - ぷりんの雑記帳
smooth-pudding.hatenablog.com/entry/2020/08/...
タグ： はてなブログ
posted at 18:40:18
Juliaは流行っているしEmacsは最も成功したテキストエディタだし、NixOpsさえあればk8sもDockerもTerraformも要らない。いいね？
タグ：
posted at 19:37:02
非公開
タグ：
posted at xx:xx:xx
#Julia言語
github.com/chualynn/Theta...
Theta.jl
A Julia package for computing the Riemann theta function and its derivatives
こういうので遊んでみると理解が進むタイプの人が結構いるかも。
タグ： Julia言語
posted at 20:58:32
辞任表明した安倍晋三首相の後継を決める自民党総裁選をめぐり、菅義偉官房長官が立候補の検討に入りました。複数の派閥幹部から推す声が出ており、出馬すれば有力候補になるのは確実です。
www.jiji.com/jc/article?k=2...
タグ：
posted at 21:06:03
#Julia言語
オブジェクト指向が嫌いなら、
Juliaタイプのmultiple dispatch
を気に入る可能性があると思う。
個人的にJuliaはJupyterとも相性がよいと思います。
Jupyterならセル1でusing Plots; plot(sin)しているあいだに、セル2でコードの入力を始められる。 twitter.com/hanatmetal/sta...
タグ： Julia言語
posted at 21:07:31
#Julia言語 でOOPをシミュレートする方法については以下のリンク先を参照。
① abstract types と composite types のtreeの階層ごとにメソッドを定義する。
② Holy traitsを使う。各型の引数を型に対応させる函数を用意しておいて、その函数の値(型になる)の違いで(multiple) dispatchする。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 21:16:34
「いちいち人格否定せず単に指摘すればよい」のなら「本当によくできる人間」を挙げる必要ないのでは？
そのタイミングでのその対比は「本当によくできる人ではない」と黒木さんの人格否定してません？
そもそもの元ツイ(アンケ)の酷さを理解してないとしか思えない。「考えの違い」ってね。 twitter.com/f_sei/status/1... pic.twitter.com/VqlY4qLsgY
タグ：
posted at 21:39:21
「大学入試で減点される！」ネタを振ることが悪いことだということを自覚さえできない人って結構いるよね。
そういう人達のそういう感覚は社会的に有害であるという認識が広まり、実際に「大学入試で減点される！」ネタを振った人が社会的に厳しく非難されるようになるとよいと思います。
タグ：
posted at 21:47:11
多くの人はCやC++やFortranを数値計算に使っており、現在Juliaも有力視されているのだと思います。
#Julia言語 がお勧め。
特に自分でコードを書いて勉強したいのであれば、素直に書いて速いJuliaを選ぶ方が合理的だと思います。
Pythonにしかないライブラリを使う特殊な仕事以外はJuliaが勝る。 twitter.com/kntoo/status/1...
タグ： Julia言語
posted at 22:00:32
自分自身で書いたコードを示さずに、実戦的にはどこまで意味のある結果なのかわからないマイクロベンチマークを引用して、CやC++の側がJuliaよりも(少し)速いことを当然視する発言をする人が結構いますが、実際に自分がコードを書いたときにどうなるかはやってみないと分かりません。
タグ：
posted at 22:07:22
SOLEIL DATA DOJO @soleildatadojo
「速くて書きやすい」次世代の言語Juliaを学ぼう! Julia勉強会#2~機械学習の実装その1~ を公開しました！ soleildatadojo.connpass.com/event/186505/?... #DeepLearning
タグ： DeepLearning
posted at 22:08:25
#Julia言語 正五角形上の自由境界条件の波動方程式の数値解のつもりの計算結果のプロット
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/FlNRY1eYbK
タグ： Julia言語
posted at 22:13:00
#Julia言語 反応拡散系によるチューリングパターン(Turing pattern)
200x200の各点での(u,v)の値の分布が右側のグラフ。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/nDXQ421Nnt
タグ： Julia言語
posted at 22:19:35
「速くて書きやすい」次世代の言語Juliaを学ぼう! Julia勉強会#2~機械学習の実装その1~ に参加を申し込みました！ soleildatadojo.connpass.com/event/186505/?... #DeepLearning
タグ： DeepLearning
posted at 22:26:38
ごまふあざらし(GomahuAzaras @MathSorcerer
cat ~/.julia/prefs/PyCall とすると
JuliaがどのPythonを呼ぶのかを知ることができるっぽい．
PyCall.jl の build.jl を眺めながら．
タグ：
posted at 22:29:12
Pythonで数値計算を勉強しようとして、「メモリの無駄遣いを気にせずにNumPyの配列計算に持ち込む」というような特殊なノウハウを学ぶようになってしまうのは非常にまずい。
そういうのは現場仕事のノウハウで時代が過ぎれば通用しなくなる。
まずは時代によらず通用する普遍的知識を学ぶべき。
タグ：
posted at 22:32:06
ごまふあざらし(GomahuAzaras @MathSorcerer
雰囲気で言わんとすることはわかるけれど，厳密に読もうとするとピィーと叫びたくなるあざらし．
#Julia言語
画像は PyCall.jl/deps/build.jl の一部より pic.twitter.com/VhumP7Dk0g
タグ： Julia言語
posted at 22:35:30
ごまふあざらし(GomahuAzaras @MathSorcerer
UseCondaPython という例外を作ってシステム上にPythonがなければその例外をキャッチして Conda.PYTHONDIR に miniconda の Python 入るのね．．．歴史を感じます．
#Julia言語
PyCall.jl を読む日記
タグ： Julia言語
posted at 22:52:04
#数楽 それ、結構あります。
固有値より固有ベクトルの方が簡単に求まる場合がある。
例えば、易しい例では、1次元格子上での離散ラプラシアン。これは古典有限型および古典アフィン型のCartan行列に一般化される。型の違いは境界条件の違いに対応。
解説↓
nbviewer.jupyter.org/github/genkuro... twitter.com/kyow_qq/status...
タグ： 数楽
posted at 22:56:55
Aegis@おもろいエンジニアの集い場 @EngineerMeetup
「速くて書きやすい」次世代の言語Juliaを学ぼう! Julia勉強会#2~機械学習の実装その1~ / soleildatadojo.connpass.com/event/186505/?... #エンジニア #zoom #相互フォロー pic.twitter.com/XvSEl9YONi
posted at 23:07:33
#数楽 アフィン型の場合も含めて、古典型Cartan行列が1次元の離散ラプラシアンの境界条件の違い(およびそのちょっとした拡張)から出て来るという易しい話の解説がググっても見当たらなかったので書いた。
単に書くだけなのですが、めんどくさかったです。
nbviewer.jupyter.org/github/genkuro...
タグ： 数楽
posted at 23:14:31
#Julia言語
/usr/local/bin/julia のリンク先が /path/to/install/folder/bin/julia になっているように見えます。
消すべきシンボリックリンクを消して、シンボリックリンクをはり直せばきっと意図した通りになると思います。
結果の確認法：
which julia
ls -l `which julia` twitter.com/misaki70533995...
タグ： Julia言語
posted at 23:30:18
ごまふあざらし(GomahuAzaras @MathSorcerer
やっぱり，ソースからビルドさせるべきですね（白目） twitter.com/genkuroki/stat...
タグ：
posted at 23:31:54
問題は歩き方のパターンが腕が違うとどのように変わるか、変わらないルールはあるか、です。実験としては、腕をつつく→歩きだすをビデオ撮影して、解析は、中心部と腕をトラッキング→統計解析というシンプルなものです。
タグ：
posted at 23:35:32
ただ、恣意的に設定した混合分布モデルや階層モデルの汎化を調べて良さを評価したいが、それにはAICは使えない。そこで広く使える情報規準WAICを使ったというわけです。
タグ：
posted at 23:35:32
#Julia言語 さすがに数値積分のライブラリは別に入れる必要があります。
1変数函数の積分では
github.com/JuliaMath/Quad...
pkg> add QuadGK
が定番で、多変数函数の積分では
github.com/JuliaMath/HCub...
pkg> add HCubature
を私は愛用しています。 twitter.com/stmtk_yu/statu...
タグ： Julia言語
posted at 23:36:56
Juliaデータサイエンス―Juliaを使って自分でゼロから作るデータサイエンス世界の探索 (日本語) 単行本
amzn.to/3ba8O0h
こんなんあるんや
タグ：
posted at 23:56:15
自分だけでなく会社での布教にも書籍は役立ちます。ネット情報は「新しいもの好きの若手がネットで見つけてきた、よくわからない人による、実績のない技術の資料」で、書籍は「勤勉な若手が購入した、執筆が出来る人間による、本が出るくらいには使われている技術の資料」です
twitter.com/genkuroki/stat...
タグ：
posted at 23:57:14
