黒木玄 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年02月23日(火)
非公開
タグ：
posted at xx:xx:xx
#数楽 #Julia言語
Remez法が「神」でした。
github.com/simonbyrne/Rem...
Remez.jl
のコード中の threshold の決め方は精度や問題ごとに自分で調節することが必要な場合もあるので(永久ループに入る場合がある)、dev をして適当に書き直してから使った方が良いかもしれません。
posted at 00:27:30
juliaでmylibのディレクトリをLOAD_PATHにpush!しただけだと動くけど、
manifest更新しろってエラーがでる。
別にUUID付与したpkgを作るほどのものじゃないから指定ディレクトリにある.jlを参照してくれるだけでいいんだけど何でだろ。
タグ：
posted at 00:31:35
#数楽 Remez.jlのメインの函数の ratfn_minimax の雑な修正の例が
gist.github.com/genkuroki/362b...
に my_ratfn_minimax にあります。
タグ： 数楽
posted at 00:33:35
#数楽 このスレッドを読んだ人の中には、自分でも sin(x) を近似する多項式を求めたい人がいると思う。
Float64に合わせて係数を作る方法の(試行錯誤そのもののナマの整理されていない)コードが、「低次の係数からFloat64での値に順次固定して求めてみる」の節にあります。
gist.github.com/genkuroki/362b...
タグ： 数楽
posted at 00:37:33
And of course, the more popular it gets, the more time it will require from you. (This is the trap!)
That's what I like to call: "The maintenance dilemma".
At some point in time, you’ll have to take a decision:
1 -You stop maintaining it.
2- Or you keep doing it.
タグ：
posted at 01:13:08
In fact, sometimes the more they want to help the more of a burden they put on your shoulders.
Now, on top of all you do, you'll have to review their pull request, understand it and potentially get into a conversation to fix that issue/feature you've never thought of.
タグ：
posted at 01:13:10
[23日01時34分頃] 岩手県 宮城県 山形県 福島県 茨城県 あたりで揺れたかも… #地震 pic.twitter.com/m0pM1UHwuF
タグ： 地震
posted at 01:34:33
Python があまりにも嫌い過ぎてサンプルコードが Python で書かれてる本に手が出なかったんだけど，Julia 書くようになってからその手の本を買うことに抵抗がなくなった，読むのは不快だけど手元では Julia の REPL で動かすので
タグ：
posted at 05:38:53
命令数4って何??って思って調べたら実際はもっと凄かった。基本的に命令は「減算して負だったら分岐」の1つのみ（なので"Sub"RISC）。それでチューリング完全だと。アセンブリ見ても人類には理解不能。 twitter.com/pc_watch/statu...
タグ：
posted at 07:39:32
作品として成立し得ないと思うほどの「ありふれた日常」が福島に来る日をボクはずっと待ち望んでいるよ twitter.com/miura_hideyuki...
タグ：
posted at 07:49:55
非公開
タグ：
posted at xx:xx:xx
あの赤い彗星も囲碁やってるんですよ！？( ͡° ͜ʖ ͡°)
#囲碁の魅力を言語化 pic.twitter.com/HKFH7XfGbf
タグ： 囲碁の魅力を言語化
posted at 09:37:17
非公開
タグ：
posted at xx:xx:xx
@6CLW77Y102 いつもながら、stさんからの重要な情報。
ご紹介ありがとうございます。
検索すると、こういうtweetがありました。
現在も日本語訳は無さそう。
#キリ番が2進数
twitter.com/genkuroki/stat...
タグ： キリ番が2進数
posted at 09:53:44
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
#数楽 実際に自分でその係数の再現に何も知らない状態から挑戦したら、「微妙な点」(Float64専用の係数を作る必要があることなど)がクソ面白かった。
Remez法だけはカンニングしてかつ既存のパッケージを使ってしまいましたが、Remez法も自分で実装すればさらに楽しめたはず。
Remez法は神。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 10:28:16
#数楽 ある程度以上の数学の勉強について「問題集はないんですか？」と質問されたら、「数学の本にあるあらゆる文が練習問題になっているので問題集いらないです」と答えています。
以下のリンク先を見れば、さらに「公開されているすべてのソースコードも数学の練習問題である」ことが分かります。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 10:38:39
出版社の方からもうすぐ発売になる『ベクトル・行列からはじめる最適化数学』をいただきました。重厚感のある書籍で、線形代数の基礎から統計・機械学習・制御・最適化まで幅広く扱っています。サポートページのコードがJuliaというのも重要なポイントですね！ pic.twitter.com/ErnVZLk4Po
タグ：
posted at 10:55:18
#数楽 こういう本は大学で線形代数を教える仕事が割り振られる可能性のある人も見ておいた方がよいと思います。
自分が教えていることが実際にどのように役に立っていて楽しいかがよく分かる。
実際には役に立つレベルの話題をする時間を取るのは難しいのですが。
目次↓
www.kspub.co.jp/book/detail/51... twitter.com/bicycle1885/st...
タグ： 数楽
posted at 11:03:12
TIL that snails hatch from eggs than eat the egg shell to help form their own shells pic.twitter.com/AllLQOpnVg
タグ：
posted at 11:21:10
#数楽
#Julia言語 での線形代数の例
①SymPy.jlとの併用で a₁₂ のような表示の成分を作って表示している。
②LU分解と行列式や逆行列や余因子
gist.github.com/genkuroki/9fa1... pic.twitter.com/U2FDAGlFoD
posted at 11:29:21
【級位者向け実戦型詰将棋part444(塚田 正夫先生作)】
こ れ は 絶 対 に
詰 ま さ な い と い け な い
◯どうしても解けない人用
⇒youtu.be/xodwXr9_gL0 pic.twitter.com/lO4tBXZuLy
タグ：
posted at 11:34:02
@konamih @nabekichi32 下のムスメは高校で
「猫が１モルいます。シッポは何モルですか。」
「１モルです。」
「足は何モルですか。」
「４モルです。」
と、繰り返し学習させられました。
タグ：
posted at 11:44:13
以下の例は #Julia言語 のグルー言語(糊言語)としての側面がよく分かる例。
SymPy.jlはPythonのsympyをJuliaから使えるようにしてくれる。
添付画像の行列計算は(Pythonではなく)Juliaで書かれたコードで実行されている。
Juliaは他言語の機能を多重ディスパッチによって自然に取り込むのが得意。 twitter.com/genkuroki/stat... pic.twitter.com/zmuT2O6IBY
タグ： Julia言語
posted at 12:03:48
a₁₂などの記号の四則演算はPython SymPy側で行われており、行列の計算はJulia側で行われている。LU分解や行列式や逆行列の計算はJuliaで書かれた線形代数の函数で行われている。
SymPyに限らず、任意の記号処理についてもJuliaに取り込むことが可能なはずである。
タグ：
posted at 12:06:50
「結局Pythonのライブラリを使うならJuliaを使う必要はない」のような発言を継続的に見かけますが、ぶっちゃけ、「こいつ、何にも分かっていないな」と思う。
ソフトウェアの世界はそういうものじゃないと思う。
すでにある道具を別の道具と組み合わせて使えるようになることのメリットは大きい。
タグ：
posted at 12:09:54
例えば、Juliaならば既存のよく使われている基本特殊函数のライブラリ(最適化が不十分！)よりも速いコードを書けます。そして、Juliaのコードの最適化の過程で数式処理の計算が必要になったとします。このとき、Julia外の数式処理系をJulia内から使えると便利です。
タグ：
posted at 12:16:23
最終的なアウトプットは「Juliaで書かれた特殊函数の爆速実装」であっても、そのコードを最適化する過程では計算速度的には全然速くない数式処理系(およびその他諸々)も非常に役に立ちます。
タグ：
posted at 12:20:00
最適化の過程では、既存のライブラリとの間の精度や速度の比較が必要になるのですが、Juliaはグルー言語として非常に優れており、比較対象のライブラリをJulia内から使えるようにできる場合が多い。Julia内でベンチマークテストまで一直線に可能になります。
タグ：
posted at 12:22:44
ググれば分かることですが、MITではJuliaを使って以上で述べたようなことがすでに教えられていたりします。
MITでの宿題の答え(指数積分函数E₁(z)のJuliaでの爆速実装)↓
nbviewer.jupyter.org/github/steveng...
最新のJuliaでも動くようにしたもの↓
nbviewer.jupyter.org/gist/genkuroki...
タグ：
posted at 12:27:03
RやPythonで解析・作図を行う学生に対して「Excelでやればいいのに」「Excelだったらもっと簡単にできる」「Excelでやればそんなことにならない」などという教員が現れる現象を”エクセルおじさん症候群”と名付けたい。
タグ：
posted at 12:48:01
結晶解析のモデリングの発表者への質疑応答で、「同じことやりたいんですけどそれってExcelでもできますか」と質問をした教員がいるという話は聞いたことがある。
タグ：
posted at 12:56:29
#Julia言語 私は LOAD_PATH にpush!とか一度もやったことがないので、こういう疑問には答えることができません。
しかし、自分専用パッケージを作ってしまえば(おそろしく簡単、後述)、LOAD_PATHを自分でいじる必要は無くなると思いました。
Juliaでは自分専用パッケージを作るのはものすごく簡単。 twitter.com/wasisama/statu...
タグ： Julia言語
posted at 13:06:47
#Julia言語 自分専用パッケージの作り方
①Juliaを起動
②作成と登録
pkg> generate MyPkg
pkg> dev MyPkg
③内容の変更と使用
julia> using Revise
julia> using MyPkg
julia> MyPkg.greet()
MyPkg/src/MyPkg.jlの内容を変更
Revise.jlのおかげで自動的に反映
julia> My.Pkg.greet() pic.twitter.com/VGMMaFAcWr
タグ： Julia言語
posted at 13:06:49
#Julia言語 もっと本格的にパッケージの開発をしたい人は github.com/invenia/PkgTem... を使えばよいと思います。
しかし、気軽なプログラミングには
pkg> generate MyPkg
pkg> dev MyPkg
だけでも十分に便利です。
参考↓
nbviewer.jupyter.org/gist/genkuroki...
PkgTemplate.jlの使い方
タグ： Julia言語
posted at 13:11:44
#Julia言語 添付画像のワークフローはJuliaでの「お気軽プログラミング」の基本中の基本だと思う。繰り返し宣伝されるべき。
julia> ]
pkg> generate MyPkg
pkg> dev MyPkg
julia> using Revise
julia> using MyPkg
MyPkg/src/MyPkg.jlを編集
編集では任意のエディタを使用可能。 pic.twitter.com/nPXV5D7n3l
タグ： Julia言語
posted at 13:23:06
ごまふあざらし(GomahuAzaras @MathSorcerer
やっぱり，実装した事ない概念は Jupyter 上で もの考えて可視化・テスト書くしかないな(・ω・｀)．
タグ：
posted at 13:27:45
自分でも気軽にパッケージを作れるようになることは、実際にやってみると、かなりうれしいです。
using 素晴らしい既存のパッケージ
using 私の作ったかわいいパッケージ
のように使えるようになる。これがうれしい。😊
そして、自分が作ったかわいいパッケージを公式登録したくなってくる。😎
タグ：
posted at 13:29:46
#Julia言語 「気軽にパッケージを作成」はJupyter notebook内からでも可能です。
]generate Foo
というセルを実行すればよい。
]dev Foo
で登録。JupyterからFoo/src/Foo.jlを編集できる。
添付画像は家庭内JupyterサーバーにiPadのブラウザでアクセスして作ったスクショです。 pic.twitter.com/5JUYupAbQ0
タグ： Julia言語
posted at 13:47:42
Re: RTs 春日氏が何かまともなことを言っているかのようにみなすような人が主要メンバーの活動には近寄りたくないし、私の友人達にも近寄らないで下さいとお願いしたいです。
これが私の正直な気持ち。
タグ：
posted at 14:34:46
須山敦志 Suyama Atsushi @sammy_suyama
「ベクトル・行列からはじめる最適化数学」拝読させていただきました。具体例や演習問題にかなりのボリュームが割かれており、独学に最適です。機械学習入門の際に基礎固めをしっかりしておきたい人に良さそうです。www.amazon.co.jp/dp/4065161967/...
タグ：
posted at 14:50:34
ベータ分布は1940年代までピアソンタイプI分布として知られていた。ガンマ分布もピアソンの研究に由来し、現代の名前を取得する前は、ピアソンタイプIII分布と呼ばれた。スチューデントのt分布は、ピアソンが1895年の論文で導入したタイプIV分布の特別なケースである:-
| en.m.wikipedia.org/wiki/Pearson_d...
タグ：
posted at 15:52:05
うちだと一年目の統計の授業の最終課題として、自分で選んだ論文と同じデータを使って分析結果を再現するというものがあります。教育効果は高いですし、再現できなかったら雑誌にコメンタリーが書けるので業績にもなります。もっと広がればいいなと思う。 twitter.com/rionaoki/statu...
タグ：
posted at 16:14:21
これとは別に、雑誌の方も出版時に分析結果を再現するためのデータとコードの提供を義務付けるべきだと思います。社会科学だと政治学が先を進んでて、社会学は遅れてます
タグ：
posted at 16:18:44
申し訳ないけど、この件に関する一連のツイートには違和感しかないです。
悪意も他意もなく、当時の事情を知らないだけなのでしょうけど。 twitter.com/tsuyomiyakawa/...
タグ：
posted at 16:27:22
#Julia言語 計算速度重視なら
github.com/Jutho/TensorOp...
TensorOperations.jl
が爆速でいいです。このパッケージの存在を知って、配列達の成分の複雑な和のコードを自分で書く気が失せた。 twitter.com/jikanbae/statu...
タグ： Julia言語
posted at 17:08:27
非公開
タグ：
posted at xx:xx:xx
春日匠氏と御用学者wiki問題について、宮川先生の「仲裁」のご努力には敬服いたしますが、以下にリンクする春日氏のブログがなんらかの回答になっているとはとても思えませんし、むしろ「回答拒否宣言」としか読めませんので、僕は春日氏の不誠実さに失望しています
blog.talktank.net/2021/02/blog-p...
タグ：
posted at 17:19:41
2011年以降、暴力的に言論を封殺しようとする風潮がありました。僕自身は暴力宣言を受けましたし、mogmemo氏は暴力的にネットから排除されました。早野先生は「ホームの端を歩かないほうがいい」と忠告されたそうです。
「御用学者wiki」はそのような暴力のひとつとして位置づけられると思います
タグ：
posted at 17:22:20
僕が春日氏に問いたかったのは、2021年から振り返って「御用学者wiki」存在そのものが暴力だったと認めるかどうかです。しかし、事前の予想通り、氏はその問題に直接答えることを避けました。それは不誠実な態度だと考えます。
御用学者wiki問題は僕自身の問題としてこれからも考えていきます
タグ：
posted at 17:27:47
「仲裁」しようとしていることは敬服できることではなく、恥ずべきことだと私は思いました。
春日氏に言い訳のない反省を促すことをきちんとやっていれば敬服できた。しかし、立場的にそれが難しくなっているのもよく分かる。
関わってしまったこと自体が運の尽きだったのかも。 twitter.com/kikumaco/statu...
タグ：
posted at 17:38:14
過去の経緯については私のツイログで検索すれば大体にところが分かると思います。
twilog.org/genkuroki/sear...
私のツイログで「skasuga OR 春日」を検索
見直すのにうんざりするような話です。
タグ：
posted at 17:43:20
#Julia言語 並列化関連資料
www.cs.unb.ca/~aubanel/Julia...
Notes on Mulithreading with Julia
Eric Aubanel
June 2020
(revised August 2020)
タグ： Julia言語
posted at 18:07:58
#Julia言語 上の資料を見つけた検索で同時に発見できたもの
Revisiting Quicksort With Julia and CUDA by Alex Ellison
link.medium.com/TQ0a0QmG6db
タグ： Julia言語
posted at 18:14:19
Qiitaで、一つ記事を修正し終わったので、次はこのプログラムについて記事を書きたいと思います。「有限要素法で、ヘリウム原子に対するHartree-Fock方程式を解いてみる（Juliaでやってみた）」という記事名にするつもりです。 twitter.com/dc1394/status/...
タグ：
posted at 18:18:36
そういえば、この記事のプログラム、CPUをCore i7-7820X（8コア16スレッド）からCore i9-10980XE（18コア36スレッド）に変えたのに、大して速くならなかった（というか若干遅くなった気もする）な…。多コアをフルに使うには、Juliaから生MKLを呼び出さないとダメなのかな…
qiita.com/dc1394/items/c...
タグ：
posted at 18:25:25
非公開
タグ：
posted at xx:xx:xx
#Julia言語 にLibxcを見つけた。これでわざわざJuliaからCを通してlibxcを呼んだりしなくても良くなるな。
github.com/JuliaMolSim/Li...
タグ： Julia言語
posted at 18:43:37
Frames Catherine Whi @oxinabox_frames
This.
I often worry that people don't realize that most(?) heavy open source contributors are (lucky enough to be) doing it "on the clock" . Whether as part of their job, or an academic project (PhD), etc.
And people not knowing this compare themselves, and feel lacking. twitter.com/IMAC2/status/1...
タグ：
posted at 19:17:41
Frames Catherine Whi @oxinabox_frames
If people don't know. Almost everything I do on #JuliaLang AD is part of my day job as a RSE at Invenia. As is NamedDims.jl, and I guess 1/2 my PRs to the Julia language itself, and a bunch more.
DataStructures.jl, and ProjectManagement.jl and ReferenceTest.jl mostly aren't.
タグ： JuliaLang
posted at 19:24:51
ご恵贈いただきました。放射能デマに与しない反原発運動をしてこられた皆さんの共著書です。科学的に冷静な検証がなされているものと思います。 pic.twitter.com/Pr7ogcARPf
タグ：
posted at 20:05:16
反原発運動家から被曝量は多くないとか甲状腺検査は中止すべきとかいう話が出ることに驚くかたもおられるかもしれませんが、この本の著者たちは科学的事実に基づいて議論してきた人たちなので、反原発に有利か不利かみたいな価値判断はされてないのですよ。そして、だからこそ信頼できるのですよ
タグ：
posted at 20:12:37
あけび書房という小さな左翼系(と言っていいんだよね)出版社からの出版ですが、実はあけび書房はつい最近社長が市民社会フォーラムを主宰する岡林さんに替わったばかりで、社長交替後初の出版だと思います。岡林さんは「まともな左翼」なので、今後の出版に期待しています。
買ってあげてね
タグ：
posted at 20:18:03
あけび書房からは吉井英勝さんの本も同時に出ています。吉井さんは東電原発事故をある意味で「予言」した元共産党議員です。まだ読んでないし、再生可能エネルギー推しの推し具合も分からないので、お勧めともなんとも言えませんけど、出ますのでご興味あれば
タグ：
posted at 20:23:50
出版社の方から話題の本を献本して頂きました🎉元々原著を読んだことがあり、素晴らしい本だと思っていたので、翻訳されて嬉しいです。すごい厚みで500ページほどありますが、思い出しながら読んでみたいと思います😃: スタンフォード ベクトル・行列からはじめる最適化数学 bit.ly/3dFzTvU pic.twitter.com/5M8bwNw0lu
タグ：
posted at 20:33:15
New post ([matplotlib animation] 99. マグヌス効果が働く斜方投射) has been published on サボテンパイソン - sabopy.com/py/matplotlib-... pic.twitter.com/uuAWaAEVk2
タグ：
posted at 21:44:50
Initial conditions:
Pendulum 1: position: (0.85, -0.57), mass: 3.75
Pendulum 2: position: (-4.04, 3.82), mass: 4.26 pic.twitter.com/dshJH60zhn
タグ：
posted at 21:45:09
「f(x)の原始関数をF(x)とする。x→∞のときF(x)→0ならf(x)→0を証明せよ」
この問題の証明が思いつきません。数列からの類推でいけそうな気がするんですが…どなたか、数学に詳しい方ご教示下さい。
タグ：
posted at 21:47:48
このシリーズを再実装しててめちゃくちゃ勉強になってるんですが，やればやるほど体系だったプログラミングの知識を持っていないことを痛感してつらくなる #Julia言語
muuuminsan.hatenablog.com/entry/2020/11/...
タグ： Julia言語
posted at 21:57:28
今日、なにげに、Python 3.8.7 で Jupyter を pip install したら、標準で、 Julia 1.3.1 が入っていて、嬉しいしすごく感動した。
タグ：
posted at 22:42:57
この本の一番のウリは教科書の内容の実際のJuliaによる実装例が公開されているところですね。(多分Juliaが数式をほぼそのまま記述できて、そのままplotできることに感動する人も多いはず) Julia Language Companion vmls-book.stanford.edu
タグ：
posted at 23:22:25
ちなみに翻訳者の玉木先生が、こちらのJuliaサンプルを日本語に翻訳して、notebook形式で公開されています🎉 github.com/tttamaki/julia...
タグ：
posted at 23:25:03
非公開
タグ：
posted at xx:xx:xx
