黒木玄 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 を見て下さい。
2018年01月26日(金)
twitter.com/tanutarou730/s...
#統計 #Julia言語 リンク先のリンク先で渡辺澄夫さん曰く【Gnが大きいときはWnが小さいということです。WnはGnを推定するために作られた確率変数ですから これは、とてもありがたくない、ということを意味】
この件について私は繰り返し重要性を指摘して来ました。続く
posted at 00:45:18
#統計 これは #Julia言語 私のJupyter notebookが渡辺澄夫さんに補足されていた件
watanabe-www.math.dis.titech.ac.jp/users/swatanab...
の「真の分布を知ることができる限界について 」のすぐ上に私の名前が!
見付かってしまったのは
nbviewer.jupyter.org/gist/genkuroki...
混合正規分布モデルと正規分布モデルの各種情報量規準の比較
posted at 00:48:58
ごまふあざらし(GomahuAzaras @MathSorcerer
互いに素の確率を求めるプロラムの件ですが,Juliaの組み込みgcdからナイーブな実装に変更すると3倍遅くなってNumbaの実装がスピードで勝ちます.なのでNumbaが遅いというわけでもないのです.Juliaが優秀なんです.ハイ.
タグ:
posted at 00:55:13
#統計 #Julia言語
n=8 nbviewer.jupyter.org/gist/genkuroki...
n=32 nbviewer.jupyter.org/gist/genkuroki...
n=128 nbviewer.jupyter.org/gist/genkuroki...
正規分布 nbviewer.jupyter.org/gist/genkuroki...
注意深く見るべきなのは添付が画像のようなプロットです。
同じ色が同じサンプルから得た数値のプロットになっています。 pic.twitter.com/04YGSBPt73
posted at 01:12:23
非公開
タグ:
posted at xx:xx:xx
#統計 #Julia言語
nbviewer.jupyter.org/gist/genkuroki...
では単純な正規分布モデルの場合に事前分布をWAICなどで選択する数値実験をやっています。やはり、黄色とシアンの部分が正しい(予測誤差が相対的に小さい)モデル選択に失敗している場合です。やはり、逆相関が失敗の主な原因になっている。
posted at 01:32:47
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
Dr. Chris Rackauckas @ChrisRackauckas
@StephenEglen @azurefriday @Python @AzureNotebooks @Jupyter @shanselman @ch9 @Azure4Research There's #julialang on Azure.
タグ: julialang
posted at 07:14:23
非公開
タグ:
posted at xx:xx:xx
#Julia言語 私の環境で Pkg.update()してから
nbviewer.jupyter.org/gist/genkuroki...
Solving single pendulums by DifferentialEquations.jl
を実行してみました。正常に最後まで動きました。
知的好奇心旺盛な人がJuliaを使うときに注意するべきことが一つあります。続く
twitter.com/adhara_mathphy...
タグ: Julia言語
posted at 09:35:18
@tdkys unicode文字が打てるのビビりますよね。ご存知かもしれませんが、jupyter notebookやREPLからだとLaTeXコマンドに似たものが使えます。βなら\beta + [TAB]で変換できますし、下添字(x₁)ならx + \_1 + [TAB]で打てます。
docs.julialang.org/en/release-0.4... pic.twitter.com/EDRywyLxI0
タグ:
posted at 10:52:29
非公開
タグ:
posted at xx:xx:xx
Ryosuke Kurihara @phys_mmc_2007
In[3]:シュレーディンガー方程式の差分化で書かれている
if dx == 0
v = 2/a^2 + V(x)
elaseif dx == 1
v = -1/a^2
elseif dx == -1
v = -1/a^2
のくだりは,漸化式を行列表示したときに対角項の前後1列目分だけ有限の行列要素を持つことを表しているのか
#Juliaメモ_Rsk
タグ: Juliaメモ_Rsk
posted at 11:17:17
GPU使ってTensorflow始めるぞ!ということで環境構築を始めるも,色々インストールしたりパス通したりcuda導入でnvidiaのドライバを死なせてブラックアウトしたりしてたところ,
RのKerasパッケージ使ったら上記のこと全部自動でやってくれて友人と大爆笑してた。たった3行。
tjo.hatenablog.com/entry/2017/06/...
タグ:
posted at 12:05:50
(´-`).。oO( 深層学習のGPU学習計算,論文の実装を確認するくらいなら,都度課金のfloydhub www.floydhub.com があまりに楽すぎ・気軽すぎて,さくらの高火力はもとよりGCPも使わなくなりそうな勢いです… )
タグ:
posted at 12:06:08
非公開
タグ:
posted at xx:xx:xx
津田和俊/急激に進行した網膜剥離と闘って @kaztsuda
わたしも色々なマスコミの記者さんとは面識があるのでフォローしますけど、細かに取材をして声を拾い上げてくれる丁寧な仕事をなさる記者さんもいれば、人の話を聞いてるのかと思いきやしゃべってもいないことを平気で書く記者さんもおり、どの新聞だからどうとも言えないように思います。
タグ:
posted at 12:21:01
津田和俊/急激に進行した網膜剥離と闘って @kaztsuda
まぁ、取材される側としては、記者さんを選べないのが困ったもので
(こちらから「指名」する制度があればいいんだべか?、でも指名料、お高いのかしら??、って、何の店の話や?!)
タグ:
posted at 12:23:40
数学の人がみんなこうだというわけではないし、奇行も基本的には無害なので、慣れたらほっといてよい。 twitter.com/genkuroki/stat...
タグ:
posted at 13:23:32
ソフトウェア工学というか、開発手法みたいなもの、テストはこうしましょうとかデザパタとかは、いわゆる典型的なソフトつまり経理システムとかOSとかそういうものを作ることが念頭にあるのであって数値計算みたいなものはほぼ考えられてない、ということが最初理解できなくて混乱の時代が長かった…。
タグ:
posted at 13:36:47
この点、OOPパラダイムも本当に混乱したなあ……。当時はなんかOOPこそが最強の開発手法なんだ!みたいに言っている人が(ネットで)多くて、数値計算を完全にOOP的にやるの色々大変だけどこんな専門家の人が自信を持って言っているのだからきっとそういうことは完璧に考慮されているはず…と悩んだ。
タグ:
posted at 13:45:21
「chopper1」という、薪を割りやすくする為にインパクトの瞬間左右にアームを展開するスプリッターレバーを搭載した『機械式の斧』が存在する事を知って、無駄にワクワクしている。
youtu.be/1AmB1rF9Mlg pic.twitter.com/kNcRrfj06Y
タグ:
posted at 15:33:01
JuliaのDistirbutions.jlにある多変量正規分布(mvnormal)が、共分散行列は半生定値なのにエラー吐いて死ぬのは固有値分解じゃなくてコレスキー分解を使っているからなのか。。。これは。。。 discourse.julialang.org/t/multivariate...
タグ:
posted at 15:33:40
すでに何時間もたってるのだけど,どうなったかな…
こういうときは,まずは志願する大学の事務に連絡するのがベスト.インフラに大きな影響が出ていることが明らかなので普通は救済的な対応をするはず.と,どっかの大学の入試委員長が言うとったで. twitter.com/plus6etc/statu...
タグ:
posted at 15:51:30
ウェブアーカイブを確認したら、共同通信社の最初の記事も無署名でした。すみません。「松吉」という人が記事の著者なのか、ツイッターの中の人なのか、わかりません。いくつか消しておきます
タグ:
posted at 15:56:55
全国の受験生とそのご家族に伝えたいのだけど,とにかく受験関係で何かあったら募集要項とか受験票とかに書いてある大学の連絡先に電話で一報入れてください.
大学の入試関係者は教員も事務職員も決して冷酷な機械ではありません.正当な理由があればそんな簡単に見捨てたりしません.
タグ:
posted at 16:00:03
非公開
タグ:
posted at xx:xx:xx
うーん、日本の家の断熱性能の低さがいろいろ指摘されているけど、研究によれば省エネの建築基準は床面積を減らすので狭い家(=貧乏人and/or日本)にはきわめて負担が大きいし不動産価値も下げるし云々。するといまの日本の狭い家は、断熱がダメなのがむしろ適応してるのかも…… twitter.com/BjornLomborg/s...
タグ:
posted at 17:45:02
ごまふあざらし(GomahuAzaras @MathSorcerer
Hello World をprintfするcのコード hello.c から
clang -emit-llvm -S -o hello.s hello.c によって
Hello.s を出した様子. pic.twitter.com/xnlqUZvF6i
タグ:
posted at 18:51:20
ごまふあざらし(GomahuAzaras @MathSorcerer
$ llvm-as -o hello.bc hello.s
$ lli hello.bc
Hello World
タグ:
posted at 18:52:08
ごまふあざらし(GomahuAzaras @MathSorcerer
以上あざらしでもできたLLVMでした.(参考文献:きつねさんでもわかるLLVM Chap3より)
タグ:
posted at 18:53:27
“これで共同通信の説明責任が果たされたと感じる読者は、果たしてどれほどいるだろうか” / “共同通信、山中氏の記事全面書き換え「編集上必要と判断」 - ITmedia NEWS” htn.to/XywCNfb
タグ:
posted at 19:09:26
tsujimotter 日曜数学者 @tsujimotter
131は素数だし、今日はプライムフライデーですね!
ちなみに p=131 は 1/p が p-1 桁の循環節を持つので、ダイヤル数を作る素数でもありますね!(画像は1/131の作るダイヤル数に 1 から 130 をかけたもの) pic.twitter.com/W7Z5ORd0Cj
タグ:
posted at 19:12:44
tsujimotter 日曜数学者 @tsujimotter
ダイヤル数であることを確認できるように,テキストファイルをGistに共有しました
gist.github.com/junpeitsuji/38...
タグ:
posted at 19:29:19
ごまふあざらし(GomahuAzaras @MathSorcerer
Nim で二つの自然数が互いに素の確率から円周率を導出するベンチマーク.
Nim(標準ライブラリのgcd 使用が遅く
自前で組んだ方がなぜか速い.それでも
Python+Numbaで組んだ方が高速に動く. pic.twitter.com/WI1syGmQlX
タグ:
posted at 20:05:34
ごまふあざらし(GomahuAzaras @MathSorcerer
あーもちろんJuliaの実装がダントツでMacBook12inchの結果 Julia 4秒ほど Python+Numbaが10秒ほど Nimが13秒ね.
Nim推し派だった私としてはとてもつらい事実. twitter.com/mathsorcerer/s...
タグ:
posted at 20:09:26
非公開
タグ:
posted at xx:xx:xx
ごまふあざらし(GomahuAzaras @MathSorcerer
そこに言語があるからです twitter.com/umikoz/status/...
タグ:
posted at 20:34:51
勝手に高専ラジオ「初社会人ゲスト!!」 youtu.be/pasKGNCw30E via @YouTube
東北で高専生のための塾を運営されているあけまつさん(@minami106)がゲストで来てくれています.
社長さんですよ!!しゃちょう!!
まさかこんなことになるとはね.
itunesの方でも同時に配信しております.
タグ:
posted at 20:52:41
補足: #Julia言語 のgcdが定義されている場所は @ which gcd(4,6) でわかります。Jupyter notebookでやるとGitHubリポジトリへのリンクが表示される。
github.com/JuliaLang/juli...
二進表示での下の方の桁の0の扱いは自明なのでガンガン削って計算する方法。
twitter.com/mathsorcerer/s...
タグ: Julia言語
posted at 21:07:00
速いとよく言われているJulia言語は gcd の類の基本的な函数の多くがJulia自身で書かれています。速さが欲しいので基本函数をCとかでゴリゴリ書きまくるのではなく、Juliaで書いている。
ソースを見て気付いたときにはびっくりしました。今は慣れたけど。
タグ:
posted at 21:10:00
#Julia言語 の基本的な数学的函数の多くがJuliaで書かれており、Juliaの文法はシンプルでごたごたした修飾が少なめなので素人でも読めます。
Jupyter notebook (juliabox.com で使える)で @ which gcd(4,6) (@ の後の空白は除く)とやってみて下さい。
タグ: Julia言語
posted at 21:13:27
今 #Julia なハッシュタグ打とうとしてたら、 #JuliaProlog って出てきたけどそんなエエもんがあるのか?
タグ: Julia JuliaProlog
posted at 21:33:04
そんなことを思って昨夜はJulia0.5.2の方で(まぁ実のところそこが一番怪しいのだけど。だって何処にどういう具合にPkg.addされてるのか完全にお任せ状態だから)やってみたけど動かなかった…シクシク。一応半年以内にメンテされてるものとかってチェックはしてるのだけど。 twitter.com/genkuroki/stat...
タグ:
posted at 21:38:29
開発中で発展途上の言語だから仕方のない部分だし諦めてるっていうかこの段階でこれほど使えることに感謝してるぐらいだけど。これとこれに類する「分からなさ」が解消されたら完全にrubyから乗り換えるのだけど twitter.com/genkuroki/stat...
タグ:
posted at 21:40:11
@genkuroki @which gcd(4,6)
としたら,
gcd(a::T, b::T) at base/intfuncs.jl:31
と返ってきました。これはintfuncs.jlをエディタで開いて見てみろってことですか?
タグ:
posted at 21:45:09
@genkuroki binary GCD algorithm ってまさに数Aの整数にドンピシャな話やん。未だ証明とかは読んでないけど。こういうのをサクサク教材にしてってのが本来の教育研究やその成果を発揮せんとあかん手やないのか。何かと言うと教育学系は手の内を秘匿する。自信の欠如と内輪受けなんやろけど
タグ:
posted at 21:48:39
@genkuroki zwxadz.hateblo.jp/entry/2015/04/...
に2進最大公約数算法の解説があるけど,「nにn-m」をってそれまんま「互除法」やないかと思うのは何か間違ってるのか?
タグ:
posted at 21:54:17
#Julia言語 binary gcd のアルゴリズムをJuliaに付け加えたのは Elliptic.jl (楕円函数パッケージ)の作者でもあります。
github.com/JuliaLang/juli...
Pull request: faster gcd #8410
github.com/nolta/Elliptic...
Elliptic Special Functions for Julia
楕円函数のパッケージにはお世話になっています。
タグ: Julia言語
posted at 21:58:23
#Julia言語 のgcdを高速化し、楕円函数のパッケージ**など**も作っているnoltaさん曰く:
nolta: Sure, but i'm lazy.
それに突っ込みが入っていた(笑)。
StefanKarpinski: You have a strange definition of lazy.
😅
こういう雰囲気は楽しくてよいです。
github.com/JuliaLang/juli...
タグ: Julia言語
posted at 22:04:00
#JuliaLang pljulia 発見。何やらPrologの方言?なSWI Prologならしい。28日前にREAD.meが更新されてる。せやけどpack_install()って何やねん。pkg.add()とどう違うねん…
タグ: JuliaLang
posted at 22:04:10
ごまふあざらし(GomahuAzaras @MathSorcerer
これは笑ってしまいましたね。 twitter.com/genkuroki/stat...
タグ:
posted at 22:18:47
#超算数
twitter.com/kanotosngk/sta...
これは、「掛け算の順序」や「正方形は長方形ではない」を凌駕する事例。
教師の採点がおかしい(採点ミス?)であって、算数教育界中枢がおかしくなっているのとは別
算数教育界中枢=超算数正規軍
おかしな指導・採点=超算数遊撃隊
タグ: 超算数
posted at 22:47:41
#超算数 この事例とどっちがすごいかな?
5人に3本ずつ、を3×5としたらバツ。
「掛け算の順序は問題文に出てきた順だから5×3が正しい」らしい。
掛け算順序強制の算数教育界中枢と、我々順序反対派、双方に喧嘩を売る怖いもの知らず^^
6828.teacup.com/amajima/bbs/229
タグ: 超算数
posted at 22:50:31
@sekibunnteisuu twitter.com/kanotosngk/sta...
この例は、凡ミスなんじゃないかなぁ。
つまり
「掛算のありがたみを教えるために、いったん面倒な足し算をさせるという意図だから、ここは足し算でやらせる。」
と思い込んで、「かけ算でやれ」と書いてあること自体を見落としたとか。 #掛算
タグ: 掛算
posted at 22:54:21
Fermat's Library @fermatslibrary
Here's a very curious property involving √2. Can you explain why this is happening? pic.twitter.com/4fAtmv9MZa
タグ:
posted at 22:58:02
@golgo_sardine 仮に採点ミスだとして、超算数擁護勢力が「採点ミスに決まっているだろうが。常識で考えたら分かるだろうが!」とか言ったら大笑いw
常識で考えられない採点をしているのだから、これも採点ミスじゃないと思われても仕方ない
タグ:
posted at 23:02:15
ごまふあざらし(GomahuAzaras @MathSorcerer
Steinのアルゴリズムで勝負したらNumba+ナイーブなgcdに追いつきました.よかった・・・.
qiita.com/SatoshiTerasak...
タグ:
posted at 23:05:36
#Julia言語 PyPlot 小ネタ。
using PyPlot
plot(sin)
とはできないが、添付画像のようにして plot 函数を拡張すればできる。
gist.github.com/genkuroki/8948... pic.twitter.com/ErAF7Jp38w
タグ: Julia言語
posted at 23:06:55
@tsujimotter Rayleighの定理より従います。差が2の倍数全体になるものに限らず、自然数nの倍数になるような無理数についても直ちに一般化ができます。
ja.wikipedia.org/wiki/%E3%83%AC...
タグ:
posted at 23:30:55
ごまふあざらし(GomahuAzaras @MathSorcerer
Rust (gcd はナイーブなバージョン)で組んだらJuliaに勝ちました.
あー満足満足.
タグ:
posted at 23:58:49