黒木玄 Gen Kuroki(@genkuroki)/2021年02月/Page 5

ちょっと息抜きにJuliaで3次元イジングモデルを作り始めた。動くんだけど、合ってるかな

タグ：

posted at 00:10:44

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
@ timeと@ code_warntypeの使い方

添付画像は目次と第5節。

docs.julialang.org/en/v1/manual/p...
を全部読んで理解できれば問題ないのですが、普通はそういうわけにも行かないので、@ timeと@ code_warntypeで確認する癖をつけておくと、経験値を貯めやすいと思います。 pic.twitter.com/0vYxCdZJLh

タグ： Julia言語

posted at 00:19:50

本屋で中学数学の参考書を立ち読み。
「倍数」の頁は「０は全ての整数の倍数」とちゃんと書いてあった。

タグ：

posted at 00:22:03

市販の参考書は高校受験を意識しているからまともなのかな？

タグ：

posted at 00:23:08

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...
@ timeと@ code_warntypeの使い方

の第4節で作った動画。1次元の波動方程式の数値解です。 pic.twitter.com/OozcaXd7Hu

タグ： Julia言語

posted at 00:23:33

#Julia言語 @ code_warntypeで強い警告色が出ないことの確認を忘れると、比較的容易に数十倍計算が遅くなる。

しかし、@ code_warntype で強い警告色の表示が出ないようにして、さらに @ time で無駄なメモリアロケーションが生じていないかを確認して改善するだけで、相当に速いコードに化けます。

タグ： Julia言語

posted at 00:26:52

#Julia言語私が使い始めたときには、

@ time の表示はメモリアロケーションに注目

しなければいけないこととか、

@ code_warntype の表示に強い警告色での警告が神

だとか、誰も教えてくれなかったので、相当な時間を無駄にしました。この2つを繰り返すだけで、様々なことに理解が進みます。

タグ： Julia言語

posted at 00:29:55

himako @himkt

requirements.txt にバージョンが書かれてないやつ、 pip の新リゾルバが「バージョンわからんけどとりあえず総当たりして組合せ見つけるわ！」っつって各ライブラリの全リリースをダウンロードすることに気づいて「オアー！」って言った（普段 poetry でバージョン管理してるので知らなかった）

タグ：

posted at 00:33:38

himako @himkt

github.com/pypa/pip/issue...
（インターンに来てくださってる学生さんの docker build が永遠に終わらなくて調査した結果見つけた

タグ：

posted at 00:40:56

#Julia言語 @ timeの表示の読み方の例

添付画像のコードでは、U, Vがそれぞれ400×8001のFloat64(8バイト)の行列なので、合わせて400×8K×8B=48.9MB程度のメモリは消費されるはずです。

ところが、@ timeが表示した数値は456MBでひと桁違っていて、アロケーションに回数も以上に多い。続き pic.twitter.com/03GJwLzYX5

タグ： Julia言語

posted at 00:48:37

#Julia言語続き。ガベージコレクションによる遅延もひどい。

こういう見方ができれば、添付画像のコードでは失敗していることは間違いないと判断できます。

失敗に気付けば色々調べて改善して経験値を増やせます。

自分の失敗に気付かないままだとかなりの時間を無駄にすることになります。 pic.twitter.com/qJhFTQXGBW

タグ： Julia言語

posted at 00:48:44

#Julia言語無駄にメモリアロケーションを増やして、ガベージコレクション地獄に陥ることは、重い計算をJuliaにやらせようとした人は一度は経験していると思う。

失敗に気付いた後に公式ドキュメントの

docs.julialang.org/en/v1/manual/p...

を読むと、

ああ、こんな所に助けの神がいたのか！

となります。

タグ： Julia言語

posted at 00:52:15

#Julia言語 @ code_warntype の出力は、最初のうちは、強い目立つ警告色で表示されている型情報のみに注目するだけで十分だと思います。

強い警告色のAnyが大量に現れていたら、それが原因で計算速度が数十倍遅くなってしまっている場合があります。

添付画像では20倍ほど遅くなっている。 pic.twitter.com/F0LF1KyfS4

タグ： Julia言語

posted at 00:56:42

#Julia言語添付では、問題を指定するデータを1つの型にまとめること(これは良いこと)をしているのですが、失敗しているせいで、@ code_warntype で赤字のAnyが出まくっています。私の環境では20倍ほど遅くなっている。 pic.twitter.com/5zxJ31eXHR

タグ： Julia言語

posted at 01:00:10

#Julia言語型を

struct Foo
a
b
end

ではなく、

struct Foo{Ta, Tb}
a::Ta
b::Tb
end

と定義して、中身の具体型はJuliaに自動的に同定させることにすれば、1つ前のツイートの20倍の速度劣化を防げます。

この方法ならJuliaの型システムについて勉強しなくても使える。

タグ： Julia言語

posted at 01:03:53

おかゆ @oka_iu_tcan

juliaってwebサーバー建てれるんだ、科学計算の言語としか思ってなかった

タグ：

posted at 01:08:17

#Julia言語 700ページ近くのJuliaのコード入りのPDFが無料。 twitter.com/sammy_suyama/s...

タグ： Julia言語

posted at 01:15:06

#Julia言語

【パフォーマンス向上のために部分的に静的型付けをすることができる】は誤りだと思います。

函数の引数の型を

function f(x::Float64)::Float64
～
end

のように指定してもパフォーマンスは何も向上しません。

function f(x)
～
end

と書いた状態でC並に速いです。 twitter.com/poyothon/statu...

タグ： Julia言語

posted at 01:19:44

#Julia言語

Juliaは函数の実行直前にネイティブコードにコンパイルして実行する仕組みなので、

function f(x)
～
end

のように型名を一切書かずに定義された函数に関する f(3.0) 実行直前のコンパイル時にコンパイラは引数3.0の型(Float64)の情報を利用できます。続く

タグ： Julia言語

posted at 01:24:12

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

#Julia言語「事前の静的コンパイル＋静的型付け」の感覚で「実行直前コンパイル＋コンパイル時に引数の具体型が分かる」のJuliaについて語ろうとするから、静的型付けができることがパフォーマンス向上に必須の事柄であるかのように誤解する。

Juliaは「静的型付け必須教カルト」を潰す側の存在。

タグ： Julia言語

posted at 01:28:54

Juliaで並列エラトステネスの篩（4並列くらいで頭打ち） pic.twitter.com/oHgN6ZtUF1

タグ：

posted at 01:39:07

JuliaHub @JuliaHub_Inc

ITensor is a software library for tensor algorithms in Julia and C++. The Julia version of ITensor now offers multithreaded contraction of block-sparse tensors. For more information click here itensor.org
#julialang #multithreading @ITensorLib pic.twitter.com/ztjLneXJFA

タグ： julialang multithreading

posted at 03:41:10

Marius Millea @cosmic_mar

Masahiro Ono 小野　昌弘 @masahirono

Just heard about this package during @zackli's talk.

Wigner 3j symbols in a couple 100 lines of #JuliaLang, 2-4X faster than the prototypical Fortran code thats been used for decades.

github.com/xzackli/Wigner... pic.twitter.com/JhIpoUnNTz

タグ： JuliaLang

posted at 05:50:41

イスラエルでワクチン接種の流行抑制への効果がでてきている可能性。イスラエルは接種率が相当先行している（60歳以上は大多数が接種済み）ので、ワクチンの流行への効果をみる試金石となる。ワクチン接種がほぼおわった６０歳以上の新規感染者数・入院数の低下が顕著。今後の経緯が注目される twitter.com/jburnmurdoch/s...

タグ：

posted at 06:48:44

medy.nim lang: @dumblepytech1

@mpyw WebはNimへ、データサイエンスはJuliaの方へ行くという尖った願望

タグ：

posted at 09:25:50

ブログ「身勝手な主張」の２０１４年の記事で、中学校では最大公約数と最小公倍数についての項目がないことが書いてあって興味深い。
ユークリッド互除法について少し触れているけど、高校では教えることになっているけど、それ以前はどうなんだろう？
＃超算数

タグ：超算数

posted at 09:27:15

ブログ「身勝手な主張」の２０１４年の記事

blog.goo.ne.jp/mh0920-yh/e/e8...

＃超算数

タグ：超算数

posted at 09:29:15

nagataka @nagataka

Pythonしか書かないマンになって久しいので、新しい言語のチュートリアルとかを研究の気晴らしにちょっとずつだらだらやりたいと思ったりする。Julia気になる（一方、C++書けるようになりたいとも思い続けて早X年、一向に手を動かさないワタシ。数学に無限に可処分時間を喰われてる）

タグ：

posted at 09:40:42

Lee Phillips @lpfeed

A list of Julia links: lee-phillips.org/julialang/

#JuliaLang

タグ： JuliaLang

posted at 09:58:00

Go Ando / PREDUCTS @goando

Webで人々の意思決定を操るデザイン「ダークパターン」

例えば退会プロセスでしつこく繰り返し慰留したり、セール価格をカウントダウンと共に表示する。

アメリカではカリフォルニア州を皮切りにこうしたデザインが禁止対象になり始めている。
@wired_jp wired.jp/2021/02/05/law...

タグ：

posted at 10:08:55

Daichi Mochihashi @daiti_m

Juliaのプロンプトはjulia> ですが、私はElkのような無駄な情報のない> という寡黙なプロンプトが好きなので、少し調べると、 ~/.julia/config/startup.jl に次のように書いておけばよいようです。 pic.twitter.com/j4z9vE95dg

タグ：

posted at 10:38:21

@biregulartree

ケフっち＠Deterministic @__cheph__

(3年前期) 機械学習
線形回帰等を julia (疑似python) で実装しようぐらいの内容
ぶっちゃけ取ってもあんまりって感じ
先生は優しいしすごくモダンな手段で課題やらせてくれるから好き
動画も資料も1回も見ずに課題できたのは草
まぁ単位稼ぎ＆軽い演習ってノリでとれば良さそう

タグ：

posted at 10:47:50

なんかな〜 JuliaでCSV.jl のreadでエラー吐かれる。昨日までは正常に動作してたんだけどなぁ〜あとから考えよう。とりあえずは、メシだ

タグ：

posted at 12:06:14

なお @nao_tuboyaki

Julia、Rustのどっちが楽しいかな？

タグ：

posted at 12:12:01

積分定数 @sekibunnteisuu

gendai.ismedia.jp/articles/-/798...

こんなくだらない分析する人はZ会の主任といこと？

だとしたら、Z会も・・・ pic.twitter.com/kYtGWzm5Mp

タグ：

posted at 12:37:28

積分定数 @sekibunnteisuu

かにになりたいね @kanikanimochi

Z会、

　共通テストにおべんちゃら言わないとならない事情があるのか？それとも本心で太郎と花子の会話の問題を素晴らしいと思っているのか？

　どちらにしても、駄目すぎる。

タグ：

posted at 12:42:38

かに pic.twitter.com/n6Z2YlOPHC

タグ：

posted at 13:01:37

おかゆ @oka_iu_tcan

juliaしててPythonの事実に気づいたけど、たしかにPythonってFor文でスコープ作らないのな（やばない？）

タグ：

posted at 13:09:48

Ukicode @opaupafz2

前にも書いたと思いますが、俺は今からPython学ぶぐらいならJuliaを学んだほうが今後のためにもなるかもだし良いと思っています。
まだJulia使ってないPython使いが言うけど、Pythonマジで遅すぎだし本当書き方に気を付けないとすぐ遅くなるからね。処理速度重視するぐらいならJuliaのほうが良い。

タグ：

posted at 13:13:16

Shoei @shoei05

"SMRは人口規模の小さい地域では変動が大きいという欠点があり、それを補うために、階層ベイズモデルを用いた縮約統計量の利活用が推奨されています。(www.kansai-u.ac.jp/step/img/publi... ;www.myschedule.jp/jpa2014/tex_ou... )" ytake2.github.io/takeblog/2018/...

タグ：

posted at 13:14:23

Ukicode @opaupafz2

PyPyやNumba、Cython使えやって話になるかもしれないけど、PyPyとNumbaはCPythonじゃ実装できなかったはずだし、CythonはPythonと書き方を変えなければまったく恩恵が得られないと言って良い。

タグ：

posted at 13:24:53

マウス @gakuseitdn

pythonだとスクラッチ実装したとき遅いから、juliaでやってみようかな。
単にコーディングがヘタクソな問題もあるけど…

タグ：

posted at 13:45:05

むっちょ @mucchyo

もうC言語禁止法出してJuliaに置き換えましょう。それがいい。

タグ：

posted at 14:56:46

#統計そこで紹介されている論文のタイトルは Beyond Newton's methodなのですが、ガンマ分布の最尤法についてそこで提案されている計算法は F(α) = α(log α - ψ(α) - c) = 0 に関するNewton法に一致しています。

全然Newton法を超えていない！(笑)

シンプルな話です。添付画像の解説を参照。続く twitter.com/nenomake/statu... pic.twitter.com/YlYY4AdU6k

タグ：統計

posted at 14:59:07

#Julia言語のDistributions.jlでのガンマ分布の最尤法も

F(α) = α(log α - ψ(α) - c) = 0

に関するNewton法で実装されています。

クイズに出したことがあった。

この話題は、ディガンマ函数ψ(α)、トリガンマ函数ψ'(α)と最尤法とNewton法に馴染むには非常に良い演習問題です。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 14:59:09

#数楽一般に x > 0 の範囲で f(x)=0 の解を求めることと、F(x):=xf(x)=0 の解を求めることは同値です。f(x)がf(x)=g(x)-cの形をしていて、xg(x)が下に凸な正値単調減少函数なら、F(x)=xg(x)-cxは局所的に直線に近くなり易くなり、F(x)=0へのNewton法の適用は効率的になります。 pic.twitter.com/jHtqM2pUWI

タグ：数楽

posted at 15:05:20

@hoxo_m #統計そこで紹介されている論文のタイトルは Beyond Newton's method ですが、ガンマ分布の最尤法について提案されている方法はNewton法です。

詳しい解説が以下のリンク先にあります。

その方法は #Julia言語のDistributions.jlでも使われています。爆速で非常に便利です。 twitter.com/genkuroki/stat...

タグ： Julia言語統計

posted at 15:11:37

相互リンク twitter.com/genkuroki/stat...

タグ：

posted at 15:12:09

函数 y = f(x) に関する方程式 f(x) = 0 について、xに関する変数変換のみを考えて、yの側の変換について考えなかったせいで、Newton法とは違う方法を見つけたと誤解してしまったのだと思います。

しかし、提案されている方法は爆速で素晴らしいです。

タグ：

posted at 15:15:07

こういう類のことはよくあって、「説明はおかしいし、本質を捉えてもいないが、結論は正しくて素晴らしい」となっている場合が結構あるので、説明や本質については自分で議論を再構成することによって理解する努力をした方が無難な場合が多いです。

タグ：

posted at 15:17:23

ややこしい数学がからむ話題では他人の議論にうしろから着いて行くと幸せになれないことが多い。

タグ：

posted at 15:20:04

FEMでヘリウム原子に対するHartree-Fock方程式を数値的に解いてみたけど、基底で展開する場合と比べて精度が低かった。具体的には、6つのGTOを使って解いた場合はエネルギーE = -2.860 (Hartree)だったけど、FEM（5000分割）で解いた場合はエネルギーE = -2.852 (Hartree)だった。SCFの収束も良くない. pic.twitter.com/bxBOPEGwuw

タグ：

posted at 15:44:55

Nimda @nimdanaoto

Hartree-FockってFEMで解けるのか twitter.com/dc1394/status/...

タグ：

posted at 15:46:55

#Julia言語

【WignerFamilies.jl is thread-safe and very fast, beating the standard Fortran routine DRC3JJ from SLATEC by a factor of 2-4.】

昔から使われている標準的Fortranのライブラリよりも速いコードをJuliaで書くといういつもの話の1つ。

常態化している。 github.com/xzackli/Wigner...

タグ： Julia言語

posted at 15:50:29

@nimdanaoto ヘリウム原子の場合には比較的簡単に解けます。一般の原子、分子に対するHartree-Fock方程式を解く場合は、Hartree-Fock交換ポテンシャルを求めるのが極めて困難なので、ヘリウム原子のようにはいきません。

タグ：

posted at 15:50:50

#Julia言語 twitter.com/dc1394/status/...

タグ： Julia言語

posted at 15:54:45

@genkuroki ソースコードはこちらで公開しています。
github.com/dc1394/helium_...

タグ：

posted at 15:56:07

教科書が #超算数

＞板橋区では読み解く力の育成を目指して
独自に「板橋 i（あい）カリキュラム」の開発
「指導書から授業をつくるのではなく、一から教科書を読んで授業を作る」
www.s4e.jp/blogs/blog_ent...

板橋区で使用している教科書
小学校用教科書　　算数大日本図書
www.city.itabashi.tokyo.jp/kyoikuiinkai/g... twitter.com/takusansu/stat...

タグ：超算数

posted at 16:46:36

ｲﾝﾔｰロボver3.9 @_HeyYeahHeyYeah

新井から、板橋区立の学校には1800人の先生がいて
一人ひとりが「教科書を使い倒す授業」を改めて考えることが大切
教科書をもう一度読み直してみてほしい
ここを試験に出すここがポイントと思わずに
まずは教科書を読むことを第一歩としてやってみてほしいとアドバイス

www.s4e.jp/blogs/blog_ent...

タグ：

posted at 16:49:02

非公開

タグ：

posted at xx:xx:xx

ミーハーだからRustとかJuliaとか触りたいンゴ

タグ：

posted at 17:07:29

TaKu @takusansu

みなみちあぷ〜( ⩌ ̫⩌)o @chiaki_at_cake

@Tei393098 「指導書」を鵜呑みにするなという認識が広がったとすれば、一歩前進したと信じたいですね。
「指導書」に批判的な声が上がれば、芋づる式に教科書も吟味する声に繋がっていって欲しいです。

タグ：

posted at 17:17:47

Juliaなる言語を遅いとは思いながら手を出したけど、これ書きやすいうえに速くて大変良い。

タグ：

posted at 17:38:46

@takusansu 板橋区は新井紀子教授のリーディングスキルテストを導入し読解力教育を推進してるのですが、新井教授の算数が超算数なんです。算数だけならまだしも
超算数的な読解法を他の教科にも取り入れてるのです。
AIに負けない読解力なので前衛的教育かも。

twitter.com/genkuroki/stat...

タグ：

posted at 17:40:34

#Julia言語複数の最適化パッケージを、いつもの

prob = OptimizationProblem(～)
sol = solve(prob, ～)

のスタイルで統一的に利用できることを目指しているパッケージ。

問題を意味する型FooProblemを定義して、solve(prob::FooProblem, ～)を定義する普遍的なスタイル。 github.com/SciML/Galactic...

タグ： Julia言語

posted at 17:46:12

@takusansu 『ハッピーになれる算数』は超算数の指導法だそうです「AI vs.教科書が読めない子どもたち」が同じ問題意識だそうなので、AIに負けない読解力は超算数指導と相性が良いということだと思います。

twitter.com/Tei393098/stat...

タグ：

posted at 17:46:28

#Julia言語問題を記述する情報を格納するグローバル変数や定数を定義して函数の中でナマで使っているコードをなぜかよく見るのですが、問題を記述するデータを保持する型FooProblemを定義して、問題を解く函数を

function solve(prob::FooProblem, ～)
～
end

と定義するとスッキリします。

タグ： Julia言語

posted at 17:50:15

Massimo @Rainmaker1973

The stack effect (aka the chimney effect) helps drive natural ventilation, air infiltration, and fires. This fire instructor demonstrates the chimney effect to trainees [read more: buff.ly/3jq2xlC] [source of the clip: buff.ly/2MY0yZA] pic.twitter.com/zLMevjpFm8

タグ：

posted at 18:05:53

#Julia言語問題を記述するデータを保持する型の定義でデフォルト値は

Base.@ kwdef struct FooProblem{Ta, Tn}
a::Ta = 1.234
n::Tn = 10^3
end

using Parameters
@ with_kw struct BarProblem{Ta, Tn}
a::Ta = 1.234
n::Tn = 10^3
end

で設定可能。

gist.github.com/genkuroki/8358... pic.twitter.com/bi2wVyfGd3

タグ： Julia言語

posted at 18:06:56

#Julia言語

const a = 1.2
const n = 10^3

function solve()
定数a,nを使ったコード
end

ではなく、

using Parameters

@ with_kw struct FooProb{T,U}
a::T=1.2
n::U=10^3
end
prob = FooProb()

function solve(prob::FooProb)
@ unpack a, n = prob
～
end

がお勧め。

タグ： Julia言語

posted at 18:12:26

#Julia言語 struct の定義自身の上書きはできないので、後で定義を変えたいstructの定義はmoduleで包んでおいて、moduleをまるごと置き換えるようにすると便利です。

Revise.jlと組み合わせると便利。

タグ： Julia言語

posted at 18:15:37

たふぇたふぇ- 大学生 @tafetafe_sns

#Julia言語

問題を記述するデータを1つにまとめる型を定義するときに、「自分が解きたい問題はどのように一般化されるのだろうか？」と考える機会が得られることも有用。

「問題と解法」という考え方は非常に普遍的。

タグ： Julia言語

posted at 18:18:44

Julia，標準で縦ベクトルなの数学と合ってて良き

#julialang #julia #Julia言語

タグ： julia julialang Julia言語

posted at 18:57:35

なお @nao_tuboyaki

天むす名古屋 Temmus @temmusu_n

Juliaと言う言語を学ぶには何からやればいいのか

タグ：

posted at 19:23:48

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

#超算数そして日本書籍の教科書にある「ひろ子のおてつだい」は、かなり露骨な掛け順こだわり強制指導を含むtwitter.com/temmusu_n/stat...。石田貞一と村田英吉はこの教科書と『算数の部』(1950)の両方に関与している。

タグ：超算数

posted at 21:14:17

非公開

タグ：

posted at xx:xx:xx