黒木玄 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 を見て下さい。
2021年03月21日(日)
私が大学1年で受けたコンピュータの教育はFORTRAN 77で、最初の講義でホーナー法の解説があった。
FORTRAN 77でも教えられることばかりだとつまらないので、数値計算におけるコードの自動生成の重要性なんかも教育に組み入れると面白いと思う。
#Julia言語 なら完全なマクロがある。 twitter.com/Atsushi_twi/st...
タグ: Julia言語
posted at 00:20:28
@genkuroki 例のランダムウォークに関して突っ込まれたとき批判に真摯に対応しなかったことを理由にT氏をhankagosaさんがブロックしたこと,数か月前からやたら主観ベイズの某氏とT氏がべったりなことが効いてそうですねえ.W先生のお弟子さんもT氏にブロックされているようですし断絶が
twitter.com/-/status/13542...
タグ:
posted at 00:27:16
@genkuroki 無用な通知がhankagosaさんやnhayashiさんに飛ばないようなリンクになっています。
スマホアプリなら大丈夫ですが、Webブラウザを使ってる場合は-を_に置き換えることで見られます
タグ:
posted at 00:32:06
#数楽 このスレッドを読んだ人の中には、自分でも sin(x) を近似する多項式を求めたい人がいると思う。
Float64に合わせて係数を作る方法の(試行錯誤そのもののナマの整理されていない)コードが、「低次の係数からFloat64での値に順次固定して求めてみる」の節にあります。
gist.github.com/genkuroki/362b...
タグ: 数楽
posted at 00:36:46
#Julia言語 Juliaでも使われているFloat64でのsin(x), |x|≤π/4をよく近似する13次多項式と同等の13次多項式をJuliaを使って自力で求めてみる話が以下のリンク先スレッドにあります。
Float64に特化した係数の計算は繊細な注意が必要な問題だった。
興味がある人は自分でやってみるとよいと思う。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 00:42:31
#Julia言語 私が求めた13次多項式の係数と伝統的に使われている13次多項式の係数は微小に違うのですが、相対誤差をほぼ同じです。グラフ的には区別ができない!
謎のレシピを謎のまま残さずに自力で求められるようにしておくことが、数学がらみの事柄では基本になります。
gist.github.com/genkuroki/362b... pic.twitter.com/oowxFdjsSV
タグ: Julia言語
posted at 00:49:10
訂正
訂正前:こういう活きた問題が時間をかけて解けば多くのことが身に付くはずです。
訂正後:こういう活きた問題が結構沢山あるので、適当に選んで時間をかけて解けば多くのことが身に付くはずです。
疑問1が解決すると派生の問題が必ず生じる。
タグ:
posted at 01:03:01
DataFrames.jl github.com/JuliaData/Data... などはまだ使っていないけど、Python + Pandasから移行出来るかも‥ twitter.com/chor0sK/status... twitter.com/kamayan1128/st...
タグ:
posted at 01:25:10
ほりたみゅ (mastodonも見てね) @Hyrodium
「FastGaussQuadrature.jlで数値積分しましょう」を公開しました!!死ぬほど分かりやすくGauss求積の解説を書いたので死なないように気をつけてください…
zenn.dev/hyrodium/artic...
タグ:
posted at 02:10:54
ほりたみゅ (mastodonも見てね) @Hyrodium
zennの記事は4本目で、文字数が単調増加してます pic.twitter.com/5Jt8SuVQcB
タグ:
posted at 02:12:51
そして…
「千葉県のシンボル/は「カタツムリ」です。八木秀次氏が〜『正論』に詳しく書いてますが、これはカタツムリが雌雄同体であることを取り上げたもの〜」
ケント・ギルバート『やっと自虐史観のアホらしさに気づいた日本人』2016
twitter.com/tr_727/status/... pic.twitter.com/tHexP5ZoWZ
タグ:
posted at 02:34:53
【19桁】第1種カニンガム鎖、1.255e18までの処理が完了。Githubのリストを更新。raw.githubusercontent.com/kazkus/Cunning... CC13は55番目の 1231790682838960529 までが正式。CC14は3番目の 1092966255671664659 までが正式。 pic.twitter.com/IdJu5MhaYY
タグ:
posted at 02:35:06
#Julia言語 10行シリーズにマンデルブロ集合を追加した。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/JdwKqbNPZA
タグ: Julia言語
posted at 02:37:21
#Julia言語 こういう動画を10行で作れます。
何百行もあるコードだと解読する気になれなくても、たったの10行ならば解読する気になれるかも。実際に解読できれば、こういうことを自分で自由にできるようになる。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/GIKHD73XsA
タグ: Julia言語
posted at 02:37:24
New post: Expand your DataFrames.jl toolbox: the flatten function - www.juliabloggers.com/expand-your-da... #julialang
タグ: julialang
posted at 02:45:18
Mark Kittisopikul ht @markkitti
OpticSim.jl: Optical Simulation software in #JuliaLang from @microsoft
github.com/microsoft/Opti...
タグ: JuliaLang
posted at 03:31:11
非公開
タグ:
posted at xx:xx:xx
@genkuroki @tsatie これ、世界を普通のトーラスでなくねじれトーラスにしたら世界がグライダーまみれになって楽しそうです。単にアスペクト比の変更でも可。
タグ:
posted at 06:07:40
#Julia言語 作り直し2 (こちらが本命、拡大率を上げた)
これももちろん10行。 pic.twitter.com/w47s7ViI5i
タグ: Julia言語
posted at 06:41:48
#Julia言語 10行でライフゲーム
初期条件は5×5
2000ステップのアニメーション pic.twitter.com/b59pxDHfwG
タグ: Julia言語
posted at 07:19:09
#Julia言語
10行シリーズ全ソースコード↓
nbviewer.jupyter.org/gist/genkuroki...
添付画像は10行ライフゲームのコード
初期条件やプログレスメーターやJuliaの文法で必要なendも含めて10行!
可読性も悪くないはず。 pic.twitter.com/5yjacEYPYe
タグ: Julia言語
posted at 07:24:05
#Julia言語 配列uに保持した状態を更新するコードを
for t in T
u = update(u)
end
の形式で行うと無駄なメモリアロケーションの嵐になる。
for t in T
update!(u)
end
のようにuの内容を更新するように書かないと損をする。
以下のライフゲームのコードもそうなっています。 pic.twitter.com/TjO2IbIeNy
タグ: Julia言語
posted at 07:27:07
#Julia言語 配列に保持した状態を時間発展させるコードはありがちなので、ライフゲームのような遊びのプログラム作成もバカにできない。10行の中に基本が詰まっている。
10行の中にプログレスメーターを入れる余裕があるとちょっとうれしい。 pic.twitter.com/zOsM2JI69Q
タグ: Julia言語
posted at 07:30:26
@VC31415 I’d add #julia for programming. It’s faster/more flexible/closer to the math. github.com/azev77/QuantEc...
タグ: julia
posted at 08:21:29
@MapleTaffyCap @VC31415 Julia job postings are growing. Besides, Julia teaches better programming habits, so many employers who list “Python” will take Julia experience seriously. Stanford engineering classes use Julia from beginning, their grads have no problems finding jobs...
タグ:
posted at 08:52:22
朝からガサガサしてるなぁと思ったら。Juliaさん脱皮っピ。なぜうちの子はイベント行く時とか脱皮するん?ちゃんとキレイに脱げました。
#レオパ #レオパードゲッコー #ヒョウモントカゲモドキ pic.twitter.com/SrlChwYHze
タグ: ヒョウモントカゲモドキ レオパ レオパードゲッコー
posted at 09:00:41
@MapleTaffyCap @VC31415 Agreed, but the fact that Stanford/MIT are increasingly using Julia in their curriculum gives an idea about trajectory. Good luck
タグ:
posted at 09:00:51
マイTips:複数言語でコードを書いてそれぞれからfigureとかtableを出すとき、wd直下データフォルダは共通にして、main_Rとmain_juliaとmain_pythonにメインのコードを入れて、julia_figuretableとR_figuretableのように言語ごとにフォルダ分けて出力を保存したほうがLatexに貼るとき混乱しない。
タグ:
posted at 09:20:39
@VC31415 @MapleTaffyCap Paul has some awesome #Julia packages, most of his recent stuff on Git is in jl: github.com/schrimpf?tab=r...
タグ: Julia
posted at 10:25:16
@VC31415 @MapleTaffyCap Mert posted R code to replicate your DML paper. I tried it on a big dataset & took forever. I have a strong prior it’d be much faster in Julia. I can try to translate when I have time. github.com/demirermert/ML...
タグ:
posted at 10:31:43
@grant_mcdermott Agree
When picking two of R/Python/Julia, also take into account that Julia is a compiled language for performance
All else being equal, this key point makes Julia a great language to compliment R or Python
- e.g. as compared to R + Python + C when you need performance
タグ:
posted at 13:04:58
ほりたみゅ (mastodonも見てね) @Hyrodium
4点を選んでLangrange補間して3次多項式を作って数値積分するGIFです
www.desmos.com/calculator/c9r... pic.twitter.com/stoFKjDnqg
タグ:
posted at 13:46:32
非公開
タグ:
posted at xx:xx:xx
津田和俊/急激に進行した網膜剥離と闘って @kaztsuda
400年前にローマ法王に謁見するため太平洋を二度往復した船を復元したレプリカ。 場所: 宮城県慶長使節船ミュージアム・サンファン館 www.instagram.com/p/CMqyWb2j0SV/...
タグ:
posted at 14:10:20
@m_ero_n seedを102に選んで自分のストーリーにとって都合の良いランダムウォークのサンプルを生成したという「不正の証拠」を自分自身で公開してしまっているのに、ああいう態度を取るのはまったく理解できませんでした。seed 101だと「ストーリー的に都合の悪いプロット」が得られることも容易に確認できる。 twitter.com/genkuroki/stat... pic.twitter.com/W68poKR1Hq
タグ:
posted at 14:37:27
@genkuroki 個人のブログなので~と庇う人が職業研究者にすらいましたが、個人のブログだからこそ謝って済むだけの話なのに変に徹底抗戦して拗らせてしまった。
ただ、ある意味かの御仁の馬脚が現れた事件でもありましたのでTJOさんの信者以外にとっては有益だったのかもしれません
タグ:
posted at 14:43:25
@m_ero_n 松浦健太郎著『StanとRでベイズ統計モデリング (Wonderful R) 』
www.amazon.co.jp/dp/4320112423
は細部に気を使って書かれている非常に良い本です。
ベイズ統計と無関係にパラパラめくってプロットの仕方を眺めるだけでも参考になる。
こういう本が売れるとよいと思います。#統計
タグ: 統計
posted at 14:46:49
松浦健太郎著『StanとRでベイズ統計モデリング (Wonderful R) 』
amazon.co.jp/dp/4320112423
が売れて、著者に興味を持った人がブログ
statmodeling.hatenablog.com/entry/bayesian...
を発見する。こういう流れができれば、みんなハッピーになれる。
この本の考え方とセンスは普遍的でStanと無関係に価値があります。
タグ:
posted at 14:58:29
#Julia言語 2次元Isingモデルのアニメーションを10行で!
10行シリーズ全ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/UVFEihNjDU
タグ: Julia言語
posted at 17:47:50
#Julia言語 10行で作った2次元Ising模型の動画その1
200×200の盤面全体を200回更新している。 pic.twitter.com/1uPE4pIznX
タグ: Julia言語
posted at 17:50:22
#Julia言語 10行で作った2次元Ising模型の動画その2
200×200の盤面全体を200回更新するごとにプロットすることを200回繰り返している。 pic.twitter.com/BvdCT5u0c4
タグ: Julia言語
posted at 17:51:30
#Julia言語 10行シリーズ
1. 信頼区間が二項分布モデルのパラメータ値を含む確率
2. 同ケースでの信頼区間の分布
3. 以上のFisher検定版
4. 3次元ランダムウォーク
5. Ising2D.jl
6. 蔵本モデル
7. ジュリア集合
8. マンデルブロ集合
9. ライフゲーム
10. 2D Ising model
↓
nbviewer.jupyter.org/gist/genkuroki...
タグ: Julia言語
posted at 17:58:11
ほりたみゅ (mastodonも見てね) @Hyrodium
@NotionHQ さん、Juliaに対応してくれると嬉しいです pic.twitter.com/z7Q2sXNmwC
タグ:
posted at 18:13:15
ごまふあざらし(GomahuAzaras @MathSorcerer
AWS の CloudShell で wget/tar で Julia インストール・展開できたでキュ.Python入ってるんで PyPlot.plot で関数をプロットして描画結果を Actions -> Donwload file で落とすことができた.
Plots.jl はアーティファクトが重いので1GB以内に多分入らない. pic.twitter.com/IXAvyUkt69
タグ:
posted at 18:54:09
@h_okumura 生協の人が保証が手厚いことを売りにしてたけど、
それに対して半額の同じスペックのパソコンを買って、
壊れたらもう半額でパソコンをもう一台買う方がお得だってぼやいてた人がいて確かにその通りだなあと思った。
もう一台買い直すときにはスペックも格段に向上しているはずだし。
タグ:
posted at 18:59:10
ごまふあざらし(GomahuAzaras @MathSorcerer
2Core, メモリ4GB
これぐらいのとこなら遊べそう.
julia> versioninfo(verbose=true)
でもっと色々わかる.
#Julia言語 pic.twitter.com/WQmPCoLtJp
タグ: Julia言語
posted at 19:11:12
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
塾も学校も課程を消して整然とした式と答えだけを残させる。真面目に有害。
矯正のために消しゴム禁止で解かせることもある。スペースの使い方やちょったした計算や図って凄い大事なのにね。 twitter.com/mimizuku_tutor...
タグ:
posted at 19:58:48
非公開
タグ:
posted at xx:xx:xx
The interactive and reactive Julia notebook system Pluto is having its own mini-conference. plutojl.org/plutocon2021
Can one of the tech giants please pour money into the development of Pluto? It is such a great tool for teaching and research communication.
#julialang pic.twitter.com/K8wx14KBiM
タグ: julialang
posted at 21:01:17
NotionをChromeでweb越しに編集すれば,GhostTextと組み合わせて好きなエディタで編集できる.
chrome.google.com/webstore/detai... pic.twitter.com/ZbOYFj2H4y
タグ:
posted at 21:04:31
That time a team of researchers recreated the famous Draupner freak wave measured in the North Sea in 1995, and found it bore an uncanny resemblance to "The Great Wave off Kanagawa' by the Japanese artist Katsushika Hokusai [read more: ow.ly/9iCt30np5oT] pic.twitter.com/17SvgqL90B
タグ:
posted at 22:32:34
📈 ピックアップ ✨
FastGaussQuadrature.jlで数値積分しましょう
by @hyrodium
zenn.dev/hyrodium/artic...
タグ:
posted at 22:38:00
#Julia言語 10行では無理なやつ
正方形上の自由境界条件の波動方程式の離散化をDifferentialEquations.jlで数値的に解いて動画を作ってみた↓
nbviewer.jupyter.org/gist/genkuroki...
これを作るために必要な行数は40行程度でした。 twitter.com/genkuroki/stat... pic.twitter.com/0QZAtKVwpd
タグ: Julia言語
posted at 22:49:30
#Julia言語 約1万個の質点がバネで繋がっている状況を記述する常微分方程式 prob を DifferentialEquations.jl に sol = solve(prob) の形式で数値的に解かせています。自分でソルバを書かなくて良いとかなり楽です。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/v1vvtPliyS
タグ: Julia言語
posted at 22:53:29
#超算数 #超国語 #十分 #星はキラキラorきらきら
公文書の読点「,」から「、」に 半世紀以上前の通知変更へ(共同通信)
#Yahooニュース
news.yahoo.co.jp/articles/88039...
タグ: Yahooニュース 十分 星はキラキラorきらきら 超国語 超算数
posted at 23:58:55
非公開
タグ:
posted at xx:xx:xx