黒木玄 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年10月31日(土)
GCCの最適化のバグの紹介。これ自体、今年の5月の報告だが、1999のバグの紹介も面白い。
これは最適化等価性問題の一つだが、等価性が正しく実装されていることを示すのは難しい。 twitter.com/hashbreaker/st...
タグ:
posted at 00:04:51
भूकंप आने के बाद भयंकर सुनामी तुर्की और ग्रीस दोनो की हालत खराब
#earthquake #Tsunami pic.twitter.com/0mJjPFPfkY
タグ: earthquake Tsunami
posted at 00:19:33
非公開
タグ:
posted at xx:xx:xx
Numerics in #JuliaLang and pure algebraic geometry don't make a good couple? Wrong! HomotopyContinuation.jl can provide rigorous proofs using interval arithmetic. Preprint is here: pbrdng.github.io/certification....
#Mathematics #Algebra #geometry
タグ: Algebra geometry JuliaLang Mathematics
posted at 01:20:27
akima@4th vaccinated @akima9936
グラフが描けた!(最初のエラーはご愛嬌) #julia言語 #julialang pic.twitter.com/LCU8yRp9mB
posted at 01:37:13
Slides from @Viral_B_Shah's recent talk. It introduces the #julialang, the size of the community, the package ecosystem, #differentiableprogramming, compiler design, and applications of scientific machine learning.
www.slideshare.net/ViralBShah1/ju...
#ML #AI
タグ: AI differentiableprogramming julialang ML
posted at 02:25:08
A Fast Data-Driven Method for Genotype Imputation, Phasing, and Local Ancestry Inference: MendelImpute.jl www.biorxiv.org/content/10.110...
タグ:
posted at 03:14:43
I saw this presented in a workshop at #ASHG20. It would take a lot for me to learn a new language at this point, but the stats here make a really compelling case for looking more at #JuliaLang. 10-100x faster than Beagle, Minimac4 for imputation. Simple API from first glance. pic.twitter.com/ihQqOyPZSi
posted at 03:14:43
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
New post: Bayesian Inference in State Space Models - Part 3 - www.juliabloggers.com/bayesian-infer... #julialang
タグ: julialang
posted at 03:25:27
@joelgrus @jeremyphoward I see. I think it depends on whether you're fully satisfied with the languages you already know.
For me Python was just too slow at some point (like coding own stuff not using c++ packages)
C++ was too complicated at that point :D
#JuliaLang fits me but has other drawbacks
タグ: JuliaLang
posted at 03:28:16
Published a new blog post today on debugging in #JuliaLang
Please check it out and post your comments:
opensourc.es/blog/basics-de...
タグ: JuliaLang
posted at 03:45:51
Bayesian Inference in State Space Models – Part 3 www.juliabloggers.com/bayesian-infer... #juliabloggers
タグ: juliabloggers
posted at 04:17:53
Huge thanks to MIT for taking over #julialang CI on short notice after our AWS Open Source credits ran out. Also thanks to AWS for hosting us this past year, we've made good use of the compute. Maybe we'll be back in the future - still a fan of those Graviton2 machines :)
タグ: julialang
posted at 04:45:11
昨夜はこれを読んだのである。
『横書き句読点の謎』
まったく決まりがないことについてだけ確信が持てるw
ri2t.kyushu-u.ac.jp/~watanabe/RESE...
タグ:
posted at 04:45:40
(And also having faster machines that allow things like PkgEval to be soooooo fast now): twitter.com/vchuravy/statu...
タグ:
posted at 04:46:17
#統計
統計関係のツイートをするようになって、尤度の定義を述べることへの人気の高さに気付くようになった。
私も学生時代に「尤度」(尤度)について講義で聴いたときに何も理解できませんでした。
ある種の人達はずっと伝統的に理解不能な解説を続けているのだと思う。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 05:10:57
Please support meta.stackexchange.com/q/355928/284632 by upvote it to get syntax highlighting for #JuliaLang
And retweet to spread this widely!
タグ: JuliaLang
posted at 05:15:25
@cameron_pfiffer Purchase price for all the hardware currently running various #julialang CI services is probably ~$300k, plus maybe $5k/month in power/bandwidth. My rough rule for cloud is that on prem purchase cost is about 1.5 years cloud cost. Not counting CDN which is sponsored by @fastly.
タグ: julialang
posted at 05:30:47
@cameron_pfiffer @fastly (Also not counting the enormous labor cost of keeping everything running of course :) ).
タグ:
posted at 05:32:07
New post: Train your brain with the Julia language - www.juliabloggers.com/train-your-bra... #julialang
タグ: julialang
posted at 06:09:33
Update: Just got an email from
@AWSOpen that they did not intend to let the #julialang credits lapse, there was just some mixup with our extension request. Since we already migrated some of the CI, my plan is to take the freed up credits to set up #Graviton2 CI.
posted at 07:20:40
かけ算順序検定
①
みかんがたくさんあります。4人に1個ずつ配ることを繰り返すと8回で配り終えました。みかんの個数を求める式はどうなりますか。
答 そんな問題は出ません。
タグ:
posted at 08:39:10
「鎖に縛られたい人だけが従えばよい鎖であれば有ってもよい」
おそらくこれは第一条であって、第二条は、
「自らに課している鎖を他人に強要してはならない」
だと思う。 twitter.com/genkuroki/stat...
タグ:
posted at 08:40:27
1次元Euler方程式のSod衝撃波管問題をHLLで解いています。コードはJuliaで書いています。波速の推定は最も単純なものにしましたがちゃんと解けますね。接触不連続はなまっています。次はHLLCやります。 pic.twitter.com/LiNkegmIdX
タグ:
posted at 09:50:43
1次元Euler方程式のSod衝撃波管問題でAUSMとHLLを比較しています。青がAUSMで赤がHLLです。コードはJuliaで書いています。AUSMのほうがシャープですがHLLも悪くありません。 pic.twitter.com/GU4U3Ldjof
タグ:
posted at 09:54:32
ちょっとしたコードを公開するにはGitHub Gistが便利。個人用電子掲示板。
例えば #Python や #Julia言語 で数学動画を作った場合には是非ともソースコードも公開して欲しいと思う。
gist.github.com/defunkt/370230
A list of Gist clients
posted at 10:15:38
Train your brain with the Julia language www.juliabloggers.com/train-your-bra... #juliabloggers
タグ: juliabloggers
posted at 10:18:22
他のプログラミング言語習得者ならば、基本的な理解は「1から〜」と公式サイトの情報が有れば、言語仕様や文法的に強いクセも少ないので短い期間で十分キャッチアップできると思う。
タグ:
posted at 10:59:12
#統計 KL情報量の-1倍を相対エントロピーと呼ぶ。
対数尤度(尤度の対数)の極限が相対エントロピーになることが、尤度の理解において決定的に重要なことは赤池弘次さんが指摘し、AICはその成果の1つである。
1980年の論説↓
www.jstage.jst.go.jp/article/butsur...
エントロピーとモデルの尤度
赤池弘次 twitter.com/genkuroki/stat...
タグ: 統計
posted at 11:08:27
#統計 訂正
❌ -(log p(Y_1|w) + … + log p(Y_n|w))/n
における "|w" を削除し忘れた。
⭕️ -(log p(Y_1)+ … + log p(Y_n))/n
確率(密度)は対数を取った後の方が扱い易い場合がある。この場合には対数をとってnで割って平均の形にすることによって大数の法則を使える。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 11:14:08
1次元Euler方程式のSod衝撃波管問題をHLLCで解いています。コードはJuliaで書いています。 pic.twitter.com/n2UJefC4iA
タグ:
posted at 11:20:50
非公開
タグ:
posted at xx:xx:xx
MonteCarloProblem って、今は EnsembleProblem に名称変更したのかな? #Julia言語
diffeq.sciml.ai/stable/feature...
タグ: Julia言語
posted at 12:15:15
#Julia言語 「動的型付け」のように見ずに、「多重ディスパッチ」(函数はデフォルトでテンプレート化)と「型情報の伝搬」の組み合わせだと見ればモヤっとせずに済むと思う。
「型情報の伝搬」については zenn.dev/aviatesk/artic... を参照。
@ code_typed optimize=false や @ code_warntype が重要。 twitter.com/skondo/status/...
タグ: Julia言語
posted at 12:45:00
非公開
タグ:
posted at xx:xx:xx
Juliaが型情報をどのように各変数に伝搬させているかは確かに分かりにくいが、Juliaで正しく函数を書けているかのチェックに @ code_warntype マクロの適用は必須に近く、それを繰り返していれば自然にどういう感じかわかって来る。
最近では平易でクリアな解説 zenn.dev/aviatesk/artic... も読める。
タグ:
posted at 12:53:56
「1つのプロジェクト」についてはGitHubの通常のリポジトリを作成した方が便利だが、ちょっとしたコードやメモを他人に見せたい場合にはGitHub Gistが非常に便利。
タグ:
posted at 13:02:41
@MathSorcerer 他のプログラミング言語に乗り換えるときってどうやって勉強すればいいんですかね。ドキュメント見ながら、作りたいものor計算したいものをとりあえず作るのがやっぱり王道でしょうか。
タグ:
posted at 13:04:52
ごまふあざらし(GomahuAzaras @MathSorcerer
@physics303 真面目な話をすると1.0.5は確かに現行のLTSなんですが、近い将来1.6がLTSになります。
有名なJuliaパッケージ(いわゆるライブラリー)は1.x(x>=2)を使わないとインストールできなかったりするのもあります。パフォーマンスも改善された1.5.2を用意しておくとイイと思います。
タグ:
posted at 13:09:15
ごまふあざらし(GomahuAzaras @MathSorcerer
@physics303 とりあえず作りたいものをJuliaで作るでどうでしょうか?それがモチベーションを保つこ
1番の策だとおもいます。
困ったことがあれば
#Julia言語 タグで呟くと答えてくれる人が多いです。
タグ: Julia言語
posted at 13:11:47
④
1 cmは10 mm、すなわち1 mm10個分です。13 cmは何mmですか。式と答を書きなさい。
答
こんな問題は出る(たぶん)。正しい式は10×13のはず。13×10で○をつける教師は超算数の素人か超超算数の使い手かのいずれか。
タグ:
posted at 13:13:12
ごまふあざらし(GomahuAzaras @MathSorcerer
@physics303 Juliaの公式ドキュメント、
パッケージのドキュメント、
JuliaのIssue、
Slack(公式のものや日本コミュニテイのも含む) 、
Discourse
discourse.julialang.org
を漁るだけでも勉強になりますね。
タグ:
posted at 13:17:13
This summarizes my experience with the new logos. twitter.com/killedbygoogle...
タグ:
posted at 13:29:30
@genkuroki コメントありがとうございます。
多重ディスパッチについては他言語での関数オーバーロード的なものを参考に理解出来ました。
抽象型や階層などまだ慣れていない部分もありますが、だいぶJuliaのキモチや設計思想が見えてきました!
新しくプログラミング言語をやはり楽しいですね。
タグ:
posted at 13:41:21
1から始めるJuliaプログラミング。
2.6 多次元配列のところまで読み進む。
後発なだけあって、後付け、別ライブラリ実装ではなく、言語仕様に入っている分、この辺りスッキリしていて良い設計な印象。
タグ:
posted at 13:53:50
1次元Euler方程式のSod衝撃波管問題をLax-Fridrichs法で解いています。やはり数値粘性マシマシでとけますね。接触不連続は見えないです。ギザギザするのは一点おきに参照するスキームだからですかね?どなたか教えて下さい。コードはJuliaで書いています。 pic.twitter.com/EIMC0KPuqV
タグ:
posted at 13:54:13
tsujimotter 日曜数学者 @tsujimotter
ぶっちさんの冒頭のスライドが濃いwww #日曜数学会 pic.twitter.com/4CW8lBmM8R
タグ: 日曜数学会
posted at 13:57:13
ごまふあざらし(GomahuAzaras @MathSorcerer
手元のMacで再現できず。Ubuntuなら動くかな?
docs.opencv.org/master/d8/da4/...
#Julia言語
タグ: Julia言語
posted at 13:58:58
MATLABコードはあまり書ていないので慣れはないけれど、Juliaのcolumn-major orderなのも、多次元配列にメモリ配置並びも、慣れればどうということはない。
タグ:
posted at 14:02:35
ほんと。税金も給食費も教材費も全て返して欲しいです。これは将来字が綺麗になるように指導しているんだそうです。この指導を3年間されたら自己肯定出来なくなり潰されます。勉強が大好きだったうちの子は2年生の後半からノイローゼになりました。 twitter.com/ma85003/status... pic.twitter.com/36GFYlUHTN
タグ:
posted at 14:05:47
#Julia言語 例
module My
struct Pt{T} x::T; y::T end
Base. show(io::IO, p::Pt) = print(io, "(", p.x, ", ", p.y, ")")
end
p = My. Pt{Float64}(1, 2)
↓
(1.0, 2.0)
Base. showの定義がないと
Main. My. Pt{Float64}(1.0, 2.0)
と煩雑に表示される。
gist.github.com/genkuroki/ea59... pic.twitter.com/LtK5lLiKIN
タグ: Julia言語
posted at 14:12:02
ごまふあざらし(GomahuAzaras @MathSorcerer
32ビットシステム”でも“使いたい人なので
嬉しい twitter.com/genkuroki/stat...
タグ:
posted at 14:19:27
2.7 モジュールまで、ざっと目を通した。
他モジュールで定義されている関数を、多重ディスパッチの仕組みを活かして拡張するスタイルがわかって、class継承やテンプレート的な機能に頼らない言語設計の考えがだんだんわかってきた。
多重ディスパッチ周りの挙動をしっかり理解するのが重要ですね。
タグ:
posted at 14:20:28
#Julia言語 実は、Diagonal型の対角行列を
using LinearAlgebra
Diagonal([1,2,3])
のように作ればお得。 twitter.com/genkuroki/stat... pic.twitter.com/6pBZ2AcPUK
タグ: Julia言語
posted at 14:26:00
@genkuroki 型を変数が持っているのが静的型付けで、オブジェクトが持っているのが動的型付けです。どちらも総称であってそれぞれ雑多なサブカテゴリーを含んでいるのですから、「静的/動的型付けだからどうした」の類はたいてい「主語が大きい」といってけっとばして困らないでしょう。
タグ:
posted at 14:39:05
tsujimotter 日曜数学者 @tsujimotter
SageMathでルービックキューブの変換が可視化されるの面白いなー!めっちゃ遊べそう! #日曜数学会 pic.twitter.com/4p9mHZPeNa
タグ: 日曜数学会
posted at 14:45:03
「1から始めるJuliaプログラミング」サラッと読了。あとは実際にコードを書きながら、精読したり公式ドキュメントを読んだりすればいける気がする。
書籍は内容もコンパクトなので、プログラミング初学者向きではないが、Juliaってどんな感じ?というのを理解し言語に触れるきっかけとして良かった。
タグ:
posted at 15:05:05
Julia自体はメタプログラミング機能のコード構文木操作とかもう少し深掘りして触ってみたい。
自動微分パッケージはすでにあるようだけれど、この辺りの事に興味がある人なら、これだけで白米ご飯三杯はいけると思う。
タグ:
posted at 15:07:33
高三 和晃 / Kazuaki Taka @takasan_san_san
素敵なQiita記事を見つけた qiita.com/Baruim/items/6...
タグ:
posted at 15:31:53
高三 和晃 / Kazuaki Taka @takasan_san_san
そもそも被積分関数を出してくるのが面倒すぎてmathematicaから移植するのを諦めていたが、数値積分がうまくいかなすぎるので、Juliaで計算することを検討中
タグ:
posted at 15:54:33
#統計 新刊の大塚淳著『統計学を哲学する』を近所の本屋で買って来ました。まだp.91にしか目を通していないのですが、
【データに基づく信念の改定というベイズ流の考え方】
とか
【ベイズ統計~ベイズ主義では確率は主観的な信念の度合いを測るもの】
と書いてあった!これはひどいと思いました。 twitter.com/genkuroki/stat... pic.twitter.com/eMFcO0Q87n
タグ: 統計
posted at 17:03:27
MATSUMIYA,Shinji @sanjyuumatsu
黒木先生のおかげで、買おうかと思っていた本を買わずに済みます笑 twitter.com/genkuroki/stat...
タグ:
posted at 17:09:54
#統計 さらに、大学でのベイズ統計のある講義では添付画像のように、「主観確率」による「ベイズ主義」によるベイズ統計は廃されており、ベイズ統計は「数学で推測が当たる道を作る」方法の1つに過ぎないという扱いになっています。
watanabe-www.math.dis.titech.ac.jp/users/swatanab...
「主義」 を心配するみなさまに
渡辺澄夫 twitter.com/genkuroki/stat...
タグ: 統計
posted at 17:18:04
#julialang でデフォルトではシングルコアしか使わないことに気づく。いつまでも計算しているなあぁって思ったらCPU使用率が20%くらい。
タグ: julialang
posted at 17:23:03
買おうか迷っていたけど、「哲学」とあったので怪しいと思っていたところ、黒木先生の評価が超参考になります🤗 twitter.com/genkuroki/stat...
タグ:
posted at 17:44:02
ごまふあざらし(GomahuAzaras @MathSorcerer
@mitsuomi_miyata Ubuntuで無事動いた
docs.opencv.org/master/d8/da4/...
タグ:
posted at 17:48:59
Togetter(トゥギャッター) @togetter_jp
いまネットで人気のまとめはこちらです。
「西松屋「わざとガラガラの店舗に」「一店一店がほどほどの売上に」コロナ禍にも繁盛した経営戦略が「柔軟」「こんなカラクリだったとは」」togetter.com/li/1615535.
タグ:
posted at 18:10:07
1次元Euler方程式のSod衝撃波管問題をLLF (Local Lax-Friedrichs) で解いています。Rusanov法とも言います。やはりcentered schemeはかなりなまりますね。コードはJuliaで書いています。 pic.twitter.com/ypWmIEiEI4
タグ:
posted at 19:04:00
#Julia言語 いまさらだけど『1から始めるJuliaプログラミング』はめっちゃいい本。 pic.twitter.com/Wb7Vx8xR3L
タグ: Julia言語
posted at 19:45:52
ごまふあざらし(GomahuAzaras @MathSorcerer
> Windows might work but is not officially tested and supported right now.
BinaryBuilder.jl で行けるか試してましたが
Windows環境のCMake が -WITH_JULIA=ON フラグを認識しないですね. twitter.com/Dsuke_KATO/sta...
タグ:
posted at 20:15:49
須山敦志 Suyama Atsushi @sammy_suyama
この認識が間違いであることは実務者全員が認識した方がよいですね。
>入力と出力の関係性が未知でもデータがあれば予測モデルを立てられる twitter.com/mass_skt/statu...
タグ:
posted at 20:32:02
@genkuroki お!お買い上げどうもありがとうございます。しかしp. 91からとはまた変わったところから読まれ始めましたね。もうちょっと前後の文脈を踏まえてコメントされたほうが良いかもしれませんね。
タグ:
posted at 20:52:51
#統計 おお!著者まで伝わった!
古典統計や頻度主義といった話題との関係とは無関係に、この21世紀に「主観確率」の「ベイズ主義」でベイズ統計について語るのはまずいです。 twitter.com/junotk_jp/stat...
タグ: 統計
posted at 22:32:18
#統計 この本での「確率モデル」という用語の使い方はp.21の脚注1で説明されています(添付画像)。
データの背後に確率法則を想定することと、純粋に数学的に確率空間を考えることは概念的に異なる。
哲学の話をしたいなら、そういう概念的な違いについて繊細な態度を取って欲しかったです。続く pic.twitter.com/AIr0sfjapK
タグ: 統計
posted at 23:41:32
JuliaでGMP使うライブラリ、Reason: Incompatible library version: hogehoge requires version 15.0.0 or later, but libgmp.10.dylib provides version 14.0.0みたいになって落ちる。GMP_jllとJuliaのGMPのバージョンがずれているのだろうか。WindowsとLinuxでは問題がなかった。
タグ:
posted at 23:44:50