黒木玄 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 を見て下さい。
2018年01月23日(火)
#数楽 あと、MITのこの講義全体の面白さはすごすぎ。
github.com/stevengj/18S09...
これは去年のやつで今年のやつは
github.com/stevengj/18S096
タグ： 数楽
posted at 00:08:34
非公開
タグ：
posted at xx:xx:xx
速度でJuliaに勝つのが難しかった話． 本日は Twitterのハッシュタグ Julia言語 で話題となっている（というかしている）Ising ModelのJuliaのコードが速いというおはなしです． Ising Model Ising… ift.tt/2DruNOX
タグ：
posted at 00:14:04
@hase3desu @sunchanuiguru 【わり算の順序を間違えるから掛け算の順序を言い出す人は…】
それは彼ら特有の錯覚から来ているのです。
彼らの筋書きでは「まず単価があって、それを数量だけ倍する、という仕組みを理解すれば、それを順序に反映させるべきと分かるだろう。だから、順序を守らないのは仕組みを（続 #掛算
タグ： 掛算
posted at 00:46:35
@ceptree 結構制約が多いんですよ...。未対応の演算が入っていてもエラーを吐かずにjitをスルーしますし。なお、オプションで（nopython=True）とするとjit出来なかった時にエラーを吐いてくれます。ご存知かもしれませんが、対応しているnumpyの型はこちらに書いてあります。
numba.pydata.org/numba-doc/dev/...
タグ：
posted at 00:52:27
なんか、こういう物言い一つ一つに、「真理は天から降って来るもの」とでもいうような、天下り的なナニかを感じるなあ twitter.com/Aki_utakata/st...
タグ：
posted at 01:00:22
もうちょっと真面目に、ローレンツ方程式を前進オイラー法から4段4次のルンゲ=クッタ法に書き換えた。最初Pythonで書いてたけど、通常のやり方やとnumbaが無名関数に使えなかったのでJuliaで書いた。 pic.twitter.com/qN3juAXGnC
タグ：
posted at 01:33:04
津田和俊／急激に進行した網膜剥離と闘って @kaztsuda
さっき除雪車が家の前を通ったので、ちょっとチェックしに表に出たら、えらい中途半端に雪寄せしてくれたなぁ。。。これじゃぁ明日の朝の子供らの登校が酷いわ、朝早起きして雪かきしないとダメか。と、もうこんな時間なのだが f(^_^;;;
タグ：
posted at 02:06:08
津田和俊／急激に進行した網膜剥離と闘って @kaztsuda
しかしホント、仙台市の除雪は下手クソやけど、他の豪雪の街のように社会インフラとして雪捨て場を確保してないので、しゃーないわなぁ。。。
タグ：
posted at 02:08:45
津田和俊／急激に進行した網膜剥離と闘って @kaztsuda
ただ雪がフカフカなので、クルマの雪払うのはさすがメッチャ楽やなぁ。
これでこの３倍くらい積もってくれれば、青葉山スキー場（謎）もパラダイスなのだが。。。
ま、あとは、明日の朝にしよ。
タグ：
posted at 02:13:21
twitter.com/k_b_ms/status/...
#Julia言語 JuliaPro とそのMKL版の違いは使っている数学ライブラリの違い。通常版はopenBLASを使っている。MKL版は
ja.wikipedia.org/wiki/Intel_Mat...
を使っている。MKLへの対応は完璧ではないのでよくわからない人は通常版選択が無難。
タグ： Julia言語
posted at 02:14:34
教員の負担軽減についてですが、外部指導員に過度な期待をするのはやめたほうがいいと思います。教員にとってのブラック状態が緩和される代わりに生徒にとってのブラック状態は加速する
タグ：
posted at 02:32:51
部活狂顧問はそれでも教員なので最低限のブレーキがついている場合が多いけど、あいつらにはそんなもん最初から無いしむしろ教員的価値観を否定するのが格好いいと思ってるような連中だったりするよ
タグ：
posted at 02:35:34
斉藤 淳 『アメリカの大学生が学んでいる @junsaito0529
日本の数学、ちょっとヤバイと思うのは、高校や学習塾の数学講師がじつは数学分かっていないのではないかと思わせることが多いこと。弊塾で数学講師を募集しつつづけているのだけど、当然分かっているであろうと前提している解析学や線形代数の問題が解けてない。
タグ：
posted at 02:38:35
斉藤 淳 『アメリカの大学生が学んでいる @junsaito0529
採用試験を課して、がちで数学できる講師が、ほとんどいない。塾業界、予備校業界で教歴が立派でも、数学で修士号持っていても、数学が分かっているとは限らないというのが、だんだん見えてきた。分かっていないのに、分かっているふりして教えているんだなー、多分。いまでどうやってたのか不思議だ。
タグ：
posted at 02:40:52
斉藤 淳 『アメリカの大学生が学んでいる @junsaito0529
例えば掛け算順序問題を眺めていると、塾や小中学校で教えている先生が、「学ぶ歓び」を忘れたまま「教える歓び」を利己的に追求した結果なんじゃないかと思えてくる。少し進んだ分野から眺めて、何が重要か分からないと、上手には教えられないと思うんだけどなー。
タグ：
posted at 02:51:35
みーくん | itmz153 @math153arclight
Burrau's problem of three bodies pic.twitter.com/2BZwZVB29z
タグ：
posted at 03:09:26
非公開
タグ：
posted at xx:xx:xx
Julia、開発は2009年に始まってるけど、2012年2月14日に例の記事が出たのが一応区切りやねんな。バレンタインデーはJuliaの日？
Why We Created Julia
julialang.org/blog/2012/02/w... pic.twitter.com/60mfO0CLw3
タグ：
posted at 03:26:00
Juliaでローレンツ方程式をRK4で解くやつ、関数の定義を関数の中に入れたら、2.016秒 → 0.013秒になった。Pythonだと関数の中に関数は対応してるけど、numbaは対応してないっぽい。関数の外で関数を定義した場合だと0.322秒。 pic.twitter.com/1jbZ8XOpHf
タグ：
posted at 03:54:56
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
単に「不要不急の外出を控えろ」と言っても何のことだか判断しかねる人が多いと思われるので、「不要不急の外出とはあなたの仕事のことだ」と言ってやる必要があると思うが、それ言われたらみんな怒りそう。
タグ：
posted at 04:55:44
World and Science @WorldAndScience
Anamorphic image of Jules Verne pic.twitter.com/xaSQpJbSNT
タグ：
posted at 05:01:08
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
本来、抽象的でかる数学を、無駄に具現化して思考能力を破壊している部分を批判しているのだと思いますが。
【教員を追い詰める世論を煽るなんて】ことはしていません。閉鎖的な教育空間で、理不尽に追い詰められているのは子供の方です。
＃初耳学 ＃掛算 ＃算数
twitter.com/shirowassaaaan...
posted at 07:23:01
@mywarstory 私とかが言っているだけなら、おっしゃる通りですが、フィールズ賞の受賞者かつ国際数学連合のトップまでもが、冒頭からそう言われたというのがショックという話。それと金がないというのは、冷暖房費がないとか、事務職員が減って雑用が増えているという部分で、研究費や給与の面ではないです。
タグ：
posted at 07:25:09
外部指導員が絡むトラブルを何件か目撃、あるいは当事者になったことがあります。ヘタな人にあたると任せるどころか言動をチェックする手間暇が増えて泥沼になります。卒業生の中から人選できればある程度はマシかもしれませんが、過去のつながりがない顧問では難しいです。 twitter.com/t2o_yama/statu...
タグ：
posted at 07:30:47
非公開
タグ：
posted at xx:xx:xx
@mywarstory 京大理学部数学科の同窓会なので、スピーチの後は、、各大学の数理系学科に在籍する大学教員で雑談していましたが、ほぼ一致した意見でした。その後、何人の研究室も訪問し、いろいろ状況を聞いて、地方大学と違ってまだましかと思ったら、そうじゃなかったということですね。
タグ：
posted at 07:52:38
非公開
タグ：
posted at xx:xx:xx
【教育現場】にいる（いた）人の中に、「小学校は学問をするところではなく、子供に従順さを身につけさせる（調教する）場所である」という感覚を持っている人が少なからずいることが怖い。
＃初耳学 ＃算数 ＃掛算
twitter.com/mayfly627/stat...
posted at 08:22:10
@mywarstory 私も別に数学が終わるとは、微塵も思っておりませんよ。そこに引っ掛かるのであれば、数学科に訂正しときます。昔は、年配の先生が事務を肩代わりしてくれて、若手は研究だけをしていればいいという雰囲気だったのが、最近はそうもいかなくなっているという傾向はあるようで。RIMSはましな方なのかな。
タグ：
posted at 08:35:37
非公開
タグ：
posted at xx:xx:xx
で、 #Julia言語 開発の本拠地のMITでの講義を発見したので見てみたら、メモリのキャッシュを意識してアルゴリズムを考えないと高速計算できないという普遍的事実をJulia言語を使って教えていた。
nbviewer.jupyter.org/github/steveng...
気楽さとこういうのが同居していて面白すぎ。
タグ： Julia言語
posted at 09:09:57
ごまふあざらし(GomahuAzaras @MathSorcerer
ミスリーディングかもしれないけれど、実際に組み込んでJuliaとその他言語とで勝負させて速度と実装の素早さ容易さの観点で痛感する経験は損にならない。 twitter.com/genkuroki/stat...
タグ：
posted at 09:20:40
Julia、試してみたい気もするが文法が気に食わんなぁ。endを書く言語好かんのよね。Pythonみたいにインデントでブロック構造を書けたらよかったのに。
タグ：
posted at 09:22:15
@mywarstory 米国の大学は、研究の面では、はるかに厳しいのでしょうが、日本の大学だと、教員が、入試の警備、校内清掃の草むしり、印刷物のホッチキス留め、不登校学生の保護者との面談、合宿のホテルの予約、自習時間の監督とかまで、持ち回り平等でやらされるので。時間を金で買いたい、という感じですかね。
タグ：
posted at 09:30:47
工夫をするとそれにこたえて高性能を発揮してくれる道具を使うことは楽しい。
MITにおける #Julia言語 を使ったコンピューターの講義はまさにそういう内容になっており、Julia言語を使う予定がない人であっても閲覧する価値があると思う。めっちゃ楽しい内容。
github.com/stevengj/18S09...
タグ： Julia言語
posted at 09:41:22
非公開
タグ：
posted at xx:xx:xx
＞同研究所は、不正防止のために研究内容を厳しく管理する体制を敷いている。３カ月に１回は全ての研究者の実験ノートを確認し、論文を発表する場合は元データや画像などの提出をルール化している。
うわっ
mainichi.jp/articles/20180...
タグ：
posted at 09:50:39
@mywarstory ええ、昨年は、出張経費の扱いについて、各大学の教員がSNSで某政治家に苦情コメントを寄せ、文科省に働きかけて改善されました。みんな大喜び。私は、大学同期でIT関係で成功した知人から、それなりの寄付をいただいているので、それで凌いでいます。が、大学全体の経費にまでは口出しできません。
タグ：
posted at 09:55:31
有期雇用どうし競争させて恐怖と不安を煽ってたら、そりゃ業績あげるために不正も出てくるだろう。おおもとの過酷な制度には手をつけずに不正を減らそうとして単純に罰を厳しくしていけば、真面目にやってる大半の研究者がただ疲弊していくだけにならへんかなあ。
タグ：
posted at 10:00:09
非公開
タグ：
posted at xx:xx:xx
「若手研究者の雇用を不安定化させることによって見事に研究不正を生産することに成功しましたね！今のお気持ちをお聞かせください：」
って文科大臣にインタビューしに行こうず(´・ω・｀)
タグ：
posted at 10:48:29
@integers_blog integers.hatenablog.com/entry/2018/01/... を拝見して今朝から5の場合がどうなってるのか考えてるけどパッと脳内では図が浮かばない。そうなると最近 #JuliaLang に嵌ってるから総当たりで答えから探すってのもアリかと思うが，どうすれば良いかで止まる。いやはや数学力が欲しい。先ずは紙に書いてか。
タグ： JuliaLang
posted at 11:03:57
@YoichiTakahashi @YoichiTakahashi オレと中川さんとは、政調会長と幹事長時代にいろいろと繋がりがあったので自民党勉強会を手伝った。発起人は中川さんと山本さん。そこに浜田先生を呼んだ。安倍さんがリフレになったは中川さんの影響だよ
タグ：
posted at 11:28:07
非公開
タグ：
posted at xx:xx:xx
#Julia言語
nbviewer.jupyter.org/gist/genkuroki...
view マクロの使い方
function test!(B, A)
B .= A.^2
end
N = 10^8
A = reshape(collect(1:2N), N, 2)
test!(@ view(A[:,2]), @ view(A[:,1]))
@ の後の空白は除く。
2次元配列の第j列を第j-1列からの計算結果にしたい場合にはこうする。
タグ： Julia言語
posted at 14:12:57
#Julia言語 応用例
nbviewer.jupyter.org/gist/genkuroki...
FFTをin-place計算で使うことによって無駄にメモリを消費せずにKdV方程式の数値解を求める方法
in-placeにする前は 1.019 s, 1.10 GiB だった計算が、 in-placeにしたら 837.614 ms, 20.48 MiB になりました。メモリ消費が50分の1に減りました！
タグ： Julia言語
posted at 14:17:31
解決した．サイトには書かれてないけど PiGPIO って StrPack 使うんだね
#RaspberryPi #julialang
タグ： julialang RaspberryPi
posted at 14:22:50
nbviewer.jupyter.org/gist/genkuroki...
にある2Dイジング模型のテストコードは、型指定を行ったり、最適化のおまじない(inboundsの類)を使ったり、アルゴリズムでさらに工夫したりすることを、意図的に一切しないようにしています。
特に擬似乱数を32ビット浮動小数で生成するのは御法度。 #Julia言語
タグ： Julia言語
posted at 14:42:54
非公開
タグ：
posted at xx:xx:xx
#Julia言語 Xorshift系の擬似乱数生成器をどう評価すればよいかは知りません。しかし、Julia言語のパッケージの文書を見るとXorshift64は使わない方がよいと書いてありました。その使わない方がよいXorshift64はデフォルトのdSFMTより確かに速いです。でもこれはダメらしい。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 14:53:35
#Julia言語 以下のリンク先はJulia言語のRandomNumbers.jlのXorshift familyに関するドキュメントです。
sunoru.github.io/RandomNumbers....
タグ： Julia言語
posted at 14:55:53
#Julia言語 個人的な意見ですが、Julia言語がデフォルトの擬似乱数発生器としてdSFMTを採用したことは慧眼だったと思います。
www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/SFMT...
タグ： Julia言語
posted at 15:02:11
@genkuroki XOR+SHIFTを64ビットでやる乱数発生機はだめでしょうね。もっと長くないと。僕はベクトル型スーパーコンピューター用に9689ビットのXOR+SHIFTを30年くらい使っていますが、特におかしな挙動には出会っていません。短いやつではひどい目にあったことがあります
タグ：
posted at 15:08:21
専門の人に聞くべきだけど、私の理解では現代的な疑似乱数は
①0～(周期-1)までの数字が一つずつ書かれた紙をシャッフルして1枚ずつ上から見る（種を指定すれば何枚目から見るか選べる）
②実際に「使う」のはカードの数字の一部の桁だけである）
③シャッフルの良さはアルゴリズム依存
という感じ
タグ：
posted at 15:29:17
その分野で学術的意義のあるものが書けるような専門性はないのでwikipeの紹介程度しかできませんが、疑似乱数発生器から出てきた数列を様々なテストにかけるほうが主流かと @genkuroki
en.wikipedia.org/wiki/Randomnes...
タグ：
posted at 15:33:28
@kyslog @genkuroki wikipeのリンクにXorshiftの開発者自身が書いたランダムネステストの論文がありますね。 www.jstatsoft.org/v07/i03/paper
タグ：
posted at 15:48:34
#数楽 続き。添付画像のように、
(-N/2～N/2-1の整数の組(a,b)で最大公約数が1のものの個数)/N^2
は6/π^2=1/ζ(2)に近くなります。
二つの整数(a,b)の最大公約数が1になる「確率」は6/π^2なのだ。 pic.twitter.com/6GevEmVccU
タグ： 数楽
posted at 15:56:03
@genkuroki ポロック予想はいかがですか。主記憶を使い切ったところでスワップを始めてカクンと遅くなるのも観測できます。d.hatena.ne.jp/wd0/20100324/a
タグ：
posted at 16:08:05
#Julia言語
-N/2からN/2-1までの整数の対(a,b)で最大公約数が1であるものの個数cを単純にカウントする。そして、sqrt(6N^2/c) を計算する。円周率に近い値が得られる。
私のパソコンではN=20000で2.8秒かかりました。4億回のループ。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/f9bnHL8s5V
タグ： Julia言語
posted at 16:30:52
さてさて、何時ものように #ruby と #LuaLaTeX に #MetaPost で「平面で3点が同一直線上にない整数距離四角形で最大距離が73」のものリストを描かせてみた。こういうのが #JuliaLang だけでサクサク描けるようになるのやろか。 pic.twitter.com/U5pLz0L4ZT
タグ： JuliaLang LuaLaTeX MetaPost ruby
posted at 16:34:28
@MathSorcerer あ、ごめんなさい。面白いので、どんどん単純速度比較をやって下さい。結構楽しみにしている人達がいるように見えます。無粋なことを言ったみたいでごめんなさい。
タグ：
posted at 16:48:30
ごまふあざらし(GomahuAzaras @MathSorcerer
@genkuroki いえいえ．無粋だなんてとんでもないです．実際にやってみるとPython+Cythonなどでも最適化の作業がバカにならないですし，他の言語で実装するのはとても勉強になります．
タグ：
posted at 18:19:28
ごまふあざらし(GomahuAzaras @MathSorcerer
@genkuroki そして自分の力で書ける言語が尽きてしまったのでQiitaで呼びかけました．
qiita.com/SatoshiTerasak...
タグ：
posted at 18:20:29
非公開
タグ：
posted at xx:xx:xx
娘中2が受けた地元大手塾の外部模試、縮小添付された回答でどう見ても合ってるのに誤答の漢字が複数あって疑問に思ってたら、塾基準のとねはね払いに合わないとバツということらしい。塾生は基準を知ってるから高得点、外部の子は誤答になる。塾生の順位を上げ、勧誘する仕組みに使ってるわけね。
タグ：
posted at 18:52:31
乗算演算子の左が掛けられる数で右が掛ける数やろうとりあえず。まぁどっちでもいいけど
#超算数 twitter.com/katsukichi1974...
タグ： 超算数
posted at 19:11:58
なんで有期雇用なんだろうね。研究なんて何年も何の芽も出なくて、ある時突然花開くことがあるなんて、誰にでも分かりそうなもんなのにね。そこにお金をかけられない今の日本の現状に、ただただ悲しくなるね。
タグ：
posted at 19:33:26
ごまふあざらし(GomahuAzaras @MathSorcerer
@genkuroki おかげさまで１７００View突破しました．いつもは３００View得られれば御の字です．
タグ：
posted at 20:18:00
#数楽 #GeoGebra
以下のリンク先で紹介されているGeoGebraのページは私がテキトーに作ったもので、Edwards曲線と呼ばれる楕円曲線(の表示における)加法構造が何であるかをプロットしたものです。つまんで動かせます。
Edwards曲線で楕円曲線暗号で有名。
twitter.com/genkuroki/stat...
posted at 20:19:34
無料で利用でいる #Mathematica な #WolframAlpha は超絶便利ではある。この程度のルート入りな計算ならルートのママやり取りして距離の関係が見事に整合していることがものの数分で確認できた。こいう事が #ruby や #Python や #JuliaLang でサクサクできる事を期待してるが出来るのか？ pic.twitter.com/NsFDJS9Efm
タグ： JuliaLang Mathematica Python ruby WolframAlpha
posted at 21:29:23
高橋将宜 Masayoshi Takah @M123Takahashi
「誤差が±3%なら、7%と思え」というタイトルのニューヨークタイムズの記事。元ネタはGelman et al.による論文で、それを一般向けにしたもの。世論調査の標本誤差が±3%の場合、無回答などの非標本誤差を考慮すると、経験的に全体の誤差は±7%はあると考えられるという話です。
nyti.ms/2dwGxUu
タグ：
posted at 21:39:15
ご安心ください。さすがに内生変数と外征変数合わせて４０００近くある大きなモデルなので被保険者数の減少などはちゃんと取り込めてます。こっちのページの方程式リストのファイルを「被保険者数」で検索してみて！
計量経済モデル及び試算関係資料
www5.cao.go.jp/keizai3/econom...
twitter.com/maseguchi/stat...
タグ：
posted at 21:42:13
サーバはスタンダードコース。ウェブコンテンツ＋データ用に100GB。独自ドメインを取らなければ最低費用515円／月でメール，ウェブ，WordPressのブログ他自作ウェブアプリ利用も。sshログインできてメンテフリー。広告の入らない個人サイトとしてとてもありがたい。
www.sakura.ne.jp
タグ：
posted at 22:46:43
@konamih 個人的に小波さんの作品群がGitHubにもあると安心感があります。
GitHubアカウントにファイルが置いてあれば、個人サイトが消えても残ります。
私も自分の作品群の一部をGitHubにうつしました。
genkuroki.github.io
ウェブサイト
github.com/genkuroki
タグ：
posted at 23:09:59
twitter.com/seiyab_/status...
#Julia言語 では特別なライブラリやパッケージを読み込まなくても、線形代数的なことは一通りできます。
FFTを使えばKdV方程式(最も有名なソリトン方程式)を容易に数値的に解けます。
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/gist/genkuroki...
タグ： Julia言語
posted at 23:22:18
Nimのマクロだいぶ書きにくい…。ASTのノードオブジェクトを逐次生成してがっちゃんこするか、文字列をparseStmtしてASTオブジェクトをつくるしかないなんて…。ほかの言語のマクロシステムはもっとマシなんだろうか…。
タグ：
posted at 23:44:41
非公開
タグ：
posted at xx:xx:xx
これがJuliaのマクロか…。NimよりはASTオブジェクトが生で書ける感がある…（NimはnimNode("hoge")とかやってノードのインスタンスを作る感じなのでちょっとつらみある）
docs.julialang.org/en/stable/manu...
タグ：
posted at 23:54:23
