黒木玄 Gen Kuroki(@genkuroki)/2021年04月16日

まとめました: julia言語でグラフを描いたりして遊びたい togetter.com/li/1698906

タグ：

posted at 23:52:38

質問者2 @shinchanchi

「FTPL風の◯経」とお聞きしたのが印象に残ってます。

要するに魔法使いの弟子か出汁

タグ：

posted at 22:56:03

質問者2 @shinchanchi

ここまで曲解してリフレ派を語れるのは、ある意味素晴らしい💦

残念な時間泥棒がネットにはいるかもしれませんのでご注意を^_^

MMTに抱きつくリフレ派｜Prof. Nemuro🏶 #note note.com/prof_nemuro/n/...

タグ： note

posted at 22:54:48

#Julia言語

* パッケージ作成の手間がほぼゼロであること

とか

* パッケージの開発にはRevise.jlが便利である

とか

* 遅延を減らすためには具体的な値の具体型の伝搬について理解しておいた方がよい

とか、この手のノウハウの普及がJuliaについてまだ十分でない。繰り返しの解説が必要。

タグ： Julia言語

posted at 22:50:55

#Julia言語続き

上書きすると、solve(prob)の実行時に再度遅延が生じます。初回よりずっと小さな遅延ですむのですが、毎回経験するには結構不快な程度の遅延が生じます。

時間発展を記述する函数fの上書きを繰り返さないように注意した方が得です。

タグ： Julia言語

posted at 22:50:54

#Julia言語続き

prob = ODEProblem(f, u0, tspan, p)

の具体型は、当たり前の話ですが、f, u0, tspan, pの型でけっていされます。

注意するべきことは、時間発展を記述する函数fの定義を(たとえ元の定義と完全に同じであったとしても)上書きしてしまわないことです。続く

タグ： Julia言語

posted at 22:50:54

#Julia言語続き

2回目以降のsolve(prob)の実行で問題probの具体型を変えなければ、すでにコンパイルは終了しているので、遅延は生じません。

だから、遅延が生じる場合を減らすには、probの具体型がどのように決まるかを知っておいた方がよい。

Juliaでは具体物も具体型がどうなるかが常に重要。

タグ： Julia言語

posted at 22:50:53

#Julia言語続き～、微分方程式の問題probを作成し、solve函数で問題probを数値的に解くことです。

このとき、solve函数は問題probの具体型に最適化されたネイティブコードのコンパイルされて実行されます。初回の実行で遅延が生じる。

続く

タグ： Julia言語

posted at 22:50:53

#Julia言語遅延対策のノウハウ追加

最近、DifferentialEquations.jl の優秀さに気付く人がかなり増えたという印象がある。典型的な使い方は

prob = ODEProblem(f, u0, tspan, p)
sol = solve(prob)

と時間発展を記述する函数fと初期値u0と時刻の幅tspanとパラメータpを与えて～続く

タグ： Julia言語

posted at 22:50:52

AＫ_TSUWO @AwoK_pleasure

超のつくほどのやっちゃダメな問題ですね。 twitter.com/mkstudying1/st...

タグ：

posted at 21:44:08

Julia Bloggers @juliabloggers

New post: Working with matrices in DataFrames.jl 1.0 - www.juliabloggers.com/working-with-m... #julialang

タグ： julialang

posted at 20:55:20

Gabriel D Weymouth @gabrielweymouth

Andhini casually rolls into the group meeting today and shows us this awesome simulation. She doesn't have a twitter, but see her github for other cool content. andhini.github.io/Flapping-foils/ pic.twitter.com/uSOULLKCBV

タグ：

posted at 20:50:35

Roland Krause*** @rrrkrause

Bun-chan(ぶんちゃん) @BunchanBamboo

Learning about a new programming language offers new perspectives on the one you’re using as your daily workhorse.

3h #julialang is a good investment for R and Python analysts and programmers. twitter.com/elixir_lu/stat...

タグ： julialang

posted at 20:45:00

Bun-chan(ぶんちゃん) @BunchanBamboo

劇場版ポケットモンスター
C++ vs julia vs fortran

タグ：

posted at 20:35:33

「今日はjuliaと一緒に心中するんだ」って言ったらmatlabユーザーの知り合いに「外人の彼女でもできたん？」って言われた

傷ついたからmatlabを煽るブログ書きたい

タグ：

posted at 20:30:52

Aww Cats @AwwCats_

pic.twitter.com/0G5UxKdbDg

タグ：

posted at 20:17:33

mathematica_command @mathematica_bot

FiniteGroupCount[n] : 位数n の有限群の数を与える

タグ：

posted at 19:35:58

@genkuroki 確かに知らなかったら躓いたと思います。ありがとうございます。

タグ：

posted at 18:48:50

萩生九 @theseaofbarley

英検もTOEFLも高額だが個人申込みができる一方、「GTEC検定版」は教員のタダ働きが文化として容認されている一部の高校の生徒しか受験できない。
こういう方式の検定を入試で使うことははっきり言って「差別的入試方式」である。関係大学の職員や学生もおかしいことに気づき内部から声を上げるべきだ。

タグ：

posted at 18:44:30

萩生九 @theseaofbarley

TOEFLやTOEICを受ける前の英語民間試験の模試扱いだったGTECを大学入試で使うと聞いた時の衝撃は大きかった。そしてあれだけ問題が明るみになった今でもこれを入試として用いている立教大学他には怒りを通り越して失笑を禁じ得ない。 twitter.com/hirotourasaki/...

タグ：

posted at 18:16:53

michio inoue l MATLA @michio_MWJ

@keisemi 【蛇足】もちろん MATLAB は MATLAB の上で使うのが一番機能が活きます。ライブスクリプトもどうぞ👍 pic.twitter.com/ZbpgjiFGV5

タグ：

posted at 18:09:00

#超算数

数学的記号法に公式ルールはありません。
数学的記号法の最高裁判所は存在しない。

8÷2(2+2)の解釈として、(8÷2)×(2+2)と8÷(2×(2+2))のどちらが正しいかは決まっていません。

こういう解答が

math.stackexchange.com/questions/3321...

で一番人気です。 twitter.com/mkstudying1/st... pic.twitter.com/8X0uUhSVnq

タグ：超算数

posted at 18:02:29

#Jupyter のipynbファイルの良い点はGitHubでサポートしてくれていること。もしもGitHubでサポートしてくれていなかったら、やはり変換したものもGitHubに置くことになるだろう。

nbviewer.jupyter.org のようにipynbファイルがどこにあっても閲覧できるサービスも便利。 twitter.com/michio_mwj/sta...

タグ： Jupyter

posted at 17:49:44

#Jupyter のように【データと使ったコード、結果とメモ全部が一つにまとまった解析ノートの様に使いたくて、必要に応じてpdfなどに出力出来る】ことへの強い需要はあるよね。

ときどき、Jupyterの類をイデオロギーとして嫌いな人がいて閉口する。 twitter.com/keiyotokei/sta...

タグ： Jupyter

posted at 17:45:48

@tune3orz 具体的には Julia v1.0 には linspace(a, b, n) がなくなっていて、代わりに range(a, b; length=n) を使うようになっており、そこで躓くと思いました。

v0.6からpre-v1.0のv0.7で仕様を大幅に変えて、2018年8月のv1.0以降はbreaking changeは無し(厳しい縛り)という約束で非常に安定しています。

タグ：

posted at 17:20:34

nope @definiteIyNOTok

Analytics India Maga @Analyticsindiam

Artwork that messes with people's reality pic.twitter.com/X0cicluXvd

タグ：

posted at 16:29:08

Julia 1.6 Released: Major Updates & Latest Features
@JuliaComputing

#BigData #Analytics #DataScience #AI #IoT #IIoT #Python #RStats #JavaScript #ReactJS #CloudComputing #Serverless #Linux #Programming #Coding #100DaysofCode #julialang

bit.ly/3aht15f

タグ： 100DaysofCode AI Analytics BigData CloudComputing Coding DataScience IIoT IoT JavaScript julialang Linux Programming Python ReactJS RStats Serverless

posted at 16:12:30

@cloez_uya

juliaってadd Arpackで疎行列対角化できるんだ

タグ：

posted at 13:25:43

まだ疲労困憊から回復できていない。

タグ：

posted at 12:40:34

山田ゴンザレス花子 @ygonzales875

#Julia言語 using CSV や using Plots などの数秒の遅延の問題を実質的になくす他の方法に、Jupyterの最初のセルで using Plots; plot(sin) を実行して、実行が終わる前に、その次のセルにコードを書き始める、というような方法もあります。

タグ： Julia言語

posted at 12:16:42

MATLAB君がライセンス通してくれないからjuliaで計算してたんだけど、「MATLABの出力結果で書け」って問題文に書かれてて意味わからんくなった

タグ：

posted at 12:10:59

#Julia言語

「パッケージを開発しなければいけないのかよ」と思ってしまいがちで、私もそうだったのですが、Juliaの場合には、たったの

pkg> generate Foo
pkg> dev ./Foo

だけで、パッケージを作って、usingで使えるようにできるので、julia foo.jl スタイルとの比較で増える手間はないです。

タグ： Julia言語

posted at 12:05:41

#Julia言語不合理な

❌ foo.jlの編集とshell上でのjulia foo.jlの繰り返し

をやめて

⭕️ パッケージFooの編集とREPL上でのFoo.f()の繰り返し

に移行すれば相当に快適になります。

もちろんREPLの代わりにJupyter, Pluto, VSCodeを使ってもよい。

タグ： Julia言語

posted at 12:00:00

#Julia言語 using Reviseの後にusing Fooしていれば、Foo/src/Foo.jl内の函数などを追加変更したときに、それが自動的にREPL側に反映されるようになります。

例えば、Foo.jl内の函数fの改訂とその実行

julia> Foo.f()

を繰り返す。これなら遅延問題を小さくできます。

タグ： Julia言語

posted at 12:00:00

#Julia言語 .jlファイルは、julia foo.jl の形式ではなく、パッケージFooのFoo/src/Foo.jlとして書き、using Fooで使うようにした方が得です。パッケージは

pkg> generate Foo
pkg> dev ./Foo
julia> using Revise
julia> using Foo
Foo/src/Foo.jlを編集

で作って使えて編集できる。

タグ： Julia言語

posted at 11:50:33

#Julia言語

JITでは実行直前のコンパイルによる遅延の問題があるので、短時間で終わる処理をjulia foo.jlのスタイルで繰り返し実行するのは不合理です。

この点で不合理なことをやって印象を悪くしてJulia入門から離れるのはもったいない。

julia foo.jl

の代わりに

julia> foo()

を使う。続く twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 11:46:44

eval parse @evalparse

@ledell @DataScienceLA Btw this is the generation code if others want to check it out

carbon.now.sh/?bg=rgba%28171... pic.twitter.com/D6GTdMMddr

タグ：

posted at 10:24:56

michio inoue l MATLA @michio_MWJ

@keiyotokei ありがとうございます。
個人的には mlx ファイル自体は MATLAB がないとそのままでは中身の確認が困難なのは不便に感じることがあります。なので markdown に変換するツールを作って github に上げる時など活用しています。

タグ：

posted at 10:22:12

scythe@八方塞がり @keiyotokei

@michio_MWJ 文字化けもあるあるですね。僕もダメでした。そうですね、データと使ったコード、結果とメモ全部が一つにまとまった解析ノートの様に使いたくて、必要に応じてpdfなどに出力出来るのが良かったです。ライブエディターで同じことが出来るので助かってます。

タグ：

posted at 10:15:21

Erin LeDell - @ledel @ledell

@evalparse @DataScienceLA Would be interesting to see how GLM.jl scales with: float columns, categorical columns, larger data. Looks very promising! 👍

twitter.com/evalparse/stat...

タグ：

posted at 10:05:12

@genkuroki そうなんですね…。
教えてくださってありがとうございます。

タグ：

posted at 09:56:36

Erin LeDell - @ledel @ledell

@evalparse @DataScienceLA I used GLMNet.jl by @skornblith during my PhD (~2013-14?). I know it's a wrapper to the R glmnet FORTRAN code, but it was really fast (at least on small data). It looks like the regularized version of #julialang GLM is Lasso.jl? Is that faster than GLMNet.jl?

タグ： julialang

posted at 09:55:11

eval parse @evalparse

@ledell @DataScienceLA Thank you! It is quite a good performer. Although in my skim-over, I am not 100% sure I turned off regularization but it does say None.

Anyway, it's quite fast at 28s on my PC so roughly on par with #julialang with compilation.

Need to have a proper assessment for features pic.twitter.com/lQNnFKkBwp

タグ： julialang

posted at 09:53:38

#Julia言語訂正: Plutoでの表示には、将来deprecatedになるTextではなく、Print、Dump, with_terminal を使った方が良さそうです。 twitter.com/genkuroki/stat... pic.twitter.com/GCpqWAEmfe

タグ： Julia言語

posted at 07:47:02

Erin LeDell - @ledel @ledell

@evalparse You might want to compare h2o::h2o.glm() in R (also available in Python). It's the fastest GLM I know on "large" datasets... (note: we do elastic net by default, like glmnet).

Also, @DataScienceLA has some similar benchmarks here: github.com/szilard/benchm...

タグ：

posted at 07:45:09

#Julia言語要するに、Plots.jl による動画作成のためには、単にプロットの仕方だけをマスターしていれば十分で、複数のプロットをforループで繰り返すコードを書いて、@ gif マクロを適用するだけで動画を作れます。

プロットの仕方さえマスターすれば誰でも簡単に動画を作れる。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 06:55:14

リフレ女子 @antitaxhike

「IMFは財政規律の重要性を今も説くが、その内容は1世代前はもちろん、10年前と比べても大きく異なる」

元はFTとはいえ、日経がこんな意見を掲載するなんて！
世界の潮流が大きく変わっていることはもはや無視できないほど明らか。日本の財務省だけが古い考えに固執している www.nikkei.com/article/DGKKZO...

タグ：

posted at 06:35:56

#Julia言語以下のリンク先の方法でPlutoの画面で

gist.github.com/genkuroki/a50f...

のURLを貼り付ければ、添付画像のようなことをできます。ただし

pkg> add Distributions DifferentialEquations Plots PlutoUI

が必要。

動画ではなく、こういうノートブックを公開するのもよいと思う。 twitter.com/genkuroki/stat... pic.twitter.com/FyM2rUthyt

タグ： Julia言語

posted at 05:45:24

#Julia言語について勉強したい人は、数学に関わる事柄の視覚化になっているような動画を作って、ソースコードと動画を公開しまくるとよいと思います。

例：蔵本モデル↓
nbviewer.jupyter.org/gist/genkuroki...

同じことを10行のコードと百数十行のコードの2種類の方法でやっている。 pic.twitter.com/VhVOf97iLu

タグ： Julia言語

posted at 05:26:46

#Julia言語 Plots.jlによるGIFアニメの作成はものすごく簡単で

@ gif for t in 1:300
t番目のフレームのグラフをプロットするコード
end

や

anim = @ animate for t in 1:300
t番目のフレームのグラフをプロットするコード
end
gif(anim, "foo.gif"; fps=20)

で行けます。

タグ： Julia言語

posted at 05:19:28

#Julia言語グラフを描くときに細かい調節をしようとすると「沼」で、Plots.jlの公式ドキュメントのattributesの表を眺めて何とかしようとすることになります。

docs.juliaplots.org/latest/attribu...
docs.juliaplots.org/latest/generat...
docs.juliaplots.org/latest/generat...
docs.juliaplots.org/latest/generat...
docs.juliaplots.org/latest/generat...

タグ： Julia言語

posted at 05:16:34

Stan (bayes.club/@mc @mcmc_stan

CmdStanR 0.4.0 is now available on Github and the Stan R packages repository.

github.com/stan-dev/cmdst...

Install instructions: github.com/stan-dev/cmdst...

The new minor release brings improvements for stability and speed, support for profiling and using GPUs with OpenCL.🧵

タグ：

posted at 05:15:34

Juliaでのプロットの仕方のサンプルコードはツイッターでのタグ #Julia言語で画像タブを見ると楽にかつ大量に得られる。

タグ： Julia言語

posted at 05:10:44

#Julia言語 Juliaでコピペでそのまま動くコードの例が欲しい場合には、Juliaがv1.0になって安定した2018年8月以降に書かれたコードを探す必要があります。

v1.0以降の安定度は非常に高いのですが、それ以前とのギャップが非常に大きい。

タグ： Julia言語

posted at 05:07:47

❌Juliaがv1.0になった2018年8月にv1.0になったより前のコードは

⭕️Juliaがv1.0になった2018年8月より前のコードは

編集していておかしくなった。

タグ：

posted at 05:04:59

#Julia言語

 youtu.be/FMVOUvWNlLE
Dan Padilha - Modern Numerical Programming with Julia for Astrodynamic Trajectory Design

スライドにJAXAのロゴが！

リポジトリ↓
github.com/dpad/OrbitalTr...

これはDifferentialEquations.jlとその周辺のパッケージ群のパワーがよくわかる講演。

タグ： Julia言語

posted at 05:04:01

#Julia言語 Juliaがv1.0になった2018年8月にv1.0になったより前のコードは現在のJuliaで動かないことが多いので要注意。

Juliaでのプロットの仕方の例を大量に欲しければ私のGitHub Gistで最近のものを検索するとよいです↓

gist.github.com/search?o=desc&...

アニメーションの例↓
gist.github.com/genkuroki/208a... twitter.com/tune3orz/statu... pic.twitter.com/JhJAb99ycC

タグ： Julia言語

posted at 04:45:16

JuliaCon 2023 @JuliaConOrg

MPI.jl: Julia meets classic HPC
Presented by Simon Byrne at JuliaCon 2020
www.youtube.com/watch?v=pV-8Yq...
#JuliaLang

タグ： JuliaLang

posted at 04:42:20

The Unexplained @Unexplained

What would you name it? pic.twitter.com/i8XBTBHaAN

タグ：

posted at 02:19:29

OpenSourcES @opensourcesblog

@DRTorresRuiz @JuliaLanguage @PhDspeaks @OpenAcademics @AcademicChatter Let us know when you have any questions. We are most active on Slack and Discourse see:
julialang.org/slack/

タグ：

posted at 00:25:12

Kim Paolo Laberinto @KimPLab

Happy #WorldArtDay ! Here's a little weird animation I made using #JuliaLang and Blender. Uses OpenSimplexNoise from Luxor.jl to bend space/time so it might also be #generativeart #proceduralart? pic.twitter.com/28RlPN7nhs

タグ： generativeart JuliaLang proceduralart WorldArtDay

posted at 00:22:48