黒木玄 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年09月09日(木)
#キャルちゃんのquantphチェック
密度行列繰り込み(Density Matrix Renormalization; DMR)群をコーディングするためのJuliaライブラリ、DMRjulia。テンソルネットワークを含む基礎演算に焦点を置いている。研究用途にできるほど高速らしい。
Juliaの波、来てるな。
arxiv.org/abs/2109.03120 pic.twitter.com/cjrTp9Gg5r
posted at 00:51:47
新しい記事がZennにアップされました！#Julia言語 #Julia日本語記事
zenn.dev/mizuto/article...👈
タグ： Julia日本語記事 Julia言語
posted at 03:59:02
そうですね。次書こうと思っていました。
twitter.com/genkuroki/stat...
このアルゴリズムはuint32_tをuint64_tにすると正しくなるなる場合があります。
github.com/herumi/blog/bl...
タグ：
posted at 06:15:50
#超算数
これは息子さんの勘違いの可能性が高いが、「順番入れ替えるのは習ってないからだめ」と思い込んでしまったとしたら、それはつまり教え方の失敗である。
仮に勘違いでないなら、それはそれでクソな教え方。 twitter.com/Ngsgr/status/1...
タグ： 超算数
posted at 06:44:25
@herumi @mkashi @hyuki #Julia言語 ジェネリック版を追加。
添付画像
①引数の型から内部で使う浮動小数点数の型を自動的に決定して計算する。非常にシンプル。
②Float64を使っている。
③④BigFloat (内部的にはGNU MPFR)を使っている。
github.com/genkuroki/publ... pic.twitter.com/SztkV8Kk2N
タグ： Julia言語
posted at 07:03:20
@ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine twitter.com/genkuroki/stat...
掛け算の順序を気にしなくてもきちんと理解できていると十分に推定できる資料はありますね。
当然ご存知だとは思いますが。
タグ：
posted at 07:18:27
#Julia言語
Juliaでの自動微分および有限差分による高速高精度の微分の近似および数式処理のパッケージ達の紹介
↓
juliadiff.org
多くの実験が行われていることがわかる。
機械学習的な複雑さを持たない場合にはFowardDiff.jlが適切な選択肢になる場合は結構多いと思います。
タグ： Julia言語
posted at 07:30:13
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 「証明しないと使えない」の機械的な適用は数学の大学院の教育においてでさえ極めて有害です。
子供相手の教育の場合には破滅的に有害だと思います。
「証明しないと使えない」のようなことは高等教育機関で数学を中途半端に習った半可通のいかにも言いそうなことなので、注意した方がよいです。続く pic.twitter.com/k9SLIMPpJ2
タグ：
posted at 07:46:40
Haven't tweeted much recently - hard at work on JuliaSPICE #julialang. Here's a nugget - a simple inverter using the ASAP7 7nm PDK - which one is JuliaSPICE, which one is Xyce ;)? Also, just put out a call for interns for this team :): jobs.juliacomputing.com/jobs/EjfgwN6sh... pic.twitter.com/PJLqOSMXWl
タグ： julialang
posted at 07:46:42
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine まるで低学年の小学生にとってかけ算の可換性が非自明であるかのように語ることも、算数教育の常識に反しています。
採択率1位の東京書籍の小2算数教科書では4の段の九九をやった直後に可換性が一般的に成立する仕組みを学ぶように編集執筆されています。
教科書くらいチェックしておくべき。 pic.twitter.com/T5htTO5HCG
タグ：
posted at 07:48:17
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 学習指導要領でも小2でかけ算の交換法則を教えることになっており、すべての小2算数教科書がそのように編集執筆されています。
算数の教科書さえ確認せずに、自分の心の中だけにある非常識な教え方を前提に、この話題に参加して来て笑い者になる人を、ここ10年間数え切れないくらい見て来ました。
タグ：
posted at 07:51:28
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine あと、現実の日本の教育では、子供は小学校外で算数に関わることを沢山学んでいます。
ある人が「碁石を長方形型に並べてかけ算を家庭内で子供に教えた」と言っていました。そのように教わった子にとって、かけ算の可換性は小2のときにかけ算について教わる前から自明になります。続く pic.twitter.com/pj2AhqoUs8
タグ：
posted at 07:54:26
@adiastra99 Make sure to drop Robert a star: github.com/mossr/Beautifu... if your life really changed : ) #JuliaLang
タグ： JuliaLang
posted at 07:55:35
If you're wondering why the output is at only half Vdd, it's because I just randomly picked a 10k Ohm load resistor, which seemed reasonable to me, but of course this is supposed to model a single (well, two) 7nm transistor, so it can't actually supply enough current.
タグ：
posted at 07:58:18
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 添付画像は現行の学習指導要領解説算数編より(警告：学習指導要領とは異なる拘束力のない単なるいち著作物であることに注意！)。これは問題のある著作物なのですが、それでも小2で添付画像のような教え方を勧めています。続く pic.twitter.com/AfGzZDctG5
タグ：
posted at 07:59:18
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 小2の算数の教科書でも、学習指導要領でも、その解説(拘束力無しの内容的に問題のある文書)でも、かけ算の交換法則は小2で習得する事柄になっています。
この程度のことさえ知らずに、勝手な思い込みを語り続ける様はものすごく恥ずかしのでそろそろやめた方がよいです。 pic.twitter.com/sKKoyj0RPx
タグ：
posted at 08:03:41
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 現実の算数教育では、かけ算の交換法則を子供に教えた後に、子供がかけ算の交換法則を自由に使えないようにするために、かけ算順序が逆ならばバツを付ける問題を出す習慣になっています。
だから、非常識だとみんな非難しているのです。
この事実を自分で確認するべき。 pic.twitter.com/TPmtSxGpaH
タグ：
posted at 08:07:17
@GreatDemon1701 @ProductsCo @kmWLdC6Gbh6tmW1 @golgo_sardine 現実に行われているかけ算順序指導の実態は、算数の教科書のマニュアル本を見ないと十分に見えて来ません。添付画像①は小6での教え方。②③は小2での教え方。
これらの教え方を全否定できない人達には「非常識」という非難が適切。
こういう教え方をする人達のお仲間にならない方がよいです。 twitter.com/genkuroki/stat... pic.twitter.com/tPHOucoRHc
タグ：
posted at 08:11:54
習ってないから使わせない。転じて子供たちは習ってないことはやらなくていいこと、とか思ってしまう
そうなると自学自習のモチベーションって絶対生まれるはずがない。(復習だけやれと言うのならば筋は通っているが、、、) twitter.com/genkuroki/stat...
タグ：
posted at 08:31:00
パターンマッチ教育の横行は、子供の算数力だけではなく、読解力も破壊していると思う。
特に算数が苦手な子ほどひどいことになっている可能性が高い。
ところが、パターンマッチ教育をやっている本人は「算数が苦手な子のためにやっている」と言い訳する。
ここまで狼藉之義也な話は珍しい。 twitter.com/genkuroki/stat...
タグ：
posted at 08:35:38
私が子供の頃嬉しかったのは、先生が言う「よくこんなこと知ってるね！」という言葉だった。
掛け算の順序逆にしてる子供には、「可換なのよく知ってるね！」くらいいってほしいものです。
タグ：
posted at 08:40:17
まお 情報通信技術（旧 ICT）支援員9 @Kyouiku_ICT
オープンソースのJuliaは「ビッグデータとアナリティクスの問題を解決するため、RやPythonなどのデータ解析に適した機能性と、JavaやC++のような生産性の高いプログラミング言語のスピードを兼ね備えている」という。japan.zdnet.com/article/351320...
タグ：
posted at 08:52:04
大久保さん、天むすさん、鰹節猫吉さん、積分定数さん、そして私などが紹介している資料を検索で発見できないまま、自分の勝手な思い込みを前提に議論に参加して来る人達は、情報収集能力にも問題を抱えている可能性が高い。
自分の思い込みに自信を持っている点もこの推測を強くサポートしている。 twitter.com/genkuroki/stat...
タグ：
posted at 08:58:29
子宮頸がんの原因はHPVというウイルスの感染でセックスの経験のある女性の80%くらいが一度は感染します。そのうちウイルスいなくなることが多いけど、い続ける人の一部が癌になるんすよ、という話は検診だったり自分の患者には伝えてるけど、あー知ってますーという人にはまだ会ってない。義務教育よろ
タグ：
posted at 09:23:11
@ProductsCo @GreatDemon1701 @kmWLdC6Gbh6tmW1 @golgo_sardine 言われて当然のことだけで嘲笑されたと感じるような見方をしているから、こんなにみっともないことになってしまったのです。
自分がどれだけ子供達の教育について有害な意見を公言してしまったかを反省した方がよい。
「自由に考え発言する機会」をあなたはこのような形で利用してしまったわけです。
タグ：
posted at 10:13:16
@ProductsCo @GreatDemon1701 @kmWLdC6Gbh6tmW1 @golgo_sardine 「証明しないと使えない」というルールをリアルな小学2年生に適用してよいと思った時点で、教育的に極めて有害なことを平気で公言する人扱いされるのは当然でしょう。
子供を殴ってよいと公言することにかなり近い。
まずは言い訳抜きの反省のひとことがあって然るべきです。
タグ：
posted at 10:17:35
非公開
タグ：
posted at xx:xx:xx
#Julia言語
function msb(x)
ε = eps(float(typeof(x)))
q = (1/ε + 1) * x
m = q - (1 - ε/2) * q
oftype(x, m)
end
これを Int64 型引数で実行すると、Juliaは添付画像①のように分析し、②のネイティブコードにコンパして実行する。
github.com/genkuroki/publ... twitter.com/genkuroki/stat... pic.twitter.com/jSlLxsNOoA
タグ： Julia言語
posted at 11:19:12
未だにイベルメクチンの在庫切れで疥癬の治療にフェノトリン外用でお茶を濁さないといけない状況、本当にどうにかして欲しいです。医師の裁量権で適応外処方でも詳記すれば大丈夫とはいえ、そのせいで本来の適応疾患に処方できないなんてどう考えても異常です。
タグ：
posted at 11:21:17
流してたけど Julia って面白い言語ですね。機械学習なんかにも向いてる。識別子に UNICODE が使えるのと、関数を数式そのままに記述可能とゆー特徴が合わさって…
長方形の面積( 縦 , 横 ) = 縦 * 横
とか出来ちゃう。これで実行速度は C++ 並なんですねぇ。
あ、正規表現が Perl 準拠だ(はーと)
タグ：
posted at 11:35:40
今回のは0.01M HCl 10mLに正露丸1粒を完全溶解したときの話。
もともと胃液は空腹時に50mLなので、この条件だと1回5粒の相当（しかも胃で完全に溶解してその濃度に達する保証はない）。正露丸の服用は成人1回3錠で、そもそも空腹時は粘膜への負担があるとして食後に水・湯で服用が推奨されてるわけで。
タグ：
posted at 12:28:22
もちろん、溶け出すときの錠剤周辺では局所的に高濃度になるだろうけど、たまたまそこにアニサキスがいるのを想定するのか、というと…。結局のところ、今回のはあくまでin vitro実験にすぎない（あと、胃なら内視鏡下で外科的に除くのがセオリーだから、腸アニサキス症に効く方がありがたいんだが）
タグ：
posted at 12:32:42
非公開
タグ：
posted at xx:xx:xx
内視鏡するときの患者の負担を考えると薬で「散らせる」方がいいのも確かだけど、アニサキス以外の食中毒を患者が自己判断で何とかしようとする方向になるとまずい。確か、以前どこかでカンピロバクター中毒の人が自分で何とかしようと正露丸飲んでた話が問題になってたと記憶してる。
タグ：
posted at 12:38:04
@GreatDemon1701 @genkuroki 小学生のころ学校が一番正しいと思って生活してた自分を思い出すと、親がこれは○でいいんだよ～とかいっても、「いや、×なんだよ。なにをいうてるん」って思ってしまったから、学校の功罪は大きいですね
タグ：
posted at 12:45:17
音響に関係する自前関数を #Julia言語 向けにいくつも作っていたんだけど、JuliaDSP github.com/JuliaDSP に同じような関数がそろって来てるので、そっちとの互換性を高める方向にしたい。
タグ： Julia言語
posted at 13:10:14
研究会の資料をGitHubにまとめておくことは、Google Driveにまとめておくよりもずっと優れた方法だと思います。
単にGitHubに資料をまとめることを超えたことをしたい研究会主催者は、以下のリンク先の門脇さんの解説を参考にすると良いかもしれません。 #Julia言語 twitter.com/kado_judo0312/...
タグ： Julia言語
posted at 13:19:03
０円公開！
０円ギャンブル！！
【嘘喰い｜ヤンジャン！】 ynjn.jp/app/title?tid=... #ヤンジャン #嘘喰い
posted at 13:31:29
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@benkyouaho 公式ドキュメントに忠実に、まず @ code_warntype の結果を確認するとよいと思います。次を実行してみてください。ただし @ の直後の空白は除く。
@ code_warntype HybridMonteCalro(S, collfunc, niter, nτ, dτ, ndim)
@ code_warntype MolecularDynamics(S, nτ, dτ, ndim, zeros(Float64, ndim))
タグ：
posted at 16:46:36
@benkyouaho あと、函数の定義での引数の型が狭過ぎる点が問題になる可能性があります。自信がなければ引数の型は書かない方が無難です。
例えば、Juliaにはfunction-like object (←公式ドキュメントの検索を参照)というのがあって、Function型ではないオブジェクトが函数のように振る舞う可能性があります。
タグ：
posted at 16:50:19
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
恥ずかしいことに配列への列追加のところですごくブサイクな書き方をしていたのでコードを修正しました
ブログを始めてから
コードは動けば良いでしょという意識を悔い改め
更生しようとしてます笑
【Julia】製造DBのデータをCSV形式で読み込み編集する方法 cae-memo.com/julia-csv/
タグ：
posted at 21:04:05
今日は「途中式を書いてなかったから全問正解なのに50点になってたんです…」というご相談があった。
内容はくり下がりのない引き算の文章題の復習プリント。その子は「読んだら答えが分かった！」と式を書かなかったそう。文中にも式を書けとの指示はなく、(式5点、答5点)の記載もない。
はぁ……。
タグ：
posted at 21:07:58
もちろん、柔軟に考えて子どもの力を伸ばして下さってたり、その為に忙しい中身を粉にして一生懸命に子どもの成長を支えて下さってる先生方もたくさんいらっしゃる事は重々承知しているんですが、最近こうした「教員の考え以外は全部✖️」みたいな話が多い……。
タグ：
posted at 21:35:51
@kado_judo0312 多分、その辺のノウハウについて知りたい人は沢山いると思います。
もしもそちらに質問が集中して迷惑をかけてしまったら、ごめんなさい。(前もって謝っておきます。)
タグ：
posted at 21:37:03
非公開
タグ：
posted at xx:xx:xx
@benkyouaho Union{Nothing, Tuple{Int64, Int64}} は無害なので無視して下さい。
あれ？ 私がやってみたら強い警告色でAnyが表示されましたが、そちらではそうならなかったですか？
Anyはとてもまずいです。
タグ：
posted at 21:39:26
@benkyouaho #Julia言語 サンプルコード
① github.com/genkuroki/publ...
② github.com/genkuroki/publ...
添付画像は②です。②はほぼ同じことをやっているつもりで書いたコードです。バグっていたら直します。しかし、解読すると(おそらく結構大変)、色々コツがつかめると思います。 pic.twitter.com/P9rSaG09TB
タグ： Julia言語
posted at 22:01:52
@benkyouaho Anyが出る理由は様々なのですが、公式ドキュメントもPerformance Tipsの最初に書いてあるグローバル変数の参照が怪しい。
S(x) = S(A, x) は引数を経由しないグローバル変数Aを含んでいます。
この辺は色々大変。私のサンプルコードのスタイルで全体を書き直せば、この問題が自然に解消されます。
タグ：
posted at 22:05:49
@benkyouaho discourse.julialang.org を見ても、多くの新参者がグローバル変数を常に引数として函数に渡すコードを書くことに苦労しているようです。
S(x) = S(A, x) のように左辺に含まれないグローバル変数Aを右辺で使うとアウト。
こういうことで困らずに済むスタイルが上のサンプルコードで示されています。 twitter.com/genkuroki/stat...
タグ：
posted at 22:15:12
@benkyouaho そのスタイルはプログラミングの世界では特別なスタイルではないです。少なくともクラスベースのオブジェクト指向言語でのプログラミングをやったことがある人にとっては普通です。
プログラミングの世界で普通のやり方をすれば多くの面倒ごとが自然消滅し、爆速のコードが出来上がります。
タグ：
posted at 22:17:45
@benkyouaho S(x)の引数はxしかありません。それを定義するコードの中で引数を経由せずにグローバル変数Aを使っているので、当たり前にアウトです。
函数の中でAを使いたければ、引数にAが含まれている必要がある。
これで困らないようにするためには全体のスタイルを変える必要があります。
タグ：
posted at 22:21:20
@benkyouaho 件のコードは
①leapfrog法でハミルトン方程式を数値的に解くこと
②HMC法でサンプルを生成
の2つの問題が含まれているので、Juliaに全く慣れていない人が「最適化ゲーム」に挑戦するのは難し過ぎるかもしれません。
「①だけ」「②だけ」を経由して「①も②も」に挑戦すると楽だと思いました。
タグ：
posted at 22:46:59
TUTORIAL IN BIOSTATISTICS
Bayesian workflow for disease transmission modeling in Stan. Léo Grinsztajn, Elizaveta Semenova, Charles C. Margossian, Julien Riou. Statistics in Medicine. onlinelibrary.wiley.com/doi/10.1002/si...
タグ：
posted at 22:47:11
この両者は一見すると互いに反しているが、実際はどちらも「決められた型に当てはめてこたえをだしましょう。それ以外のやり方は認めません」という思想がベースになっている。いずれにせよマジキチなのには変わりないが twitter.com/moritatsu/stat...
タグ：
posted at 22:52:20
#Julia言語
低次元空間内での常微分方程式の数値解法では、StaticArrays.jl を使うことが、Juliaでは定跡化されていると思います。
ハミルトン方程式の数値解法のコードを、SVectorの初期値の型が伝搬するように書くのは易しい。solve函数の素朴なコードを見て下さい。
github.com/genkuroki/publ... twitter.com/genkuroki/stat... pic.twitter.com/asn2jjby5y
タグ： Julia言語
posted at 23:00:46
#Julia言語 しかし、ハミルトン方程式の数値解法のコードを、ポテンシャル函数をφ(x)とするときの分布 p(x) = exp(-φ(x))/Z のサンプルをHMC法で作るために利用する場合には、運動量ベクトルを計算途中で多次元正規分布でランダムに与えるので工夫が必要になります。 pic.twitter.com/8FVXdnQX4l
タグ： Julia言語
posted at 23:00:49
Rational Tangles should be called "Magic with Numbers"
(Thanks @PBSInfinite and @math3ma for such wonderful video!!)
bit.ly/2X5F0KZ
#math #science #iteachmath #mtbos #visualization #elearning #numbers #geometry pic.twitter.com/rVFPK5PSuM
タグ： elearning geometry iteachmath math mtbos numbers science visualization
posted at 23:01:26
#Julia言語 SVectorを型安定に作るためには、値が確定しているdimについて SVector{dim}(randn!(rng, vtmp)) のようにする必要があります。
そのようなコードを含む函数は dim の値の違いでディスパッチするようにしておく必要がある。
これはJulia新参者には無理に近い難しいことだと思います。 pic.twitter.com/PqQMpCboCb
タグ： Julia言語
posted at 23:04:11
#Julia言語 leapfrog法をどのように行うかの情報が含まれるstruct LFProblemの最初のパラメータをdimとしておき、函数達を
function f(lf::LFProblem{dim}, ～) where dim
～SVector{dim}を含むコード～
end
の形式で書く必要がありました。
github.com/genkuroki/publ... pic.twitter.com/nKWdSBBw0F
タグ： Julia言語
posted at 23:07:47
コロナのデータを用いて、Stanの微分方程式の機能を使ってmodeling (SIR model + ロックダウン)を改善していく例が論文になっています。チュートリアル論文なのでコードや図も多くあって親切でした。 twitter.com/triadsou/statu...
タグ：
posted at 23:12:32
式の意味云々を言ってる筑波大学附属小学校の田中博史氏は、みはじではなくて４マスなのですが、
まぁ、みはじも４マスも大差ないですな。
#超算数 twitter.com/moritatsu/stat...
タグ： 超算数
posted at 23:30:38
@benkyouaho 以下のリンク先の反応は間違っているので注意してください。
「～のはず」じゃダメで、実際にコードを書いて検証することが科学的な態度です。
S(x) = S(x, A) はAnyを発生させ、添付画像のような害が小さな場合にさえ、大量のメモリアロケーションを発生させます。
github.com/genkuroki/publ... twitter.com/cometscome_phy... pic.twitter.com/hZG94DmteA
タグ：
posted at 23:42:27
@benkyouaho 私もおかしなことを言う場合があります。
そういう場合には、小さなコードとベンチマークを突きつけて「間違っています！」と教えて頂ければ助かります。
タグ：
posted at 23:47:12
非公開
タグ：
posted at xx:xx:xx
