黒木玄 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年05月30日(日)
#Julia言語 では、a.f(x,y) スタイルを積極的に排除した方が良い理由があると思う。
例えば、パッケージAで型Sが定義されていたとします。
私は、パッケージAの作者に無断で、自作のパッケージBで函数 f(s::A.S, x, y) を定義して便利に使用できます。続く
discourse.julialang.org/t/allowing-the...
タグ: Julia言語
posted at 23:56:19
#Julia言語 を使えば10行でガウス過程回帰のアニメーションを作れちゃう。
これもシンプルな線形代数。
ただし、線形写像としての行列ではなく、二次形式としての行列の線形代数。
行列を線形写像と決め付ける感覚になった人は悪しき過学習に陥っているので要注意。そう単純なものではない。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 22:55:04
@mat_der_D 添付画像を見れば誤解せずに済むのですが、念のために補足。
上で【1.2~1.3倍程度速くなりました】と述べているのは、Juliaのコードどうしの比較です。もっと速くできるかどうかは知りません。
タグ:
posted at 22:45:09
@mat_der_D #Julia言語 Rust版はu32で、Julia側ではすべてInt64で計算しているので、そのままだとRust版とJulia版の比較は無意味だと思いました。
個人的な意見では、同じLLVM仲間でここまで単純な計算で速さを比較しても無意味だと思いました。双方頑張って最速にするとおそらく同じ程度の速さになる。
タグ: Julia言語
posted at 22:38:51
@mat_der_D #Julia言語 Rustのコードのis_primeをJuliaに翻訳したら、私が試した範囲内で1.2~1.3倍程度速くなりました。
module AのコードはオリジナルのJuliaのコードで、Bのis_primeはRust版の翻訳で、CではBase.srem_intを使っています。
これ実質的に剰余計算の速さの計測です😅
gist.github.com/genkuroki/5cb3... pic.twitter.com/ZfXyNc0Q2i
タグ: Julia言語
posted at 22:33:37
ごまふあざらし(GomahuAzaras @MathSorcerer
#Julia言語
以前は htmlpreview.github.io で共有できていたがそれができなくなってた漁ったら下記のIssueを見つけたからです.
github.com/fonsp/Pluto.jl...
タグ: Julia言語
posted at 22:10:20
ごまふあざらし(GomahuAzaras @MathSorcerer
使い方:
1. Pluto.jl でノートブックを適当に書いて HTML に保存する(Plutoノートブックの右上のアイコンからHTMLに変換するボタンがある.)
2. それを gist などにアップする.
htmlview.glitch.me
に行って HTML がある場所URLを指定する.
それだけ.
#Julia言語
続く
タグ: Julia言語
posted at 22:10:20
ごまふあざらし(GomahuAzaras @MathSorcerer
#Julia言語
Pluto.jl から変換したHTMLを他の人にも見てもらうためのサンプルを作りました.
gist.github.com/terasakisatosh...
続く
タグ: Julia言語
posted at 22:10:20
非公開
タグ:
posted at xx:xx:xx
OokuboTact 大久保中二病中年 @OokuboTact
文科省の「指導要領」と「指導要領・解説」の法的強制力の違いを知らない教師が多いというのはヤバい。
「指導要領」も昔は法的強制力が無かったはずだけど。
数学(算数を含む)教育で「法的強制力」を気にするのは変だと思う。
文科省が「1+2=5」と命令したら従うのか?
#超算数
タグ: 超算数
posted at 21:46:32
OokuboTact 大久保中二病中年 @OokuboTact
#超算数
今の算数教育は社会主義国家のマルクス主義理論みたいなもんだな。
大学で文科省の「指導要領・解説」をテキストにしているところ多いけど、批判せずに教条的に教えるんだろうな。
タグ: 超算数
posted at 21:42:48
@zobomi #Julia言語 タプルのベクトル ys = g.(xs) を添付画像の方法で y に変換して、plot(xs, y[1]) とかすればよいと思います。
t[i] = getindex(t, i)
なので、
getindex.(ys, i)
は
(t -> t[i]).(ys)
とも書けます。タプルの第i成分を集めた配列を作れます。 pic.twitter.com/TfZcy6CtLA
タグ: Julia言語
posted at 21:39:01
非公開
タグ:
posted at xx:xx:xx
今朝、新しい記事を書いています。文字コードについての説明がメインです
muuumin.net/julia-intro-pr... twitter.com/muuumin20/stat...
タグ:
posted at 20:17:31
Σ(1/n) が30になるまでどのくらいかかるかの計算.Juliaで実行中だが,2時間半経っても終わらないw スピードアップのため,途中経過も何も表示せず,とにかく30を超えたら終了というループを走らせている.
タグ:
posted at 20:12:47
#Julia言語 Juliaで10行シリーズを修正更新した。
nbviewer.jupyter.org/gist/genkuroki...
たとえば
⪅(x, y) = x < y || x ≈ y
を
x ⪅ y = x < y || x ≈ y
に書き直した。
タグ: Julia言語
posted at 19:00:45
スウェーデンの人たちは、誰もこんなにたくさん人が死ぬとは想像もしてなかったんだよね。今は政府に対する不信感から政権交代も起きそうなほど。
そして不思議なのは、日本のマスコミがスウェーデンの実情を報道していないこと。 twitter.com/siamese1125/st...
タグ:
posted at 18:59:47
#Julia言語 自画自賛していた。
こも「Juliaで10行シリーズ」を結構気に入っている。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 18:08:39
#Julia言語 Juliaならたったの10行で以下が全部可能!
1. 信頼区間が二項分布モデルのパラメータ値を含む確率
2. 同ケースでの信頼区間の分布
3. 以上のFisher検定版
4. 3次元ランダムウォーク
5. Ising2D.jl
6. 蔵本モデル
7. ジュリア集合
8. マンデルブロ集合
9. ライフゲーム
10. 2D Ising model twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 18:05:15
#Julia言語 私のツイログを
twilog.org/genkuroki/sear...
P値 信頼区間 実装 Julia
で検索すると、P値や信頼区間を2~3行で実装する方法が見つかります。
nbviewer.jupyter.org/gist/genkuroki...
のpvalとciがそれぞれP値と信頼区間を求める函数です。 twitter.com/mathworld4/sta... pic.twitter.com/bQoBT23BkQ
タグ: Julia言語
posted at 18:00:27
大失敗です。
新型コロナの感染による集団免疫を目指した結果、世界でも人口比で日本の十倍以上の被害を出しました。
また、変異ウイルスには既感染者も普通に感染することが分かり、国王まで「失敗だった」と異例のコメントを出しました。
また経済的ダメージも極めて大きかったです。 twitter.com/siamese1125/st...
タグ:
posted at 18:00:01
#Julia言語
nbviewer.jupyter.org/gist/genkuroki...
波動方程式を質点がバネで繋がった系で近似しているのですが、自由境界条件と固定境界条件のコードの違いは、添付画像中の赤線で囲んだ部分しかありません。
どうしてこれだけの違いでよいかがノータイムで分かることはとても基本的。イメージすれば終わり。 pic.twitter.com/Vj1avC8gO6
タグ: Julia言語
posted at 17:50:12
#Julia言語 修正後のソースコード
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/gist/genkuroki...
前者にはついでに固定境界条件(ディリクレ境界条件)の場合も追加した。境界での反射するとき赤と青が反転する。 pic.twitter.com/Ss1zYWH6Sx
タグ: Julia言語
posted at 17:37:33
#Julia言語 明らかにおかしいのは添付画像の部分。赤枠の1つ上の行をコピー&ペーストしたまま、必要な変更を施さずに放置してしまったミス。
直すのは簡単。
これで、失敗例の証拠を明瞭に残せたので、gistの方のコードも後で直しておきます。
失敗の例が公開されていることはとても大事。 pic.twitter.com/u3Cyhcjm69
タグ: Julia言語
posted at 17:20:47
スウェーデン、調べたら普通に子供亡くなってるじゃん。
この亡くなった子供たちって、本当に全員、基礎疾患あったの? twitter.com/MIKITO_777/sta... pic.twitter.com/GLoN8ODdpW
タグ:
posted at 17:07:46
ほりたみゅ (mastodonも見てね) @Hyrodium
@genkuroki 境界部分のdv[end,1], dv[1,end], dv[end,end]を編集してそれっぽい出力になりました! pic.twitter.com/pYtQGwjwKm
タグ:
posted at 16:34:18
#Julia言語 微妙に失敗しているっぽいが、面白いのでそのままにしておきます。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 16:16:37
ほりたみゅ (mastodonも見てね) @Hyrodium
t=0.75付近で、右上と左上から発生している波は何なのでしょう twitter.com/genkuroki/stat... pic.twitter.com/GCqslx7TP5
タグ:
posted at 16:07:35
ごまふあざらし(GomahuAzaras @MathSorcerer
workaround
```
x(t) = cos(t)
y(t) = sin(t)
ẋ(t) = -sin(t)
ẏ(t) = cos(t)
p(t) = [x(t), y(t)]
# ṗ(t) = jacobian(p, t)[begin]
ṗ(t) = [ẋ(t), ẏ(t)]
``` pic.twitter.com/5pOPwgDej0
タグ:
posted at 16:01:12
一般人へのインターネットの世界的普及によって、世界中の「かわいい」を紹介してもらえるようになった。
科学技術は「かわいい」の楽しみを増やすために使われるべき。
タグ:
posted at 16:00:11
#Julia言語
ノイマン境界条件(固定境界条件)では、外が赤(正)で内が青(負)の波がぶつかると、波の正負が反転して、外が青で内が赤の波になって跳ね返って来る。 twitter.com/genkuroki/stat... pic.twitter.com/WQlFZOcOuZ
タグ: Julia言語
posted at 15:53:29
ごまふあざらし(GomahuAzaras @MathSorcerer
#JuliaLang
#Julia言語
Does not work... :(
```julia
using QuadGK, Zygote
x(t) = cos(t)
y(t) = sin(t)
p(t) = [x(t), y(t)]
ṗ(t) = jacobian(p, t)[begin]
s(t) = quadgk(
t -> sqrt(ṗ(t)[begin]^2 + ṗ(t)[end]^2)
, 0, t
)[begin]
s'(3.14) # :sob: :sob: :sob:
``` pic.twitter.com/lGrXPk4IDH
posted at 15:47:03
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/gist/genkuroki...
コードはちょっと複雑になっていますが、やっていることは単純で、質点がバネで繋がった系で波動方程式を近似している。境界条件は境界でのバネの繋がり方の違いで表現されている。
任意の領域を扱えます。
タグ: Julia言語
posted at 15:38:21
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
ノイマン境界条件
びよんびよん😊 pic.twitter.com/SIidCvgBYo
タグ: Julia言語
posted at 15:35:32
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
ノイマン境界条件 pic.twitter.com/6vujsRBhbN
タグ: Julia言語
posted at 15:33:56
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
ノイマン境界条件 pic.twitter.com/suaeAAtLWd
タグ: Julia言語
posted at 15:32:35
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
自由境界条件 pic.twitter.com/nHPejLpqXC
タグ: Julia言語
posted at 15:30:55
#Julia言語 正五角形上での波動方程式
nbviewer.jupyter.org/gist/genkuroki...
自由境界条件 pic.twitter.com/n7izrD707Y
タグ: Julia言語
posted at 15:29:13
#Julia言語 自由境界条件であることはこちらの動画の方が見易い。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/XmIKKgbdZb
タグ: Julia言語
posted at 15:26:28
#Julia言語 上の動画とこのツイートの動画は自由境界条件。正方形の境界が自由に動ける。 pic.twitter.com/Nz9VRdiWpv
タグ: Julia言語
posted at 15:25:04
#Julia言語 私が書いた波動方程式のコード↓
nbviewer.jupyter.org/gist/genkuroki...
DifferentialEquations.jlを使って楽をしている。 twitter.com/HmaimesAlice/s... pic.twitter.com/epfoxoTboV
タグ: Julia言語
posted at 15:22:21
vscode上でjupyter notebook使ってjuliaしてたらエラーが出てなぜかセルの追加ができなくて(ファイルを読み込みなおすとセルが追加された状態で開かれる)、ググった情報をもとにvscode insidersで試してみたらエラーが出なくなった
タグ:
posted at 13:55:07
前に貼ったヒゲのおっちゃんは本垢復活したらしい。やっぱ中国系ではこの人の料理が一番うまそう
vt.tiktok.com/ZSJmE6tPA/ pic.twitter.com/iLZP1y4Pa8
タグ:
posted at 13:53:06
非公開
タグ:
posted at xx:xx:xx
@physics303 ついでに
pkg> add PyCall PyPlot SymPy
julia> using PyPlot
julia> using SymPy
として、~/.julia/conda 以下にmatplotlibやsympyも入れてもらった方がよいと思います。Pythonさん御一行はJuliaのライブラリとしても非常に便利です。
タグ:
posted at 12:53:53
速い!!!!
#Peing #質問箱 peing.net/ja/qs/1006562650
posted at 12:52:25
#Julia言語 私と同じインストールの仕方を再現したい人は
nbviewer.jupyter.org/github/genkuro...
を見て下さい。ちょっと時代遅れになっている部分がありますが、現在もほぼ有効な解説になっています。
'zmq.h':No such file or directory 問題に出会ったら、以下のリンク先を参照。
twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 12:30:34
私はJulia用のPythonさん御一行のインストールはJuliaさんに全部まかせています。まっさらな状態で
julia> ]
pkg> add IJulia
pkg> バックスペース
julia> using IJulia
julia> notebook()
とすればJuliaさんがJupyter関連のインストールを自動的に行ってくれます。
運が良ければこれだけで使える。 twitter.com/physics303/sta...
タグ:
posted at 12:27:29
ごまふあざらし(GomahuAzaras @MathSorcerer
それすらめんどくさい場合は Docker と winget を導入してMyWorkflow.jl で整えた環境をカスタムしてもろて。。。
初めてJulia触る人は
有名な
1から学ぶJulia
か
特集のInterface2020年6月号を買ってください。(電子版あるらしいです)
タグ:
posted at 12:14:49
ごまふあざらし(GomahuAzaras @MathSorcerer
既存のPythonの環境と連動したい場合はENVにPythonとJuliaのパスを突っ込んで IJulia PyCall をadd build してリンクさせる。
何を言ってるのかわからない場合はPyCall 側でよしなに新しいMiniconda環境がHOME/.julia/conda に入るのでそれに身を任せる
続く
#Julia言語 twitter.com/physics303/sta...
タグ: Julia言語
posted at 12:14:48
#Julia言語 関連情報
github.com/fonsp/Pluto.jl...
Cannot use Revise.jl's includet() in Pluto #1101
タグ: Julia言語
posted at 11:56:19
#Julia言語 Pluto.jlでRevise.jlを使用できれば問題解決だと思うのですが、その方法がさっぱり分からない。
github.com/fonsp/Pluto.jl...
Revise support! #416
なにをどうサポートしたのか、理解できませんでした。
理解していてかつ実際に動く例を提供できる人がいたら教えて下さい。 twitter.com/glakneshoe/sta...
タグ: Julia言語
posted at 11:49:00
#Julia言語 あと、Juliaは数学記号や絵文字のユニコード文字の入力方法もサポートしており、Julia対応環境であれば、変数名などに σ², x₁, x₂, Δt などを自由に使えるので、アルゴリズム説明用の擬似コードを使う代わりに、Juliaで実際に動くコードを使えます。
bicycle1885.hatenablog.com/entry/2014/12/...
タグ: Julia言語
posted at 10:35:23
#Julia言語
sciml.ai/news/2018/04/3...
【It has been mentioned that "every great open source math library is built on the ashes of someone's academic career". What we needed was a new way to connect software development ...】
これへの答えの1つがJulia言語のcomposability.
タグ: Julia言語
posted at 10:35:23
#Julia言語 2018年8月にv1.0になった瞬間にヘビサイド函数のようにGitHubにおけるJuliaのスターが増えました。v1.xの間はbreaking changeを行わないと約束した点が大きかったのだと思います。
Juliaでは所謂expression problemが解消していている点も大きいと思う。
www.juliabloggers.com/newsletter-jan... twitter.com/tanki_nobu443/... pic.twitter.com/rBjPF3ekEc
タグ: Julia言語
posted at 10:35:20
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
#JuliaLang
julia> a … b = a > b ? (a:-1:b) : a:b
… (generic function with 1 method)
julia> a = collect(4 … 1)
4-element Vector{Int64}:
4
3
2
1
julia> a = collect(1:4); a[end…begin]
4-element Vector{Int64}:
4
3
2
1 pic.twitter.com/ej2zdMKS4u
タグ: JuliaLang
posted at 02:04:10
非公開
タグ:
posted at xx:xx:xx
もしかしたらこの伊能によるデータは、1804年以前の鳥海山の標高の参考データになるのかもしれない。もう少し詳しいデータを伊能が残していてくれたらとは思うが、仕方ない。 参考 www.data.jma.go.jp/svd/vois/data/...
タグ:
posted at 00:34:11
1800~1804年の鳥海山の噴火活動で、現在の新山付近に溶岩ドームが形成。1804年の地震では象潟が隆起して消滅。伊能の測量隊も、この時期に出羽沿岸を通過しているので、もし地震や噴火に巻き込まれていたとしたら、あの偉業は成し遂げられなかったということにもなるわけか。
タグ:
posted at 00:34:10
伊能忠敬による測量データと現在のデータの照合をしていたら、どうしても鳥海山の標高が合わない。鳥海山から18kmの近距離で測量していたはずなのになぜか角度が1度近くも狂っている。確認したら、伊能が象潟付近を通過した前後、鳥海山は活発に噴火活動していたのでそのせいだと予想が付いた。
タグ:
posted at 00:34:10