黒木玄 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年06月29日(火)
#Julia言語 最大最小値を求めるときには全順序が必要です。
< は全順序ではないので使うとNaNなどでまずいことがおこる。
Juliaではislessを使うことが定番なのですが、findmin, argminには適さず、isgreaterを別に定義することになりました。
github.com/JuliaLang/juli...
github.com/JuliaLang/juli...
タグ： Julia言語
posted at 23:51:08
#Julia言語 このスレッドの主役の v1.7.0-beta2 における argmax(f, X) については
discourse.julialang.org/t/add-a-2-arg-...
↓
github.com/JuliaLang/juli...
を参照。findmin(f, X), argmin(f, X)の実装にisgreater(x, y)という新函数が使われている所が非自明です。
Colin Caineさんに感謝。
タグ： Julia言語
posted at 23:45:25
#Julia言語 訂正
❌最大値f(x)とされを与えるxの組
⭕️最小値f(x)とされを与えるxの組
こういう書き間違いはよくある。
* 二択である。
* 自明な事柄。
* 面倒な入力に疲れている。
こういう場合につまらないミスが起き易い。
多発しても訂正は自明なので害はないと思う。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 23:31:44
Juliaの外に良いライブラリがあればその存在についても教えてもらえます。いい人が集まり過ぎ。
数値計算に関する質問への非常によく見るアドバイスは「逆行列の計算は避けた方がよい」です。かなり定番になっている。
タグ：
posted at 23:27:11
Juliaのdiscourseでは「そんなに〇〇したいなら、Juliaではなく△△を使ったら？」というアドバイスが結構普通に出て来る。
基礎的な事柄であっても、数値計算についてJuliaに関係があるということにしてdiscourseで質問すると、親切な回答が得られ易いように見える。
あの雰囲気は独特。
タグ：
posted at 23:24:43
Juliaのdiscourseを見ていれば分かることですが、#Julia言語 を使っている人達は(Juliaについて技術的に変なことを言われると反対しますが)、Juliaは便利な道具の1つに過ぎず、便利な道具達を可能な限り小さな手間で組み合わせて使うことが好きなようです。 twitter.com/mtknnktm/statu...
タグ： Julia言語
posted at 23:16:44
フィッシャーが NP とそんなに仲が悪くなかったっていうの、事実ならそれはそれでって感じではあるんだけど、もっと性格悪いのがフィッシャーではないのか…って感じで俺の中のフィッシャー像はだいぶ崩れる。
タグ：
posted at 22:53:28
#Julia言語 argmax(f, X)を(f(X[i]), i)を返す仕様にしたがっている人達がいるのですが、将来breaking changeが許されたタイミングで名前を変えるとしたら、Xの中で最大のX[i]を与えるiを返すargmax(X)の方だと思う。新名称はindmaxとか。
Jeff Bezansonさんは正しい。
↓
discourse.julialang.org/t/findmax-and-...
タグ： Julia言語
posted at 22:32:26
#Julia言語 f(x, y)のグラフを直接プロットせずに、その対数 log(f(x, y)) のグラフをプロットしています。
X = Y = range(-5, 5; length=1001) は X, Y を区間 [-5, 5] を千等分したものにすることを意味し、Iterators.product(X, Y) は大雑把に直積集合 X×Y に対応しています。 pic.twitter.com/gwJORtdghB
タグ： Julia言語
posted at 22:17:37
#Julia言語
f(x, y) = (1 - x)^2 + 100(y - x^2)^2
の直積 [-5, 5]×[-5, 5] 上での最小値mとそれを与える (x, y) が欲しければ、以上で紹介した函数を使って
f(v) = f(v...)
X = Y = range(-5, 5; length=1001)
m, (x, y) = valargmin(f, Iterators.product(X, Y)) pic.twitter.com/JKqOWKwUep
タグ： Julia言語
posted at 22:13:39
#Julia言語
以上のスクショのソースコード
↓
github.com/genkuroki/publ...
スクショはコピペで使えないので、必ずコピペできるように別に公開するべき。
タグ： Julia言語
posted at 22:09:30
#Julia言語 この話のどこが面白いか。
以上の構成は、函数fの配列X上での最大値f(x)を与えるx∈Xを返す函数しかないときに、x = X[i] となるインデックス i を得る方法が pairs(X) という函数を使えばシンプルなコードで得られることを示しています。
これはちょっと面白いと思いました。 pic.twitter.com/rDUHQr9fFE
タグ： Julia言語
posted at 22:07:43
#Julia言語 v1.7.0-beta2に限らず、上の仕様の函数 valargmin(f, X) が与えられているとき、さらに
valindargmin(f, X) = valargmin(f∘last, pairs(X))
とおくと、配列Xについて、最大値f(X[i])とそれを与えるインデックスiとx=X[i]のPair i => x の組を返す函数が得られます。 pic.twitter.com/9R19PPKwEw
タグ： Julia言語
posted at 22:04:12
#Julia言語
Julia v1.7.0-beta2 では argmin(f, X) という函数が使えます。これは函数fのX上での最小値f(x)を与えるx∈Xを返します。これを使って、
valargmin(f, X) = (x = argmin(f, X); (f(x), x))
とおくと、最大値f(x)とそれを与えるxの組(f(x), x)を返す函数が得られます。 pic.twitter.com/b3AN63xLzL
タグ： Julia言語
posted at 21:59:15
おお、Juliaでigraph 使えr… "run(`$(PyCall.python) -m pip install python-igraph`)" A tutorial on igraph for Julia www.juliabloggers.com/a-tutorial-on-...
タグ：
posted at 21:12:05
突発性難聴やっちまいました。ツイッタで「突発性難聴なったら一刻も早く耳鼻科行ってステロイド処方してもらえ」と予習してたので、すぐ行動に移せました。昨冬の帯状疱疹もフォロワのお蔭ですぐ動けた。だから皆さんも覚えておいて下さい。「帯状疱疹と突発性難聴はすぐ病院」な。
タグ：
posted at 21:08:37
Even soap bubble freeze and when they do it, it's an amazing show [full video: buff.ly/2MqeNRH] pic.twitter.com/PyrQvvuCB4
タグ：
posted at 21:00:00
#Julia言語
Arb github.com/fredrik-johans... の作者とPython mpmathの作者は同じだったのか。作者曰く
discourse.julialang.org/t/hypergeometr...
【Nemo.jl (Arb)はここではmpmathよりも厳密に優れているため、mpmathを比較する必要はありません。】
これからはJuliaからArbを使うことにしよう。
タグ： Julia言語
posted at 20:39:06
#Julia言語 添付画像に転載したプロット結果は
discourse.julialang.org/t/ann-graphmak...
[ANN] GraphMakie.jl - Plotting graphs… with Makie!
より。これはJuliaのExprオブジェクトのツリーなのですが、こういう表示は非常に教育的でよいと思います。誰か、これで色々遊んで見せてくれないものか。疲労状態。 pic.twitter.com/hT2Son3tgE
タグ： Julia言語
posted at 19:09:34
The ability for real-time animations is pretty cool! pic.twitter.com/I8u33BdhCW
タグ：
posted at 18:46:30
「○○はパロディの許可を取っててえらいみたいなツイのバズりをいくつか見てるが、あんまりえらくない。「パロディをやるときには許可を取るべき」という不寛容を助長するのではないか。許可を出した側はふつう。
タグ：
posted at 18:23:10
Graph plotting in Makie finally got a huge update:
discourse.julialang.org/t/ann-graphmak... pic.twitter.com/D4PIOjmSnL
タグ：
posted at 18:18:40
「妊娠女性における mRNA Covid-19 ワクチンの安全性に関する予備的知見」 www.nejm.jp/abstract/vol38... は米国ワクチン有害事象報告システム（VAERS）の 20 年 12 月 14 日～21 年 2 月 28 日のデータを用いたもの。「今回の予備的知見からは…妊娠女性における明らかな安全性シグナルは示されなかった」
タグ：
posted at 18:17:52
「全国規模の集団接種という状況下での mRNA Covid-19 ワクチン BNT162b2」は、イスラエルの研究で、治験のRCTと一致する結果が得られたという論文の紹介。www.nejm.jp/abstract/vol38...
タグ：
posted at 18:15:22
坪野さんは、ぼくが疫学ってなんぞやと気になり始めた2000年代の始めの頃に、いろいろ教えてくださった方です。これからどんなふうに展開していくのか楽しみです。ちなみに、坪野さんが監修しているNEJMの日本版のウェブサイトはこちら（だと思う）。www.nejm.jp
タグ：
posted at 18:12:23
坪野吉孝さん（東北大＆早稲田大）が、COVID-19の研究を例とり疫学を語る新連載が、勁草書房のウェブサイトで始まりました。NEJMの日本版監修を20年！も続けている坪野さんが主にNEJMに出た論文をとりあげて語るそうです──コロナ時代の疫学レビュー 第1回 感染と情報の爆発 keisobiblio.com/2021/06/29/tsu...
タグ：
posted at 18:08:15
【可愛い】藤井聡太王位が注文「ぴよりん」が話題に、サイトは一時ダウン
news.livedoor.com/article/detail...
藤井王位(棋聖)が王位戦七番勝負の第1局で、午後のおやつに「ぴよりん(ひよこ型バニラアイスクリーム)」を注文。ネット上で大きな話題となった。 pic.twitter.com/jE74kxsrHz
タグ：
posted at 17:57:04
Juliaもめちゃくちゃ雑に書ける。単に遅くなるだけ。
Pythonがものすごく普及してしまったことが原因で生じているPythonへの悪しき過学習(過剰適合)の治療が今後問題になると思う。
タグ：
posted at 17:15:20
algebraic geometer @BarbaraFantechi
@TamasGorbe I typed for the first time \def\eps{\varepsilon} in 1988 and never went back.
タグ：
posted at 17:15:11
@togekiss1234 #Julia言語 最近書いたin-place計算でメモリアロケーションを削減して高速化した例が以下のリンク先にあります。アロケーションの削減をした後に並列化すると結果的に全体で10倍以上速くなりました。
多分、みんなやりたいのはこれ。
github.com/genkuroki/publ...
↓
nbviewer.jupyter.org/github/genkuro... pic.twitter.com/VXBv84RH4K
タグ： Julia言語
posted at 17:01:51
言っている意味が分からない。採点の画像との関係も不明。 twitter.com/chikoriroom/st... pic.twitter.com/yvTXAoF7Pd
タグ：
posted at 17:00:11
実際には計算速度の劣化が気にならない場面は非常に多くて、型の伝搬が一意に決まってなくても動くことは、Juliaが気楽に使えるプログラミング言語であることを支えまくっている。
しかし、計算速度が重要な場合には型の伝搬がうまく行くように書く必要がある。それが @ code_warntype が必要な場面。
タグ：
posted at 14:08:43
Juliaは実行時に引数の値の型を見てから、それに合わせて最適化されたネイティブコードが得られるようにコンパイルを実行する。
だから、型の伝搬がうまく行かずに、型が一意に決まらない場合が出てくると可能な型すべてで動くようにコンパイルを実行して、計算速度は大幅に劣化する。
タグ：
posted at 14:08:43
重要なポイントは、例えば、
f(x::Real) = 2x
というコードを書いたときに、
❌引数xの型Realがどのように伝搬するか
を考えるのではなく、
⭕️Real型のxの値3や5.7(それぞれInt型とFloat64型)がどのように伝搬するか
を考えなければいけないこと。
タグ：
posted at 14:08:42
Juliaだと、引数の型注釈などが無用になって、結果的にジェネリックなアルゴリズムを素直に書いただけに見えるコードになり、すっきり見易くなる。
そのときに、支払われるトレードオフは @ code_typed や @ code_warntype で函数の引数の値の型がどのように伝搬するかの確認が重要になること。
タグ：
posted at 14:00:34
discourse.julialang.org/t/cas-benchmar...
【～a macro to add type declarations to every subexpression,～then a good lisp implementation would generate fast code.】
【Maybe julia stands for “Jeff’s uncommon lisp is automated”?】
だそうです。
J→Jeff’s
u→uncommon
l→lisp
i→is
a→automated
😊 twitter.com/trop1000/statu...
タグ：
posted at 12:57:35
入試問題(現代文)が長文化すればするほど、子どもたちは文章とじっくり対峙する機会を奪われ、世の中には、「いかに効率よく拾い読みするか」という、情報処理に特化した技術ばかり溢れてしまう。社会を生きていくためにはもちろん必要な技能だけれど、そればかりじゃ、粗い読み手しか育たんのですよ。
タグ：
posted at 12:28:43
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
まお 情報通信技術（旧 ICT）支援員9 @Kyouiku_ICT
英語に自信がある方は、Julia言語について基礎から学べるオンライン講座がありますので、そちらを参考にしても良いかもしれません。
費用は調べてみましたが、今のところ、それらしい記述はありませんでした。
juliaacademy.com/courses twitter.com/techcampinfo/s... pic.twitter.com/hwghF5oRn2
タグ：
posted at 12:14:46
@kado_judo0312 X -= U/2
Y -= V/2
としている理由は、矢線の中点を対応する座標に来るようにプロットするためです。これは好みの問題。
タグ：
posted at 12:13:13
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Amazonはギグワーカーの配達夫、Flexをアルゴリズムで首にしていて既に問題が多発しているが、代わりはいくらでもいるという態度を取っている。
Amazonは糞 twitter.com/arstechnica/st...
タグ：
posted at 04:42:25
まとめ: NASAは月面着陸船の開発予算に$3.3Bを要求。議会は$850Mしか出さず。NASAは最も安い選択肢としてSpace Xを選択。議会が抗議。NASAはもっと金を出さない限りSpace Xの選択は変わらずと回答 twitter.com/arstechnica/st...
タグ：
posted at 04:38:52
Re: RT 「何を言っているかがさっぱりわからないので質問してみた。しかし、回答のお陰で何を言っていたかが分かった。想像以上に馬鹿げたことを言っていてびっくりしたが、回答してくれてありがとう」という意味で「理解した」と言っているのに、「同意した」という意味だと誤解されることがある(笑)
タグ：
posted at 02:29:19
Julia Inclusive @juliainclusive
Because gender is multifaceted, anyone who sees their gender as under-represented within the Julia community (incl. women, non-binary people, trans people of all genders, people exploring or questioning their gender, ...) is welcome to join!
medium.com/@quinncrossley...
タグ：
posted at 02:14:25
Julia Inclusive @juliainclusive
Julia Gender Inclusive is a new initiative to improve gender diversity and inclusion in Julia through community building, outreach, and support!
Sign up here for an upcoming social meet-up: forms.gle/dqc4EHHxB9LpEz...
discourse.julialang.org/t/announcing-j...
タグ：
posted at 02:11:30
学術論文は公的言説です。卒論も同じ。あなたが「何を理解し」「何に同意するか」を述べるところ。Aを理解したが同意しない、なぜなら、とその論拠を述べる。「納得」とか「腑に落ちる」はその前段階の舞台裏の出来事にすぎません。エッセイならぬ学術論文の仕事は舞台裏を見せることではないのです
タグ：
posted at 01:57:43
Julia programming language tackles differential equation challenges: Scientists in Germany are using it to explore the challenges involved in transitioning to all-renewable power generation. tinyurl.com/32pn6c9v
#JuliaLang #energy #ML
posted at 01:50:52
New post: Coloring in Scientific Publications - www.juliabloggers.com/coloring-in-sc... #julialang pic.twitter.com/DyAIDLwzs4
タグ： julialang
posted at 01:45:00
New post: Optimizing Julia code: Improving the performance of entropy-conservative DG methods in Trixi.jl - www.juliabloggers.com/optimizing-jul... #julialang pic.twitter.com/injlcFkZjn
タグ： julialang
posted at 01:44:55
The mesmerizing data visualization of 24 hours of marine traffic on the Baltic Sea [source and credits: bit.ly/2iDeZjM] pic.twitter.com/bnRkGNwVo8
タグ：
posted at 01:30:00
Mark Kittisopikul ht @markkitti
It's not a race, but I'm quite impressed with the efficiency of the #JuliaLang ecosystem. At the time of this tweet, the Conda binaries are more than a year old and the PyPi wheels are a week old.
anaconda.org/conda-forge/c-...
pypi.org/project/blosc2/
タグ： JuliaLang
posted at 01:13:22
Mark Kittisopikul ht @markkitti
For my current work, the exciting potential about Blosc2 is easily compressing lightsheet microscopy data larger than 2 GB faster than I can copy the same data in memory.
タグ：
posted at 00:59:47
Mark Kittisopikul ht @markkitti
Bindings for #JuliaLang against the above binaries are available here:
github.com/waralex/Blosc2...
タグ： JuliaLang
posted at 00:59:47
Mark Kittisopikul ht @markkitti
The Blosc2_jll binary distribution has been updated to reflect the C-Blosc 2.0.0 release. Download shared libraries for Windows, Mac, and Linux as available for x86, ARM, and PowerPC platforms.
github.com/JuliaBinaryWra... twitter.com/Blosc2/status/...
タグ：
posted at 00:59:47
新井紀子氏が、いかに駄目駄目であるか、あらためてあきらかになりました。
「書き順を理解」って、なんだそりゃ！
合理的な根拠のない、学校の内部でしか通用しない人為的なきまりを「理解」できるものなのか？
正常な思考力をもつ人間には「書き順を理解」することは不可能だ。
#超算数 twitter.com/RochejacMonmo/...
タグ： 超算数
posted at 00:28:14
まあしかし、俺の歳になると、あんまり痩せてると老けるんだよね。めっっちゃシェイプアップしてるひと、シワ多いひと多いし。やっぱりもうこの歳になると、適度に肉がついてたほうがいいんだよね。
と言いながらいまこの時間にカステラ食ってます。カステラうめー
タグ：
posted at 00:22:03
Mark Kittisopikul ht @markkitti
@HenriquesLab @napari_imaging Or Java tools like View5D from @HeintzmannLab :
github.com/RainerHeintzma...
タグ：
posted at 00:15:26
Mark Kittisopikul ht @markkitti
@HenriquesLab You can still use Python tools like @napari_imaging :
github.com/mkitti/Napari.jl
タグ：
posted at 00:14:41
Mark Kittisopikul ht @markkitti
@HenriquesLab Here's a demo of seam carving image processing being taught at MIT:
youtu.be/rpB6zQNsbQU
タグ：
posted at 00:10:11
Mark Kittisopikul ht @markkitti
@HenriquesLab #JuliaLang has no GIL and can compile code to your favorite GPU framework if you want to do 2021 computing.
juliagpu.org
タグ： JuliaLang
posted at 00:07:20
