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

#数楽これはもしかしたらその手の仕事に向いている人にとってはチャンスかもしれず、徹底的に最適化された汎用特殊函数ライブラリを作って配布すれば有名になれるかも。

そういう仕事は連分数による函数近似を理解できないと無理。プログラミングにおける最適化以前に数学に詳しくないと無理。

タグ：数楽

posted at 23:57:21

#数楽 Python scipyのようなメジャーなパッケージで採用されている特殊函数ライブラリは徹底的に最適化されていると思っていたのですが、Julia言語を使ってあっさり5倍程度速くできるのを見ると、必ずしもそうではないようです。おそらくユーザーがあまり多くない函数は十分に最適化されていない。

タグ：数楽

posted at 23:55:16

#数楽 Python scipyではFortranで書かれた特殊函数の数値計算用のライブラリを使えるのですが、MITの講義ではJulia言語で指数積分函数を実装して、そのライブラリよりも5倍程度の速くできることを実演しています。最適化は偉大。Julia言語は将来的にはFortranの代わりに利用されるかもしれません。

タグ：数楽

posted at 23:53:08

#数楽リンク先のJupyter notebookは In [134] にバグがあって、一ヶ所 E₁_cf を E₁_cf64 に直さないとダメ。2017年度のMITでの講義の一部分です。

nbviewer.jupyter.org/github/steveng...

これ滅茶苦茶面白いです。特殊函数の数値計算のコードの最適化の仕方がよくわかります。

タグ：数楽

posted at 23:49:30

三乗根 @cubic_root3

これ吹奏楽の強豪で有名な市船高校の顧問だ。
まあここに限らず強豪校はアホみたいに長時間練習させてますよ。
吹奏楽部はかなり闇が深い。 twitter.com/soverro/status...

タグ：

posted at 23:47:20

#数楽 #連分数 #数値計算 #Julia言語

特殊函数の数値計算の仕方をググればわかることですが、連分数は最も基本的な手段になっています。

連分数とテイラー展開の併用で指数積分函数の数値計算をJulia言語でする話はMITの数値計算の講義でも扱われています。

nbviewer.jupyter.org/github/steveng...

タグ： Julia言語数値計算数楽連分数

posted at 23:47:02

@kuri_kurita

ゴルゴ・サーディーン @golgo_sardine

「かけ算の順序」をマスターさせないと子供は「割り算の順序を間違える」だとかその変種を言う人が良くいるけど、全く子供をバカにした話だし、第一、（恐ろしい事に）何の証拠もないだけでなく事実に反するわけで、単なる「信仰の告白」とか「教義の暗唱」の類。

タグ：

posted at 23:46:52

「#掛算の順序はどうでも良い」というのは、俗世間の大人なら当たり前のこと。
彼ら順序派にとって勉強とは単元ごとにブツ切りな物なので、教育学部で順序固定を教えられた時、「え？　自分の日常の中ではこんな事しないぞ？　子ども達にこんな事を守らせる必要があるのか？」と振り返らない。（続

タグ：

posted at 23:30:09

koh_t @koh_t

俺達のメディアアート(言ってみたかっただけ)
koh_t
とある近接勾配法アルゴリズムの学習曲線 / a learning curve of a proximal gradient descent method.
Mac, Julia, Jupyter notebook
2017
#締切からの逃避 pic.twitter.com/twwLTJP10M

タグ：締切からの逃避

posted at 23:24:37

天むす名古屋 Temmus @temmusu_n

やべーダウンロードおせー

タグ：

posted at 23:24:23

#掛算の順序を入れ替える癖を矯正しないと #割算の順序を間違えるという #超算数言説がまた聞かれるのでいつものやつを
twitter.com/temmusu_n/stat...。割算の順序はみんなそんなに間違えません。

タグ：割算掛算超算数

posted at 23:24:01

Julia ProとMKLは何が違うんだろう？
Pro⊂MKL？

タグ：

posted at 23:21:03

オースティン @keen339

@tmnghryk @ochiyan5 なんて事！！

タグ：

posted at 23:20:39

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

#Julia言語のチュートリアル：

nbviewer.jupyter.org/github/bicycle...

nbviewer.jupyter.org/github/bicycle...

プロットについては

goropikari.github.io/JuliaPlotsGall...

nbviewer.jupyter.org/github/goropik...

「あれ？遅い？」と思ったら

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

myenigma.hatenablog.com/entry/2017/08/...

を見れば何とかなるはず。

タグ： Julia言語

posted at 23:19:51

BinaryBuilder.jlマスターした

タグ：

posted at 23:17:37

てらモス @termoshtt

連分数のカオス解析での有効性は
散逸力学系カオスの統計力学 (非線形科学シリーズ)
でも触れられています(/・ω・)/
www.amazon.co.jp/dp/4563023450

タグ：

posted at 23:17:05

せっかくなのでJuliaProにしてみようかと。

タグ：

posted at 23:16:07

koh_t @koh_t

すごい勢いでねじれる学習曲線その2。こういうのメディアアートって言っても良いんじゃないんですかね？（炎上案件 #SorkScrewLearningCurve part2 pic.twitter.com/0I8L10F6lR

タグ： SorkScrewLearningCurve

posted at 23:14:53

非公開

タグ：

posted at xx:xx:xx

koh_t @koh_t

今日のバグ学習曲線はコークスクリューのようになった pic.twitter.com/yZZcaTPzWH

タグ：

posted at 23:11:38

#Julia言語で

Pkg.add("DifferentialEquations")

を実行すると微分方程式関係の巨大パッケージ群がインストールされる。Julia言語の微分方程式パッケージはあまりにも巨大でびっくりしまくった。誰かが網羅的な紹介をしたら尊敬する。とにかくでかい。

github.com/JuliaDiffEq/Di...

タグ： Julia言語

posted at 23:07:38

まわしてみた。途中で脳がバグってくる。 pic.twitter.com/2TLKuRqHae

タグ：

posted at 23:03:09

#Julia言語～書いてあります。私の環境では

C:\Users\username\AppData\Roaming\jupyter\kernels

のようなディレクトリに Jupyter notebook のカーネルの設定が保管されている。その中身をよく観察して、必要ならフォルダごとコピーを作って内容を適切に書き換えれば終了です。

タグ： Julia言語

posted at 23:01:46

#Julia言語続き。その上書きをもったいない(例：v0.6.0とv0.6.2の違いをJupyter notebook上で確認したいなどの)場合には手動でv0.6.0のJupyter notebookへの登録を復活させることができます。その方法は(Windowsのケースなのでフォルダ名が全然違うでしょうが)、

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

に～続く

タグ： Julia言語

posted at 23:00:21

@genkuroki ありがとうございます。いつも参考にさせてもらっています。

タグ：

posted at 22:59:26

twitter.com/adhara_mathphy...

#Julia言語 Julia v0.6.*がJupyter notebookのkernelとして登録されているときに、Julia v0.6.2を入れて、さらにJulia v0.6.2の方でPkg.add("IJulia")すると、Jupyter notebookの方の Julia v0.6.*のkernle登録をv0.6.2で上書きする仕様になっています。続く

タグ： Julia言語

posted at 22:57:56

twitter.com/k_b_ms/status/...

#Julia言語

Windows上の話に限定しますが、Julia言語関係のファイルは邪魔になったら、単に消せば大丈夫な感じです。

全部消して最初からやり直して、クリーンな環境を作りたい場合にはJulia関係のファイルは全部消してよいと思います。

タグ： Julia言語

posted at 22:50:46

shin @Matz_Shin

@genkuroki ご教授ありがとうございます。前職で熱方程式を、陽解法、陰解法、QR分解と勉強しながら解いた事があるのですが、FFTは使わなかったので、思わず聞いてしまいました。

タグ：

posted at 22:50:39

@Matz_Shin FFTの単純な利用は周期境界条件以外の境界条件の場合には苦しいです。

タグ：

posted at 22:45:12

twitter.com/k_b_ms/status/...

#Julia言語私はいつも Jupyter notebook 上で Pkg.add などをしています。Jupyter notebookに作業の記録を残せるので便利です。

タグ： Julia言語

posted at 22:43:46

shin @Matz_Shin

@genkuroki この方法って、領域の形が複雑なときでも使えますか？

タグ：

posted at 22:43:28

面倒なのでJuliaだけアンインストールしたものの、残ったフォルダはゴミなような。これ全部消していいんだっけ？

タグ：

posted at 22:40:38

#Julia言語もとの函数がガウシアンのケースではFFTを使った高階の導函数の数値計算が極めて優秀に見えた。区間を16個で区切っただけで目が超絶荒いのに誤差が非常に小さい。 #数楽 #数値解析

タグ： Julia言語数値解析数楽

posted at 22:39:49

#Julia言語目が粗い場合の画像 4/4

g(x) = exp(-x^2/2) の3階の導函数の2通りの数値計算。
左側はFFTを使って計算した導函数で、右側は

(g(x+3h/2) - 3g(x+h/2) + 3g(x-h/2) - g(x-3h/2))/h^3

のプロット。 pic.twitter.com/DDH3Em6osO

タグ： Julia言語

posted at 22:37:52

#Julia言語目が粗い場合の画像 3/4 (訂正一つ前のやつは2/4)

g(x) = exp(-x^2/2) の2階の導函数の2通りの数値計算
左側はFFTを使って計算した導函数で、右側は

(g(x+h) + g(x-h) - 2g(x))/h^2

のプロット。 pic.twitter.com/ecOnewVcY3

タグ： Julia言語

posted at 22:35:40

#Julia言語目が粗い場合の画像 1/4

g(x) = exp(-x^2/2) の1階の導函数の2通りの数値計算
左側はFFTを使って計算した導函数で、右側は

(g(x+h/2) - g(x-h/2))/h

のプロット。 pic.twitter.com/WdzgzoKLfB

タグ： Julia言語

posted at 22:34:28

#Julia言語目が粗い場合の画像 1/4

函数 g(x) = exp(-x^2/2) の荒い離散化のグラフ pic.twitter.com/qk33k2ubtV

タグ： Julia言語

posted at 22:32:45

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
FFTを用いた熱方程式やKdV方程式などを数値解法

に「数式を使った解説」と「目の粗い場合」を追加しました。

FFTを使った偏微分方程式の解法の基礎はFFTして掛算してIFFTでもとに戻すと微分の数値計算になるようにできることです。その辺の仕組み関係の追加。

タグ： Julia言語

posted at 22:26:52

twitterの制限が15 MBなので、ImageMagickで「convert input.gif -fuzz 10% -layers Optimize output.gif」とやると20.5 MBが355 KBになった

How can I compress the size of gif images with ImageMagick?
stackoverflow.com/questions/8578...

タグ：

posted at 22:18:37

Paul Painlevé @Paul_Painleve

@subarusatosi @neet2go ガウス　数論論文集
www.chikumashobo.co.jp/product/978448...
で翻訳されてますが、Werke2巻は全部ここで読めます（翻訳5論文はラテン語です）
archive.org/details/werkec...

タグ：

posted at 22:10:29

となることです。

タグ：

posted at 22:06:07

となり、sagemath8.0を起動すると自動的にjupyter notebookが呼ばれ、 pic.twitter.com/7P6iSGvv9M

タグ：

posted at 22:04:35

両眼視差を用いたステレオグラム(平行法)。せんしちさん(@tdkys)の下記のツイート参考にさせて頂きました。公開する許可も頂いてます。上２つがずらしなしで、下２つが5°ずらしてます。

せんしちさんのツイート↓
twitter.com/tdkys/status/8...

ステレオグラムはこちら↓
ja.wikipedia.org/wiki/%E3%82%B9... pic.twitter.com/4e8s08bYLQ

タグ：

posted at 22:03:41

となって、
シェルからjupyter notebookの時は、 pic.twitter.com/HcLW44ZLt0

タグ：

posted at 22:02:01

鰹節猫吉 @sunchanuiguru

RT 左派は、成熟社会だからもう経済成長しなくていいなんて言ってる人は排除して、貧困をなくすために経済成長をめざす政策を提示しないといけないのですが。そうでなければ、恵まれていない人に目を向けているとはいえない。

タグ：

posted at 22:00:10

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

jupyter notebookについては（私が）わかっていないことがまだあって、
julia内部からusing IJulia→notebook()の時は、 pic.twitter.com/W81pu5C3XE

タグ：

posted at 21:59:51

今年、働き方改革も部活問題も、いよいよ大詰めを迎える。

スポーツ庁は、先日のガイドラインで部活顧問辞退を阻止しようとしている。
文科省は、給特法改正に踏み込まれないように、あらゆる防波堤を築くだろう。

教員に対する人権侵害に、どう立ち向かうか。

悔いのないように、挑みたい。

タグ：

posted at 21:49:16

旧juliaとjupyterの縁が切れた仕組みはよく分かっていません。

タグ：

posted at 21:47:48

JuliaPro（0.6.2）をインストールしたらJupyterからはこちらはがすぐ使えるようになりました。
一方、旧Julia(gitからコンパイルしたもの、0.6.0)はKernel一覧に出て来なくなりました。
旧旧Julia（homebrewからコンパイル、0.5.1）は一覧に残っています。 pic.twitter.com/tQzfnoeOfh

タグ：

posted at 21:37:03

@genya0407

しまじろうよう社会常識身に付け中 @you_s1025

Julia良さそうやな。

タグ：

posted at 21:07:49

ガウシアンカーネルを使ったニ値分類の実装を Julia で試してみた。先週は結果がおかしくてストレスフルなまま終了したけど今回はちゃんと再考して思い通りに動くようになったのでかなり気持ち良い。平面(二次元)で試して大体の感覚が分かったので一般次元に拡張しても特に問題は無いはず。いいね

タグ：

posted at 20:28:39

春 @sunbluesome

いつだったかPythonとJuliaでODEソルバの速度比較をした話
hotoke-x.hatenablog.com/entry/2017/12/...

タグ：

posted at 20:13:44

Akinori Ito @akinori_ito

仙台行きが来ないので、もう金沢か越後湯沢に帰りたい

タグ：

posted at 19:57:36

いま「Julia(Jupyter notebook)で微分方程式の数値解法」って教科書だしたら俺に売れるで

タグ：

posted at 19:57:29

tmnghryk @tmnghryk

@ochiyan5 研究費が足りないというより、事務職員が粛清され、間接経費がないので、機器がが交換できない、本が買えない、膨大な書類を教員が書かないといけない、冷暖房を切られる、出張がしづらい、授業担当が増える、とか時間が奪われることが大きいです。研究費からは一定額までしか流用できないし。

タグ：

posted at 19:38:53

ochiyan @ochiyan5

@tmnghryk 批判ではなく純粋な疑問なのですが、数学はお金が無いと研究できないのですか？哲学や数学はお金がなくとも研究できる印象があります。

タグ：

posted at 19:33:33

#Julia言語あと、繰り返し紹介している例ですが、以下のリンク先で公開されているMITでの講義のJupyter notebookは、内容的には普遍的なのですが、Julia言語抜きで講義するのはちょっと難しいと思います。

github.com/stevengj/18S09...

タグ： Julia言語

posted at 19:29:55

#Julia言語数学的に基本的な部品まで戻って計算の仕方をJulia言語のコードを示すことによって説明している素晴らしい例に

machine-learning.hatenablog.com/entry/2018/01/...
→ github.com/sammy-suyama/M...

と

github.com/sammy-suyama/B...

があります。「素直な書き方で速い」というJulia言語の特徴がないとこういうのは無理。

タグ： Julia言語

posted at 19:25:47

goropikari @goropikari_

@genkuroki Webサイト版もあるのでよかったらどうぞ
goropikari.github.io/JuliaPlotsGall...

タグ：

posted at 19:01:05

非公開

タグ：

posted at xx:xx:xx

#Julia言語

既成のパッケージをブラックボックスとして使うことに心理的抵抗がある人や、数学的に基本的な部品まで戻って計算の仕方を説明したい場合には、Julia言語 + Jupyter notebook の方がよいと思う。

ブログ記事を書くよりJupyter notebookを書く方が楽、という人は結構多いはず。

タグ： Julia言語

posted at 17:14:23

#Julia言語

Rの世界には統計処理に関するありとあらゆるパッケージが揃っているので、通常のケースでは統計処理をしたいときにはRを使うのが合理的だと思います。

続く

タグ： Julia言語

posted at 17:09:57

すげー・・・． twitter.com/genkuroki/stat...

タグ：

posted at 17:07:23

twitter.com/goropikari_/st...

#Julia言語素晴らしいです。こういうのが欲しかった！

タグ： Julia言語

posted at 17:05:00

ʇɥƃıluooɯ ǝıʇɐs @tsatie

#JuliaPro 無事にインドールもできた。昨夜の問題の @genkuroki さんのコードも無事に動かせて確認できた♬ さてはて Sam Loyd のディオファントスなパズルは解けるやろか。いやどないして解くのかか。総当たりは辛い。

タグ： JuliaPro

posted at 16:22:00

Haruhiko Okumura @h_okumura

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

@kikumaco Juliaまだよくわからないので，ぜひ例題のJulia版をネット公開してください

タグ：

posted at 15:54:10

奥村さんの本の例題を片っ端からJuliaで試してみればいいのかな。しかし、さまざまな検定などはRのほうが便利なんだろうなあ

タグ：

posted at 15:45:09

Yuki Nagai @cometscome_phys

須山敦志 Suyama Atsushi @sammy_suyama

Juliaでトポロジカル物質のハミルトニアン書いてそれに壁を入れて対角化して壁状態を見ることができるのって楽だ

タグ：

posted at 15:00:42

須山敦志 Suyama Atsushi @sammy_suyama

これも好き。
It is better to be vaguely right than exactly wrong. (Carveth Read)
厳密に間違えるよりおおよそ合っている方が良い。

タグ：

posted at 13:46:33

これはぜひ多くの人に知っていただきたい。ある意味、純粋数学とエンジニアリングの相違点とも見れそうです。 twitter.com/TJO_datasci/st...

タグ：

posted at 13:32:34

goropikari @goropikari_

Plots.jl を使ったプロット例

nbviewer.jupyter.org/github/goropik...

#julialang

タグ： julialang

posted at 13:25:18

TJO @TJO_datasci

"All models are wrong, but some are useful"教を打ち建てて布教して回ろうかと最近考え始めてゐる en.wikipedia.org/wiki/George_E....

タグ：

posted at 13:05:32

ななみ @aheuv

@genkuroki @GoITO 仕事で文科省の方と話す事がありますがとんでもない人か相手の話しを聞かない人が多くて嫌になります

タグ：

posted at 12:18:38

ʇɥƃıluooɯ ǝıʇɐs @tsatie

いいなぁ。こういう風に使うのか。ふむふむ。for だけなんか浮いてる感じがするなぁ。仕方ないのか。 twitter.com/genkuroki/stat...

タグ：

posted at 12:03:05

らじうむ小山_PPPMP @Ra_koyama

UFO教授 (藤木文彦 Fumihiko @UFOprofessor

@sekibunnteisuu @genkuroki 掛算順序は信仰になってるんだなぁ…

タグ：

posted at 11:45:42

@genkuroki 人の無能に付け込む人には腹が立ちますね。自分の利益しか考えていない。人間社会がより良くなるためにはそうした無能者を減らすとか救済するとか考えるべきと思います。ここが、教師や研究者とお役人や富裕層との違いなんでしょうね。

タグ：

posted at 11:30:30

どんなに勉強しても知らないことの方がたくさんあります。ある種の人達はそういう無知につけこんでやりたい放題やって来ます。

タグ：

posted at 11:05:47

「文科省けしからん」と言いつつ、財務省が言っていることを根拠に「国は借金で大変だ」と思っているとしたら、その知的無能さはちょっと救われないレベルかも。

twitter.com/fly15f/status/...

タグ：

posted at 11:02:32

非公開

タグ：

posted at xx:xx:xx

#Julia言語でのテストコードは

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

のIn [2]。見ればわかることですが、型宣言は一切していないし、fastmath inbounds simd などの最適化系マクロも使用していません。(アルゴリズムの工夫も十分にしていないのでコードは短い。)

twitter.com/MathSorcerer/s...

タグ： Julia言語

posted at 10:56:58

大学関係者で「文科省はけしからん」と「国も借金で大変だから仕方がない」の組み合わせで考えている人達はできるだけ早く自分自身が知的に無能であったことを知るべきだと思う。国が借金で大変？根拠は？もしかしてテレビで「国の借金が1000兆円」と言っているのが根拠？まさに知的無能の典型例。

タグ：

posted at 10:51:02

#Julia言語初等整数論的もしくは組み合わせ論的な実験をコンピューターで気軽にやりたい人にはJulia言語は向いていると思います。なぜならば、基本文法が素直でかつ速いから。

タグ： Julia言語

posted at 10:33:42

#Julia言語での解答例

lst = [n for n in 1:2017 if gcd(n,30) == 1]
length(lst), sum(lst)
→ (538, 542708)

twitter.com/tsatie/status/...

タグ： Julia言語

posted at 10:02:27

Haruhiko Okumura @h_okumura

ゴルゴ・サーディーン @golgo_sardine

brew cask install julia でバイナリを簡単にインストールできた。ただし/Applicationsに入れるので一般ユーザではなくsudoerであることが必須 twitter.com/h_okumura/stat...

タグ：

posted at 09:41:45

非公開

タグ：

posted at xx:xx:xx

こうして #掛算の順序問題が多くの人の知るところとなっても、順序派は、「ＴＶで怪しい医療情報を扱った翌日の医師の立場」と同じつもりでいるのだろうな。

タグ：掛算

posted at 08:37:34

砂___の___女 @vecchio_ciao

@LimgTW 【テクニックとして「順番はどうでもいい」と教えるのは違うだろ】

この解釈は斬新ですね。

＃掛算

 twitter.com/busido0531/sta...

タグ：掛算

posted at 08:03:16

かなか @kanakanaka_a

数学者に義務教育を語らせるのがお門違いって言ってる人多いけど、デタラメな「算数のきまり」を何の疑問もなく受け入れてる「すうがくなにそれ」な層に #超算数の問題を伝えるのに「えらいせんせいのことば」は効果的なんだと思う。 #初耳学

タグ：初耳学超算数

posted at 08:00:21

かなか @kanakanaka_a

うちの子の担任ははっきり言った。「算数力よりも人間力をつけたい」と算数の授業のなかで。そして根拠のないルールを押し付け、従わない子を否定した。思考停止して指示に従うのがこの教師のいう「人間力」。 #超算数

タグ：超算数

posted at 07:40:07

非公開

タグ：

posted at xx:xx:xx

春 @sunbluesome

そこでJuliaですよ twitter.com/ktrmnm/status/...

タグ：

posted at 07:16:12

きゅい　#Julia言語 twitter.com/MathSorcerer/s...

タグ： Julia言語

posted at 05:19:26

Ising Model で niter = 1e9 での実行時間を計測しましたので共有します． pic.twitter.com/LloDgLUwQM

タグ：

posted at 05:18:53

君のRustは20倍遅い - 簡潔なQ qnighy.hatenablog.com/entry/2017/05/...

タグ：

posted at 04:34:47

Dr. Chris Rackauckas @ChrisRackauckas

IsingModelでNim Go の実装が遅いと言ってしまってたけれども，アルゴリズム内でのループの設定が間違っていたため遅いという結論を出してしまいました・・・．しかしながら修正した実装での計測の結果，結論は不変でした・・・．#Julia言語　ぱやい．

タグ：

posted at 04:10:02

DifferentialEquations.jl in #julialang will be making some big breaking changes. Blog post coming soon, but here's a sneak peek. Last call for options:

github.com/JuliaDiffEq/Di...

タグ： julialang

posted at 03:49:40

Julia Bloggers @juliabloggers

Thread Parallelism in Julia www.juliabloggers.com/thread-paralle... #julialang (Re-post)

タグ： julialang

posted at 02:41:43

Yamanami @yamanamitakeshi