黒木玄 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月12日(金)
@ham144430 長期ではなく短期の職業訓練としては Python でもいいとは思います。Pythonも非常に便利です。
しかし、C, C++, Fortran 並の計算速度が欲しい場合もある分野の勉強をしたければ、PythonよりもJuliaの方が明らかに勝ります。
タグ：
posted at 23:51:36
@ham144430 @h_okumura 私の #Julia言語 版のコードも函数化しているだけ(プログラミングを教えるなら函数化した形で教えた方が良いと思う)です。
添付画像のJulia版とC版を比較して見て下さい。
Julia版の方が圧倒的にシンプル！
科学技術計算や最近の新しい計算技術を学びたいなら #Julia言語 がよいと思います。 pic.twitter.com/0OJLh466g3
タグ： Julia言語
posted at 23:46:47
「日本語LaTeXをバリバリ使っている💪、でもTeX言語はアンマリ／サッパリ／ゼンゼン😓」
という人がpLaTeXカーネル／LuaTeX-jaパッケージ開発への「貢献」を何かしたいと考えているなら、まずは「開発版のテスト」への参加を検討してほしい。
#TeX #TeXLaTeX twitter.com/wtsnjp/status/...
posted at 23:39:46
#Julia言語 同名が存在する野良パッケージ
github.com/genkuroki/MyUt...
MyUtils.jl
を使った「e のモンテカルロ計算の並列化・分散化の例」
gist.github.com/genkuroki/fa10...
この例を知っているだけで、並列化・分散化を気軽に楽しめる人が増えると思う。 pic.twitter.com/tElToWFgXu
タグ： Julia言語
posted at 23:19:35
この半年くらい「色覚」のワードでずっとウォッチしているのだけれど、色覚検査を受けて、医師、警官、自衛官、消防士etc...になれないと言われたする当事者（or親）が多くてびっくりしている。いずれも無条件になれないものではないのに、心痛む。
タグ：
posted at 23:09:31
#Julia言語 例の e のモンテカルロ計算の @ threads を使ったコードが不適切だった。添付画像は訂正後の結果。
もとの話題ではJuliaで28秒かかったと書いていたが、普通に書けば私の環境だと1.4〜1.5秒程度でその20倍くらい速いし、並列化・分散化すると、0.27～0.29秒！
gist.github.com/genkuroki/5212... pic.twitter.com/DeATYzRVS8
タグ： Julia言語
posted at 22:43:45
@genkuroki @h_okumura juliaが速いとは聞いていましたが、ここまで驚愕の速度が得られるとは知りませんでした。
ここで相談なんですが、プログラミング初心者(CとFortranを少しかじった程度の機械工学科の学生)がこれから覚えるとしたらpythonとjuliaどちらがお薦めですか？
タグ：
posted at 22:43:25
@genkuroki @h_okumura アドバイスありがとうございます。
コードの後半部分を関数化せず、オリジナルのforループをそのまま使っています。プログラミング初心者の息子に見せるために、読みやすさを優先し最適化や拡張ライブラリは使っていません(詳しくないので使えてませんｗ)。
タグ：
posted at 22:29:11
なんやいまどきの大学生はこんなに読ませるんか。しかも今期のB1とは、ええもん見たなあ
眼鏡フェチの二次元オタクが将棋にドハマりした話｜加藤葵 #note #スキしてみて note.com/saneomi/n/n610...
posted at 22:01:09
NATROMさんの過剰診断の解説は分かりやすいので、一読を勧めます |
（主に）過剰診断の定義に関する記録 - NATROMのブログ natrom.hatenablog.com/entry/0021/02/...
タグ：
posted at 21:59:15
NASA Uses Julia to Model Spacecraft Separation Dynamics: NASA’s Jonathan Diegelman presented Modeling Spacecraft Separation Dynamics in Julia at the SIAM 2021 Conference on Computational and Scientific Engineering - youtu.be/tQpqsmwlfY0
#julialang #aerospace #simulations
タグ： aerospace julialang simulations
posted at 21:58:45
#Julia言語 DifferentialEquations.jl
この動画が弾性衝突を扱っているとは普通思わないよな。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/XghEhHFsgC
タグ： Julia言語
posted at 20:40:26
UNSCEARのレポートに過剰診断の可能性があるとまで書かれて、それでも甲状腺検査を続けますかね、福島県は。
もう充分なんじゃないでしょうか。充分すぎると思いますよ
タグ：
posted at 20:24:06
開発版 pLaTeX で問題を生じるケースを発見された場合の報告先はこちらの GitHub リポジトリです．最悪私宛てのリプライでも構いません．よろしくお願いいたします．
github.com/texjporg/plate...
タグ：
posted at 20:07:28
日本語 LaTeX が窮地なのは冗談ではなく本当です．日本語 LaTeX メンテナンスの難しい点の1つは「どんな不具合が生じるか予見できない」点にあります．ご協力いただける方は，ぜひ開発版 pLaTeX のテストにご参加ください．その方法は下記で解説しています．
#pTeX #TeXLaTeX
qiita.com/wtsnjp/items/7...
posted at 20:06:21
式と答えばかりで数学の問題を解くわけではないのだから、
「説明 結果 」みたいに自由の利く方が良い気がします。 twitter.com/genkuroki/stat...
タグ：
posted at 19:44:56
@nyanko_box ```
for row in eachrow(A)
println(row)
end
```
が一番スマートかな、と。
インデックス取りたい場合は `for r in axes(A, 1) ～ end` の方が少しだけスマート。
タグ：
posted at 19:39:14
@OokuboTact #超算数 座談会号、猫が運送してきました！
ここに出てくる人は、矢口新さんをのぞきすべて1951年指導要領関係者です。 pic.twitter.com/VrFUEYIQHN
タグ： 超算数
posted at 19:37:47
いいなぁ
モフモフなうさぎさんが最高にモッフモフ！ 最大限に味わう飼い主さんがうらやましい nlab.itmedia.co.jp/nl/articles/21... pic.twitter.com/8r5FApObEg
タグ：
posted at 19:05:00
#Julia言語 これらのマクロ達(@ my_threads, @ my_distributed)はものすごく頻繁に使いそうになるように思えたので、野良パッケージ化しました。
github.com/genkuroki/MyUt...
この野良パッケージは今度私の落書き帳のようになる可能性があります。
タグ： Julia言語
posted at 19:02:13
#Julia言語 @ my_threads マクロの仕様を変更した。
複数のスレッドの分けて実行するforループの前だけではなく、後に実行するコードも指定できるようにしておかないと困るという当たり前のことに気付いた。
gist.github.com/genkuroki/1efd...
添付画像の左右のコードは実質的に等価だと考えてよい。 pic.twitter.com/LN8sy0CKXY
タグ： Julia言語
posted at 18:57:49
Carlos Balsalobre @balsalobrephd
How to explain Normal Distribution to a bro in the gym pic.twitter.com/WZN9JHUm9b
タグ：
posted at 18:47:20
Ligun (りぐん) ☠️ Nenio @choonggeun
@antimon2 ありがとうございます！
まだJulia知らないこと多いのでとても助かります。
タグ：
posted at 18:22:39
ごまふあざらし(GomahuAzaras @MathSorcerer
はいはいはいはい！
#Julia言語
github.com/JuliaLang/juli...
タグ： Julia言語
posted at 18:21:29
Pythonだとforにnumpy.arrayを渡すと1行ずつ処理できるけど、#Julia言語 で同じことやるときはどうするんじゃろ。 pic.twitter.com/VwKg9jLiv3
タグ： Julia言語
posted at 17:12:32
@tsatie @sekibunnteisuu 試行錯誤の結果捨て去られるかもしれないが、集合と写像n〇言葉だけで作れる定式化(抽象化)を行なってみて、白黒を付け易い方向に持って行くとよいと思う。
その定式化内で「やりたいこと」(←これも自分で決める)ができるかどうかが最初のチェックポイント。
すでに私が色々やって見せています。
タグ：
posted at 16:31:53
@tsatie @sekibunnteisuu 【月曜日だと云われた時に表と裏の確率は半々なのだから1/2,1/2と考えるべきではと思う】
そのように思うことを、集合と写像の言葉だけで作れる数学的構造で定式化してみるとよいと思います。そういう試行錯誤が必要。
曖昧なことを曖昧なまま考えることをどんなに続けても堂々巡りになるだけ。
タグ：
posted at 16:28:10
@OokuboTact #超算数
等分除・包含除でこの誘導尋問を思い出した。子供は最初から答えが分かっている。手間暇かけて一体何を教えているのやら・・・ pic.twitter.com/mEjtA0AkXj
タグ： 超算数
posted at 16:27:47
@OokuboTact #超算数
これ、子供の感想を肯定的にとらえて、自分の授業を自画自賛しているんですかね？
馬鹿ですね。 pic.twitter.com/96MqwE9DI7
タグ： 超算数
posted at 16:17:03
@choonggeun `abs(z)<Inf` は `isfinite(z)` で良いと思います。複素数にも対応してるはずだし多少なりともより効率良くなるはず。
タグ：
posted at 15:47:36
#Julia言語
@ my_threads - threads version of @ my_distributed
gist.github.com/genkuroki/1efd... pic.twitter.com/n4pGlAsEX2
タグ： Julia言語
posted at 15:13:51
留学生寮での友達(今はUKにいる)が久々に連絡くれてよもやま話をしてたんだけど、彼がJET(Japan Teaching and Edu) programmeで頑張ろうとしてるという話になり、ちょうどこっちもJETで羽ばたこうとしてるから、ただのダジャレだけど、「やっぱ俺ら繋がってんな！」って謎に熱い感じになったという話✈️
タグ：
posted at 14:31:39
@sekibunnteisuu 「具体的な状況や場面を忠実に翻訳した数学的構造が一意に決まる」という考え方はチョー算数の中心的教義の発展版になっています(笑)。
「我々」の立場的には当然否定しておきたいところ。
タグ：
posted at 14:05:50
@sekibunnteisuu P(表月)=P(裏月)=P(裏火)=1/3のスタイル(もしくはその一般化)でも適切な数学的構造を補ってやれば、π:{表月,裏月,裏火}→{表,裏}経由と同等のことを常にできます。
数学的に頑張りたい人は適切な「同型」の概念を定義すればよいでしょう。
タグ：
posted at 14:03:18
@sekibunnteisuu P(表月)=P(裏月)=P(裏火)=1/3
をコインで表が出る確率が p の場合に一般化すると、
P(表月) = p/(p+2(1-p))
P(裏月) = P(裏火) = (1-p)/(p+2(1-p))
となります。
この公式を使えばこのスタイルでの確率の値からpを逆算できる。
ギャンブラー的には裏月と裏火が一塊という情報も忘れてはいけない。
タグ：
posted at 13:59:14
@sekibunnteisuu トゥーンさんには誤解されているような気がします。
私は「筋の通った考え方は複数通りある」とはっきり述べているし、表月, 裏月, 裏火のそれぞれの確率の値を 1/3, 1/3, 1/3 とする考え方をしてもよいと述べています。
ただし、それだとギャンブラー的には情報が落ちているので補う必要がある。 twitter.com/movmovology/st...
タグ：
posted at 13:42:06
「福島県民がん増える可能性低い」 被曝線量を下方修正：朝日新聞デジタル www.asahi.com/articles/ASP39...
【科学委は報告書で、被曝の影響ではなく、高感度の超音波検査によって「生涯発症しないがんを見つけた過剰診断の可能性がある」と指摘した。県の評価部会の専門家も同様の指摘をしている】
タグ：
posted at 13:19:30
Ligun (りぐん) ☠️ Nenio @choonggeun
昨日の #Julia言語 でのマンデルブロ集合のプログラム、少し修正してみた。
関数で分けてシンプルにかけたかなと思う。数学的な定義を素直に書けてるかなと
ちょっとはJuliaっぽい感じになったかな？ pic.twitter.com/pbKM5YSeN4
タグ： Julia言語
posted at 12:41:14
OokuboTact 大久保中二病中年 @OokuboTact
＃超算数
算数教育の専門雑誌を読んでいたら、掛け順の話があった pic.twitter.com/NlIbGuYBwH
タグ： 超算数
posted at 12:27:16
#Julia言語 例の e のモンテカルロ計算については
gist.github.com/genkuroki/5212...
を参照。色々試してみたが、私の環境では @ my_distributed 版がシンプルでかつ一番速かったです。 pic.twitter.com/LKvfTB3d7C
タグ： Julia言語
posted at 12:25:31
#Julia言語 で @ distributed を使ってモンテカルロシミュレーションの足し上げを行うときに、forループの内側で rand() を使うと遅くなります。固定された rng について rand(rng) の形式で実行する必要がある。
解決策：@ my_distributed マクロを自作して使う！
gist.github.com/genkuroki/a6dc... twitter.com/genkuroki/stat... pic.twitter.com/QvwSVYuaV7
タグ： Julia言語
posted at 12:21:03
@sekibunnteisuu 個人的には、「複数の筋の通った考え方がある」「パラドックスは存在しない」という感じで納得できているように思うのですが、まだ考えが足りないかも。
タグ：
posted at 10:57:30
@sekibunnteisuu このギャンブルでは、コイントスで裏が出たときの裏月と裏火を分割できないので、期待値を考えるときに、期待値を取られる側の確率変数({表月,裏月,裏火}上の函数)の集合{裏月,裏火}上への制限は定数函数でなければいけない。
これがπ:{表月,裏月,裏火}→{表,裏}経由での考え方。
タグ：
posted at 10:54:52
@sekibunnteisuu このとき、π:{表月,裏月,裏火}→{表,裏}経由で考えて、さらに手で確率を裏月と裏火の単体に拡張した立場での、常に「はい」と答えた場合の期待値は
(1/2)a + (1/4)(b-c) + (1/4)(b-c) = (1/2)a + (1/2)(b-c)
になります。これを(1/2)a + (1/4)b + (1/4)(-c)とするのは間違った考え方になります。
タグ：
posted at 10:45:30
@sekibunnteisuu 表月, 裏月, 裏火のそれぞれで「今日は月曜ですか？」と質問して「はい」と答えたら、正解の表月, 裏月のときにはそれぞれa, bの賞金を支払い、不正解の裏火にはcの罰金を取り立てることにします。続く
タグ：
posted at 10:43:17
@sekibunnteisuu 確率の表(≒確率測度)を例えば
表月 1/2
裏月 1/4
裏火 1/4
に設定した場合につい考える場合には、「曜日に関係なくどの時点であっても表である確率は1/2であるはずである」という「直観」を無断で使っちゃダメで、設定した確率の表と数学的定義に忠実に確率を計算しないと奇妙な結果が導かれます。
タグ：
posted at 10:31:23
@sekibunnteisuu 表月, 裏月, 裏火の確率がそれぞれ1/2, 1/4, 1/4のとき、定義に忠実に計算すると
月曜である確率=P(表月)+P(裏月)=3/4
月曜という条件の下での表の条件付き確率
=月曜でかつ表である確率/月曜である確率
=(1/2)/(3/4)=2/3
月曜である確率×月曜という条件の下での表の条件付き確率
=3/4×2/3=1/2 twitter.com/sekibunnteisuu...
タグ：
posted at 10:21:38
Soraさんの意見にはいつも賛成です。
ただし、「立式」という用語は一般的な用語ではなく、算数教育では「場面や考え方を忠実に式に翻訳すること」という意味で使われており、たし算やかけ算の悪しき順序指導の根拠としても使われる言葉です。意味を代わりに簡単に説明させて下さい。 #超算数 twitter.com/musorami/statu...
タグ： 超算数
posted at 10:03:45
非公開
タグ：
posted at xx:xx:xx
Check out "Natural Language Processing in #JuliaLang" presented by Ayush Kaushal @ JuliaCon 2020
youtu.be/hHCi8ojazqk
タグ： JuliaLang
posted at 08:15:00
@TrevorABranch Oh yes, I like it! #JuliaLang The first to use that in a peer-reviewed paper gets a beer 🍻 pic.twitter.com/YIIBkVezpQ
タグ： JuliaLang
posted at 03:37:52
@StefanKarpinski Thanks for sharing the full story. The biggest harm is the effort spent, it could have been worst if the whole deceiving identities hadn't been discovered
タグ：
posted at 03:17:54
Stefan Karpinski @StefanKarpinski
@matbesancon In short, I think it was not something they originally intended to do. I think it was petty, personal and done because of a bruised very fragile ego in a misguided attempt to “punish” the project. And of course, in the end, comically revealed by their own incompetence
6/n
タグ：
posted at 02:03:20
Stefan Karpinski @StefanKarpinski
@matbesancon They also either recruited their friend to create yet another fake identify to sow FUD in the same thread. The friend may also just have shared login info with the original person.
5/n
タグ：
posted at 02:00:58
Stefan Karpinski @StefanKarpinski
@matbesancon To do so without looking as petty and personal as it actually was, they created fake identities, including a fake “free software activist” to write a really bad hit piece and an innocuous looking curious academic type to bring up a different (real) critique of the project
4/n
タグ：
posted at 01:58:24
Stefan Karpinski @StefanKarpinski
@matbesancon They decided after that discussion that they hate the project, or at least the moderators and leadership (ie me), and wanted to bring negative attention to it as punishment.
3/n
タグ：
posted at 01:56:01
Stefan Karpinski @StefanKarpinski
@matbesancon b) that the thread was slowed down to prevent overheating despite that action explicitly not being directed at them.
They then recruited a real friend to help drum up outrage about the (entirely reasonable) moderation action
2/n
タグ：
posted at 01:53:15
Stefan Karpinski @StefanKarpinski
@matbesancon My impression is that this was a person who was initially quite pro-Julia, gave some mixed feedback, due to a fragile ego and was offended that a) it wasn’t all simply accepted without any discussion and
1/n
タグ：
posted at 01:51:43
効果音とかBGMで全然印象変わるから報道ではやめてほしいよね。
nico.ms/sm663696 twitter.com/genkuroki/stat...
タグ：
posted at 01:43:54
Ligun (りぐん) ☠️ Nenio @choonggeun
@genkuroki なるほど、Juliaだと一般的なテクニックなのですね
色々なサンプルありがとうございます
まだ一月くらいしか触ってなくJulia力が足りないため追いきれない場所も多いですが、ゆっくり確認してみます！
タグ：
posted at 01:39:17
@LirimyDh #Julia言語 私もその辺をよく理解していないのですが、なぜか Random.default_rng() がかなり遅いので、それを呼ぶ回数を減らすとさらに速くなると思います。
gist.github.com/genkuroki/5212...
gist.github.com/genkuroki/771f...
タグ： Julia言語
posted at 01:16:08
村）Ａ級昇級を決めた永瀬王座にインタビューしました。「６連勝の時は安全かなと思ったが、それから崩れてしまった。後半かなり苦しんだ」「鈴木先生（＝鈴木大介九段）には２０代のうちに上がるように、と言われていた」「（名人）挑戦にからめるようにしないといけない」 pic.twitter.com/MVVG8ybFHx
タグ：
posted at 01:12:38
僕はいろいろなことをほったらかしてきてしまったし、事実上なにも後世に残すものは作れていないのだけれども、それでも「いちから聞きたい放射線のほんとう」はいい仕事をしたと自分でも思っています。
小峰さんとおかざきさんのおかげです。読んでくれたみなさんにもほんとうにありがとう
タグ：
posted at 01:02:48
耕）Ｂ級１組順位戦、最後に残った対局は、午前０時37分、木村一基九段が深浦康市九段に131手で勝ちました。阿久津主税八段が勝っているため、深浦九段は降級が決まっています。
タグ：
posted at 00:40:13
Juliaの文字列検索で2-10倍くらいの性能向上ができそうだ。github.com/bicycle1885/St... pic.twitter.com/zuDO2N9g8H
タグ：
posted at 00:38:53
村）永瀬拓矢王座と近藤誠也七段の感想戦が続いています。近藤七段の玉は入玉した後、自陣に戻って粘りましたが、永瀬王座の正確な指し手の前に屈しました。永瀬王座は１３、１４日に島根県で王将戦第６局を戦います。 pic.twitter.com/Jm3TeNYvOT
タグ：
posted at 00:34:40
@shg #Julia言語 絶対誤差 < 1e-16 の部分を相対誤差にすれば大丈夫です↓
gist.github.com/genkuroki/e69a...
タグ： Julia言語
posted at 00:05:06
