黒木玄 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 を見て下さい。
2017年12月22日(金)
「これが分かると何が嬉しいかというと」という表現は医学系ではあまり聞いたことがないけど、理学的、工学的、あるいは知的好奇心の意味でこの言葉を使う時に垣間見える、各々が学問をする動機とそこから覗く学問の景色が魅力的で、積極的に使っていきたいし色んな人が使ってくれるといいなと思ってる
タグ：
posted at 00:09:42
@happy__owl いきなりで失礼します。
【場合の数でかけ算使ったら減点されるらしい】
そうなのですよね。
何故か、掛算 でやらずに樹形図を書いて全部数えるのですよね。
hokkaido-sannsuu.com/pdf_sidouan/06... #掛算
タグ： 掛算
posted at 00:15:36
@happy__owl さきほどの学校の「指導案」ですが、この問題の解き方が「ちょっと変わって」いますね。
高校生以上だったらこの問題は普通、４×３ という式にするところなのですが、なぜ逆なのでしょうね。（まあ、「逆は間違い」という物ではないでしょうけれど。） #掛算 pic.twitter.com/i1xYXSHCSV
タグ： 掛算
posted at 00:22:13
#JuliaLang
DifferentialEquations.jl で単振子を数値的に解くときに複数のシンプレクティック積分法を比較してみました。
nbviewer.jupyter.org/gist/genkuroki...
単振子をDifferentialEquations.jlで解いてみる
twitter.com/kikumaco/statu...
タグ： JuliaLang
posted at 00:48:34
#JuliaLang Yoshida6()にすると精度が上がる。
Yoshida6() pic.twitter.com/qZryEm9q1r
タグ： JuliaLang
posted at 00:54:57
#JuliaLang 次に単振子は楕円函数で解けるので、厳密解と数値解の差の大きさを見てみました。VelocityVerlet()だと誤差が非常に大きい。ずれまくり。
VelocityVerlet() pic.twitter.com/mXSMfiAZze
タグ： JuliaLang
posted at 00:56:27
#JuliaLang あまりにずれっぷりに「うっひゃあ」と思ったのですが、Yoshida6()にしたら結構な精度で計算できました。Δt=0.1です。他にも
nbviewer.jupyter.org/gist/genkuroki...
では色々試しています。使用できる常微分方程式の数値解法のリストは
docs.juliadiffeq.org/latest/solvers...
にあります。
Yoshida6() pic.twitter.com/Nrc7XcXTXh
タグ： JuliaLang
posted at 00:59:27
Dr. Chris Rackauckas @ChrisRackauckas
@genkuroki That's a cool notebook. Could that be added to DiffEqBenchmarks.jl? One interesting thing to check as well is how low of a timestep you need for VelocityVerlet to perform well here. It will converge, it just needs a lower timestep than the other methods.
タグ：
posted at 01:00:04
#JuliaLang exact solutionがわかっている場合を全然知らないとかなり困りますね。数値計算の誤差がどの程度かの感覚がつかめなくなる。
タグ： JuliaLang
posted at 01:01:23
非公開
タグ：
posted at xx:xx:xx
#JuliaLang
Hamiltonian だけを与えて常微分方程式は自動的に決定させているのですが、
github.com/JuliaDiffEq/Di...
のVector case testの方のコードを参考にしました。
Scalar case test の方は{at}test_brokenになっていて、「失敗するべき」ということになっているらしく使えません。
タグ： JuliaLang
posted at 01:06:44
#JuliaLang Julia言語のDifferentialEquationsパッケージは超巨大パッケージになっていて、滅茶苦茶たくさんの種類の微分方程式を数値的に解けるようです。
github.com/JuliaDiffEq/Di...
これの詳しいレビューがあるなら見てみたいものです。
タグ： JuliaLang
posted at 01:10:40
#JuliaLang DifferentialEquations.jlパッケージを使う場合には注意するべきことがありました。他に古いパッケージを入れていると自動的にバージョンダウンさせられてしまっていることが原因で、プレコンパイルできなくなった！
結局、私はすべてのパッケージを入れ直しました。
タグ： JuliaLang
posted at 01:13:50
非公開
タグ：
posted at xx:xx:xx
トランプが1日10本コーラ飲んでハンバーガーを何個も食べてあの歳ですげぇなと言われてるが。日本にも1日に板チョコ10枚食べてカルピスやコーラを3リットル以上飲む78歳がいてだな…。 pic.twitter.com/I8ffgBqbDF
タグ：
posted at 08:17:11
反ワクチン批判の村中先生のマドックス賞受賞を渋々ながら新聞が書き始めたのはネット医師達の地道な活動があったからだと思うけど、これは稀有なこと。一般の国民にはここまで強い力はない。
タグ：
posted at 08:49:19
ぐるぐるうづまき 東V30b・2日目・C @guruguruuzumaki
今日は冬至ってことで、いつもの図。
冬至の日は、一年で昼の時間が一番短い日ですが、日の出が一番遅い日でも日の入りが一番早い日でもないのです。
twitpic.com/bnhvxm
タグ：
posted at 08:56:40
#JuliaLang
gist.github.com/genkuroki/563e...
Solving single pendulums by DifferentialEquations.jl
MIT license pic.twitter.com/YC1zN05ojs
タグ： JuliaLang
posted at 09:40:35
#JuliaLang
Very precise numerical solution of the single pendulum solved by SofSpa10, 10th order explicit symplectic integrator
docs.juliadiffeq.org/latest/solvers... pic.twitter.com/DmgFk3kpBu
タグ： JuliaLang
posted at 09:45:18
adhara_mathphys @adhara_mathphys
ハミルトン系の数値解法としてよく使われるシンプレクティック積分法の講義ノートです。絵が可愛いです。
www.unige.ch/~hairer/poly_g...
タグ：
posted at 09:58:12
adhara_mathphys @adhara_mathphys
@phykm 実際のところ色々な講義ノートはありますが、このノートでは色々あるシンプレクティック積分法(例えば吉田公式と鈴木公式)の比較をしており珍しいです。
タグ：
posted at 10:04:06
adhara_mathphys @adhara_mathphys
@phykm アルゴリズムを無心で眺めると、シンプレクティック積分法は時間を進めたり戻ったりするとなぜか精度の良い数値解を得られるように見えます。
リー群・リー代数の公式であるところのBCH公式のなせる業なのですが、何度見ても不思議ではあります。
タグ：
posted at 10:08:51
adhara_mathphys @adhara_mathphys
もう少し新し目の話題は、多様体上で方程式の数値解を求めるというものでしょう。
同著者の講義ノートで
“Solving Differential Equations on Manifolds”
というものがあります。
www.unige.ch/~hairer/poly-s...
タグ：
posted at 10:15:55
adhara_mathphys @adhara_mathphys
ロボット制御などでは平坦な空間というよりは多様体上で物事を考える必要が出てきます。物理系だけではなく機械系でも微分幾何学が必要になってきている背景はここにあります。
タグ：
posted at 10:20:44
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Zenに起きた日本ルールの問題。このまま終局なら左上は白地。黒は劫を仕掛けて白に2手手入れさせれば1目得。この碁は半目勝負なので黒は打たなければならない。東洋囲碁の死活判定ミスで結果は黒勝ちになったけど() pic.twitter.com/jS5Pi9SSw2
タグ：
posted at 13:05:38
財政再建病が酷い
／／…「予定通り財政健全化を中心に使う」と「財政健全化に充てる額を減らし、その分を幼児教育の無償化や高等教育の負担軽減に用いる」…衆院選の当選者は「財政健全化」が１１％、「教育無償化」は８４％だったが、有権者は「財政健全化」が６７％…逆転digital.asahi.com/sp/articles/DA...
タグ：
posted at 13:06:54
事前に親に体罰の許可をもらう「合理的スパルタ塾」の高実績（NEWS ポストセブン） - Yahoo!ニュース headlines.yahoo.co.jp/article?a=2017... @YahooNewsTopics
タグ：
posted at 13:33:31
@YahooNewsTopics これ、合法なの？
当人の同意があっても公序良俗に反する契約は無効だし、そもそも同意しているのは当人じゃなくて親。
暴行罪は親告罪じゃない。
タグ：
posted at 13:35:29
@tajiriyuto まず、終局は両者が着手放棄した時点でできるので（画像1枚目参照）、このまま終局はできないという局面は存在しません。
そして、終局した時点で次は死活確認となるのですが、それは画像2枚目の通りに従うことになります。
という解釈をしているのですが違いますかね… pic.twitter.com/Bdjs791SyP
タグ：
posted at 14:01:23
@spicyspyigo 本コウは手入れ要、一手ヨセコウは手入れ不要で、
あの形で終局しても白は手入れしませんというのは許されなかったような…
終局はできるんですね。
タグ：
posted at 14:13:10
ですよねえ。何だろう、そのような文を書く背景に「何ら根拠や裏打ちや検討が無いにも関わらず何かしら尤もらしい結果を上げる／でっち上げる／まるでそれしかないように教え込み無垢な人々を愚衆に引き摺り込む」という意図や思惑や馬鹿さ加減が見える「気がする」のだけど単なる妄想であって欲しい。 twitter.com/genkuroki/stat...
タグ：
posted at 14:15:51
@sekibunnteisuu @YahooNewsTopics 何故親に許可が出せるのか？ #意味不明 #愚鈍
posted at 14:23:01
出生数 明治32-平成27のグラフ（出展→ www.mhlw.go.jp/toukei/list/dl... ）
twitter.com/hayano/status/... pic.twitter.com/eZcAGwxMX7
タグ：
posted at 15:38:37
岡山県と福井県が実施していたEM菌による水質浄化実験の詳しい報告書を入手。結果はやはり、どちらも「EM菌なし」の対照と比べて効果は確認できず。岡山県の方は、対照の1つとしてオートクレーヴで滅菌処理したEM培養液を使っており、微生物の活動としての効果が確認できないという念の入れ様でした。
タグ：
posted at 16:24:50
ちょっと真面目な、医薬品の適正使用に関するお話です。
最近、浅田飴を水やお湯に溶いて服用してよいのか？という質問を多く頂くのですが、解説させていただきますね。
タグ：
posted at 16:30:33
固形浅田飴は医薬品であり、有効成分を規定の量含有することで、せきやのどの痛みに効果をあらわします。生薬由来の独特の風味を服用しやすくするために、水飴等に溶かしドロップ化しております。
一度にとる量（用法・用量）も決まっており、水には溶かさず口中でゆっくり溶かして服用してください。
タグ：
posted at 16:33:46
一方、食品ののど飴（キャンディー）はハーブの香料などで風味を出したりすることはございますが、医薬品のような有効成分は含みませんし、もちろん効能効果を謳うこともできません。
舐めることでのどを潤すので、これを水に溶くのは問題ございませんが、衛生面にはご注意いただければと存じます。
タグ：
posted at 16:36:08
@t2o_yama @YahooNewsTopics こんなのがまかり通っていいわけがない。刑事事件として立件して欲しい。
sks.nobody.jp/skshp%20situmo...
タグ：
posted at 17:09:02
とりあえず、保存系の問題は6次と10次(あるいは3次と6次)のSymplectic Integratorで解いて、結果が一致すればオッケー、みたいな緩いのりでよさそう
タグ：
posted at 17:47:28
右の画像の式が合ってるのに
左の画像の式はなして違うの？しかも答えは合ってるが式が間違ってるって式は合ってるべや！
教えてエロい人 pic.twitter.com/o1YfEeqoMr
タグ：
posted at 18:06:44
@TJO_datasci あの本，数学をターゲットにしているはずなのに，数学の内容がかなり微妙でしたし，君そもそもPythonでコード書けないでしょ，まともな深層学習モデル書いたこと無いでしょ，でも自分でできないのに本を書かなければならないのか的な事情察してそっと閉じて書棚に戻しました…
タグ：
posted at 19:14:09
ケプラー問題の軌道。GM=1にセット。6次のsymplectic integratorだと、ちゃんと軌道が閉じる pic.twitter.com/WusNMERPdf
タグ：
posted at 19:20:20
Kaggleの良い所は、分析系とかで集めると間違いなく喧嘩する人達を、『まぁKaggleだし』の免罪符でそれなりに建設的な議論が出来る優しい世界を作れる所だと思う
分析官・基盤・Python派・R派を集めて血を見ないのはKaggleだけ説
タグ：
posted at 19:28:09
東京メトロのベイジアンモデリングによる地下鉄トンネル外壁の健全性推定。見事な取り組みだと思います www.datascientist.or.jp/symp/2017/pdf/...
タグ：
posted at 20:26:53
#JuliaLang 添付画像解説
*すべて単振子の数値計算
*左上と右上はq(t)=θ(t)とp(t)=q’(t)のプロット
*左下はエネルギーのプロット。一定になって欲しい。
*右下はy(t)=sin(θ(t)/2)のプロット。数値解と厳密解の比較。ぴったり一致して欲しい。厳密解はy(t)=k sn(t,k).
タグ： JuliaLang
posted at 21:30:39
#JuliaLang
nbviewer.jupyter.org/gist/genkuroki...
単振子をDifferentialEquations.jlで解いてみる
MITライセンスを付けておいたので、誰でもコピー、改変、マージ、配布、販売できます。
タグ： JuliaLang
posted at 21:32:44
#JuliaLang
単振子をDifferentialEquations.jlで解いてみる
nbviewer.jupyter.org/gist/genkuroki...
を見ればJulia言語で自励系のハミルトン系の数値解を求めてプロットする方法がわかります(ハミルトニアンを書くだけでよい)。積分法も複数選択できる。
タグ： JuliaLang
posted at 21:35:47
どんな形であれ、財政緊縮に向かう流れには反対する。金融緩和しても財政出動絞り続けたら意味が無い事くらいそろそろ国民も理解をする時期だよ｜東京新聞:６７％で生活保護費減額、厚労省 来秋、単身高齢者は７６％:話題のニュース(TOKYO Web) www.tokyo-np.co.jp/s/article/2017...
タグ：
posted at 21:58:40
#数楽 シンプレクティック積分法の比較については次のリンク先の講義録を参照。図がかわいい。添付画像はその講義録より。単振子H(q,p)=p^2/2-cos(q)の相図。
単振子ブーム？(笑)
www.unige.ch/~hairer/poly_g... pic.twitter.com/FsBPvcbs2v
タグ： 数楽
posted at 22:59:21
〝この検討会では、HPVワクチン接種後に生じた痛みやけいれんなど多彩な症状の多くは、ワクチンの成分が原因ではなく、「機能性身体症状」であると整理している〟
ここまで言っておいて、勧奨再開の検討しないのは何故？
HPVワクチン 国民に理解してもらえる情報提供を www.buzzfeed.com/jp/naokoiwanag...
タグ：
posted at 23:00:18
もしかして、この検討会には勧奨再開どうこうする権限がないのかな。「積極的接種勧奨中止」という訳の分からん自体は今までなかったので、誰も勧奨再開どうすれば出来るのか分かっておらず、敢えて火中の栗を拾う人もいないと。
タグ：
posted at 23:04:50
非公開
タグ：
posted at xx:xx:xx
