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

Ultra super coooooooooooooooooooooool!!!!! 🎉🎉🎉🎉🎉🎉🎉🎉🎉🥳🥳🥳🥳🥳🥳🥳🥳🥳🥳🍾🍾🍾🍾🍾🍾🍾🍾🍾🎈🎈🎈🎈🎈🎈🎈🎈🎈🎊🎊🎊🎊🎊🎊🎊🎊🙏🙏🙏🙏🙏🙏🙏 > Juliaで数値計算　その1：コードサンプル〜基本的計算編〜 qiita.com/cometscome_phy...

タグ：

posted at 23:14:13

Shisato @4310sy

Julia For Data Scien @JuliaForDataSci

Julia版詳解確率ロボティクスの解説記事
今回はシンプルなガウス分布を例に、確率分布の計算、その確率密度関数や累積分布関数を扱うJuliaサンプルコードを紹介します。
www.eureka-moments-blog.com/entry/2021/04/...

タグ：

posted at 22:53:01

📢 New Post! First Steps #4: Digging Into DataFrames. www.juliafordatascience.com/first-steps-4-... #JuliaLang

タグ： JuliaLang

posted at 22:44:21

高梨陣平 @jingbay

Windows Insider版のWSLでLinuxのGUIアプリ対応が提供された。Waylandを用いている。WSLがマイクとオーディオにも対応しているので動画再生等も可能に。

なお、タスクマネージャーも拡張され、特定のアプリを絞ることが可能なecoモードが追加されたとの話 twitter.com/verge/status/1...

タグ：

posted at 21:17:59

みょうが @mrkn

感染拡大防止を阻止するのだけはやめて

タグ：

posted at 21:16:01

大' @satodainu

問題の論文: github.com/QiushiWu/Qiush...

ミネソタ大の声明: cse.umn.edu/cs/statement-c...

筆頭著者が学生で共著者が指導教官らしい。その指導教官のツイート。
スレッドの前後を辿ると、いろんな人によるいろんな立場からの意見が読める。
twitter.com/kengiter/statu...

タグ：

posted at 19:17:59

大' @satodainu

なんともひどい話。こんなことへの警戒が必要になると、オープンソースコミュニティは存続できなくなるんじゃないだろうか。
ミネソタ大学は事態を深刻に受け止めているようなので、対応に期待。

＞Linuxカーネルに意図的にバグを混入したとして大学にコミュニティ出禁措置
gigazine.net/news/20210422-...

タグ：

posted at 18:54:07

tamaspapp.eu/post/common-li...
Switching from Common Lisp to Julia

数値計算のエコシステム構築には、完全なマクロによるメタプログラミングが使えた方がよく、さらに単に多重ディスパッチなだけではなく、型情報の伝搬を助けるパタメータ付きの型が使えることが重要らしい。

タグ：

posted at 18:22:56

#Julia言語 Juliaの場合には、多重ディスパッチがパラメータ付きの型で行われる点も重要。

Common LispからJuliaに来た人によれば、パラメータ付きの型で多重ディスパッチできる点は数値計算のエコシステム構築において、Common Lisp Object Systemより、Juliaの方が優れているらしいです。

タグ： Julia言語

posted at 18:19:35

#Julia言語「多重ディスパッチ」自体は「単なるシンプルな仕様」に過ぎないので、それ自体を眺めても何が良いのか理解できるとは思えない。

エコシステム全体でどのように使われているかを見ないと何が良いのか理解できないと思う。

タグ： Julia言語

posted at 18:14:22

#Julia言語 Juliaのリポジトリの昔の議論を見ると、クラスベースのオブジェクト指向デザインパターン(最も典型的には「継承」)を採用し易くなるような仕様が提案されている場合もあって(結果的に却下されている)、相当な試行錯誤の結果、現在のJuliaの仕様が決まっていることが分かります。

タグ： Julia言語

posted at 18:14:21

#Julia言語多重ディスパッチでも擬似的にクラスベースのオブジェクト指向デザインパターンをシミュレートできます。

しかし、そういうスタイルがエコシステム全体もしくは中核部分の無視できない領域にはびこると、クラスベースのオブジェクト指向の不自由さが復活してしまう。

タグ： Julia言語

posted at 18:14:20

#Julia言語「多重ディスパッチ」自体は「単なるシンプルな仕様」に過ぎない。

多分、重要なのは、多重ディスパッチに適したスタイルがパッケージ開発者達で共有されていて、エコシステム全体が複数のパッケージが想定外の組み合わせで利用される可能性を大事にしていること。

タグ： Julia言語

posted at 18:14:20

iwmn @iwmn_bs

VScodeでTexとjulia書けるようにしたけどgitと連携しやすそうだしもしや最強なのでは

タグ：

posted at 17:52:11

非公開

タグ：

posted at xx:xx:xx

Mうら @tchaikovsky1026

積分区間も「g_1(x)からg_2(x)」とかにしたら、いい練習問題になる。 twitter.com/genkuroki/stat...

タグ：

posted at 16:18:07

#数楽

強制力付き古典調和振動子および強制力付き1+1次元波動方程式に関する解説↓

genkuroki.github.io/documents/Diff...

タグ：数楽

posted at 15:40:24

#数楽以上のような計算の価値は非常に高い。

たまに数学科卒らしいのに残念な数学半可通になってしまった人が「微分と積分を勝手に交換しちゃダメ」とか言って、その手の計算を抑制するような雰囲気を出そうとするのが困りもの。

数学科卒ならノータイムで正しい設定&証明が心に浮かんで欲しい。

タグ：数楽

posted at 15:31:34

ごまふあざらし(GomahuAzaras @MathSorcerer

#数楽強制力付き1+1次元波動方程式

u_{tt}(t,x) - u_{xx}(t,x) = f(t,x)

は両辺をxについてFourier変換すると

U_{tt}(t,ω) + ω²U(t,ω) = F(t,ω)

の形になり、1つ前のツイートの結果に帰着できる。

要するに、chain ruleとFourier解析を知っていれば強制力付き1+1次元波動方程式を解ける！

タグ：数楽

posted at 15:26:59

そういうことか！(目からアザラシ)！ twitter.com/genkuroki/stat...

タグ：

posted at 15:20:50

#数楽

∫_0^t sin(ω(t-s))/ω f(s) ds

をtで微分すると

∫_0^t cos(ω(t-s)) f(s) ds

になり、もう一回微分すると

f(t) - ω²∫_0^t sin(ω(t-s))/ω f(s) ds

になる。

タグ：数楽

posted at 15:20:46

マナマの生マナマコ @chor0sK

あとjuliaならsomefunc.(arr1, arr2)で済ませられることをpythonだとわざわざ[some_func(x, y) for x, y in zip(arr1, arr2)]とか書かないといけないのつらすぎる。

タグ：

posted at 15:18:51

#数楽例えば、ω>0とし、

u(t) = a cos(ωt) + b sin(ωt)/ω + ∫_0^t sin(ω(t-s))/ω f(s) ds

は強制力付き古典調和振動子

d²u(t)/dt² + ω²u(t) = f(t)

の解になる。これの確認で上の「公式」にあたる計算が出て来る。

これを知っていると、1+1次元波動方程式の場合もFourier変換経由で解ける。

タグ：数楽

posted at 15:18:19

chain ruleの特別な場合に過ぎず、習っていなくてもその場で考えてそういう計算をできないとまずい場合なので、習わないのが当然の公式。 twitter.com/wa1t_sush1/sta... pic.twitter.com/DrzcOWmohW

タグ：

posted at 15:04:12

イナバノクロウサギさんの言う通りで、chain ruleの特別な場合を特別に「公式」扱いするのはやめた方がよいと思います。 twitter.com/kunne_isepo/st...

タグ：

posted at 15:00:17

ごまふあざらし(GomahuAzaras @MathSorcerer

using StatsPlots すれば plot(Normal(3, 5)) できます。

個人的にはPackageCompiler.jlの使用例にはPlots.jlだけではなく、StatsPlots.jlも入れて欲しいでしゅ。

教育目的で確率分布のプロットを扱わない選択肢はないので、StatsPlots.jlはほぼ必須。 twitter.com/mathsorcerer/s...

タグ：

posted at 14:40:11

イナバノクロウサギ @kunne_isepo

Plots.jl レシピ仕様変更された？？
julia> using Plots, Distributions
julia> plot(Normal(3, 5), lw = 3)
ERROR: Cannot convert Normal{Float64} to series data for plotting pic.twitter.com/34b6eyWvY9

タグ：

posted at 14:19:02

れいな@底なし沼の魔女 @MegaBlackLabel

連鎖公式からすぐ従うからでは twitter.com/wa1t_sush1/sta...

タグ：

posted at 14:06:22

アカン、Linuxカーネルにバグ投げ込んだ学生の担当助教授がこの程度の認識しか持ってない。
そりゃ、ミネソタ大学の倫理委員会が機能してないってなるわな。 twitter.com/kengiter/statu...

タグ：

posted at 13:54:28

PG_MANA@TYO @PG_MANA_

「意図的にバグぶち込んだら出禁にする」という文言を盛り込んだ
GPLv3.1/2.1やApache Licence V2.1が爆誕する回

タグ：

posted at 13:39:05

Hiroyasu Kamo @kamo_hiroyasu

@EzoeRyou 被害がLinuxだけなのかが気になります。

タグ：

posted at 13:25:04

#Julia言語「多重ディスパッチ」は「単なるシンプルな仕様」に過ぎず、それ自体だけで嬉しくなるようなものではないと思う。

みんながそのスタイルに合わせてパッケージを開発していれば、複数のパッケージ達を想定外の組み合わせで使える御利益が生まれて嬉しい。Juliaではそうなっている。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 13:11:02

@tikani_nemuru_M C++で有名なStroustrupさんはx.f(y)のスタイルの採用について "__extremely__ fashionable" と述べているのですが、Sokal&Bricmontの本のタイトル "Fashionable Nonsense" のfashionableと同じような意味だと思いました。😆 pic.twitter.com/p2rnTE8FNU

タグ：

posted at 11:49:54

@genkuroki なるほど。元ネタの増田もその辺りの文脈なのですね。

タグ：

posted at 11:44:37

@genkuroki 「ファッショナブルなオブジェクト指向」と聞いて、日本の現代美術の一派である「もの派」というのを連想しました。
「もの派」は大きな流れとしてはミニマリズムの一種と言えるかな。ファッショナブルとは逆方向に見えて、アンチファッションというファッションだったと思っています。どうでもええ話w

タグ：

posted at 11:42:23

論文が出たあとですらなお「初心者への嫌がらせをやめろ」と白々しく言っているようではなぁ。

タグ：

posted at 11:38:34

多重ディスパッチって、生物の認識に近いような気がしますね。 twitter.com/genkuroki/stat...

タグ：

posted at 11:36:45

@tikani_nemuru_M 「オブジェクト指向」(鉤括弧付き)に関連する言説はソーカルさんが馬鹿にした難解な現代思想のスタイルに近く、常に冗談のネタにもなって来たと思います。

「意識の高い哲学的なプログラマーが語る難解な思想に満ちたオブジェクト指向」は現代的には冗談のネタそのもの。

タグ：

posted at 11:36:23

論文のことについて問い詰められたら、「ちがうよースタティックアナライザーによる警告をみて修正したコードだよーそんなに厳しくまさかり投げるならもう貢献してやんねーよー」
Greg KH「ミネソタ大BANするわ」
lore.kernel.org/linux-nfs/YH%2...

タグ：

posted at 11:33:43

@tikani_nemuru_M これ、初めて聞くと「どんな魔法だ！」と感じるのですが、多重ディスパッチという名前は難しいのですが、内容はシンプルで簡単なので、自分自身で書いたコードで「どんな魔法だ！」と感じていたことを簡単にできるようになります！

この実用性と比較すると、難解な現代思想的な言説は無用の長物。

タグ：

posted at 11:30:18

@tikani_nemuru_M 「任意の数値」として「幅を持つ数値」を使うだけで、シミュレーションに使用するパラメータを摂動したときに宇宙船の分離がどうなるかをシミュレートできる。

宇宙船の分離のシミュレーションのコードを変えることなく、入力として「幅を持った数値」を与えればよい。

タグ：

posted at 11:26:48

@tikani_nemuru_M Juliaの微分方程式のパッケージは多重ディスパッチに合わせた「型安定」なコードで書かれているので、通常の浮動小数点数で微分方程式を数値的に解くだけではなく、「任意の数値」で微分方程式を数値的に解けます。続く

タグ：

posted at 11:25:02

問題の論文。研究倫理として勝手にOSSコミュニティを被験者にして意図的に脆弱性を紛れ込ませ、さらにレビューでツッコミが入ると当初嘘の回答をしたことが問題視されている。そりゃ当然だ。
github.com/QiushiWu/Qiush...

タグ：

posted at 11:24:49

ミネソタ大学がLinuxに意図的に脆弱性を含んだ貢献を試みた上で論文を書いたので、信頼が毀損されたとしてLinuxカーネルからミネソタ大由来の貢献をすべて除去する作業が進められている。マジでやるらしい。
lore.kernel.org/lkml/202104211...

タグ：

posted at 11:22:59

@tikani_nemuru_M 多重ディスパッチの御利益はNASAでのJuliaの使い方の1つを見るだけでも分かります。

NASAでは宇宙船の分離をJuliaの微分方程式のパッケージで計算しています。それとは無関係に開発された「幅を持った数値」を扱うパッケージを使うと、「幅を持った数値」で微分方程式を数値的に解けます。続く twitter.com/genkuroki/stat... pic.twitter.com/wprRgj6VCM

タグ：

posted at 11:22:10

@tikani_nemuru_M こんな感じで、まるで難解な現代思想のごとく語られている「オブジェクト指向」(鉤括弧付き、とてもファッショナブル)のスタイルへのこだわりを相対化して行くと、「なんてアホらしい議論を何十年も続けてやっているんだ！」という気持ちになります。

続く

タグ：

posted at 11:16:53

@tikani_nemuru_M 多重ディスパッチでは、メソッドが単独のオブジェクトの保有物だと考えないだけの話で、x.f(y, z)のスタイルにこだわる必然性はなく、C++の開発で知られるBjarne Stroustrupさんもその点を強調しています。

続く

タグ：

posted at 11:14:28

@tikani_nemuru_M 続き。

x.f(y,z)のスタイルでは「xが保有するメソッドf」という発想になる。

多重ディスパッチでは、函数f(x, y, z)で実行されるメソッドはx,y,zの型の組み合わせで決まるだけで、メソッドはx,y,zのどれの所有物でもないと考えます。

続く

タグ：

posted at 11:12:27

@tikani_nemuru_M シングルディスパッチでは、同名のf(x,y,z)で何が実行されるかはxの型で決まり、x.f(y,z)と書くことが普通になっている。

多重ディスパッチでは、同名のf(x,y,z)で何が実行されるかはx,y,zの型の組み合わせで決まる。

こういう言い方をすると多重ディスパッチの方が自然に見える。続く

タグ：

posted at 11:08:57

@genkuroki 一連の処理にあたって、用意された適切な道具をその都度呼び出していくというイメージなのですか？＞多重ディスパッチ

いあ、プログラミング「も」まるでぜんぜんワカランチンなのですが。

タグ：

posted at 10:55:20

まるで、まるで現代思想のごとく語られる「オブジェクト指向」(鉤括弧付き、すごくファッショナブル)の雰囲気の悪さが気になっていたのだが、「オブジェクト指向」にこだわる人達が #Julia言語の開発者達にみっともなく喧嘩を売っているのを何回も見て、印象の悪さが誤解ではないことが分かった。

タグ： Julia言語

posted at 10:54:18

定義域が異なる同名の函数を作れるようにして(多重ディスパッチ)、定義域にあたるもの(型)をユーザーが定義できるようにするだけの仕様はものすごく自然でシンプルだが、その仕様のもとでパッケージ化されたコード群を組み合わせて仕事するのが非常に楽になることはかなり非自明。

タグ：

posted at 10:49:52

「多重ディスパッチ」という言い方をすると意味不明になるが、その実態は、同名の函数について引数達の型の組み合わせで何が実行されるかが変わるだけ。

定義域が違う同名の函数を作れるだけのシンプルな仕様。

タグ：

posted at 10:46:30

多重ディスパッチ(引数達の型の組み合わせで実行されるメソッドが変わる)の簡単な例(#Julia言語)

f(x, y) = x*y
f(x::Int, y::String) = y^x
f(x::String, y::Int) = x^y

のとき

f(3, 4)
→12

f("pine", "apple")
→"pineapple"

f(3, "meow")
→"meowmeowmeow"

f("meow", 2)
→"meowmeow" pic.twitter.com/Gx5kTGgzsS

タグ：

posted at 10:42:03

ずかし♡手作りツイート @ZukashiT

多重ディスパッチでコードを書くようになったら、「オブジェクト指向」(鉤括弧付き、すごくファッショナブルなやつ)なるものが浅はかなアイデアに見えて仕方がない状態になった。

C++の開発で知られるBjarne Stroustrupさんも多重メソッド=多重ディスパッチ支持の立場。

discourse.julialang.org/t/is-julias-wa... twitter.com/tikani_nemuru_... pic.twitter.com/MrgKoNCZQY

タグ：

posted at 10:29:03

IntroとMethodsはすでに書き上げてある、というのがプロポーザルの前提であり、IntroとMethodsというのはそういう前提で書かれたものというトレーニングの一環ですからね。実験したあとにMethodsを書くのは（かなり広義の）不正っていうのを知る機会になりますし。 twitter.com/ykamit/status/...

タグ：

posted at 09:58:38

Hiroo Yamagata @hiyori13

東日本震災のときのこれもすごかった。自分は何やら裏の情報ルートを持っているんだぜ、という低級なめくばせ。あのとき馬脚が出た人は実に多い。
twitter.com/miyadai/status...

タグ：

posted at 09:55:06

Hiroo Yamagata @hiyori13

れいな@底なし沼の魔女 @MegaBlackLabel

宮台信司はブルセラ援交支持で、変なこと言うともてはやされる刷り込みのせいでどんどんおかしくなったのか、もともとのおかしさがどんどん露出してきただけなのか、過去の人になった焦りで変になったのか。ともあれ、もうどうしようもない。 twitter.com/Kontan_Bigcat/...

タグ：

posted at 09:49:10

ミネソタ大学の研究者、マジモンのサイコパスでLinuxコミュニティもキレるわというか、これはミネソタ大学のパッチが他のOSSでもrevertされる流れにならない？

タグ：

posted at 09:34:37

じょい @joy1192

ミネソタ大学、マジでゴミな研究やったなあ……
水質調査がちゃんと行われてるか調べるために、井戸に毒を投げ込んでみましたみたいなやつじゃん
迷惑系YouTuberかよ
（※わざとLinuxカーネルにバグコード混入させようとした件）

タグ：

posted at 09:26:26

Hiroo Yamagata @hiyori13

楽しい。ドーキンスが「男が『心は女だから女扱いしろ』と言うのはトランスと賞賛されるのに、白人が黒人として活動したら石を投げられるって変だよねー」と言ったら、全米ヒューマニスト協会が過去の賞を剥奪決定したそうな。過去の受賞者たちがいっせいに抗議/嘲笑してる。
www.mailplus.co.uk/edition/commen...

タグ：

posted at 09:10:22

天むす名古屋 Temmus @temmusu_n

#JuliaLang

Gapjm.jl
github.com/jmichel7/Gapjm...
Experimental port of some GAP functionality to Julia

①②Kazhdan-Lusztig polynomial (Gapjm.jl is faster than GAP3)

③invariants of Coxeter group of type F₄ pic.twitter.com/QrL4UH4XxO

タグ： JuliaLang

posted at 09:01:15

#超算数 It's not maths, but #schmaths !

タグ： schmaths 超算数

posted at 08:29:37

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

藤井　恒 @my_flute

天むす名古屋 Temmus @temmusu_n

昨年の4月に愛知県であったコンタミ事例。

www.pref.aichi.jp/site/covid19-a... twitter.com/concert_1750/s...

タグ：

posted at 07:58:39

#超算数 This has been going on for >= 150 years in Germany. Blame in on Diesterweg e.g. for defining multiplication by the lang form A mal B, written A × B, or the *situation* it denotes twitter.com/temmusu_n/stat.... The total *number* of fingers is 20, 4 × 5, 5 × 4, or 2 log₂1024. twitter.com/CarstenPaul198...

タグ：超算数

posted at 07:31:25

非公開

タグ：

posted at xx:xx:xx

Jacob Quinn @quinn_jacobd

Well, here it is folks, at long last. DataFrames.jl 1.0 release. So proud of all the work, discussion, and effort that went into making a fast, user-friendly, stable API for #JuliaLang. pic.twitter.com/bIWOrLPale

タグ： JuliaLang

posted at 07:03:07

非公開

タグ：

posted at xx:xx:xx

Mark 1x Engineer @virullius

I just found a reason to try #JuliaLang
docs.juliaplots.org/latest/example...

タグ： JuliaLang

posted at 06:03:52

ミネソタ大学が見解を示した模様
twitter.com/UMNComputerSci...

タグ：

posted at 05:43:24

Dean Markwick @DeanMarkwick

I accidentally wrote a very inefficient function. dm13450.github.io/2021/04/21/Acc... #JuliaLang

タグ： JuliaLang

posted at 04:20:32

We will report our findings back to the community as soon as practical.

Sincerely,

Mats Heimdahl, Department Head
Loren Terveen, Associate Department Head

タグ：

posted at 04:15:21

We take this situation extremely seriously. We have immediately suspended this line of research. We will investigate the research method & the process by which this research method was approved, determine appropriate remedial action, & safeguard against future issues, if needed.

タグ：

posted at 04:15:02

The research method used raised serious concerns in the Linux Kernel community and, as of today, this has resulted in the University being banned from contributing to the Linux Kernel.

タグ：

posted at 04:14:37

Leadership in the University of Minnesota Department of Computer Science & Engineering learned today about the details of research being conducted by one of its faculty members and graduate students into the security of the Linux Kernel. pic.twitter.com/QE9rrAyyMX

タグ：

posted at 04:14:31

簡単に言うと、バグがあるパッチ送って Linux Kernel コミュニティを研究対象として論文を書いた、って感じです。倫理的にどうなのっていう。

タグ：

posted at 02:10:24

#Julia言語体裁が整ったパッケージ公開を補助するツールが

github.com/invenia/PkgTem...

github.com/JuliaDocs/Docu...

現時点では、手動で色々やらなけれないけないところがあり、かなり面倒。可能な限り自動化されることが望ましい。

こういうところにも面白そうな問題がある。

タグ： Julia言語

posted at 02:04:47

#Julia言語自分専用ローカルパッケージは

pkg> generate Foo
pkg> dev ./Foo

で作って使えるようになる。Foo/src/Foo.jlなどを編集する。

generate Fooで作ったFoo/以下をGitHubでFoo.jlというリポジトリ名で公開すれば、公開された野良パッケージになる。

これだと体裁は全然整っていない。

タグ： Julia言語

posted at 01:59:17

まずは体裁が全然整っていない自分専用野良パッケージをGitHubで公開することから始めるという手もあると思う。

タグ：

posted at 01:55:58

GitHubでmasterブランチの名前がmainに変わったことや、PkgTemplates.jlの開発速度が結構速いことに注意が必要かも。

Juliaでのパッケージの公式登録方法に関する最新の事情に基く解説が必要かも。 #Julia言語

自分専用野良パッケージ作成は易しい。

体裁を整えた公式パッケージ登録は面倒。 twitter.com/taki__taki__/s...

タグ： Julia言語

posted at 01:54:09

問題の論文はこちら
github.com/QiushiWu/qiush...

タグ：

posted at 01:48:14

ミネソタ大学がLinuxから追放された模様...
“I will now have to ban all future contributions from
your University and rip out your previous contributions”
lore.kernel.org/linux-nfs/YH%2...

タグ：

posted at 01:24:20

Stefan Karpinski @StefanKarpinski

@0x00B1 @QiuhaoLi @gregkh That means the UMN IRB cannot be trusted to prevent this kind of research from being done on open source project, so a university-wide ban is warranted

タグ：

posted at 00:54:51

Stefan Karpinski @StefanKarpinski

@0x00B1 @QiuhaoLi @gregkh Initially they didn't get IRB review but then when there was public pushback, they did and "The IRB of UMN reviewed the study and determined that this is not human research (a formal IRB exempt letter was obtained)." (from that document)

タグ：

posted at 00:52:50

Stefan Karpinski @StefanKarpinski