黒木玄 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月25日(木)
#Julia言語 もう何度目になるか分かりませんが、私もdual numbersを実装してみた。
添付画像2の @ code_warntype と @ btime の結果に注目。この2つの確認は必ずやった方がよい。続く
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... twitter.com/abap34/status/... pic.twitter.com/4DUnCw4MXJ
タグ: Julia言語
posted at 00:00:49
#Julia言語 逆行列を求めさせるにはabsとislessも実装する必要があった。Real dual numberの場合にしか実装していない。
absは少なくともComplex dual numberの場合には実装しておくべきなのだが、そこはサボっている。
興味を持った人はやってみるとよいと思う。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 00:24:50
理想ボーズ気体のノートを作っているので、ζ(3/2)の級数がどれくらいで収束するか馬鹿正直に足し算してみたら、1億項くらいでなんとか4桁収束した。これくらいなら許容範囲
タグ:
posted at 00:28:32
@rickasaurus Which videos are we talking about? Also all hip quantum computing physicists are using #JuliaLang
タグ: JuliaLang
posted at 02:00:15
#Julia言語 (1-2^(1-s))ζ(s) は交代級数になりEuler変換が使える。Euler変換を使えば数十項でFloat64の精度を使い切った計算が可能。
①オイラー変換のウェイトのグラフ
②64項の和の計算結果
③正確な値との相対誤差の常用対数、横軸Lは項数
Euler変換はかなり強力。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/kikumaco/statu... pic.twitter.com/uD8yVsYAPb
タグ: Julia言語
posted at 02:28:11
Julia and the Satellite Amazonia-1: Dr. Ronan Arraes shares his experience of using Julia for the Amazonia-1 project and how #Julialang has become part of the Brazilian space program. discourse.julialang.org/t/julia-and-th...
#Aerospace #Satellite #space
タグ: Aerospace Julialang Satellite space
posted at 02:29:52
#Julia言語 これ、非常にいい話。
ブラジルでの宇宙開発でJuliaをv0.3のときからずっと使い続けて来たという話。スクリーンショットも載せてくれている。
宇宙開発の仕事で実際に使っているコンピュータの画面のスクリーンショットって見たいですよね。
NASAでもJuliaを使っている。 twitter.com/juliacomputing...
タグ: Julia言語
posted at 02:45:52
Is there a nullstellensatz for trigonometric polynomials?
mathoverflow.net/questions/3870...
タグ:
posted at 03:01:47
#Julia言語 以下のリンク先の言い方は正確でなかった。
@ timeで測ったメモリアロケーションの回数が異様に多い原因は
① 型安定性が崩れている。特にグローバル変数やstructに注意。
② ループの内側で無駄にメモリ割当が生じるコードを書いている。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 03:04:30
When you lose to your twin brother by 0.001 sec... pic.twitter.com/gMkYJ6L8us
タグ:
posted at 07:42:07
2022年度から施行される新カリキュラムの数学1では、データの分析に仮説検定が入るらしい。
ん、確率を扱う前に仮説検定をやるの?
理屈を説明しないで単にやり方だけ教えるってこと???
タグ:
posted at 08:11:48
#Julia言語 リポジトリの方の公式マニュアルに以下が追加されましたね。
* 引数の型の過剰な制限はよくある間違いです。疑わしいなら引数の型を書くのをやめましょう。
* Juliaでは戻り値の型宣言はほとんど使われません。一般に「型安定」な函数を書くべきです。
github.com/JuliaLang/juli... pic.twitter.com/oGKK9pBZYS
タグ: Julia言語
posted at 10:21:09
#Julia言語 無用に引数の型を制限してしまううようだと、NASAでは仕事をできなくなります。(NASAでもJuliaを使っている) twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 10:30:04
これまで、カリキュラム設計側に対し「高校数学への統計の導入の仕方が強引、数学的な順序が不適切」と沢山批判してきたが、「確率の前に仮説検定を扱う」というカリキュラムを作る時点で本当にやばい。設計側に数学の力が無いと言わざるを得ない。内部に止める人もいないのか。何とかしないと…。 twitter.com/HirokazuOHSAWA...
タグ:
posted at 10:30:42
#Julia言語
function f(x::型名) ~ end
と引数の型を書く場合には型名として最も広い抽象型を選ぶ。よく分からなかったら、型名は省略するべき。
不適切に返り値の型を
function f(x)::型名 ~ end
と指定してしまう問題については以下のリンク先を参照。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 10:33:52
Release v1.6.0 · JuliaLang/julia github.com/JuliaLang/juli... ひっそりとJulia 1.6がリリースされてる。#julialang
タグ: julialang
posted at 10:40:43
#Julia言語 典型的な過学習は以下のリンク先で確認できる。
f(x::Int)::Int = x/2 のように::Intを付けてもJuliaはコンパイル時に返り値がIntになることをチェックしてエラーを出してくれません。f(2)はコンパイルされて正常に実行されてしまう。f(1)は実行時エラー。
discourse.julialang.org/t/why-specify-...
タグ: Julia言語
posted at 11:39:09
#Julia言語 Juliaでは、予想外のバグを取るために型宣言は役に立たず、@ code_warntype や @ code_typed で型情報の伝搬の仕方を見ることが役に立ちます。
その自動化の方向での素晴らしい成果(開発進行中)が
github.com/aviatesk/JET.jl
です。こういう新しい試みをみんなで応援して行くべき。
タグ: Julia言語
posted at 11:59:35
大学生相手であっても、必要な数学の実力は結構高いので、統計学を教えるのに苦労しています。
高校生相手に検定が「お墨付きが得られる道具」であるかのように教えられてしまうようになったら最悪。
あと、信頼区間がモデル依存であることも(大学生と同様に)教えることにならないと思う。 twitter.com/hirokazuohsawa...
タグ:
posted at 12:07:49
現実には世界的にかなり悲惨なことになっていて、論文を日常的に書いている研究者であっても、統計的検定を「お墨付きが得られる道具」扱いしている人達が沢山いるんじゃないか?
そういう現状は若くて優れた研究者が育つことを妨害していると思う。
こういう問題を維持固定しないような教育が必要。
タグ:
posted at 12:10:39
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
【19桁】新しい長さ14の第1種カニンガム鎖、5037591577013475479 が見つかった。進みの遅い 3.4e18からの区間を除き、4.0e18までの処理が完了。現在は 5.4e18までを処理中。19桁はようやく1/3が完了することになるけど、まだ1/2には達していない。ぼちぼちwheel sieveを改良してスピードアップしたい‥
タグ:
posted at 12:21:21
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
#Julia言語 @ code_warntype や @ code_typed の出力は見難いので、
通常のJuliaのコード + 型の伝搬の情報
を表示してくれるとよいような気がする。
誰か、そういう函数を書いてマクロとして使えるようにしてくれないかな? twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 12:38:38
非公開
タグ:
posted at xx:xx:xx
Julia v1.6.0がリリース。 / “Release v1.6.0 · JuliaLang/julia · GitHub” htn.to/42ye7jFkDN
タグ:
posted at 14:23:41
#統計 簡単な確率を計算できる能力抜きに、統計学を理解できるはずがない。
例えば、表と裏が等確率で出るコインを10回投げたとき、表の回数が2回以下になる確率を求められないなら完全にアウト。
高校で統計学を無理に教えずに、統計学を見据えた確率の話をしっかり教えた方がうまく行くと思う。 twitter.com/hirokazuohsawa...
タグ: 統計
posted at 15:00:33
これ、大学の先生でコメントした方が良い人達がツイッターにも結構いるんじゃないか?
直接制度設計に関わったことがない人であっても、大学で統計学の講義をする仕事もしたことがある人なら(多分苦労しているはず)、高校数学への統計の導入の仕方が強引すぎることのまずさは明らかだと思う。 twitter.com/hirokazuohsawa...
タグ:
posted at 15:08:21
#統計 現実の硬い金属製の精巧なコインで試すと、「表と裏の確率はそれぞれ1/2で、10回のコイン投げは独立」というモデルは極めて適切なモデルであることを確認できます。
サイコロでも同様のことをできる。
それらの場合は、数学的モデルの現実への適用が妥当な場合。
妥当性が不明の場合もある。 twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:08:48
気温が上がってきたので、チゴガニのダンスが活発になってきました。
このダンスはウェイビングと呼ばれる行動で、オスが縄張りを主張したり、メスにアピールする目的があると考えられています。 pic.twitter.com/LRfl8Gcx6C
タグ:
posted at 16:15:55
#統計 関連スレッド
accept, rejectという用語選択は不幸な失敗だった。(統計的検定でお墨付きが得られるという考え方は誤り)
フィッシャーさんとネイマン&ピアソンのピアソンさんによれば、統計的検定は「学習の手段」。
「統計的に有意!」は「ひかえおろう!」という意味ではない(笑) twitter.com/genkuroki/stat... pic.twitter.com/AIpl48Dn7k
タグ: 統計
posted at 16:17:26
#統計 新型コロナ関連の推定や予測という困難な仕事でなくても、モデルの妥当性について考えるべき場合は普通にあります。
例えば共通テストの得点の分布の分析では、統計学入門の教科書によくある「正規分布の仮定」は妥当ではありません。
現役生と浪人生では分布が違う。続く twitter.com/genkuroki/stat...
タグ: 統計
posted at 16:26:55
#統計 続き。得点分布を平均と分散の値で代表することは、実質的に正規分布モデルで考えているのと同じです。
現役生と浪人生の特典分布が違っていることが分かっているのですから、単純な正規分布モデルをやめた方が生産的でしょう。
www.ct.u-tokyo.ac.jp/images/koudai-...
twitter.com/genkuroki/stat... pic.twitter.com/QBIMIacjhz
タグ: 統計
posted at 16:26:57
juliaでPlotsを使ってグラフを表示したいのだが、相変わらずvscode上では表示できない。なんだよもう。vscodeでREPL開いてplot(sin)とかやるとvscodeのペインに表示されるんだけど。
タグ:
posted at 17:18:12
『Juliaで学ぶ計算論的神経科学』の書籍化企画が講談社サイエンティフィクで承認されました!出版予定は4年後(卒後3年目)です.体系的かつコード付きの本を目指します.若輩者ゆえ荷が重くはありますが,神経科学・Julia界隈の方々には何卒応援いただきたく存じます.よろしくお願いいたします. twitter.com/tak_yamm/statu...
タグ:
posted at 17:44:06
#数楽 その辺については、中学校の学習指導要領解説(≠学主指導要領、厳密に区別すること!)にもひどく不明瞭な説明があってひどいことになっています。
おそらく数学が苦手な人達の数学教育業界ではそのひどいやつが標準化している。 twitter.com/genkuroki/stat...
タグ: 数楽
posted at 18:46:01
#超算数 「操作の方法」「操作の結果」という説明の仕方がひどく曖昧で悲惨なことになっている。
数学の言葉を使って明瞭に説明せずに、独特の説明で全国の中学校数学教師に悪影響を与える方針。
ひど過ぎ。 twitter.com/genkuroki/stat...
タグ: 超算数
posted at 18:50:27
非公開
タグ:
posted at xx:xx:xx
私も高校の数学Iで仮説検定を扱うことは無理筋との意見に同意だが、前提となる事実認識の確認として、確率の考え方自体は、中学1年および2年で学習するので、確率の考え方を全く知らないまま仮説検定を学ぶことになるわけではない。(またあくまでも数Iで高校数学の学習を終える人向けの話題。)
タグ:
posted at 19:23:05
中1で「多数の観察や多数回の試行によって得られる確率の必要性と意味」が中2から低学年化して、ペットボトルのふたを投げる例とかがここで扱われ、中2では場合の数を基にして得られる確率が扱われている@新学習指導要領 中学数学 pic.twitter.com/gCeLeRYmB2
タグ:
posted at 19:24:42
#統計 中1での確率の話は、某所でのケースでは、教科書本体には掲載されておらず、別の薄い冊子に5ページだけ載っていました。
ペットボトルを投げる話が書いてありますが、全部で5ページなので薄っぺらい内容。
さすがにこれで中1で確率について習ったという扱いは無理。
中2については知らない。 twitter.com/rochejacmonmo/...
タグ: 統計
posted at 19:41:39
@genkuroki 今回「ペットボトルのふたを投げる」話が中1に移行されたことを初めて知りました。以前見たことのある中学の教科書は中2のもので20頁弱くらい。コインを投げて表裏の回数、ペットボトルのふたを投げて表裏横、サイコロをなげてどの目が出るかを記録といった話から「同様に確からしい」を導入し、
タグ:
posted at 19:49:59
@genkuroki 確率を求める例題の中に「2枚の硬貨を同時に投げるとき、1枚は表で1枚は裏となる確率」とか「3枚の硬貨を馬ゲルとき、少なくとも2枚は表となる確率」を求めるものが取り上げられていました。発展的な内容のところにくじびきの話題があり、
タグ:
posted at 19:52:22
@genkuroki 5本のうち、あたりが2本のくじが入っているとき、A,Bがこの順に1本ずつ、ひいたくじを元に戻さずに引くとき、当たりやすさに違いがあるかという問題を樹形図で考えるものなどが紹介されていました。
タグ:
posted at 19:53:26
#数楽 定番のEuler-Maclaurinの公式でζ(2), ζ1.5)などを計算してみました。18~20項でFloat64なら十分な精度でそれらを計算できます。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/MklPU3Q1cm
タグ: 数楽
posted at 20:00:29
ほりたみゅ (mastodonも見てね) @Hyrodium
Julia v1.6.0入れた🎉🎉🎉🎉🎉🎉 pic.twitter.com/DyEsQjIJxq
タグ:
posted at 20:24:42
ほりたみゅ (mastodonも見てね) @Hyrodium
パッケージインストール進捗かっこよくなってるな pic.twitter.com/NCrkZCYdFs
タグ:
posted at 20:27:53
reading: Plots/GR: グラフ package のおすすめ · julia について
www.cas.cmc.osaka-u.ac.jp/~paoon/misc/ju...
タグ:
posted at 21:28:29
Juliaのbroadcastで、ある引数は固定したい(broadcastさせたくない)ときはRefで囲えばいいんですね。。知らなかった。恥ずかしい。。:Avoid broadcasting on an argument in Julia - Stack Overflow stackoverflow.com/questions/4734... pic.twitter.com/zRaQR68JFc
タグ:
posted at 22:29:39
@Atsushi_twi foo.(([1,2],), [3,4,5]) のように1要素タプルのほうが構文的に軽くて好みです。性能はどうか分かりませんが。
タグ:
posted at 22:35:11
#Julia言語 VSCode で Plots.jl のグラフを表示するやつ。
.jl ファイルを実行するのではなく、 REPL を開きっぱなしにして、入力したコードを Ctrl + Enter などで REPL に流し込むことが多い。 pic.twitter.com/SmgDLjqmBL
タグ: Julia言語
posted at 22:38:46
@bicycle1885 おお、ありがとうございます。これでもいけるんですね。知らなかったです😁。自分の手元のベンチマークでは、性能はほぼ同じのように見えます。 pic.twitter.com/yjrwakvOfy
タグ:
posted at 22:50:12
#Julia言語 計算結果とプロット結果を知るために、コードの再実行が不要であることは、#Jupyter の大きな利点です。
計算結果とプロット結果を見たいのに、コードの再実行が必要なせいで、やめてしまうことがある。
diffはコード部分だけに適用したいので、
github.com/mwouts/jupytext
を使う。
posted at 23:02:17
本日3月25日(木)の各都道府県の新型コロナウイルス新規感染者数(18:40時点)を、東京都の人口に換算した値。
※換算値ですので注意して下さい。
※色付きは300以上。
※実数はNHK特設サイト参照
www3.nhk.or.jp/news/special/c... pic.twitter.com/LfTJRNB2Pd
タグ:
posted at 23:04:26
“Julia and the satellite Amazonia-1 - Community - JuliaLang”
discourse.julialang.org/t/julia-and-th...
タグ:
posted at 23:21:17
#超算数 「多項式函数」ブームの元ネタはこれか!
有理函数を函数とみなすときの定義域は有理函数ごとに変わる。
例えば、有理函数1/xを函数とみなすときには0を除外しなければいけないが、有理函数x²/xは多項式xになるので函数とみなすときに0を除外しなくてよい。
本質的にこれだけの話。 twitter.com/mike_twinkle/s...
タグ: 超算数
posted at 23:23:57
非公開
タグ:
posted at xx:xx:xx
shared-mime-infoの開発者がRubyのmimemagicの開発者に苦情。前者はGPLv2であり後者はMITだが、ソースを利用しているのだからGPLv2にすべしと。後者の開発者は直ぐに現行のリリースをGPLv2にし、過去のリリースを削除、GitHubのrepoをarchiveにした。つまり開発は継続しないとのこと。Ruby on Railsを twitter.com/TheRegister/st...
タグ:
posted at 23:37:35
始め、数多くのプロジェクトで利用されており、今直ぐでなくても影響は大きい。Railsチームは現状、回避策を模索しているが、何にしても簡単に実現はしないという。
タグ:
posted at 23:40:37
非公開
タグ:
posted at xx:xx:xx
New post: Julia 1.6 Highlights - www.juliabloggers.com/julia-1-6-high... #julialang
タグ: julialang
posted at 23:52:11
日本語に訳すときに失敗するパターンもあって、行列も行列式もどちらも式なのですが、行列と行列式は違う式になります(笑)
今さらdeterminantを「行列式」と呼ばずに済ますのは難しい。
タグ:
posted at 23:58:41