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

jl.devもくもく会（6月23日水朝）に参加を申し込みました！何やるかまだ決まってません！ jldev.connpass.com/event/215987/?...

タグ：

posted at 23:39:53

もふもふ動画 @tyomateee2

子猫を絶対に見捨てない大人の猫達 pic.twitter.com/djnk9e4KHc

タグ：

posted at 23:35:14

Sheehan Olver @dlfivefifty

InfiniteLinearAlgebra.jl+ForwardDiff.jl:

julia> J̃ = z -> (SymTridiagonal(-2*(0:∞)/z, ones(∞)) \ [1; zeros(∞)])[1];

julia> J̃(0.5) ≈ besselj(0,0.5)/besselj(1,0.5)
true

julia> derivative(J̃,0.5) ≈ derivative(x -> besselj(0,x)/besselj(1,x), 0.5)
true

タグ：

posted at 23:04:03

#Julia言語このスレッドで見せた動画に関連した勉強をしたいなら、今なら Julia をマスターしてしまうと非常に具合がよいと思います。

DifferentialEquations.jl の周辺の現在進行中の仕事を理解すれば最新の考え方が身に付く可能性が高いです。

タグ： Julia言語

posted at 22:43:44

#Julia言語全部、質点がバネで繋がった系に変換してから数値的に解いています。

境界条件について「斜めの境界の法線」のようなものは一切考えずに、もっとも素朴な発想でコードを書いているのですが、結構もっともらしい計算結果が出ている。

タグ： Julia言語

posted at 22:41:47

#Julia言語初期値だけを変えたもの

nbviewer.jupyter.org/github/genkuro...

吸収境界条件の波動方程式

最後に波が消えてなくなる。

吸収境界条件のプログラムがバグっていると波が発生しないはずのところから波が発生しまくったりする。 pic.twitter.com/TU1dIPXbzG

タグ： Julia言語

posted at 22:37:58

#Julia言語初期値だけを変えたもの

nbviewer.jupyter.org/github/genkuro...

ノイマン境界条件の波動方程式

自由境界条件なので端がびらびらする。楽しい。 pic.twitter.com/nVEfSB3KNL

タグ： Julia言語

posted at 22:37:55

#Julia言語初期値だけを変えたもの

nbviewer.jupyter.org/github/genkuro...

ディリクレ境界条件の波動方程式

波動方程式はこうやってゆっくりな動画にすると結構きもい。 pic.twitter.com/8hvoYgAttm

タグ： Julia言語

posted at 22:37:53

早川尚男 @hhayakawa

全く同意。ウチも。。。 twitter.com/takahirohatano...

タグ：

posted at 22:23:47

内村直之 @Historyoflife

東北大の河田さんの元現代新書で、ウェブにあがってるのはいいと思います。 twitter.com/qumaiwao/statu...

タグ：

posted at 22:08:55

波多野恭弘 @takahiroHATANO

政府は、留学生の入国禁止措置を解除せよ - 田中駿介　webronza.asahi.com/national/artic...
ジンさんのことも書かれている。
うちの研究室でも、秋入学の大学院生2人と交換留学生1人、更にポスドク1人が入国できずに研究が開始できない。一律禁止なら諦めもつくが五輪ではザル入国なので憤りを禁じ得ない。

タグ：

posted at 22:04:19

#Julia言語まだバグが残っていそう。

この手のことを扱うコードを書いてみて、細かいところを人間の手で書くと、バグが滅茶苦茶生じやすいです。

何らかの自動化は非常に重要だと思いました。

nbviewer.jupyter.org/github/genkuro...

タグ： Julia言語

posted at 21:51:21

#Julia言語「斜めの辺」なるものを一切考慮せずに、吸収境界条件を実装してどれだけうまく行くかに興味があってのですが、思っていたよりうまく行きました。

数値計算の動画は失敗した場合も面白くなりがちなので、たとえ失敗していても面白い話になると思ってやりました。

タグ： Julia言語

posted at 21:50:11

#Julia言語誰もが思うことは「斜めの辺での境界条件をどのように扱っているか？」だと思います。

コードを見ればわかるように、任意形状領域対応になっていてかつ、「斜めの辺」なるものは考慮されていません。格子状に縦横にバネで繋がっている系として扱っている。

nbviewer.jupyter.org/github/genkuro...

タグ： Julia言語

posted at 21:48:37

#Julia言語正五角形版

nbviewer.jupyter.org/github/genkuro...

お待たせしました！

吸収境界条件の場合！

微小に吸収に失敗しているが、概ねうまく行っている。 pic.twitter.com/vuZCxj10M0

タグ： Julia言語

posted at 21:47:03

#Julia言語正五角形版

nbviewer.jupyter.org/github/genkuro...

ノイマン境界条件の場合 pic.twitter.com/RrwwZuykx4

タグ： Julia言語

posted at 21:45:31

#Julia言語正五角形版

以前のコードのバグを発見した。私が書いた類似のノートブックを見たら、それ間違ってますから。バネが縦横に走っているつもりが斜め45度になっていた。まあ意味のある計算ではあるが意図とは違う。

nbviewer.jupyter.org/github/genkuro...

ディリクレ境界条件の場合 pic.twitter.com/5QgH2tb5xx

タグ： Julia言語

posted at 21:44:43

Viral B. Shah @Viral_B_Shah

The work that @KenoFischer has been driving on Diffractor along with compiler support from @vtjnash and @JeffBezanson is going to be a huge step forward for #DifferentiableProgramming in general, with #julialang at the cutting edge! twitter.com/tp2750/status/...

タグ： DifferentiableProgramming julialang

posted at 21:42:52

Szymon Bęczkowski @owiecc

知念実希人　物語り @MIKITO_777

Without looking, what is the result of this code in #julialang

!isnothing(Nothing)

タグ： julialang

posted at 20:12:15

つらい……

「観客上限数、私たちの提言はほとんど反映されなかった」「心を病んでしまった専門家、距離を置くようになった専門家も」“専門家有志の会”メンバーが明かす政治との“距離”news.yahoo.co.jp/articles/8fba6...

タグ：

posted at 18:38:23

稲葉可奈子 @kana_in_a_bar

『ワクチンは打っても打たなくてもいいんですけど』てアナウンサーがサラっと言っていて、

コロナに打ち勝つために世界中が英知を結集させた、これからどれだけの命を救うかわからないワクチンが、
あってもなくてもみたいな軽々しい扱いかと思うと、もうなんか涙が出てきました。

タグ：

posted at 18:25:43

でえもん @GreatDemon1701

@musorami そこから5進法の考え方も紹介できるくらいの余力は欲しいかな。

タグ：

posted at 18:19:33

でえもん @GreatDemon1701

Dr. Chris Rackauckas @ChrisRackauckas

@musorami すると、なおさら5個(円)のまとまりを作ることに合理性があることを思い知ることになりますね。

タグ：

posted at 18:15:38

How hard can it be to do GPU-accelerated 3D PDE solving with high stability time stepping methods? Pretty simple in #julialang curtesy of #sciml DifferentialEquations.jl with ParallelStencil.jl. Code:

gist.github.com/ChrisRackaucka... pic.twitter.com/Y16DGRAcco

タグ： julialang sciml

posted at 17:16:52

さわさわ @heinriches0000

「ちょっと詳しくないのでお伺いしたいのですが」の時点で俺だったらもうちびってる twitter.com/muscle_penguin...

タグ：

posted at 16:52:41

霞＠ @ljyubjya

おお、これが噂のめっちゃ詳しい人がいう「ちょっと詳しくないのでお伺いしたいのですが」…！ twitter.com/muscle_penguin...

タグ：

posted at 15:59:25

にょの人@ @ayano9747

@rnk_examer わかる、『はじめまして、ちょっと詳しくないのでお伺いしたいのですが』って最高じゃん？ってなったｗ
峰先生ってたしかガッキーの結婚で仕事できないって落ち込んでた人だよね？笑

タグ：

posted at 14:40:49

Yusuke Mazda MD PhD @yusuke_mazda

「ちょっと詳しくないのでお伺いしたいのですが」on Twitter

「素人質問で恐縮ですが」@JCIRC_IPR

一見腰の低い態度で入ってきて，バッサリ斬る⚔️
→ 虎伏絶刀勢

二段構えの技がないと論破できない😅 twitter.com/muscle_penguin... pic.twitter.com/kHlxt7fcqo

タグ：

posted at 14:25:15

藤色 @erico763409641

峰先生「ちょっと詳しくないのでお伺いしたいのですが」と腰低い状態からの一刀でバッサリが、ついったらんど様式でカッコいい♡

タグ：

posted at 14:11:50

@antimon2 @kuri_ohkawa @akinori_ito ああ、ごめんなさい。元ネタがよく分かってなかった。

タグ：

posted at 13:54:37

あんちもん2 @antimon2

@genkuroki @kuri_ohkawa @akinori_ito でデリミタを指定した場合のkeepemptyのデフォルト値はtrueで、つまり length(split(str,delim))はnum of delims+1になります。デリミタを指定しない場合（暗黙にホワイトスペースをデリミタとする場合）はkeepenmpty=falseがデフォルトという仕様です。

タグ：

posted at 13:54:00

あんちもん2 @antimon2

@genkuroki @kuri_ohkawa @akinori_ito スミマセン勘違いしてました（最初の @genkuroki さんの例が間違ってます）元々の話は「空文字列のsplitでデリミタ（区切子）を指定した時にどうなるか」なのですね。だったらJuliaもsplit("",",")==[""]です。split()は関数で、第1引数がsplit対称の文字列、第2引数がデリミタ（省略可）です。

タグ：

posted at 13:48:51

キャットアパートメントコーヒー @CatApartment

@kuri_ohkawa @akinori_ito より正確な #Julia言語の仕様

デリミタを指定するか否かで、""のsplitの結果が変わる。

split("") → 空配列
split("", ",") → ""のみを成分とする配列

だから、JuliaとSwiftでは違います。

しかし、Juliaのsplitではkeepemptyオプションの指定でどちらにもできる。

さらに錯綜してきた！😅 twitter.com/antimon2/statu... pic.twitter.com/khQsTrUK9p

タグ： Julia言語

posted at 13:35:40

姫は世界一美しいラグドールなんじゃないかしら pic.twitter.com/K2KVRDnN29

タグ：

posted at 13:23:26

Ｋ @SPMK35

静岡の皆さんはこの青山衆議院議員には気をつけた方がいいと思うの🥺

医学的に間違ってると何度指摘されてもきにしてない。

政治家が命や健康に関わる医療についての姿勢がこれでは怖い。 twitter.com/io302/status/1...

タグ：

posted at 13:20:09

るり @ruri_psychiatry

「ちょっと詳しくないのでお伺いしたいのですが」から「まちがいです」の流れほんと好き twitter.com/muscle_penguin...

タグ：

posted at 12:54:32

Akinori Ito @akinori_ito

改訂版
空文字列をsplitした結果→
空リスト派: Perl, Ruby, R, D, Go, Tcl, Julia
空文字列のリスト派: JavaScript, Python, Java, Kotlin, C#, Rust, Crystal

タグ：

posted at 12:50:05

Juliaを使っていると、「何かを犠牲にして、別の何かを得る」という経験をすることが実に多い。

トレードオフに常に直面する。

教育的にはこれは非常によいことだと思います。

タグ：

posted at 12:50:03

#Julia言語引数を経由せずに函数の中でグローバル変数達をナマで使ってしまうことは、「悪い癖」「悪い習慣」「悪い教育の産物」に過ぎません。やめてしまっても不便なことはないので、やめてしまった方がお得。

高等教育の側がそういう悪い癖を維持させて来た疑いが強い。

タグ： Julia言語

posted at 12:27:18

#Julia言語 v1.7以上では

a = 1
b = 2.0
c = 10^3
d = 4.567
p = (; a, b, c, d)

function f(x, p)
(; d, a) = p
～
end

のようにして、pからd,aのみをシンプルに取り出せて便利です。v1.6でも

using Parameters

function f(x, p)
@ unpack d, a = p
～
end

とできる。

タグ： Julia言語

posted at 12:21:06

#Julia言語函数に引数として渡したいグローバル変数がa,b,c,dのとき、

function f(x, a, b, c, d) ～ end

と毎回書くのはしんどいので、

p = (; a, b, c, d)

のように1つの変数にまとめて

function f(x, p)
a, b, c, d = p
～
end

のように書くと楽です。

タグ： Julia言語

posted at 12:17:47

#Julia言語

②constを使えば、計算速度的に函数内で引数を経由させずに使って大丈夫になるが、型が固定されてしまうせいで、インタラクティブに型の変更を伴う値の変更ができなくなる。

①に慣れた方がお得ですが、すでに書いてしまったコードを必要最小限の手間で直したい場合には②が便利です。

タグ： Julia言語

posted at 12:13:27

#Julia言語トレードオフの問題になる。

①グローバル変数の型が固定されておらず、インタラクティブに値を自由に変更して試行錯誤可能だが、計算速度的には、函数内でそれらの値を使うときには常に引数として渡す必要がある。

続く

タグ： Julia言語

posted at 12:13:26

#Julia言語「グローバル変数aを常に f(x, a) のように函数に引数として与え、f(x)=hoge(x, a)のように使わない」は計算速度が重要な場合には徹底する必要があります。

それがしんどい人は、不便になりますが、const a = 1.0 と書いて
aの型を固定すれば、f(x)=hoge(x,a)が型不安定なくなります。

タグ： Julia言語

posted at 12:04:04

こたろう@坐骨神経痛のサランさん @tatatamaoioi

#Julia言語 aがグローバル変数のとき、f(x)=hoge(x, a)は引数xとグローバル変数aを函数hogeに与えて呼び出す函数になります。

グローバル変数aが定義されていなければ、f(1.0)は実行時エラーになる。

f(x)=hoge(x,a)は何か難しいことをやっているわけではないです。 twitter.com/choonggeun/sta... pic.twitter.com/WyLejcQnzZ

タグ： Julia言語

posted at 12:00:30

峰先生の「ちょっと詳しくないのでお伺いしたいのですが」に（出たっ💕）ってなってしまい次の「まちがいです👶」に悶絶しました、ごめんなさい twitter.com/muscle_penguin...

タグ：

posted at 11:50:12

やもす @ShIma_Megagauss

ちょっと詳しくないのでお伺いしたいのですが……… twitter.com/muscle_penguin...

タグ：

posted at 11:40:20

Shoichi Koyama @sh01

Published MeshRIR - room impulse response dataset on meshed grid points💥 You can visualize wave motion🌊
- website: sh01k.github.io/MeshRIR/
- preprint: arxiv.org/abs/2106.10801 pic.twitter.com/dxVpBOmZ0h

タグ：

posted at 10:54:18

内村直之 @Historyoflife

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

ピントがずれているように思います。先送りというよりは根拠なく方向性を打ち出したことの方が問題でしょう。見極めることが可能でしょうか？ twitter.com/asahi_kantei/s...

タグ：

posted at 10:48:12

@genkuroki ありがとうございます。
以下の記事で混乱していたのですが、改めて引用リツイートみたら黒木さんが回答されてて参考になりました汗
twitter.com/MathSorcerer/s...

タグ：

posted at 10:22:01

非公開

タグ：

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

nope @definiteIyNOTok

Improbawall pic.twitter.com/a9nAa4nw0X

タグ：

posted at 07:42:03

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

Thomas Poulsen @tp2750

So much to look forward to, including what @KenoFischer has been up to lately: twitter.com/juliaconorg/st... pic.twitter.com/NYTExRNxxe

タグ：

posted at 03:56:36

JuliaのPlots.jlの主な仕様はずっと変わっていないと思います。 twitter.com/bamboo_nova/st...

タグ：

posted at 03:40:07

【コンパイラ言語何やろうかなぁ】

Julia twitter.com/dynyanko/statu... pic.twitter.com/dCryyCflA3

タグ：

posted at 03:38:54

#Julia言語グローバルスコープで

a = 1.0
hoge(x, a) = exp(-x^2/(2a))
f(x) = hoge(x, a)

と書くと、函数f(x)の中で引数を経由せずにグローバル変数aを使ったことになります。f(x)をループの内側で使いまくりたい場合には避けた方が無難です(ちょっとした計算なら問題ない)。 twitter.com/choonggeun/sta... pic.twitter.com/d8P45BpOdX

タグ： Julia言語

posted at 03:35:31

#Julia言語自分が入力したコードがどのようにコンパイルされるかを自分で確認したい人は

julialang.org/downloads/

から公式バイナリのインストーラをダウンロードして実行すればよい。そしてJuliaを実行してこのスレッドで紹介したコードを自分で入力してみる。

タグ： Julia言語

posted at 02:08:45

Devtalk @dev_talk

Julia v1.7.0-beta2 released! forum.devtalk.com/t/10714 #JuliaLang #devtalk

タグ： devtalk JuliaLang

posted at 02:06:44

非公開

タグ：

posted at xx:xx:xx

アセンブラのコードも確認して最適化を行うことは、#Julia言語界では結構よく行われている。 #Julia言語はそういう最適化の作業を楽にできるインタラクティブな環境の一種だともみなされる。

タグ： Julia言語

posted at 01:57:16

#Julia言語の公式文書のスタイルガイドには、函数の引数に狭すぎる型を書くことはやめた方が良いと書いてあって、実際、引数の型を何も書かずに使った方が合理的な場合が多いのだが、コンパイルの様子をインタラクティブに確認でき、即ベンチマーク可能なので、コードの詳細な最適化も楽にやり易い。

タグ： Julia言語

posted at 01:51:30

非公開

タグ：

posted at xx:xx:xx

文字列の #Julia言語の函数の定義のコードがパースされて、実行されて函数が定義され、その函数がどのような中間段階を経てネイティブコードにコンパイルされるかを見た様子。途中の様子が興味深い。
↓ pic.twitter.com/Uo4swHiVwY

タグ： Julia言語

posted at 01:46:36

Togetter（トゥギャッター） @togetter_jp

コンパイラがループを削除して計算結果が埋め込まれたコンパイル結果を生成することの確認は、添付画像のスクショのように #Julia言語を使っても確かめられるよね。

コンパイルに必要な中間段階のコードもインタラクティブに確認できる。

添付画像ではllvmに渡す時点でループが削除されている。 twitter.com/shigerufujita/... pic.twitter.com/SKh2oQjHA5

タグ： Julia言語

posted at 01:33:28

賃貸を引き払う時に原状回復費用11万円請求されていたが消費者センターに電話したら7500円返還に変わった話
togetter.com/li/1733524.

タグ：

posted at 01:20:07

Julia Bloggers @juliabloggers

New post: Analyzing international football results with Julia - www.juliabloggers.com/analyzing-inte... #julialang

タグ： julialang

posted at 01:02:23

Diego Javier Zea @diegojavierzea

This VS code extension will make it easier to write and modify #JuliaLang docstrings without worrying about line length 😃 twitter.com/DynamicWebPaig...

タグ： JuliaLang

posted at 01:00:55

筋肉博士 @muscle_penguin_

全く勉強もしていないのに鼻息荒くデマを流す議員と真摯に付き合う専門医のリスバトルにおけるいいね数の違いについて pic.twitter.com/WItH6g3tJj

タグ：

posted at 00:44:36

こちらが一月くらい前に公開されたプレプリントです。あくまで小規模の研究ですが、mRNAワクチンを接種したヒトの血中には、ウイルスを中和する抗体はバリバリにあったけど、胎盤のタンパク質であるSyncytin-1と結合するような抗体は「なかった」ということが明らかです。
www.medrxiv.org/content/10.110...

タグ：

posted at 00:04:32

こちらもNY Times
www.nytimes.com/2021/01/26/opi...

タグ：

posted at 00:04:32

少なくとも昨年12月の時点で否定されている記事がありますね。ファイザーの元従業員が言い出した根拠のない主張だと。
NY Times
www.nytimes.com/2020/12/10/tec...

タグ：

posted at 00:04:31

いまだに「新型コロナワクチンの接種で胎盤を攻撃する抗体ができて不妊が〜〜」と言っている人や、なんと国会議員までいるらしいんですが、半年前の時点で専門家から否定されている上、「そんな抗体はできん」という研究まで出てきてもはやデマでしかありません。

タグ：

posted at 00:04:31

Hiromitsu Takagi @HiromitsuTakagi