Twitter APIの仕様変更のため、「いいね」の新規取得を終了いたしました

黒木玄 Gen Kuroki

@genkuroki

  • いいね数 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 を見て下さい。
Favolog ホーム » @genkuroki » 2021年06月05日
並び順 : 新→古 | 古→新

2021年06月05日(土)

マナマの生マナマコ @chor0sK

21年6月5日

正確にいうと、
- jupyterコマンドは使えててほしい
- でもpythonのコマンドは汚したくない
- 環境ごとにカーネルを変えることができるようにしたい
- org-modeからjupyter使いたい
- Juliaからも使いたい
- 全部NixOS上で宣言的に設定したい
を全部叶えた。

タグ:

posted at 23:24:35

天むす名古屋 Temmus @temmusu_n

21年6月5日

10月までの再延期なら、尾身氏の述べた懸念事項はほとんど回避することができそう。政府にとり、極小開催と延期とどちらが受け入れやすいだろうか。

タグ:

posted at 23:17:24

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 2つに整数が互いに素(最小公倍数が1)になる確率が 6/π² = 1/ζ(2) になることは大昔から知っていたが、コンピュータを使ったモンテカルロ法で確認したのは #Julia言語 を使い始めてから。

コードを書くのが楽でかつ速いというのは非常に快適で素晴らしい。

タグ: Julia言語 数楽

posted at 23:13:02

天むす名古屋 Temmus @temmusu_n

21年6月5日

日本にも新型インフルエンザのために作られたワクチン接種計画があるはず(すくなくとも議論はあった)だが、新型コロナの接種実態は必ずしもその計画通りではない。首相らが訪米に際して接種しているのもどういう基準だったのかわからないまま。何度改定しても良いので方針は示してもらいたいのだ。

タグ:

posted at 23:12:33

Murawaki @murawaki

21年6月5日

漢語の借用を一切考慮せずに、このデータで日本語と朝鮮語の関係を議論するのは無謀もいいところだ。

タグ:

posted at 23:11:06

Murawaki @murawaki

21年6月5日

LingPyでアラインメントして適当な閾値で同源性を判定しているみたい。結果を直接見てみたい。まあ無茶苦茶だと思う。外群としてドラビダ語族を使うという判断はわからなくもないが、それだけでは駄目。中古漢語も考慮できないと厳しい。NorthEuraLexに載っていないけど。

タグ:

posted at 23:09:23

Murawaki @murawaki

21年6月5日

私がアルタイ系の議論に手を出す気になれないのは、データが最大の問題であって計算機の出る幕はないと思っているから。著者はなぜかNorthEuraLexから単語を集めてきている。どう見てもアルタイ方面は手薄なんだけど。

タグ:

posted at 23:07:48

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 ℤ×ℤには「2つの整数が互いに素になる確率」が得られるような自然な確率測度は入らない。しかし、ℤのある種の「完備化」である

ℤ̂ = proj lim ℤ/n!ℤ = ℤ₂×ℤ₃×ℤ₅×ℤ₇×ℤ₁₁×…

はコンパクトAbel群なのでHaar速度で確率測度になるものが本質的に一意に定まる。

タグ: 数楽

posted at 23:07:43

Murawaki @murawaki

21年6月5日

アルタイ語族の仮説の系統樹に借用の枝を追加したネットワークの尤度を計算。仮説同士の尤度を比較したところ、マクロアルタイ説に対して否定的な結果を得た。著者は生物系の計算機屋さんらしい。hal.archives-ouvertes.fr/hal-03241668/

タグ:

posted at 23:07:00

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 このスレッドの数学ネタ解説

2つに整数が互いに素(最小公倍数が1)になる確率が 6/π² = 1/ζ(2) になることを使ったモンテカルロ法で円周率を近似計算している。

これは rand(Int) と gcd(x, y) の速さを測っているに近い。

タグ: 数楽

posted at 23:02:35

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

スレッド並列のマクロを使うと begin ~ end でも速い。その理由は並列処理のために呼ばれる函数内でほとんどすべての処理が行われるから。

注意:自前の便利なマクロを使っている。 pic.twitter.com/eJIAr1qHzV

タグ: Julia言語

posted at 23:02:35

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

函数にしなくても速い場合。

let ~ end で囲めばその内側で作られた変数はローカル変数になり速い。

begin ~ end で囲んでも無意味。

github.com/genkuroki/publ... pic.twitter.com/cNvBClaERr

タグ: Julia言語

posted at 23:02:33

Akinori Ito @akinori_ito

21年6月5日

これと同じ問題をDで実装して時間を測ってみたが,25まで計算したところで指数っぽい時間になってたのであきらめて外挿してみたら30まで計算するのに7時間半ぐらいの予測.Julia速いな

タグ:

posted at 22:13:47

次郎作 @zilousaku

21年6月5日

「今は流行しているんだから、流行を止めるんでしょうが」
「今はとにかく流行を止めるぞ」
と一喝して、目に涙を溜めてみんなをいさめた、

という出来事が書かれていて、読んで目頭が熱くなりました。
このタイミングで尾身先生が日本にいてくださってるのは、奇跡だよなぁ
www.amazon.co.jp/%E7%90%86%E8%A...

タグ:

posted at 22:11:30

次郎作 @zilousaku

21年6月5日

尾身先生の凄さを垣間見たのは、西浦先生の本でした。

専門家会議の意見が厚労省側で書き換えられるようなことが続き、若い専門家が不満をためている時に、
「厚労省がちびちび書き換えるとか、そんなしょうもない話はどうだっていいんだ。責任取れと言われるんだったら俺がとるぞ」

続く

タグ:

posted at 22:11:29

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan @y_nya_kc #Julia言語 での典型的なワークフローの解説は公式ドキュメントにもあって、その3つ目が、

using Revise
using MyPkg

のスタイルの解説になっています。

パッケージの編集結果の自動反映は魔法のような感じで面白いです。

docs.julialang.org/en/v1/manual/w...

タグ: Julia言語

posted at 21:48:07

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan @y_nya_kc #Julia言語 structの仕様がほぼ確定したパッケージ Foos を書いているときには、Jupyterで

using Revise
using Foos

の形式でFoos.jlを読み込んでおけば、Foos/src/Foos.jl を変更すると、変更結果がJupyter側に自動反映されて便利です。魔法のように便利。

タグ: Julia言語

posted at 21:44:06

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan @y_nya_kc #Julia言語 何を言われたのか見えないのですが、私はJupyterを頻繁に利用していますが、試行錯誤中に頻繁にカーネルを再起動する必要は生じていません。

struct Fooの定義を試行錯誤で決めたい時には、定義を

module O ~ end

で囲んで、O. Foo を使い、定義の変更はモジュール丸ごと上書きで。 twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 21:40:58

三毛山 @XUFd6G7AkfvCAl3

21年6月5日

@genkuroki 勝手に行く息子でよかった。
電話かかってきて何故か怒られたけど🙄自然現象ですからねーと流しました。トイレの話なだけにw
小2の時だったなぁ。膀胱の発達が遅い子はやたらとトイレ近いしそのあたりも知ってほしいと思います。

タグ:

posted at 21:36:38

黒木玄 Gen Kuroki @genkuroki

21年6月5日

教員採用試験で十分に難しい問題を出して、大学で学問をきちんと勉強した人から優先的に採用される仕組みになっていればよいのですが、そうなっていない。

ものすごい弊害が生じていると思います。

教員採用を学力で決めるようにしないとまずい。 twitter.com/sekibunnteisuu...

タグ:

posted at 21:24:07

QDくん⚡️AI関連の無料教材紹介 @developer_quant

21年6月5日

Haskellは一部の金融システムベンダーで使われてますがレアケースかと。
QuantLibにはHaskellに移植中のものが2つあり、一部で好まれているようです
HQuantLib
github.com/paulrzcz/hquan...
QuantHas
github.com/simoncourtenag...
しかしJulia版の方がかなり進捗してます
QuantLib.jl
github.com/pazzo83/QuantL... twitter.com/hanatmetal/sta...

タグ:

posted at 21:23:27

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 NamedTuple bar を bar.a, bar.b もしくは (; b, a) = bar の形式で使っているコードは、bar を対応する struct で置き換えてもそのまま動きます。

後で場合によってはstruct化するかもしれないが、ひとまずそうする必要がない場合にはNamedTupleを使っておけばよい。

タグ: Julia言語

posted at 21:09:54

でんボ @Dmb_fish

21年6月5日

ハコフグの歯から滲み出る巨神兵感。
こんなんなってんのや…びっくり〜🦷
見れば見るほど変やなぁ、、本当に魚なんかな。 pic.twitter.com/UUPt5G5zy8

タグ:

posted at 21:03:22

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia入門

NamedTupleなら抽象型フィールドによる型不安定性の問題が生じない。

structは後で定義を変えるのが面倒になる。

多重ディスパッチを積極的に使うつもりがないなら、NamedTupleで十分。

NamedTuple bar は、bar[1], bar[2] と bar.a, bar.b の両方の使い方が可能。

タグ: Julia入門

posted at 20:59:39

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 structを使って

Base.@ kwdef struct Bar{A, B}
a::A = 1
b::B = 2.0
end
bar = Bar()

のように書くよりも、

bar = (a = 1, b = 2.0)

とNamedTupleを使った方が楽な場合は多いと思う。

タイピングの量が段違いで、structだと初心者は抽象型フィールドを作ってしまいがち。 twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 20:55:30

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 お手軽プログラミングでは「structは重い」と感じることが多く、上のように struct の bar ではなく、NamedTupleを

bar = (a = 1, b = 2.0, c = "three")

と作って、bar.a, bar.b, bar.c を利用したり、v1.7(nightly build)以上なら

(;c, a) = bar

も使ったりすると楽です。

タグ: Julia言語

posted at 20:50:22

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

天むす名古屋 Temmus @temmusu_n

21年6月5日

名誉と実益があったら、ためらわず実益をとれる位置付けだということ。

タグ:

posted at 20:40:58

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 はv1.0になって仕様が固定されたのは2018年8月なので、ググって見つけたそれより前の情報の多くが時代遅れになっています。

おまけ。v1.7(nightly build)であれば、上に続けて

(; c, a) = bar
@ show c a;

c = "three"
a = 1

とできて非常に便利です。 pic.twitter.com/fwv3yLgSOq

タグ: Julia言語

posted at 20:26:09

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 別ツイでは疑問を誤解‼️ごめんなさい‼️

正しい回答は以下の通り。

Base.@ kwdef struct Bar{A, B, C}
a::A = 1
b::B = 2.0
c::C = "three"
end
bar = Bar()

Bar{Int64, Float64, String}(1, 2.0, "three")

mutableでも同様ですが、mutableはあまり使いません。 twitter.com/aeba_chan/stat... pic.twitter.com/6VHg5Ts0Rx

タグ: Julia言語

posted at 20:22:21

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 おまけ

* mutable structはほとんどの場合に不要。例えば

struct Foo{A} a::A end
foo = Foo(Int[])
push!(foo.a, 123)
foo

Foo{Vector{Int64}}([123])

push!(foo.a, 456)
foo.a[1] = 999
foo

Foo{Vector{Int64}}([999, 456])

* お手軽プログラミングにはNamedTupleが便利。 pic.twitter.com/Mry6yktbtM

タグ: Julia言語

posted at 20:05:20

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 だから、問題を記述するパラメータを1つの変数にまとめて、その変数を函数の引数として常に渡すようにするのがよい、というのが結論になります。

そうした方がコードがスッキリして見易くなることが多いです。

サンプルコード↓ twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 19:59:12

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 しかし、実際に使ってみると、定数は不便です。

* すでに定義されているグローバルと同じ名前で定数を作れない。

* 別の型を持つ値に変えられなくなる。たとえば、Int64型の定数をBigInt型の定数に変えて、巨大な整数値でコードを動かしてみると言うようなことができない。

続く

タグ: Julia言語

posted at 19:55:21

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語

Juliaの定数(const)は型が固定されます。例えば

const g = 9.80665

とすると、gの型はFloat64に固定されます。

だから、このgを含むコードを「gはFloat64型である」という仮定を使ってコンパイルできるようになります。だから、グローバル変数を定数に変えると計算が速くなる。続く

タグ: Julia言語

posted at 19:52:28

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 解説

Juliaのグローバル変数は型を固定できません。

だから、グローバル変数を含むコードはその変数がどんな型でも動くようにコンパイルされる。だから計算が遅くなります。

続く

タグ: Julia言語

posted at 19:50:17

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 具体的なコードを表に出せば、適切なアドバイスが得られる可能性が高いと思います。 twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 19:41:47

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan #Julia言語 グローバル変数をナマで使うことはJuliaの計算速度劣化の原因1位だと思います。

公式ドキュメントのPerformance Tipsを見ると多くの疑問が解消し、まずいコードを書く可能性を大幅に減らせます。

docs.julialang.org/en/v1/manual/p...

私は「これ↑の存在を最初に教えてくれよ」と嘆きました。 twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 19:39:00

Takashi SASAKI 佐々木隆志 @TakashiSasaki

21年6月5日

民事上の示談で済ますことなく、きっちりと暴行罪または傷害罪に問うてほしいと思う。 twitter.com/genkuroki/stat...

タグ:

posted at 19:35:01

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 統計モデル D とパラメータ θ を結びつける方で、Distributions.jl とのコンパチビリティを保つことができれば(できるかどうかはやってみないとわからない)、スムーズに置換できる可能性がワンチャンあるかもしれません。

タグ: Julia言語

posted at 19:14:29

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 固定されたパラメータ θ に関する統計モデル D のサンプルを生成したい場合には、Dをθと結合させる方が自然です。

Distributions.jlは後者のDとθの結合を強制してくるのでつらい。

統計モデルDとデータYの方を結合させて欲しい。

誰か、究極も確率分布パッケージを書いて欲しいです。

タグ: Julia言語

posted at 19:11:03

でんボ @Dmb_fish

21年6月5日

いつか拾いたいと思ってたハコフグ!!
この状態のやつを自分で拾って見てみたかったのでめっちゃ嬉しい😆
骨板かっこよすぎるなぁ…。 pic.twitter.com/xQBuvzKnoK

タグ:

posted at 19:10:10

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 確率分布 D と結合させるオブジェクトをデータ Y とパターン θのどちらにした方がよいかは場合によります。

観測データ Y が得られていて、モデルのパラメータが未知の場合には、DとYを結びつけて「統計モデル D 内で観測データ Y と同じ値が偶然生成された状況」と解釈するのが自然です。

タグ: Julia言語

posted at 19:06:38

Ninja DAO | CryptoNi @CryptoNlnjaNFT

21年6月5日

「東京が失敗し翌年の北京が成功した場合にはアジアにおける覇権の移動を強烈に印象づけてしまう」←今時、東京にアジアの覇権があると思ってる人なんてどこにもいないって。時代錯誤。www.sponichi.co.jp/entertainment/...

タグ:

posted at 19:03:20

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 #数楽

* 確率分布 D とデータ Y とパラメータ θ = (μ, σ) から尤度 L(D, Y, θ) の値を計算する。

* 尤度を計算するメソッド L は、D, Y, θ のどれの保有物でもない。

とするのがJulia流だと思うのですが、Distributions.jl では尤度メソッドが実質D(θ)の保有物になってしまっている。

タグ: Julia言語 数楽

posted at 18:53:02

@kankichi57301 @kankichi57301

21年6月5日

何この情報量ゼロのまとめ。 togetter.com/li/1720716#c92...

「#超算数【かけ算に順序は不要って言ってると、こんな大人になる】※発言に確認が..」togetter.com/li/1720716 にコメントしました。

タグ:

posted at 18:49:37

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 もっと正確に言うと、Distributions.jl では、尤度函数を計算するときに常にパラメータが固定された確率分布オブジェクトのコンストラクタを経由する。例えば、

L(μ, σ, y) = logpdf(Normal(μ, σ), y)

yはデータ(n=1)で、Normal(μ, σ)は正規分布オブジェクトのコンストラクタ。 twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 18:45:09

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 #数楽

横に長すぎて、一目で見にくいので、括弧を使って改行をいれてみました。 pic.twitter.com/DnXFShnH2N

タグ: Julia言語 数楽

posted at 18:35:41

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 計算効率的に優れたパッケージ群の共有という観点から見た場合には、クラスベースのOOP的な発想は捨てて、多重ディスパッチをうまく使って計算効率とコードの共有の調和を目指してくれた方が、ユーザー側としてはありがたいと思います。

タグ: Julia言語

posted at 18:27:50

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 しかし、そのせいでパフォーマンスを落として、足を引っ張っている疑いを個人的には持っています。

尤度函数が常に確率分布オブジェクトのコンストラクタを経由して計算される仕様はいかにもJulia的では__ない!__

尤度函数は常に最大速度で計算されて欲しい。

タグ: Julia言語

posted at 18:27:47

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 Juliaの基礎パッケージの中にも、ほぼそのままクラスベースのオブジェクト指向デザインになっているとみなせるものがあります。具体的には確率分布を扱う Distributions.jl のことです。

JuliaでOOPしたい人はDistributions.jlを研究すればよいと思う。

タグ: Julia言語

posted at 18:27:46

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 私もずっとabstract typeやstruct は「クラス」の直接的対応物だと誤解していました。

Juliaにはクラスなるものはないし、むしろないことが積極的にJuliaのエコシステムを支えているというようなストーリーを理解するまで相当な時間が必要だった。

タグ: Julia言語

posted at 18:27:45

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 の struct 内の field は「メンバ変数」ではないです。

abstract type や struct そのものを「クラス」の直接的対応物だと思ってしまうと失敗する。

discourse.julialang.org では定期的に話題になります(笑)

クラスベースのOOP的な考え方を捨てると楽になれます。 twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 18:27:43

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

21年6月5日

授業中に、自由に水分補給できなかったり、トイレにも自由に行けないようにしている教員が処罰されるような世の中になるとよいと思います。

教員側が禁止しているせいで、水分補給できずに子供がぶっ倒れたり、トイレに行けずに漏らしてしまう子供が出た場合には刑事罰がくだされた方が良いと思う。 twitter.com/nksmgr/status/...

タグ:

posted at 16:50:07

さのたけと @taketo1024

21年6月5日

@genkuroki こちらこそ!😆

タグ:

posted at 16:47:51

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@taketo1024 おもちゃであるがゆえに、シンプルで面白く、教育的だと思ったので、直訳にトライしてみました。もともとやっていることがシンプルなので一発で動きました!

楽しい話題、いつもありがとうございます。

タグ:

posted at 16:44:35

黒木玄 Gen Kuroki @genkuroki

21年6月5日

@aeba_chan すでにC, C++, Fortranで書き直して高速化できているなら、それと同じ程度まではJulia版も速くなるはずです。

グローバル変数をナマで使った計算を #Julia言語 でしてしまうという初歩的なミスをよく見かけます。 twitter.com/genkuroki/stat...

タグ: Julia言語

posted at 16:41:39

さのたけと @taketo1024

21年6月5日

@genkuroki ありがとうございます😆

(forAnySmall などのオモチャ実装をそのまま直訳させてしまって恐縮です😂)

タグ:

posted at 16:41:19

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 続き

④Pluto.jlを使えば(制限はあるが)、structの再定義も自由に可能になる。

⑤この辺のノウハウをまとめた文書の存在は私は知らない。しかし、discourse.julialang.org を適切に検索すれば何か出て来るだろう。現時点ではdiscourseが最良の情報源。

タグ: Julia言語

posted at 16:33:20

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

①変数名は定数名や函数名や型名やモジュール名として使えなくなる。他の組み合わせも同様。

②module My ~ end 内で定義されたexportされていないものは module まるごと自由に上書きできる。私は一時的に使うモジュール名としてOをよく使っている。

③Revise.jlも便利。

続く twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 16:33:18

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 JuliaではコードがS式になっていないので、通常のコードとそのパースの結果であるExprオブジェクトの対応の仕方を理解していないと、メタプログラミングは困難です。その学習用に作ったパッケージがこれ↓

github.com/genkuroki/Meta...

タグ: Julia言語

posted at 16:07:13

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 JuliaはLispのような完全なマクロを持つ現代的なプログラミング言語としても有名です。Juliaを

 S式を捨てたLisp(笑)

とみなす人もいるし、

 Jeff’s uncommon lisp is automated の略

という説もあるようです。

discourse.julialang.org/t/cas-benchmar...

タグ: Julia言語

posted at 16:04:32

鐘の音(除夜の鐘)ダイエット中−20kg @kanenooto7248

21年6月5日

尾身先生を見ていると、勇者が人間を見捨ててスローライフを送りたくなる気持ちがわかる。

タグ:

posted at 16:04:03

鐘の音(除夜の鐘)ダイエット中−20kg @kanenooto7248

21年6月5日

尾身先生「酒を飲むなと言ってるわけではない。短時間喋る時はマスクをして」と、感染防御と経済の両立で苦心してきた人なのだが、「そんなこと無理ゲー」でいうこと聞かず感染が広がったら「イエスマン」とバッシングを受け、挙げ句の果てには「アリバイ作り」とか、勇者でも人間見捨てるレベルである

タグ:

posted at 16:01:41

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 マクロで構文糖衣

こういう遊び方もできる。

github.com/genkuroki/publ...

も更新した。 pic.twitter.com/9IpqF9jxmI

タグ: Julia言語

posted at 16:00:00

天むす名古屋 Temmus @temmusu_n

21年6月5日

#超算数 応用問題もよく解けるとか、くもわで割合を指導すると子供に理解しやすいとかいえばいいんです。これらは反証可能性のある主張です。例えばくもわ指導をしても、割合がよくわかる子供にはなりませんtwitter.com/temmusu_n/stat...。他の主張にも、私は根拠をもって反論しています(ツイログ参照)。

タグ: 超算数

posted at 15:59:58

S (ツイートはスレッド全体をご確認く @esumii

21年6月5日

(「順序を区別しないと、単価と数量を明記せず暗黙に仮定しているときに間違える!」みたいな主張の論理的倒錯がすごい twitter.com/esumii/status/...

タグ:

posted at 15:49:48

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

SwiftからJuliaへの直訳

x -> x^2 は x を x^2 に対応させる無銘函数(ラムダ式)。

github.com/genkuroki/publ...

julia isContinuous.jl でも実行可能 → github.com/genkuroki/publ... twitter.com/taketo1024/sta... pic.twitter.com/KLoA7gJXGE

タグ: Julia言語

posted at 15:21:48

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

iterm2 の結果も

```juila-repl
julia> run(`imgcat yourimg.png`)
```

でできるんだ.ちょっと感動

タグ:

posted at 15:19:11

非公開

タグ:

posted at xx:xx:xx

非公開

タグ:

posted at xx:xx:xx

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

READMEのサンプルは Linux と Windows なら動くようです.
github.com/hyrodium/Image...

タグ:

posted at 14:27:16

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

ImageClipboard.jlパッケージの中身見て Mac でできないかなーって思って screencapture コマンドというものを初めて触って動かしてるけれど何故かマウスで囲んだところのデスクトップ背景がキャプチャーされる...

```
screencapture -i goma.png
imgcat goma.png
```

なんでやねん.
#Julia言語 pic.twitter.com/hhvWoueVeP

タグ: Julia言語

posted at 14:27:16

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

1 + 1/2 + … + 1/n



log(n) + γ + 1/(2n) - 1/(12n^2)

でさらによく近似されます。

この近似は、階乗のStirling近似の補正項 1/(12n) まで考慮したものに対応しています。

シンプルで相当に実用的な近似式になっています。 pic.twitter.com/hdbGoNezEj

タグ: Julia言語

posted at 14:05:58

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 以上に書いたことを、ガンマ函数については無限積表示を認めて使った手計算と、#Julia言語 のように気楽に使えるプログラミング言語による数値計算を行なっていれば、ガンマ函数やディガンマ函数が怖くなくなり、実用的に使うには十分な理解に達すると思います。

タグ: Julia言語 数楽

posted at 14:00:02

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

ψ(x+1) = log(x) + 1/(2x) + …

を使えば

H_n = 1 + 1/2 + … + 1/n = ψ(x+1) + γ



log(x) + γ + 1/(2x)

で近似できることがわかります。

プロットするとぴったり重なります。 pic.twitter.com/tQLJsx5r23

タグ: Julia言語

posted at 13:57:30

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 訂正版

harmonic numbers の素朴な計算と digamma 函数を使った計算を比較。

非常によく一致しています。 pic.twitter.com/idRU5BJmrx

タグ: Julia言語

posted at 13:53:51

シンカー/釣り/PT @shinka_fs_ch

21年6月5日

ワクチン接種の日本ペースが世界ペースを上回った!この調子や!
このペースなら秋には集団免疫獲得までいけるかもしれん。マジで凄い pic.twitter.com/1FVQJMymDk

タグ:

posted at 13:50:46

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 ガンマ分布がらみの統計の計算では、ガンマ函数やディガンマ函数が必須であり、ガンマ函数は階乗の連続拡張になっており、ディガンマ函数 + γ はharmonic number 1+1/2+…+1/n のnを連続的に拡張したものになっています。

数学的に基本的であることは明らかだと思います。

タグ: 数楽

posted at 13:38:08

Fabio Sodre @fabiorsodre

21年6月5日

Apesar da boa proposta de JuliaFinance para se trabalhar com dados financeiros , me caiu bem operações básicas com números fracionados e fechar com uma conversão e formatação dos resultados em string usando nativamente recursos da própria linguagem. #JuliaLang pic.twitter.com/4vZLY00iNm

タグ: JuliaLang

posted at 13:34:58

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 ガンマ函数もディガンマ函数も統計がらみの計算で必要なので、コンピュータで計算するための基本特殊函数のライブラリに含まれており、ユーザー的には容易に計算できます。

だから、harmonic number H_n = ψ(n+1) - ψ(1) も容易に計算できます。

タグ: 数楽

posted at 13:34:48

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 ガンマ函数 Γ(x) は正弦函数 sin(x) より一段基本的な函数なのですが、大学1年まで教わらないようです。

ガンマ函数が基本的なので、ディガンマ函数=ガンマ函数の対数微分も基本的なのですが、残念ながら教わる機会はほとんどないと思う。続く

タグ: Julia言語

posted at 13:31:45

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 H_n = ψ(n+1) - ψ(1) = ψ(n+1) + γ はガンマ函数 Γ(x) の無限積表示を対数微分すれば容易に示せます。Stirlingの公式

log Γ(x+1) = x log x - x + (1/2)log n + (1/2)log 2π + 1/(12n) + …

を微分すれば、ディガンマ函数 ψ(x+1)の近似公式も得られます。続く

タグ: 数楽

posted at 13:31:44

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#数楽 H_n = 1 + 1/2 + 1/3 + … + 1/n をharmonic numberと呼びます。これの数値計算は、公式

H_n = ψ(n+1) - ψ(1) = ψ(n+1) + γ

を使うとよいです。ここで、γはEuler定数 0.5772… で、ψ(x) は

ψ(x) = Γ'(x)/Γ(x) = d(log Γ(x))/dx

で定義されるディガンマ函数です。続く

タグ: 数楽

posted at 13:31:43

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

ラズパイ買ってMPIするのが社会常識だと思ってた節があります。 twitter.com/hika_ruriruri/...

タグ:

posted at 13:05:46

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 数学がらみのパッケージで、数学なので内容が古くなりようがないものは、Project.tomoの[compat]の設定を更新していないものによく出会います。

そのようなパッケージはローカルにダウンロードして展開して、Project.tomlの[compat]節を更新してから、dev で使うようにするとよいです。

タグ: Julia言語

posted at 12:53:45

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 現時点では既知の多くの試みによって「JuliaがC, C++, Fortranなみに速いこと」は確定した事実と言ってよく、疑っている人達の知識は周回遅れになっています。

しかし、遅くなった原因を調べるためのベンチマークテストは結構面白いことが多いので、みんな遊んでみるべきだと思います。

タグ: Julia言語

posted at 12:50:06

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 のPerformance Tipsの最初の節は

docs.julialang.org/en/v1/manual/p...
Avoid global variables

グローバル変数を使った計算をして失敗するのが基本パターン。

JuliaでC, C++, Fortran並の計算速度が出ていない場合にはどこかで失敗しています。

あと、BLASが行う計算速度はどこでも同じ。 twitter.com/aeba_chan/stat...

タグ: Julia言語

posted at 12:45:39

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 警告

❌pkg> add Polylogarithms

すると、SpecialFunctions.jl がv0.10.3 にバージョンダウン!

私は以下のようにしました。

pkg> dev Polylogarithms

~/.julia/dev/Polylogarithms/Project.tomlを編集
SpecialFunctions = "0.10.3" の右辺を "1" に変更

pkg> up

タグ: Julia言語

posted at 12:36:49

非公開

タグ:

posted at xx:xx:xx

@PKAnzug

21年6月5日

ところでこれ、ELECOMさんの1TBのSSDなんですけど、ここまで小さくなると取り回しも普通のUSBメモリと変わらなくなりますね。1万5千円くらいなんで、容量を考えたら特別高いわけでもない。うちは大容量の外付けストレージを多用してますけど、これはインパクトがすごいです。 pic.twitter.com/ny830I8i59

タグ:

posted at 12:17:04

積分定数 @sekibunnteisuu

21年6月5日

@OokuboTact @uKi2wQXyG7rx3gL @garapago03 @sekai_tankyu sumatome.com/su/11655587132...

 これを見ると日本のおかしな算数教育に嵌ってしまって、アジア諸国の真っ当な算数教育をdisっていることが分かります。

 ということは、掛け算順序指導は必然ではなく、日本の算数教育だとたまたまそうなっただけ、と推測できます。

タグ:

posted at 12:16:01

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 パッケージを使ったチート(ずる!)

浮動小数点数(Float64)の 1/k を足して 30 を超えるまで足し上げるより、小さなnで30を超えるみたいですね。

github.com/genkuroki/publ... twitter.com/dannchu/status... pic.twitter.com/mVouIzXW2i

タグ: Julia言語

posted at 12:15:35

Erde @Erde_st7

21年6月5日

TeX撃ってRubyとJuliaいじって疲れたのでqk

タグ:

posted at 11:31:21

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

あとifは文なんだっけ?

タグ:

posted at 11:01:00

ごまふあざらし(GomahuAzaras @MathSorcerer

21年6月5日

実質なんも書いてなくてマジゴミだった件 twitter.com/bicycle1885/st...

タグ:

posted at 10:58:28

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 うけた!

昔:1日10時間くらいで5年間

今:M1 MacBook で8時間弱 twitter.com/dannchu/status...

タグ: Julia言語

posted at 10:58:22

非公開

タグ:

posted at xx:xx:xx

Julia Bloggers @juliabloggers

21年6月5日

Why Julia – Meet Dispatch (Re-post) www.juliabloggers.com/why-julia-meet...

タグ:

posted at 10:37:20

abap34 @abap34

21年6月5日

@genkuroki これで治りました、ありがとうございます😭(まともに作業できない状況だったので、助かりました…)

タグ:

posted at 10:22:38

iwata kenichi @kenic

21年6月5日

kahan summationでやるとかなり遅く、M1 MacBook Airで7時間6分、n=6000022499693までで Σ1/n が30を超えました ^^; pic.twitter.com/cIqlkHStyV

タグ:

posted at 10:08:41

Yoshi Yonemitsu MD P @yoshipatho

21年6月5日

尾身先生を軽々しく批判する方に読んで欲しい本。

公衆衛生を適切にやるには、政府の介入が必要。だけど政治は公衆衛生に都合良く動いてくれない。

だからあえて100点を狙わず政治との対話を続け、少しでも政治が受け入れる様に60点を目指す。

この意識が、尾身先生の発想の根底にあると思う。 twitter.com/md_green/statu...

タグ:

posted at 09:48:03

とーちか @touchica

21年6月5日

これ詰むのか。
4手先の局面が脳内将棋盤では、詰みセンサーに引っかからなかった。
一間竜つえー。
「詰めチャレ」で難易度2105の問題が解けませんでした..。あなたは解けますか? #詰めチャレ #将棋クエスト pic.twitter.com/sUNc1116sC

タグ: 将棋クエスト 詰めチャレ

posted at 09:41:23

sekkai @sekkai

21年6月5日

尾身先生の「今の状況で五輪をやるのは普通ない」という科学者の立場からの発言に対して、自民党議員から「餅は餅屋だ」と不快感を表す発言があったって、自民党議員はどれだけ不遜なんだろうか。餅は餅屋なら、感染対策は当然尾身先生達分科会に任せるのが当然でしょう。 news.yahoo.co.jp/articles/caad5...

タグ:

posted at 09:16:54

CG GEEKS@CGニュースサイト @cg_geeks

21年6月5日

地球で一番いい月のデータ!!NASAがクリエイターに向けて「月周回無人衛星から撮影された"月"の3Dテクスチャ素材」を公開中!! (12位 / CG GEEKS 2020年 ベスト30) cg-geeks.com/archives/29295 pic.twitter.com/fG53gojI9C

タグ:

posted at 08:59:50

Yuki Nagai @cometscome_phys

21年6月5日

なぜ小学二年生男児は筆箱は持って帰ってくるのに中身の鉛筆は一本もないのか。鉛筆はどこに消えたのか

タグ:

posted at 08:53:35

Keno Fischer @KenoFischer

21年6月5日

I am genuinely super excited about this. Proper #julialang support is still in the works. As always if your organization buys large amounts of compute from a vendor, tell them you want better #julialang support :). twitter.com/jekbradbury/st...

タグ: julialang

posted at 08:45:54

@kankichi57301 @kankichi57301

21年6月5日

@sekibunnteisuu □や返り点を計算機が食える形に編集して「読む順を出力するプログラムを作れ」ならやる気がする。プログラミング教育の演習にもなって一石二鳥(てへ)

タグ:

posted at 07:50:10

清水 団 Dan Shimizu @dannchu

21年6月5日

やっとJuliaで計算終わった!M1 MacBook でn= 6000125006292でした。7時間46分49秒でした。もっと速くできると思うのですが,とりあえず終わってよかった! twitter.com/con_integral17... pic.twitter.com/h5MhKYFZZO

タグ:

posted at 06:37:39

積分定数 @sekibunnteisuu

21年6月5日

自分がこれで教わったらもしかしたら漢文躓かなかったかもしれない。私、何で躓いたのかあまり覚えていないけど、最初から「なんだこりゃ?」となって拒絶しちゃったような気がする。 twitter.com/hosen_school/s...

タグ:

posted at 02:50:06

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 他人が書いたパッケージの内容を理解するときには、「ほんのちょっと書き変えて、どうなるかを確認する」が結構効果的な場合があります。

devしてソースに手を入れずに調べたいこともよくあり、そういう場合には @ eval ModuleName begin ~ end が非常に便利だと思います。

タグ: Julia言語

posted at 01:44:59

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 あと、パッケージのリポジトリで「更新された最終日時」などを確認して、「健康を保っているかどうか」を確認して、ずっと更新されていないパッケージは入れない、または、dev で入れて問題が生じたら自分で手を入れるという方針にした方が無難だと思います。

dev はよく使う。

タグ: Julia言語

posted at 01:33:42

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 Juliaのパッケージマネージャーは使い易いので、最悪の場合には、~/.julia を ~/.julia-bak などにリネームしてすべてのパッケージを入れ直すという手も使えます。

実際にやってみると本当に大したことはないです。

タグ: Julia言語

posted at 01:30:42

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語

パッケージのaddやupに伴うGeneral registryの更新には時間がかかる場合があって、その途中でjuliaを強制終了させることは、パソコンのOSのアップデート中に電源を切る行為に近く、当然の災厄が降って来ます。

しかし、1つ前のツイートの処方箋を知っていると大丈夫。

タグ: Julia言語

posted at 01:27:20

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 General registry が壊れた可能性があるなら

pkg> registry rm General
pkg> registry add https ://github.com/JuliaRegistries/General
pkg> registry up General

を試してみるといいかも。httpsの直後の空白は除く。

以上は自分用のメモでもある。 twitter.com/abap34/status/...

タグ: Julia言語

posted at 01:24:12

Hideki Kawahara: WAS @hidekikawahara

21年6月5日

CAPRICEPの付属情報を使って、失敗した測定点の値を近傍から補間して修正。少し修正し残しがあるけれど見やすくなった。MATLABに用意されている実時間処理のテストツールを使って、測定プログラムをしぶとい版に修正します。DOBOT用にスピーカースタンドを用意する必要があるが、部屋があふれそう。 pic.twitter.com/N7ljiinTna

タグ:

posted at 00:52:09

台風 @taifu21

21年6月5日

@sekai_tankyu @golgo_sardine 【あれ?これって5年生の問題では?
順序解除されるんじゃなかったん?】

twitter.com/genkuroki/stat...

順序解除どころか、
教師用指導書に従うならば
たとえ6年生でも「順序間違い」は児童の「つまづき」扱いですよ?

5年の教師用指導書には何が書いてあるのでしょうね?

タグ:

posted at 00:43:26

黒木玄 Gen Kuroki @genkuroki

21年6月5日

#Julia言語 Plots.jl の中身をほんのちょっと上書きするだけで済ませたいときに、exportされていない変数や函数のすべてに Plots. を付けて書くか、1つずつ同定して import するかを考えるより、 @ eval Plots begin ~ end で囲んでしまったほうがずっと楽。

タグ: Julia言語

posted at 00:35:13

ほりたみゅ (mastodonも見てね) @Hyrodium

21年6月5日

一部ちゃんと動かないけどWindowsのサポート追加した
github.com/hyrodium/Image...

タグ:

posted at 00:32:23

左巻健男(サマキタケオ) @samakikaku

21年6月5日

次亜塩素酸水溶液普及促進会議の記者会見のアーカイブを見ると、三重大学教授の福崎氏は「もし人がいたら」ということで、実験結果を示して「人がいる空間ではこのようなうすい次亜塩素酸水では効かない」と述べた。「人間は汚染源ワーストワン」と。

タグ:

posted at 00:24:11

@genkurokiホーム
スポンサーリンク
▲ページの先頭に戻る
ツイート  タグ  ユーザー

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

掛算 統計 超算数 Julia言語 数楽 JuliaLang 十分 と教 モルグリコ 掛け算

※タグはスペースで区切ってください

送信中

送信に失敗しました

タグを編集しました