黒木玄 Gen Kuroki(@genkuroki)/2018年01月31日

ちょうどjuliaでカルマン渦試していたところだった

タグ：

posted at 23:13:24

正直な話、二度と使いたくない。

タグ：

posted at 23:00:49

ceptree @ceptree

非情なコード twitter.com/genkuroki/stat...

タグ：

posted at 23:00:35

FORTRAN77で教育を受けた世代。

タグ：

posted at 23:00:04

#Juliag言語小ネタ：Juliaでgotoを使ってみた。@ の後の空白な除く。

begin
@ goto a
println("とばさないで！")
@ label a
println("とばしました。")
end
→とばしました。

gotoを使ったはもしかしたら大学1年のときにFORTRAN77を使ったとき以来かもしれない。 pic.twitter.com/TGaJOpk3Hv

タグ： Juliag言語

posted at 22:59:22

非公開

タグ：

posted at xx:xx:xx

Daigo Akane (赤根大吾) @dgakane

@genkuroki 黒木さんありがとうございます。まだちゃんと理解したとは到底言えない状態ですが、教えていただいた筋道を自分なりのペースでたどってみます！ ^_^

タグ：

posted at 22:45:36

juliaで列を取り出したいときもlist[:,n]でええんやな

タグ：

posted at 22:43:47

須山敦志 Suyama Atsushi @sammy_suyama

#数楽ここ数十年のあいだに、数学がとても役に立つこと(普通の意味で役に立つこと)がばれてしまったのですが、それはわけのわからないことをやっていた人達の昔の仕事も含めて普通に役に立っているのであり、役に立つ数学研究をしているから役に立っているのではありません。そこを誤解するとまずい。

タグ：数楽

posted at 22:40:58

ベイズでは欠損値補間は特別興味の対象にならないことを示しています．記事ではロジスティック回帰の入力欠損値をMCMCで推論していますが，代わりに変分推論を使えばVAEと（ほぼ）同じになります．
machine-learning.hatenablog.com/entry/2017/08/...

タグ：

posted at 22:39:50

#数楽ヤコビさんの仕事の真の価値を理解していた人達は200年前にはほんの少ししかいなかったと思う。通常の意味で「役に立つ」ものだとはきっと誰も考えなかったと思います。

数学の研究はそういうものなので、わざわざ数学を研究してくれる稀有な人達には好きなようにやらせるのが正解だと思う。

タグ：数楽

posted at 22:38:54

須山敦志 Suyama Atsushi @sammy_suyama

#数楽ヤコビさんの楕円函数に関する研究を読もうと挑戦した人は結構いると思いますが、あれ、私には正直読めません。元祖ヤコビさんの記号法は現代でも結構通用していて、雰囲気はわかるのですが、複雑な計算がたくさん書いてある感じでちょっと読めない。

タグ：数楽

posted at 22:37:19

コメント欄では激励を受けていますが，個人的にはお気に入りです．かなり素直な式展開で変分推論が導けた例です．
machine-learning.hatenablog.com/entry/2017/08/...

タグ：

posted at 22:37:15

#数楽こういう話を本当に理解しようとすると大変ですが、私が示した円、楕円、双曲線、放物線、Edwards曲線の画像を見れば、全部連続的につながっていることが誰にでもわかると思います。円なら小学生でも知っている。そこから新しいタイプの数学テクノロジーの話題は連続的に繋がっているのです。

タグ：数楽

posted at 22:35:05

#数楽 Jacobiの楕円函数のきれいでシンプルな公式がコンピューターによる高速計算に向いており、しかも例外的な場合分けを含まないので、セキュリティー的にも優れている。こういう形で200年前のヤコビさんによる楕円函数の研究は現代の楕円曲線暗号で役に立っているわけです。

タグ：数楽

posted at 22:33:03

#数楽楕円函数にも色々あるのですが、Jacobiの楕円函数の加法公式は相当にきれいな形をしています。コンピューター社会においてきれいな公式は高速計算の源泉になります。そして例外を場合分けしないきれいな公式はセキュリティー的に重要です。例外的な処理が起こった瞬間は外部から測定可能。続く

タグ：数楽

posted at 22:30:14

#数楽「円の群構造＝三角函数の加法公式」の一般化として「楕円曲線の群構造＝楕円函数の加法公式」があります。楕円曲線をEdwards形式で実現している場合には「Edwards曲線の群構造＝Jacobiの楕円函数の加法公式」となります。200年前のヤコビさんの仕事(研究)。

タグ：数楽

posted at 22:28:25

#数楽 Edwards曲線は楕円曲線のある表示のこと。楕円曲線の同値な定義がたくさんあって、「楕円曲線とはEdwards曲線のことである」という方法で楕円曲線を定義することもできます。そして、Edwards曲線はパラメーターが特別な場合として円(二次曲線)を含みます。続く

twitter.com/dgakane/status...

タグ：数楽

posted at 22:26:19

Daigo Akane (赤根大吾) @dgakane

そんなことを呟いていると、黒木玄 @genkuroki さんが何やら関係ありそうなツイートをされている。楕円曲線とエドワーズ曲線というのは別々のものではなく、何か繋がっている話？ん？三角関数も？ twitter.com/genkuroki/stat...

タグ：

posted at 22:23:09

#数楽 Edwards formでの楕円曲線の群構造は200年前のヤコビさんの楕円函数の仕事までさかのぼれます。ヤコビさんは自分が作った楕円函数(三角函数の一般化)の加法公式が暗号の基礎理論になっていてかつたくさんのコンピューターで繰り返し計算される社会になっていることを想像できなかったはずです。

タグ：数楽

posted at 22:19:51

#数楽高校数学を正しく理解している人達は「ああ、なるほど。あの易しい三角函数の加法公式の話の一般化ね！」と気楽に理解する道がひらけるわけです。

タグ：数楽

posted at 22:17:37

#数楽そして、私が示した図の方法で、楕円曲線のEdwards formに足し算の構造を与えることができます。その足し算の構造は三角函数の加法公式の直接的な一般化になっています。そして、その加法公式の一般化が楕円曲線暗号の形式で我々の社会の中で役に立っているわけです。

タグ：数楽

posted at 22:16:32

#数楽楕円曲線暗号の普及のおかげで、「楕円曲線ってなんかすごそう」のように思っている人は増えています。しかし、現実の楕円曲線は「すごそうなあちら側の世界の住人」ではなくて、円 x^2 + y^2 = 1 を x^2 + y^2 = 1 + d x^2 y^2 に一般化すれば得られ、図もコンピューターで簡単に描けます。

タグ：数楽

posted at 22:15:08

#数楽一つ前のツイートの添付画像の曲線の式は

x^2 + y^2 = 1 + d x^2 y^2 (楕円曲線のEdwards form)

の d = -10 の場合です。この式のd=0の場合が単位円で、高校で習う三角函数の話になっています。それがそのままd≠0の場合に一般化されるのです。

図で形を認識することは誰にでもできる。

タグ：数楽

posted at 22:12:15

#数楽一つ前のツイートの添付画像もこのツイートの添付画像も楕円曲線のEdwards formです。

geogebra.org/m/hnvMuRtb pic.twitter.com/H6obDEubLZ

タグ：数楽

posted at 22:08:53

#数楽そして、じっと眺めると、点-(A+B)および点A+Bを赤色の双曲線を用いて構成できることもわかります。この双曲線を用いた二次曲線の群構造の特徴付けは楕円曲線のEdwards formにそのまま完全に同じスタイルで一般化されます。

www.geogebra.org/m/V3DbvDqC
pic.twitter.com/GeBICI095q

タグ：数楽

posted at 22:06:33

#数楽以上の添付画像には直線ABに平行な点Oを通る直線だけではなく、点Oを通る接線とそれに平行な点A+Bを通る直線も描かれています。それだけではなく、赤色の双曲線も描かれている。

タグ：数楽

posted at 22:04:31

#数楽放物線の場合

この場合は単純で、x座標の和が A,B→A+B に対応しています。

www.geogebra.org/m/eqJerCxh pic.twitter.com/wUgTFs8C1z

タグ：数楽

posted at 22:01:05

#数楽双曲線の場合

geogebra.org/m/FYuGc48z pic.twitter.com/bl0sSdOUyL

タグ：数楽

posted at 21:59:56

#数楽そして、その作図法には角度の概念が一切使用されていません！
だから、即、楕円を含む任意の平面二次曲線に一般化されます。
添付画像は楕円の場合です。geogebraのリンクをクリックすればマウスでぐりぐり動かして遊べる世界に突入！

www.geogebra.org/m/P6MpBarH
pic.twitter.com/GEaRNZQTqU

タグ：数楽

posted at 21:58:42

#数楽点Aと点Bから点A+Bを作図するのは簡単で、直線ABと平行な点O=(1,0)を通る直線と単位円の交点をA+Bとするだけです。この作図を座標を使って計算し直せばcosとsinの加法公式が得られます。計算が必要な点はちょっと面倒になりますが、図は圧倒的にシンプルになっています！続く

タグ：数楽

posted at 21:56:25

#数楽三角函数の加法公式ネタがこんなに拡散されるとはびっくりなのですが、加法公式を添付画像のようなこれまたシンプルな図で理解することも可能です。角度の意味で点Aと点Bの和が点A+Bになっている。続く

www.geogebra.org/m/xbrMSVt5 pic.twitter.com/Q9iiJnr08j

タグ：数楽

posted at 21:54:30

Haruhiko Okumura @h_okumura

アンケート結果
1366 （右利き）右耳に当てる
1025 （右利き）左耳に当てる
91 （左利き）右耳に当てる
174 （左利き）左耳に当てる
twitter.com/h_okumura/stat...

タグ：

posted at 21:26:09

HIROPON @hiropon_matsu

「共形場理論入門」「結び目の数学と物理」両方とも出版社は培風館初版部数少ないし、増刷しないので古本でもリーゾナブルな価格で出回らない。

タグ：

posted at 21:21:37

#数楽 BPZはすでに34年前の論文になってしまったのですが、純粋数学的には今このタイミングで読む価値がありそうな論文の一つになっていると思います。なぜならばBPZの論文に書いてある内容は単独二階の線形常微分方程式のモノドロミー保存変形の理論の正準量子化になっているからです。

タグ：数楽

posted at 21:17:13

当時のロシア人数学者達は

a+b+c+
+d+e

のように行の終わりと始めに+を書く習慣があって、しかも当時は手描きだったので、

a+b+d+d+e

の意味ではなく、a+b+c+?+d+e の ? の部分が欠けているという疑いを晴らすために大変な計算をしなければいけなかった。私は結局REDUCEを使ったことがある。

タグ：

posted at 21:12:40

#数楽 BPZというのはこれ↓

scholar.google.co.jp/scholar?cluste...

Feigin-FuchsのVirasoro代数の表現論に関する仕事を知らないと大変かもしれないが、今なら情報がたくさんある。

タグ：数楽

posted at 21:11:13

optical_frog @optical_frog

レンルイスせんせが紹介してるこの例みたいな「ネオリベラル」の論じ方もあるよね: goo.gl/FYPwH9

タグ：

posted at 21:04:05

#数楽で、z-wについて負べきの方は有限で止まっていて、正べきの方は形式べき級数として無限に伸びている場合を考えるのですが、形式べき級数が出て来ない本質的に同じやり方があって、それがBeilinson-Drinfeldのchiral algebraです。

タグ：数楽

posted at 21:02:23

#数楽続き～、vertex algebra の複雑に見える公理系がそのまんま出て来ます。そして、昔の共形場理論の基本文献はそういうスタイルで書いてあります。理論物理学者が書いた文献なので、公理系の形で整理されていない**だけ**で、よく見ればそう理解できることが書いてあります。

タグ：数楽

posted at 21:00:16

#数楽～だから、BPZまで戻って勉強した方が楽ということはあると思う。アイデアは単純で、a(z)とb(w)のあいだの

a(z)b(w) = Σ (z-w)^{-k-h} (ab)_k(w) (z-wに関するLaurent展開型)

のような形の積構造を適切に定義するにはどうすればよいかという話。可換性や結合律を書き下せば～続く

タグ：数楽

posted at 20:57:44

#数楽今だと2次元量子共形場理論を勉強しようとすると、vertex algebraの公理系を勉強させられるはめに陥るパターンが多いような気がします。たぶん、そういう経路で勉強すると公理系の見た目の非本質的な複雑さに騙されて苦しい思いをしてしまう。

だから、Belavin-Polyakov-Zamolodchikovまで～続

タグ：数楽

posted at 20:53:13

Juliaびっくりするくらい簡単に使えるから全能感がパない

タグ：

posted at 20:43:18

Rあまり使いたくないので努力を積んでJuliaガチ勢になろうと思う

タグ：

posted at 20:41:46

[入門Julia] よく使うパッケージまとめ on @Qiita qiita.com/saicologic/ite...

タグ：

posted at 20:40:34

＃超算数　「どうして日本の小学校ではトンデモ算数が堂々と行われているか？　（仮）」みたいなタイトルで特集して欲しい。

タグ：超算数

posted at 20:35:24

＃超算数　　この理由は、「算数教育界の複雑な権力構造」を知らないと理解できない。NHKの教育テレビで解明して欲しいと思う。あるいは小学校教育に関心のあるライターが本にしたりしないのかなあ？

タグ：超算数

posted at 20:34:00

＃超算数　小学校の掛け算順序はネットではかなり有名になり、理論的なことは出尽くした感じがある。でも、そのようなトンデモ算数がどうして日本全国の小学校算数授業で教師に支持されているか？・・・という理由を知らない人は多い。

タグ：超算数

posted at 20:31:25

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

＃超算数　『新しい算数研究』の２月号を購入。今月号は特別号なので高い。２８００円。しかし、この号は算数教育の権威リストにもなっているので、持っとくと何かと便利

タグ：超算数

posted at 20:28:09

ついにJuliaが1万スター！

タグ：

posted at 20:27:45

@MathSorcerer @Qiita #Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
PyPlot と Plots の比較

に

nbviewer.jupyter.org/gist/genkuroki...
JuliaカーネルのJupyterでjavascriptアニメーションを動かす

の方法を追加した。

一つの道具(Julia言語)で複数の選択肢があるのはうれしい。

Plots.gr()が相対的に爆速。

タグ： Julia言語

posted at 20:27:10

Viral B. Shah @Viral_B_Shah

The #julialang repo on GitHub now has 10,000 stars. Nice milestone! github.com/JuliaLang/julia pic.twitter.com/ymLVbs0bKu

タグ： julialang

posted at 20:14:14

@MathSorcerer @Qiita qiita.com/yoku_001/items...
Jupyterでmatplotlibのjavascriptアニメーションを動かす

おお！素晴らしい！速い！

ことを #Julia言語でも簡単にできました！

Julia言語からPythonライブラリを簡単に使えるので、Python関係の情報の多さはかなりありがたい。

nbviewer.jupyter.org/gist/genkuroki...

タグ： Julia言語

posted at 20:07:58

#Julia言語

 github.com/JuliaPlots/Plo...

の変更(L557-L564)のせいで、Plots.gr()を使った場合にのみ、pngおよびgifのドットサイズが6倍になるという問題が生じているという話。

github.com/JuliaPlots/Plo...

も参照。

タグ： Julia言語

posted at 19:27:49

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
PyPlot と Plots の比較

を更新した。Plots.gr()によるGIFアニメを現在のPlots 0.15.0で作るとサイズが6倍になってしまう問題対策の部分。

タグ： Julia言語

posted at 19:24:33

optical_frog @optical_frog

訳したYO: ジョセフ・ヒース「『批判的』研究の問題」 wp.me/p1ydPE-7Zn

タグ：

posted at 19:14:20

@goropikari_ おお、さんくす！

タグ：

posted at 18:47:48

goropikari @goropikari_

@genkuroki 私のプロット例集の Pkg.pin の部分です。
最新に戻す時は Pkg.free を使ってください。
nbviewer.jupyter.org/github/goropik...

タグ：

posted at 18:20:37

@goropikari_ 1行ですませる方法がわからないので無理。

タグ：

posted at 18:18:12

しまリアル @shimarealoki

この記者は去年11月にも、地元をきちんと取材した客観的な記事を書いてますね。
ジャーナリストらしい人は、毎日の本社では嫌われて地方に飛ばされちゃうのでしょうかね？
mainichi.jp/articles/20171...

タグ：

posted at 18:05:49

しまリアル @shimarealoki

ごまふあざらし(GomahuAzaras @MathSorcerer

@RedsHiromitsu ほんとうにそうですね。客観的な報道する記者は、本社から嫌われているんですかねぇ

タグ：

posted at 17:56:53

いいね．Jupyterでmatplotlibのjavascriptアニメーションを動かす on @Qiita qiita.com/yoku_001/items...

タグ：

posted at 17:46:58

goropikari @goropikari_

@genkuroki 6倍サイズ問題はPlots.jlのバージョンを一つ下げるほうが対策としては楽です。1行で済みます。最新版にこだわりがあるなら別ですが。
既知の問題なのでそのうち直ると思います。
github.com/JuliaPlots/Plo...

タグ：

posted at 17:21:46

菅野洋充 @KannoHiromitsu

@shimarealoki これ、記事自体は客観的でなんの偏りもないんだよね。それを名前も出さないツイッター担当にこんな使われ方してかわいそうすぎる…

タグ：

posted at 17:07:42

斉藤ひでみ・現職教師（西村祐二） @kimamanigo0815

@kikumaco @mainichi_dmnd twitter.com/kumappus/statu...

(1) 「これでいいのか」と書いたことの真の意図は何か？

(2) 真の意図がどのように誤解されたと考えているのか？

タグ：

posted at 16:50:30

「校長の苦悩」。すごく分かる。。
中教審の中間まとめで「部活は学校の業務」とされた。とても複雑な気分だった。この位置づけに一番頭を抱えるのは校長だろう。そして校長は、ヒラの教員に部活顧問を押し付けてくるだろう。

教員の部活指導　校長の苦悩　(内田良) news.yahoo.co.jp/byline/ryouchi...

タグ：

posted at 16:48:34

ceptree @ceptree

@tkf 検索かけたらJulia出て来て笑いました

twitter.com/ceptree/status...

タグ：

posted at 16:48:32

非公開

タグ：

posted at xx:xx:xx

岸政彦 @sociologbook

こんにちは、表記ゆれを絶対に直さない市民の会です。

タグ：

posted at 16:23:35

竜田一人 @TatsutaKazuto

今日も「誤解を招いた」で済ましたメディアが出たそうで、再掲 twitter.com/TatsutaKazuto/...

タグ：

posted at 16:12:13

しまリアル @shimarealoki

この元記事書いた松倉展人記者は毎日新聞地方版で加計学園誘致に関する記事を好意的に書いてきた人で、マスコミが加計批判一色になってから今まで干されてた人だったのか。松倉記者がまた干されませんように。
ksl-live.com/blog13542 twitter.com/mainichi_dmnd/...

タグ：

posted at 15:55:58

#Julia言語 Plots.histogramによる表示では各バーの境界に黒い線が入ります。それを消したければ、linealpha=0とすればよいみたいです。

タグ： Julia言語

posted at 15:50:20

#Julia言語 Jupyter notebookを公開するときには

nbviewer.jupyter.org/gist/genkuroki...

の先頭にある diaplayfile函数の方法でノートブック内に作成した画像や動画を埋め込んでから公開した方がよさそうですね。そうしておかないとせっかく公開したのに動画などを見てもらえなくなる。あとでそうしておきます。

タグ： Julia言語

posted at 15:47:42

#Julia言語

Plotsを使う利点：Julia言語における標準になりつつある。本質的に異なる複数のバックエンドを切り換えて利用できる。

Plots.gr()を使う利点：速い。

Plots.gr()の欠点：画像を出力させると6倍サイズになる(対処可能)。日本語フォントを使用できない。

タグ： Julia言語

posted at 15:39:14

#Julia言語

PyPlotを使う利点：細かな調節がやり易い。Python関係の情報がたくさん見付かる。Juliaにおけるプロット法だけではなく、Pythonでのプロット法もマスターできる。

PyPlotがPlotsより劣る点：JuliaではPlotsが標準になりつつある。

PyPlotがPlots gr()より劣る点：遅い。

タグ： Julia言語

posted at 15:36:07

#Julia言語 PyPlot と Plots の両方を同時に試してみたい場合には

using PyPlot
using Plots

などと**せずに**、

import PyPlot
import Plots

として、PyPlot.plotやPlots.plotを使うのがよいと思います。

タグ： Julia言語

posted at 15:33:00

#Julia言語 Plots & gr() における「6倍サイズ問題解決のために gr_display を上書き」の方法は次のリンク先の「Plots gr() backend によるアニメ作成」の節にあります。

nbviewer.jupyter.org/gist/genkuroki...
PyPlot と Plots の比較

PyPlotとPlotsは一長一短なのでどちらを使ってもよいと思います。

タグ： Julia言語

posted at 15:31:35

あ〜る菊池誠(反緊縮)公式 @kikumaco

#Julia言語でGIF動画を作成するときには、

* Plotsパッケージを使う
* gr() バックエンドを使う
* 6倍サイズ問題解決のために gr_display を上書き

とするのが良さそうです。速いです。(pyplog()バックエンドを使うとものすごく遅くなります。10秒未満ですむ仕事が5分以上かかるようになる。)

タグ： Julia言語

posted at 15:29:24

あ〜る菊池誠(反緊縮)公式 @kikumaco

@mainichi_dmnd 誤解ということでしたら、本来の意図をご説明ください

タグ：

posted at 15:27:07

じゃあ、どんな意図だったんだろうねえ。何についてどういう意図でおわびしているんだろう twitter.com/mainichi_dmnd/...

タグ：

posted at 15:25:35

#Julia言語 PlotsパッケージによるGIF動画作成法は、PNGファイルをたくさん作ってそれをまとめてGIFファイルにするという素朴な方法でした。

タグ： Julia言語

posted at 15:25:29

#Julia言語 PyPlotとmatplotlib.animationで1分半程度作成にかかる動画が、Plots GRバックエンドの6倍サイズ問題解決版なら10秒未満で作成できてしまいます。1分半待つのは結構つらいので、10秒ですむことはとてもありがたいことです。

添付動画はその10秒未満で作った動画です。 pic.twitter.com/xU57o1Poc6

タグ： Julia言語

posted at 15:16:52

#Julia言語次のリンク先も更新しておきました。

nbviewer.jupyter.org/gist/genkuroki...
PyPlot と Plots の比較

そこでは Plots.gr_displayの上書きを行って、6倍サイズ問題を解決しています。

Plots gr()バックエンドによる動画作成はとても速いです。
PyPlotを使った場合よりも10倍以上速い。

タグ： Julia言語

posted at 15:14:15

毎日新聞デジタル報道グループ @mainichi_dmnd

#Julia言語 Plotsパッケージをgr()バックエンドで使ってGIF動画を作るとサイズが6倍の動画ができてしまう問題が解決しました。PNGで保存するときに常に6倍サイズになるようになっていた！

github.com/JuliaPlots/Plo...

その部分を削れば問題解決。続く

twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 15:10:20

昨晩投稿したこの記事に関するツイートは、受験生を傷つける意図はまったくありませんでしたが、誤解を招く表現であり、おわびして削除します。

加計学園：獣医学部、一般入試に１０００人超応募 - 毎日新聞 mainichi.jp/articles/20180...

タグ：

posted at 13:55:24

Norio Nakatsuji @norionakatsuji

UFO教授 (藤木文彦 Fumihiko @UFOprofessor

最近世界で注目され流行ってるオーガノイド研究。従来の培養平面に接着した2次元細胞集団ではなく、3次元の細胞集団と組織構造を持つ「小さな擬似臓器」を使って人間臓器に近い状態で調べようとする分野。繰言になるけど、故笹井芳樹氏が世界のパイオニアだったんだよね。誠に惜しい研究者を亡くした。 twitter.com/regmednet/stat...

タグ：

posted at 13:26:46