黒木玄 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年03月22日(月)
#超算数 #超国語 #十分 #星はキラキラorきらきら
文科省は「大臣の記者会見や、」
www.mext.go.jp/b_menu/b002.htm
傘下の文化庁は「政府は,」www.bunka.go.jp/seisaku/bunka_...
タグ: 十分 星はキラキラorきらきら 超国語 超算数
posted at 00:02:41
自分の作った TranscodingStreams.jl がどれくらい使われてるのかJuliaHubで見てみたら、841パッケージもこのパッケージに直接・間接的に依存しているらしい。全Juliaパッケージの16%にもなる。
タグ:
posted at 00:24:44
村上真雄 MURAKAMI Shinyu @MurakamiShinyu
PDFの漢字が康熙部首に化ける問題、Chromeの最新Canary版(Chrome 91)のPDF出力では直っていること確認しました🌟
PDFから元のテキストを取り出せるようになり、アクセシビリティがだいぶ改善します。
ChromeのPDF出力を利用する #CSS組版 @Vivliostyle で問題だったことがやっと解決へ。 twitter.com/MurakamiShinyu...
タグ: CSS組版
posted at 00:35:59
乾いたパチンという音はハサミがぶつかって出ている音ではなく、キャビテーションの泡がはじけるときに出る音です。その撮影には数万フレーム毎秒のハイスピード撮影が必要です。このテッポウエビはハサミも大きく、ほんとにパッチンするのかなと疑問でしたが、実際に聞いて納得したときの動画です。
タグ:
posted at 01:13:10
キャビテーションが音の原因であることを実証したのは上の記事でも紹介されている次の論文です。
Versluis, M., Schmitz, B., Von der Heydt, A., & Lohse, D. (2000). How snapping shrimp snap: through cavitating bubbles. Science, 289(5487), 2114-2117.
science.sciencemag.org/content/289/54...
タグ:
posted at 01:25:17
ただし、ハサミの形からも分かるかもしれませんが、異なる種のテッポウエビです。
どのような形ゆえにキャビテーションがハサミからある程度離れたところではじけるのか、関心ある謎でまだよく分かってません。
タグ:
posted at 01:28:32
#Julia言語 正五角形上の自由境界条件のもとでの波動方程式の数値解の動画。値がNaNだと作画しないことを使っています。
実際には23729個の質点がバネで繋がっている状況の運動方程式をDifferentialEquations.jlで数値的に解いています。
こういうことがシンプルなコードで可能。 pic.twitter.com/chZvlwq2zB
タグ: Julia言語
posted at 01:29:54
この論文のすごいところは、方程式から導出される圧力の変化と、水中マイクの時系列データとがほぼ一致しているというところです。理論と現実がほぼ一致しているのです。
タグ:
posted at 01:38:43
#Julia言語 ソースコード↓
nbviewer.jupyter.org/gist/genkuroki...
In[4]より前の部分は準備です。
以下の添付画像中の解説を見てから、In[4]より前の部分を見た方が多分楽に解読できる。
これと似たような時間で重力で相互作用している300体問題の数値解を求められると思います。誰かやらないかな? pic.twitter.com/mwnjUp1fPh
タグ: Julia言語
posted at 01:57:34
#Julia言語 #数楽 こういう動画を作るときには、境界での値が固定されるディリクレ境界条件よりも、境界での値もびらびら動くノイマン境界条件の方が見た目的に面白い動画が出来上がります。
境界条件は私のコードでは力(もしくは加速度)を記述する函数 g! の定義の中に自然に組み入れられている。 pic.twitter.com/DH2bwGEMoy
posted at 02:04:38
Tried Julia for first time...🤩🤩🤩🤩. Cool one, give it a try.
@JuliaLanguage
#julialang
#DEVCommunity
#javascript
タグ: DEVCommunity javascript julialang
posted at 02:06:59
amazing technique by @zozuar for creating an infinite mountainscape, by adding ridged details at increasing frequencies twitter.com/zozuar/status/... pic.twitter.com/LxtIbFc5h8
タグ:
posted at 02:15:11
In 1930 the Indiana Bell building was rotated 90°. Over a month, the structure was moved 15 inch/hr, all while 600 employees still worked there. There was no interruption to gas, heat, electricity, water, sewage, or the telephone service they provided. No one inside felt it move. pic.twitter.com/07lf2hsgGT
タグ:
posted at 03:11:37
People could still enter/exit the building thanks to an entryway which was connected to a special curved sidewalk.
Bell had planned to demolish it but that would've interrupted phone service for a big chunk of Indiana, which they didn’t want to do.
タグ:
posted at 03:11:37
They lifted the whole building with steam-powered hydraulic lifts, then set it on enormous pine logs. It was moved via hand-operated jacks, which pushed it over the logs 3/8" at a time. Once the building rolled far enough forward, the last log would be moved to the front.
タグ:
posted at 03:11:37
The rotation plan was conceived & executed by famous architect Kurt Vonnegut Sr (father of the famous author).
The building was demolished in 1963.
タグ:
posted at 03:11:38
第1種カニンガム鎖の桁数と長さの表を3つ作ってみた。2つ目は各長さについて、その桁数までの累計。3つ目はさらに各桁数について、長さ以上の合計。2つ目は twitter.com/nekomath271828... と一致。3つ目は arxiv.org/abs/2103.04483 の推定値と比較するため。 twitter.com/kazmuzik/statu... pic.twitter.com/uEMp12I3hc
タグ:
posted at 04:45:01
小菅 信子 Nobuko M KOSUG @nobuko_kosuge
今回の妖怪大戦争で「人文系は信頼をなくした」とツイートされているかたがた、
人 文 系 は こ れ が デ フ ォ で す。
タグ:
posted at 08:20:53
arxiv.org/pdf/2103.04483... のConjecture 3.7の式で、Table 7にあるadjustment factor Bkの値を使って、N=10^9のpredicted ratioを求めてみた。k>=8についてはBkの値が乗っていないので、 Hardy-Littlewood constantsを使って計算しなければならない。 pic.twitter.com/ujSJ9uZbJg
タグ:
posted at 08:32:25
Julia版詳解確率ロボティクスの第3章をスタート。
まずは2輪モデルのロボットが別々の動きをするアニメーション実装。徐々にロボティクスらしくなってきた。
アニメーション実装はPythonよりもJuliaの方が書くのが楽です。
github.com/ShisatoYano/Ju... pic.twitter.com/boymQE18fd
タグ:
posted at 08:56:35
工学系だと「X氏の提案したAという手法はここが不十分なので新たにBという手法を提案します」みたいな流れだけど、ある種の人文系だと「X氏の主張は不十分なので俺がより正しい主張をします」みたいな言い方になるのかなあ
タグ:
posted at 09:00:55
ハーディ・リトルウッド定数のリストがどこかにないか探してみたけど見つからない。とりあえず、ハーディ・リトルウッド予想 ja.wikipedia.org/wiki/%E3%83%8F... について理解を深めるところから。日本語のブログでは、integers.hatenablog.com/entry/2019/04/... や python.atelierkobato.com/twin/ など‥
タグ:
posted at 10:27:57
Ninja DAO | CryptoNi @CryptoNlnjaNFT
嫌韓的な発言の内容は、自分の専門からして困ったなぁ、やめて欲しいなぁ、と思うけど、今の日本では余りに多すぎて、この点をつつきだしたら、かなりの割合の人がアウトになると思うんだよね。それくらい広がってしまってる。
タグ:
posted at 10:42:52
Ninja DAO | CryptoNi @CryptoNlnjaNFT
自分がブロック中心の対応になっているのも、こういう状況に対して強い無力感があるからなんだよね。だって我々が袋叩きになってても皆さん、助けてなんかくれないもんね。
タグ:
posted at 10:46:26
福島の甲状腺検査については、手術実施が10年間で200例という圧倒的な多さを問題にすべきだと思います。これが全て「必要な手術」だったとするなら、2000人にひとり程度の子どもが手術を必要としていることになります。全国では膨大な数の子どもが放置されているはずというその結論は正しいのか
タグ:
posted at 11:33:03
常識で考えて正しくない。スクリーニングをしなければ小児甲状腺がんは年間100万人あたり数人ですから、これが「手術が必要な子ども」の上限です。
必要を遥かに超える多数の子どもたちの甲状腺を手術してしまったことは認めるしかないのではないでしょうか
タグ:
posted at 11:36:03
#Julia言語 普通にJuliaの(野良または公式)パッケージ化した方がよいかも。(Juliaなら野良パッケージ作成は相当に簡単)
ただしその場合にはパッケージ名に"Julia"も文字列は入れない方がよい。
pkgdocs.julialang.org/v1/creating-pa...
野良及び公式パッケージ作成の基本→ github.com/invenia/PkgTem... twitter.com/4310sy/status/...
タグ: Julia言語
posted at 11:49:54
#Julia言語
Juliaでは入力の値の具体的な型を見てからそれに合わせてネイティブコードにコンパイルする仕組みで、支払わなければいけないトレードオフは当然発生する。
しかし、入力の値の型ごとに最適化されたネイティブコードに自動的にコンパイルされる仕組みは多くの場合に便利です。 twitter.com/ray_model2k00/...
タグ: Julia言語
posted at 11:59:34
年間100万人あたり数人を超える多数の甲状腺がんを治療するのがなんらかの意味で有効というエビデンスはありません。
逆に、もし有効なのだとしたら、全国の子どもたちに甲状腺検査をするべきだという結論になるはずですが、それは相当に常軌を逸しています
タグ:
posted at 12:00:43
#Julia言語 私はよく使う自分専用の道具を、自分専用のパッケージMyUtils.jlに入れて使っています。野良パッケージとして公開もしています。
github.com/genkuroki/MyUt...
nbviewer.jupyter.org/github/genkuro...
①前処理と後処理付きの@ my_threadsおよび@ my_distributed
②printfの函数版
③インデックスを1つ削除
タグ: Julia言語
posted at 12:29:26
#Julia言語 では函数fの型は「fの型」になります(笑)
fの型がA→Bになったり、a→Φ(a) (∀a)になったりしない。
f(x)のコードを書くときには、xの値の具体型からf(x)の値の具体型が適切に決まるように書くことが基本になるのですが、必ずしもそう書く必要がないお陰で色々助かっている面もある。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 13:10:34
#Julia言語
❌静的型宣言による事前の厳密な型検査がない言語は劣っている
❌配列の1始まりインデックスの仕様は0始まりより劣っている
❌クラスがないとべたがきするしかない
のような馬鹿げた考え方は潰す必要がありますが、そのせいでJuliaがお気楽言語であるという事実が伝わらないのも困る。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 13:23:10
An Amazing Prime Heuristic arxiv.org/abs/2103.04483 の著者 Chris K. Caldwellの UTM(University Tennesee at Martin)には PrimePages primes.utm.edu があり、Cunningham chainのページ primes.utm.edu/glossary/page.... に、B8=677.356 があった。
タグ:
posted at 13:39:08
#Julia言語 では任意の型のオブジェクトfにf(x), f(x,y), …のように使える函数の機能を追加できます。
だから、函数を引数とする函数φを function φ(f::Function) ~ end と定義すると、見た目的に動いて欲しいコードが動かなくなるリスクが増える。
function φ(f) ~ end と書くべき。 twitter.com/genkuroki/stat...
タグ: Julia言語
posted at 13:40:18
.@n_petton さんのコメント「最初のグラフ、強姦事件の認知件数(用は警察が確認した数)を認知率に変えるというデータ改竄を行ってるのはさすが社会学者と言わざるを得ない。」にいいね!しました。 togetter.com/li/1685531#c89...
タグ:
posted at 13:44:46
k=8も含めて、10^9, 10^10, 10^12までの predicted ratioを求めてみた。Predicted integralの方が精度が良いので、今度は Julia言語で数値積分をやってみる‥ pic.twitter.com/5xqvf4kJTW
タグ:
posted at 13:59:00
そういえば、ある点をある点のまわりに回転するときは、前の課程では行列を使っていたけれど、今は、複素数平面になった。問題を解いているときに、回転のためだけに一度中断して、複素数平面を出すのも変な話と思っているのだけれど、国際標準はわからないけれど、国内標準になってしまった。
タグ:
posted at 15:03:01
餌という目標に対して数ミリ秒から数十ミリ秒程度の弾道的運動を<身体ばね>を活用して実現していてすごい。どうやってるんでしょうね?
この点シャコパンチも通じるところがあります。身体に情報処理をアウトソースしないと生理学的に実現不可能と思われます。
twitter.com/uta_31/status/...
タグ:
posted at 15:17:47
@kamo_hiroyasu @kale_aojiru Killing formは私は結構よく使う曰く付きの数学用語。
Killingさんが実名な点は小事で、そもそもKilling formの概念を考えたのはKillingさんではなく、E. Cartanさんで、Killing formという数学用語を作ったA. Borelさんは「どうしてその名前にしたのか覚えていない」と言っている。🤣🤣🤣🤣🤣
タグ:
posted at 16:56:13
@kamo_hiroyasu @kale_aojiru 今さら名前を変えられない専門用語であっても、この手のひどい話(笑)があったりするので、人の名前が付いている専門用語を見てもその概念を提唱したのがその名前の人だと安易に信じちゃうのもまずい(笑)
個人的な意見ではこの点に関しては特に数学者達がひどい。😊
en.wikipedia.org/wiki/Killing_f... pic.twitter.com/AyrsxBF52p
タグ:
posted at 17:00:07
鼻血は否定するのが当たり前なんですよね。否定できない人はおかしい。
甲状腺がんについても、甲状腺等価線量が何ベクレルと分かったところでそれだけでがんの発生リスクなど分からないわけで、チェルノブイリ事故と比較するんですよ。だから甲状腺等価線量で比較して全く問題ない twitter.com/hiroma_misumi/...
タグ:
posted at 17:15:18
#Julia言語 27+826=853パッケージに増えていた!
juliahub.com/ui/Packages/Tr... twitter.com/bicycle1885/st... pic.twitter.com/zwRbso5Bn9
タグ: Julia言語
posted at 17:16:19
ホモトピー論だと、Eilenberg-MacLane空間を用いたファイブレーションでホモトピー群を計算するkilling methodがあるので、さらに混乱する。 twitter.com/genkuroki/stat...
タグ:
posted at 17:33:00
#数楽 用意した配列の数十%を無駄に使っている点がちょっと気になったので、五角形版のコードを書き直しました。
数値解の計算が12~13秒程度から7~9秒程度に高速化。
2万数千個の質点がそれぞれ隣の質点とバネで接続されている系
ゆっくり版動画も作った。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/BVFMIg2lsu
タグ: 数楽
posted at 17:33:14
非公開
タグ:
posted at xx:xx:xx
#Julia言語 整理したので微分方程式の記述(特に g!)がものすごくシンプルになりました。
複雑なことはモデルのパラメータの情報が全て格納してあるMy.Param型変数 p にすべて押し付けられています。
「モデル・問題→数値解」の流れはよくある典型的パターンの1つ。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/81AFuD6lfw
タグ: Julia言語
posted at 17:49:11
#Julia言語 添付画像のコード
⓪ p = My.Param(~)の行でモデル記述に使われる情報が1変数にまとめられている。
① prob = DynamicalODEProblem(~)の行でモデル全体の情報が1変数にまとめられている。
② sol = solve(prob)の行で問題を解いている。
これがJuliaでのよくあるパターンです。 pic.twitter.com/ZNajRFkClm
タグ: Julia言語
posted at 18:00:32
非公開
タグ:
posted at xx:xx:xx
「恐そうな数学用語」で、取り上げてました。
twitter.com/toyo9/status/1... twitter.com/genkuroki/stat... pic.twitter.com/aUKQy7jYZU
タグ:
posted at 19:09:00
#Julia言語 ノイマン境界条件(端の動きが自由)の場合から、ディリクレ境界条件(端が固定)の場合に変更するためには、微分方程式を記述している函数 g! をほんの少し書き変えるだけで済みます。
添付画像①②を比較してみて下さい。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/VBcDaNNuV5
タグ: Julia言語
posted at 19:41:03
#Julia言語 これで(かなり手抜きになってしまいますが)、任意の形状の領域上での波動方程式の数値解をDifferentialEquations.jlで得る方法が分かりました。
添付画像はノイマン境界条件の正五角形の場合。しっかり境界が固定されています。 pic.twitter.com/R5QfSd6OMk
タグ: Julia言語
posted at 19:42:55
#Julia言語 色を変えたノイマン境界条件の正五角形の場合。
ソースコード↓
nbviewer.jupyter.org/gist/genkuroki...
解読し易さに配慮しているつもりです。 pic.twitter.com/DeE2UTGs59
タグ: Julia言語
posted at 19:44:00
#Julia言語 以上では私が好きだという理由で正五角形上の波動方程式を扱いましたが、正五角形は任意の領域に容易に置き換えられます。
添付画像はアニュラスの場合(ノイマン境界条件)
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/6VcDkvdMza
タグ: Julia言語
posted at 20:39:09
#Julia言語
アニュラスの場合(ノイマン境界条件)
色だけ変えた
結構キモイ
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/DNYseuRcLe
タグ: Julia言語
posted at 20:40:15
#Julia言語
アニュラスの場合(ディリクレ境界条件)
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/v685gh89Ae
タグ: Julia言語
posted at 20:41:07
#Julia言語
アニュラスの場合(ディリクレ境界条件)
色だけを変えた
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/upbGBC5qad
タグ: Julia言語
posted at 20:41:56
#Julia言語 アニュラスの場合には94232個の質点がバネで繋がった系で離散化を行っています。離散化で得られた常微分方程式をDifferentialEquations.jlでさくっと数値的に解いている。
DifferentialEquations.jlを使う場合にはJulia v1.6以上を推奨。私はnightly build
nbviewer.jupyter.org/gist/genkuroki...
タグ: Julia言語
posted at 20:44:03
非公開
タグ:
posted at xx:xx:xx
#Julia言語 【ネタ】
function f(N)
s = zero(N)
for k in N
s += k
end
s
end
は合法的なコードで f(10^17) も一瞬で計算してくれる!(笑)
f(10^17)
→100000000000000000
たまにこれをやって凹む。 twitter.com/wananananabe/s...
タグ: Julia言語
posted at 22:02:44
#Julia言語 Base.summarysizeというのがあります。
julia> using Primes
julia> Base.summarysize(Set(primes(10^8)))
150995136
あとvarinfo()もたまに使う。 twitter.com/ppp3141592ppp/... pic.twitter.com/TPbelXN96H
タグ: Julia言語
posted at 22:10:56
#Julialang and the satellite Amazonia-1 discourse.julialang.org/t/julia-and-th...
タグ: Julialang
posted at 22:32:15