黒木玄 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 を見て下さい。
2019年12月24日(火)
#Julia言語
y₁ + y₂ のような書き方もできるバージョンを作ってみました。次の1行でSymPyの変数達 y₁,...,y₁₀ を作れる。
y = [eval(Meta.parse(numtosub("y$k = symbols(\"y$k\")"))) for k in 1:10]
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/ARUWN8wAsH
タグ： Julia言語
posted at 00:51:57
もぐら叩きできりがないが、
『被乗数も乗数もごっちゃにしてよければ、ますます読解力が育たないし、場面のイメージが持てない子が育つ』
が典型的誤謬。数式に場面のイメージをそのまま表す機能はない。場面を捨象するから数式が普遍的に立てられる。
togetter.com/li/1446316
タグ：
posted at 02:05:16
#Julia言語 あれっ！！！
Julia v1.3.0は速度的に劣化していますね！！！
gist.github.com/genkuroki/e320...
twitter.com/genkuroki/stat... pic.twitter.com/5HRxPhKody
タグ： Julia言語
posted at 03:40:06
#Julia言語 v1.3.1で以上の問題が解決していると良いですね。
discourse.julialang.org/t/julia-1-3-1-...
タグ： Julia言語
posted at 04:13:14
非公開
タグ：
posted at xx:xx:xx
ごまふあざらし(GomahuAzaras @MathSorcerer
不安定というかBinaryBuilderが提供している shared object が armv7l 用の物なので Illegal Instruction 出してるだけでした．
apt で導入した /usr/lib/arm-linux-gnueabihf/ 直下にある so たちを
~/.julia/packages/パッケージ/.../deps/usr/lib に適切な形で置き換えると動作確認できました pic.twitter.com/9l4Z6m0lWq
タグ：
posted at 08:25:02
維新には文化的なもの自体への憎悪のような感情があるんだろう。
大阪から出てこないでほしい。
大阪には大阪の事情があるので、大阪が維新をかついで滅ぶのを外からとやかく言ってはいかんのだそうだし。
タグ：
posted at 08:56:19
昨日のR言語やPythonよりも実行速度が速い Julia
// 式の記述がプログラミング的ではなく数学的に書けるのすきかも
println("HELLO WORLD")
#毎日HELLOWORLD
タグ： 毎日HELLOWORLD
posted at 09:00:00
ごまふあざらし(GomahuAzaras @MathSorcerer
Since BinaryBuilder.jl does not support shared library for armv6l 😅, We should resolve dependencies by ourself to run IJulia on RPiZero.
By applying prescription to RPiZero, it will gain the result like an attached image.
#julialang
Yay. 👍 pic.twitter.com/Ur0Lj3UMhv
タグ： julialang
posted at 09:21:02
ごまふあざらし(GomahuAzaras @MathSorcerer
＃julialang
I’d like to know whether there is a way to switch shared objects to default ones (not <PkgName>/.../deps/usr/lib)
adding LD_LIBRARY_PATH ?? (not tested)
タグ： julialang
posted at 09:25:07
ごまふあざらし(GomahuAzaras @MathSorcerer
apt で導入した so を Julia の適当な場所にコピーして依存解決するという手動黒魔術を使ってしまった．
タグ：
posted at 09:38:45
ごまふあざらし(GomahuAzaras @MathSorcerer
地味にアップデートされてて好き
qiita.com/mt08/items/e8e...
タグ：
posted at 09:56:50
World and Science @WorldAndScience
Relative orbital periods of the planets orbital periods rather than orbital velocities
(Credit: Dr. James O'Donoghue/@physicsJ) pic.twitter.com/IEPFPQv2eu
タグ：
posted at 10:15:19
#数楽 現在の
nbviewer.jupyter.org/github/genkuro...
高校数学の話題
にはベルヌイ多項式を経由しない
3.3 べき乗和の直接的な取り扱い
の解説もある。
タグ： 数楽
posted at 10:19:28
エネルギー転換を論じるなら、何年スパンの話なのかを言わないとまったく意味がありません。30年？50年？それまでの繋ぎはどうする？そういう現実的な話がだいじなんですよ。お花畑のような話しかしない政治家はいらない twitter.com/naotokan/statu...
タグ：
posted at 10:26:14
大学無償化と言っているのは、実質的には無償範囲の縮小で、低所得層の範囲が厳しくなって、これまでなら大学独自の授業料免除を受けられていた層の相当数が授業料免除から外れるということのようですね。これは制度改悪でしょう。なんのために消費税増税したんでしたっけ
タグ：
posted at 11:12:22
@UUQjfllQ68CKPia 内科医です。熱がないインフルエンザがありうるという情報提供ありがとうございます。ただ、必ずしも検査は必要ありません。他の症状が軽ければご自宅で療養するという選択肢がお勧めです。ほとんどの場合、薬を飲まなくてもインフルエンザは自然治癒する疾患です。
タグ：
posted at 11:48:26
@UUQjfllQ68CKPia 「他人にうつすかもしれないので検査が必要」というお考えもあるかもしれません。しかし、インフルエンザであっても検査で陰性に出ることはままあります。検査で陰性だからと大手を振って出勤したほうがかえって感染を拡大するかもしれません。検査の有無に関わらず自宅で安静にするのが望ましいです。
タグ：
posted at 11:49:10
@UUQjfllQ68CKPia リプライありがとうございます。「仰るように検査で陰性で実は罹患していた場合と、普通の風邪だと思って症状が落ち着いたら出歩くのとさして変わらない」。その通りでしょう。
タグ：
posted at 13:58:43
@UUQjfllQ68CKPia ただですね、検査で陰性であった時点で「インフルエンザじゃないというお墨付きをもらった」として、症状があっても出歩く人たちも中にはいるんじゃないかと思います。これは検査のデメリットです。
タグ：
posted at 13:58:58
@UUQjfllQ68CKPia さらに、インフルエンザではない人が「念のための検査」目的に受診して、そこでヤブヘビで感染することもあるのではないでしょうか。流行期の病院の待合室は感染のリスクが高いと考えられます。「念のための検査」目的の受診者が増えれば待合室が込み、待ち時間も長くなります。
タグ：
posted at 13:59:11
@UUQjfllQ68CKPia なので、流行期にインフルエンザを疑う症状があれば、原則として（あくまで原則です）、インフルエンザとみなして自宅待機、というのが理想だと考えます。受診の必要性が高い人は、高リスク群（高齢者や糖尿病患者等）およびred flag sign（呼吸困難や意識混濁等）のある人です。
タグ：
posted at 13:59:30
@UUQjfllQ68CKPia いまだに「検査陽性でない限り出社せよ」という無理解な会社もあります。個別の社会的事情がある人は受診してもかまいません。red flag signがあるのかどうか判断困難な場合も受診してかまいません。
タグ：
posted at 14:00:16
@UUQjfllQ68CKPia 「発熱のないインフルエンザもある」という情報提供は素晴らしかったです。同時に、「インフルエンザについて、検査や薬は必須とは言えない」という原則についても十分に知られるべきだと考えますが、いかがでしょうか。
タグ：
posted at 14:00:54
#Julia言語 Julia v1.3.0で遅くなってしまったのはrand()のようです。4倍近く遅くなっている。
gist.github.com/genkuroki/9b07...
添付画像
1. πをモンテカルロ法で求める函数が4倍近く遅くなっている。
2. rand()を含まないforループならば問題無し。
3. rand()の10回実行が4倍近く遅くなっている。 pic.twitter.com/6kSbUrNy2t
タグ： Julia言語
posted at 15:01:39
#つぶやきProcessing #p5js
C=0
draw=_=>{createCanvas(W=720,W)
D=x=y=W/2
for(i=C++;i<C+150;i++){
x+=cos(noise(i*.01)*PI*9)*5
D*=.965
line(x+D,y+D,x-D,y+D)
if(i%32==0){
for(j of[-D,D]){
line(X=x+j*1.5,y-D,X,y+D)
line(X,y-D,X-j/2,y-D*1.2)
}}}} pic.twitter.com/i5orW34n26
タグ： p5js つぶやきProcessing
posted at 15:14:53
#数楽 コンピューターがらみの例を示すときにも、数学的にちょっと面白い話を混ぜるように努力しています。
sinc(x) = sin(x)/x (x≠0), sinc(0) = 1
のとき、
∫_{-∞}^∞ sinc(x) dx = π
は有名なディリクレ積分。以下のリンク先の添付画像は
Σ_{k=-∞}^∞ sinc(k) = π
の数値的確認。 twitter.com/genkuroki/stat...
タグ： 数楽
posted at 15:24:11
#Julia言語 以上で紹介したように
rand()
がJulia v1.3.0で遅くなっています。
しかし、a = Array{Float64, 1}(undef, 10^7) のときの
Random.rand!(a)
はJulia v1.3.0で遅くなっていません。
モンテカルロ法ではv1.3.0で遅くなった rand() を結構使う。 pic.twitter.com/ALe3xCvpS2
タグ： Julia言語
posted at 16:34:39
これについてるリプを見ていても、まぁ、世間一般の人は「お金を刷るのは禁じ手であってそんなことをしたらハイパーインフレになる」ってな具合に素朴に信じてる人が多いというか、ただ多いだけじゃなくてどうやらそれが多数派らしいんだ。 twitter.com/nikokosan/stat...
タグ：
posted at 16:45:59
@MathSorcerer Are you talking about binaries provided through BinaryProvider.jl? You can remove those in that directory and replace them with a symlink to what you want. However this is a hacky solution that won't work with JLL packages. If you have problems, better report to BinaryBuilder.jl
タグ：
posted at 17:20:15
ごまふあざらし(GomahuAzaras @MathSorcerer
@MoseGiordano Thank you for your advice :+1 . This is, almost same as you said, an example of a workaround that avoids getting an error Illegal instruction on RPiZero devices.
github.com/Julia-Embedded...
May I report my issue to BinaryBuilder.jl ?
タグ：
posted at 17:36:35
@MathSorcerer Yes, I recommend you to open an issue to github.com/JuliaPackaging... since you have problems using binaries created with this package.
タグ：
posted at 17:39:55
ごまふあざらし(GomahuAzaras @MathSorcerer
あざらし，サンタになったでキューーーー．
github.com/Julia-Embedded...
タグ：
posted at 17:40:19
ごまふあざらし(GomahuAzaras @MathSorcerer
@MoseGiordano Alright, Thank you so much!!!
タグ：
posted at 17:40:43
この 300 ページの解説（というか理論体系）を著したグレン・パケットさんが『熱力学：現代的な視点から』の英語版の共著者です。
というわけで、もともとしつこい『田崎熱力学』がさらにしつこい『田崎パケット熱力学』になって世界デビューします。 twitter.com/yubais/status/...
タグ：
posted at 18:09:01
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@UUQjfllQ68CKPia 発熱のないインフルエンザの存在をご理解していただくことはとてもいいことだと思います。ありがとうございます。コストがかかり、他人や自分に感染を広げるリスクを伴う、「念のための検査」の推奨以外の方法で、より多くの人に知ってもらえればいいですね。
タグ：
posted at 18:28:20
「組体操がそんなに嫌ならなんでしてるか調べて、国に訴えたらどうですか？」多くの人がそうしましたけど。その結果、スポーツ庁が中止の方向で勧告を出したのですけど。 twitter.com/uu01203/status...
タグ：
posted at 18:30:38
『21世紀の統計科学』の無償配布。シェアされたサイトがService unavailable になっていますが，こちらの日本統計学会のサイトもダウンロードできます。
www.jss.gr.jp/book/books/ twitter.com/learn_learning...
タグ：
posted at 21:06:17
Julia言語でJulia集合を描いた。色付けがイマイチだけど、ようやくやってやろうと思っていたことができた。 pic.twitter.com/P4tnr2sU4n
タグ：
posted at 21:22:20
How to profile julia code?
#Julia言語
www.juliabloggers.com/how-to-profile... @juliabloggersさんから
タグ： Julia言語
posted at 21:28:22
掛け算に順序の決まりがあるとか世の中の数は全てかける数とかけられる数に分けられると思ってる人は、まず事実誤認しており嘘を教えてることを自覚して欲しい。何のために教えてるのか省みて欲しい / “掛け算順序祭りのフィナーレは、馬鹿教師のお焚き上げ - Togetter” htn.to/45MQDtcLXU
タグ：
posted at 21:46:25
正しい答えを得るための考え方・計算方法は無数にあり、そのどれもが正しい。なぜ機械的に特定の書き方以外を誤りとするのか。まず数的な事実ではないし、世間で一般的な慣習でもない
タグ：
posted at 21:50:20
サンドイッチだのずつだのはじきだの、理解を放棄したわけのわからない謎ルールを機械的に当てはめさせようとする一方で、様々な答えへ至る考え方を理解することを放棄するのは、正当性がある指導態度とは到底思えない
タグ：
posted at 22:00:27
#Julia言語 Julia v1.3.0では rand() を繰り返し呼ぶコードを書くと非常に遅くなります。(v1.2.0以前よりrand()の繰り返しの呼び出しが数倍遅くなる。)
その理由はv1.3.0ではrand()を呼ぶごとにRandom.default_rng()を呼ぶからです。だから、rand(rng)の繰り返しならばv1.2.0以前の速さに戻ります。 pic.twitter.com/hcLQq9XxS9
タグ： Julia言語
posted at 22:12:44
#Julia言語
Julia v1.3.0においては using Random: default_rng として、
rng = default_rng()
for k in 1:N
rand(rng)を使った計算
end
のように書いた方が速い。
nbviewer.jupyter.org/gist/genkuroki...
Distributed and Threads examples
もそれに合わせて書き直しました。 pic.twitter.com/CxjkOveOtB
タグ： Julia言語
posted at 22:24:38
#Julia言語
Random.default_rng() の最新の定義↓
github.com/JuliaLang/juli...
おそらく、スレッド並列のために、これを rand() を呼ぶごとに実行する仕様に変更されたので、rand()が非常に遅くなった。
これってどのように解決可能なんですかね？
自分でコードを書くときには困らないのですが。
タグ： Julia言語
posted at 22:31:07
Particle Life: A Game of Life Made of Particles
video youtu.be/Z_zmZ23grXE
code github.com/fnky/particle-...
demo fnky.github.io/particle-life/ pic.twitter.com/ZBKkrzVytU
タグ：
posted at 22:39:49
