黒木玄 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 を見て下さい。
2020年10月05日(月)
@SGThr7 #Julia言語
using Plots
plot(sin, 0, π; label="sin(x)")
plot!([π/2, π/2], [sin(π/2), 0]; label="")
hline!([0]; label="", color=:black)
xtickstr = ["0", "π/4", "π/2", "3π/4", "π"]
s2f(s) = float(eval(Meta.parse(s)))
plot!(xtick=(s2f.(xtickstr), xtickstr))
gist.github.com/genkuroki/084f... pic.twitter.com/U26VZFEeo9
タグ: Julia言語
posted at 00:13:06
@SGThr7 #Julia言語 三角函数のプロットではtanのプロットがちょっと面倒で「NaNがプロットされない」という性質を使わないと添付画像1のように失敗します。
添付画像2は修正版。
gist.github.com/genkuroki/084f... pic.twitter.com/jv87AfeDEs
タグ: Julia言語
posted at 01:01:10
誘導法則を研究するにつれて無限の応用や問題がとり上げられることを諸君は見るだろう。電気機械の設計の研究はそれ自体としても生涯の仕事になる。その方向にあまり立入れないが、誘導法則を発見した時に、突然我々の理論が巨大な実用上の発展と結ばれたことをよく認識すべきである。
タグ:
posted at 01:31:41
@marksaroufim @tomkwong @staticfloat @MoseGiordano built BinaryBuilder which allows you to reliably and easily cross compile binaries for Windows, MacOS,... as well as targeting different architectures.
So instead of having to build things locally, we can just download and unpack a tarball.
タグ:
posted at 02:47:55
@marksaroufim @tomkwong @staticfloat @MoseGiordano BinaryBuilder.jl also takes care of these pesky little things like C++ and Fortran ABI's. It is serious magic!
タグ:
posted at 02:49:27
@vchuravy @tomkwong @staticfloat @MoseGiordano The benefits from the users perspective are so obvious. I wonder why pip or conda or brew never got around to doing this
タグ:
posted at 02:49:49
@ChadScherrer @tomkwong We only really need a nice command line toolkit that does this and exposes it to the rest of the system.
タグ:
posted at 02:50:25
@vchuravy @ChadScherrer @tomkwong can be done with github.com/Roger-luo/Como... easily probably
タグ:
posted at 02:59:37
私が「まるでshell scriptのように仕事をしている場合」を多用しているものと完全に同じものが、スーパーコンピューターの応用先の典型例である気候モデルを1から作り直すためのプログラミング言語として採用されているという事実は非常に面白い。
github.com/CliMA
タグ:
posted at 03:02:45
コンパイルのタイミングをうしろにずらしているお陰で、即時コンパイルの分のオーバーヘッドは生じるが、ライブラリAとライブラリBを合わせて使用した場合には、それぞれのコードが全体で最適化された形で常にコンパイルされることになる。
タグ:
posted at 03:22:21
実行直前にコンパイルする仕組みのお陰で、バイナリが対応しなければいけない函数の引数の型の組み合わせの指数函数的爆発の問題を一切気にすることなく、各分野の専門家が作成した型に基くパッケージどうしを滑らかに貼り合わせることができる。
タグ:
posted at 03:33:49
実行直前のコンパイルによるオーバーヘッドが気にならないような重たい計算をする場合には、実行直前にコンパイルする仕組みは計算速度の面でかなり有利になるものと思われます。
気候モデルのようなスーパーコンピューティングでJuliaが採用されることは不思議なことではありません。
タグ:
posted at 03:37:46
最近ツイッターでも話題になっていた「#R言語 においてGPUで加速された微分方程式のソルバを使用できるようになった」というニュースも #Julia言語 ネタ。
Juliaは多くのモノを貼り合わせてできているだけではなく、他のグルー言語による貼り合わせの対象も提供している。
www.stochasticlifestyle.com/gpu-accelerate...
タグ: Julia言語
posted at 03:44:24
Stefan Karpinski @StefanKarpinski
#JuliaLang’s BinaryBuilder is by far the easiest and most reliable way to install libraries in #macOS twitter.com/tomkwong/statu...
posted at 03:52:31
Patricia Schwarz, Ph @TessEractica
@phalpern And I testify to this from personal experience.
タグ:
posted at 06:33:38
Shozaburo Nakamura @shownakamura
1から始めるJuliaプログラミング、Juliaプログラミングクックブック
iiyu.asablo.jp/blog/2020/10/0...
タグ:
posted at 08:19:21
@NeXTSTEP2OSX もちろんです。サブ機の一つである Core m3 モデルです。Apple Silicon 版 Mac を買うのは見送ろうかなと思い始めてます(笑)
タグ:
posted at 08:53:31
加藤公一, 가토우 기미카즈(はむかず) @hamukazu
オラオラ高給で雇えと言ってる博士は(いるかもしれんが)ほぼいなくて、みんなが言ってるのはただ「仕事をください」ですね。 twitter.com/hiroki_f/statu...
タグ:
posted at 09:56:05
Birthday "paradox" states that in a room of 23, 50% chance 2 or more have same birthday
Making explainer video on some Julia programming technique by simulating it
Down to 78ms/million simulations
Possible to do better? Other langs welcome!
discourse.julialang.org/t/fastest-way-...
#julialang pic.twitter.com/RNy08xf3Ds
タグ: julialang
posted at 10:54:28
小学校の掛け算順序に関する問題って、順序を問わず正解にしている教師と、順序によって不正解にしている教師がいる、という状況なのだろうか。(とつぜんの超算数)
タグ:
posted at 11:47:19
回帰なんかをStanで実際に推定しようとすると、分散揃える謎処理した方がかなり速くなるんですよね。
NUTSのパラメータ調整が入るので、スケールが違いすぎると限られたwarmup期間では上手くいかない。
タグ:
posted at 12:03:54
須山敦志 Suyama Atsushi @sammy_suyama
あー、確かにこれはありますね。計算アルゴリズムが幾何的な情報をうまく扱えていない。 twitter.com/nan_bayesstat/...
タグ:
posted at 12:17:10
石原 俊/ISHIHARA Shun @ishihara_shun
この写真の男性がその後の人生で隣の女性と連れ添い、今も生きていて、しかも名乗りをあげたというニュースに、衝撃を受けている。
被爆翌年撮影の恋人は「私と妻」 広島の男性、写真出版で名乗り | 2020/10/5 - 共同通信 this.kiji.is/68566465115992...
タグ:
posted at 12:27:38
非公開
タグ:
posted at xx:xx:xx
Am I doing CUDA GPU properly?
It's _relatively_ easy to write a CUDA kernel in Julia for Bday paradox!
The performance can reach a scary 2ms on my RTX2080 but average out to abt 68ms which I think is gc time. But 2ms for 1million sims!
See discourse.julialang.org/t/fastest-way-...
#julialang twitter.com/evalparse/stat... pic.twitter.com/QlElLYOtLV
タグ: julialang
posted at 13:06:33
境界条件処理を追加したついでに、Juliaのソースコードも載せました。
水素原子に対するSchrödinger方程式を有限要素法で数値的に解いてみる(C++とJuliaのソースコード付き) qiita.com/dc1394/items/c... #Qiita
タグ: Qiita
posted at 13:48:27
#Julia言語 これは天文学におけるスーパーコンピューティング
juliacomputing.com/media/2017/09/...
juliacomputing.com/case-studies/c...
【Achieved peak performance of 1.54 petaflops using 1.3 million threads on 9,300 Knights Landing (KNL) nodes】 twitter.com/juliacomputing...
タグ: Julia言語
posted at 13:51:23
#Julia言語 これはすでに古くなっているCUDA.jl v1.0の情報
現在のCUDA.jl は v2.0.0
ものすごいスピードで開発が進んでいる印象あり
Julia+CUDAで気軽にGPU計算 qiita.com/Shoichiro-Tsut... #Qiita
posted at 13:58:04
#Julia言語 CUDA.jl v2.0 情報
juliagpu.org/2020-10-02-cud...
新しい技術を取り込んで、ユーザーにとって使い易くしてくれるのはとてもありがたい。 twitter.com/juliabloggers/...
タグ: Julia言語
posted at 14:03:33
加藤公一, 가토우 기미카즈(はむかず) @hamukazu
lets-get-arrestedにリトアニア語翻訳のプルリクエストが来ていた。ありがたい。
github.com/hamukazu/lets-...
タグ:
posted at 14:25:47
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
#Julialang 最小二乗法な回帰分析の対象としたかった関数が2種類:本質的なところは乗り越えたと思う。あとは単純作業の繰り返しでいけるはず。ご助言感謝!#JuliaNLSolvers
posted at 16:57:50
絶望して死を選ぶ前の最後の望みで来る人も居るわけで、この対応は殺人と変わらない
//…仕事はコロナで休業になってしまって、うつっぽい状態でもう働けません』…『一人暮らしです』って答えたら、『地元はどこなの?実家に帰れないの?』news.livedoor.com/article/detail...
タグ:
posted at 18:51:51
イギリスで大量の感染者の報告漏れ。検査結果をまとめる EXCEL の行数が限界に達していたのに1週間気づかなかったのが原因。…気を付けよう!← twitter.com/MaxCRoser/stat...
タグ:
posted at 19:04:45
短時間のもくもく会やります。今回は開発に限らずjuliaに関するワークならなんでもOKって感じにしてます。
#julia言語 twitter.com/takuizum/statu...
タグ: julia言語
posted at 19:22:17
Instead of "one" problem, they now have "several".
#rstats #julialang to the rescue! twitter.com/MaxCRoser/stat...
posted at 19:40:44
#Julia言語 Re:RTs
【We install all the things for you】(笑)
github.com/JuliaPackaging...
Juliaのパッケージで使う各プラットフォームでのバイナリを供給するための仕組み。ユーザー側はこれのお陰で滅茶苦茶楽をできる。
例 github.com/JuliaBinaryWra...
例 github.com/JuliaBinaryWra...
タグ: Julia言語
posted at 20:50:14
もじゃもじゃのほう (Koki Tsuy @antiplastics
Rの
library("パッケージ名")
は、pythonで言うところの
from パッケージ名 import *
みたいなものだと思うが、別にlibraryしないでも
パッケージ名::関数名
で使えてしまうので、だったらずっと::だけで書いててもいいのか?(よく知らない
タグ:
posted at 21:11:54
#Julia言語 公式バイナリのインストーラーを
julialang.org/downloads/
からダウンロードして、インストールすれば、数学なんかの勉強にすぐに利用できるようになります。
julia> ]
pkg> add Plots
pkg> precompile
pkg> バックスペース
jukia> using Plots
julia> plot(sin, -pi, pi)
とか。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 21:48:33
ごまふあざらし(GomahuAzaras @MathSorcerer
3乗したら1が連続して2020個並んでいるような数はあるか?って言われたんだけれど,頭働かないのでとりあえず腕力で実験.
#Julia言語 pic.twitter.com/J1a0yvoOPO
タグ: Julia言語
posted at 21:53:18
はてなブログに投稿しました #はてなブログ
水素原子の基底状態解は原点での境界条件を陽に考慮せずに得られる話 - わかばめにっき wakabame.hatenablog.com/entry/2020/10/...
タグ: はてなブログ
posted at 21:55:37
ごまふあざらし(GomahuAzaras @MathSorcerer
@hirako_zonepre キュッキュ😃(
こちらにインストールの手順がありますので
お時間がある時に動かしてみてください(上のスクリーンショットは Jupyter で動かしているところです
julialang.org/downloads/plat...
)
タグ:
posted at 22:00:43
#Julia言語
本当の話をすると、Juliaでは最小二乗法の計算を本質的に
b = X\y
の1行で可能です。Juliaは線形代数と相性がよく、最小二乗法は線形代数で理解できます。
sin(x)+ノイズを3次式で近似するのもこんな風に簡単。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat... pic.twitter.com/0mu0P4t9zw
タグ: Julia言語
posted at 22:10:00
非公開
タグ:
posted at xx:xx:xx
@MathSorcerer @hirako_zonepre #Julia言語 【Pythonそっくり】と言われたので紹介。
添付画像1はPythonでmatplotlibを使っている様子。
添付画像2はJuliaでmatplotlibを使っている様子。
Pythonのライブラリの多くが容易にJuliaでも使えます。
PythonのおかげでJuliaも便利になっている。
nbviewer.jupyter.org/gist/genkuroki... twitter.com/genkuroki/stat... pic.twitter.com/Ef5qJ4c9rs
タグ: Julia言語
posted at 22:25:42
https://mathstodon.x @abel_siqueira
As aulas de hoje de Otimização em @JuliaLanguage começam a tratar dos problemas com restrições.
- youtu.be/NH8mufxxZNw
- youtu.be/l-5HUcILnjA
- youtu.be/71W-meyPak8
#otimização #orms #JuliaLang pic.twitter.com/BOkJJzZa74
posted at 22:38:06
#Julia言語
Juliaで配列Aの各要素ごとに函数 exp(-x^2)*x を作用させるには
exp.(.-A.^2).*A
と演算子や函数に . をつければよいのですが、. の付け方が足りなくなるバグによく悩まされます。
f(x) = exp(-x^2)*x
f.(A)
と書いた方が失敗が減らせる。 twitter.com/genkuroki/stat... pic.twitter.com/Lj3giJuKVI
タグ: Julia言語
posted at 22:43:08
非公開
タグ:
posted at xx:xx:xx
#Julia言語 パラメータ a = [a[1], a[2]] を何らかのパッケージを利用して最適化する場合には、a[1], a[2] が動ける範囲に合わせて変数変換を行った方が無難です。
例えば、a[1]の部分が正ならば、a = [exp(b[1]), b[2]] の形でaを使い、パラメータ b = [b[1], b[2]] を最適化するようにする。 pic.twitter.com/E1loxIW4Hw
タグ: Julia言語
posted at 22:47:03
しかし、MCMC法のパッケージを徹底的に使いこなそうとするときには、内部で何をやっているかに関する予備知識が必要になるので、以上の話は知っておいて損がありません。
タグ:
posted at 22:52:04
@hoxo_m @antiplastics 横入ですが,私の記憶だとkfasで両社の挙動が変わったことがあったような記憶があります.詳細は覚えてないのでアレですがまあまあ有名なパッケージでも起きるという
タグ:
posted at 22:53:08
@MathSorcerer こんなんでどうでしょう?
#Julia言語
gist.github.com/antimon2/19fb6... pic.twitter.com/1jwjrEjJWw
タグ: Julia言語
posted at 23:01:18
#Julia言語
配列A内のx以下の成分の割合は
ecdf(A, x) = count(a ≤ x for a in A)/length(A)
で求まるのですが、
xs = range(-5, 5, length=1001)
の各成分に作用させるために
ecdf.(A, xs)
と書くと、配列Aにもブロードキャストが効いて失敗します。代わりに
ecdf.(Ref(A), xs)
と書く! pic.twitter.com/k1bhx8gOKS
タグ: Julia言語
posted at 23:15:51
#Julia言語
?Ref
でヘルプを読むと添付画像のようにRefはブロードキャスト内でスカラーのように解釈されることも書いてあります。 pic.twitter.com/TrifdMPtMN
タグ: Julia言語
posted at 23:19:55
こういうのも提供してくれる人がいました。
版が違ってもこれなのだそうです。
本気で信じているとしか思えません。。。 twitter.com/f_sei/status/1... pic.twitter.com/E7TEAzBG3L
タグ:
posted at 23:46:03
笑い事ではないけれど「3つのベクトルのどの2つも平行でなければ1次独立」って笑っちゃうな。
世の中にはデタラメの書いてある参考書や問題集も少なくないけど,何とか規制できないのかな。
タグ:
posted at 23:54:01