Thread by @genkuroki

#Julia言語 値の集合から確率を指定してサンプルを1つ抽出する方法

using StatsBase
a = sample([2, 3, 5], ProbabilityWeights([0.2, 0.3, 0.5]))

using Distributions
a = rand(DiscreteNonParametric([2, 3, 5], [0.2, 0.3, 0.5]))

複数取り出す場合は添付画像。

github.com/genkuroki/publ…

https://twitter.com/qman63933857/status/1438687818016980993

#Julia言語 Juliaでappend!は配列に配列を追加するために使うことが普通で、配列に要素を追加するために普通使うのはpush!です。REPLで

julia> ?
help?> append!

julia> ?
help?> push!

を確認してみてください。REPLに ? で質問すると沢山のことが分かります。続く

https://twitter.com/qman63933857/status/1438755234034057221

#Julia言語 Pythonから来た人達がよくやっている失敗は

a = []
for i in 1:n
xを計算
push!(a, x)
end

というコードで長さnの1次元配列を作ること。

このaのeltype(a)はAnyになるので、aを使う計算が数十倍遅くなる。

具体的に他のどのような書き方が適切であるかはケースバイケース。続く

#Julia言語 多くの人がよくやってしまう失敗は、計算速度向上を目指しているのに、グローバル変数を函数の引数を経由せずに使ってしまうことです。計算が数十倍遅くなります。

函数として書いていないコードを実行すると大抵そうなります。(let ～ end で囲んでいる場合は例外的)

#Julia言語 Fortranから来た人は主にグローバル変数の問題を気にしていれば何とかなり易いという印象あり。もともと数値計算に慣れている点も大きい。

しかし、Pythonから来た人は、Pythonが特殊で風変わりであることに気付いていないと、Juliaで(実際には他でも)失敗するパターンが非常に多い。

#Julia言語 特にPythonから来た人は

docs.julialang.org/en/v1/manual/p…
Performance Tips

を最初に読んでおくと、計算を遅くする多くの失敗を事前に防ぐことができて、時間の無駄を省けると思う。

私自身はこの文書の存在を気付いたとき「最初にこれの存在を教えてくれよ😭」と泣きが入りました。

多分、需要が大きいのは、

　Pythonから #Julia言語 への翻訳例

で

docs.julialang.org/en/v1/manual/p…
Performance Tips

の要求をすべて守っているもので注意するべきことの解説付きのものの集まり。

クラスを使っている部分の翻訳の仕方も含まれる。

誰か頑張って書くと良いと思う。かなり大変そう！

「Pythonから #Julia言語 への翻訳で注意するべきこと」はものすごく多いので、個人的には普通の人間には大変すぎる仕事になるような気がします。

あと、ユーザー数が巨大でお客さんの立場で論評し易いPythonの世界と、全然そうではないJuliaの世界の違いを最初に認識することも大事だと思う。

#Julia言語 っぽい書き方については

v1.6以下
github.com/genkuroki/publ…

v1.7以上
github.com/genkuroki/publ…

を参照。

扱っている題材は一様重力下での質点の投げ上げなので中学校レベル。地球と月の場合を計算している。

#Julia言語 できるだけFloat64とかInt64のような具体的な型名を書かない方針になっていることにも注目。

Float64とかInt64のような具体的な型名を書かない代わりに、Problemタイプの定義で型パラメータG, Y0, V0, TSを使っている。

v1.7以上でのコード
github.com/genkuroki/publ…