黒木玄 Gen Kuroki(@genkuroki)/2021年03月/Page 16

なるほど、プレゼンテーショングラフィックの極意ですな（意味不明）。
#ビジネス数学 twitter.com/bee_hmdx/statu...

タグ：ビジネス数学

posted at 00:07:47

#Julia言語 julia foo.jl の繰り返しの効率が悪い場合には

github.com/dmolina/Daemon...
DaemonMode.jl

を使うという選択肢もあるかも。

ANN: DaemonMode.jl, a package to run faster scripts in Julia discourse.julialang.org/t/ann-daemonmo...

タグ： Julia言語

posted at 00:48:14

朝桜さく @sacasacla

Julia日本語記事まとめ @julia_kizi

Juliaはコーディングをあまり勉強せずに雑に書いて雑に速いプログラムが欲しいというfortranユーザーのニーズにちゃんと答えているのが偉いですよね

タグ：

posted at 01:25:43

新しい記事がZennにアップされました！
zenn.dev/ligun/articles...👈

タグ：

posted at 01:35:11

#Julia言語で間違って

A = zeros(10000, 100000)

とかしてしまったときには

A = nothing
GC.gc()

でメモリを解放できます。

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

タグ： Julia言語

posted at 02:50:37

Kaz Muzik @kazmuzik

Hideki Kawahara: WAS @hidekikawahara

【17桁】17桁の第1種カニンガム鎖の処理は 2.73e16あたりで、Mac Pro 2009から、EPYCワークステーションへ移行。今回は JULIA_NUM_THREADS=30として、@ threadsのforループでは10倍の300を使ってみた。Juliaがうまくスレッドにdispatchしてくれている感じ。これで物理的なCPUを意識せずにCPUを有効に‥ pic.twitter.com/9Tpugo4DgK

タグ：

posted at 02:55:55

MATLAB R2021aをインストール。とりあえず動作することを確認した。まだ幾つか締切が残っているので，しばらく封印しておく。と言いながらjulia+Plutoで遊ぼうとしている。自己矛盾。

タグ：

posted at 04:18:54

barsalmocto @barsalmocto

ケフっち＠Deterministic @__cheph__

「詐欺グラフ」は最近「ビジネス数学」という立派な名前をもらって出世した。「詐欺」なんて言葉はポリコレじゃないから、肯定的な表現に直されたんだね twitter.com/bee_hmdx/statu...

タグ：

posted at 06:48:19

(｢・ω・)｢ｶﾞｵｰ @bicycle1885

うちの元ボスのプログラミングの講義、Fortranだから今度Juliaを布教しておこう。最近、なんかディープラーニングとかいう言葉を覚えたらしいし。笑

タグ：

posted at 07:17:01

富岳でJuliaが動く（予定）➝ Juliaでは俺のコードも動いてる ➝ 富岳で俺のコードが動く！

タグ：

posted at 07:38:28

Takeshi Nishikawa @tnishikawa

アンポンタン・ポカン博士 @tikiwa3756

@kaityo256 @gachacomplete 富岳クラウドプラットホームでJuliaが使えるようにして、若手が気軽に使えるようにします。

タグ：

posted at 07:40:24

このためにプログラムを書くこともあるけれど、あくまでも「モデル立てとシミュレーション、結果の評価」をすることが目的なので、あらかじめコードは配布されるし任意の言語を用いて計算しても良い。(僕はJuliaを使っていた)課題内容に至るまでプログラミング能力は殆ど問われない設計になっていた。

タグ：

posted at 09:50:16

#Julia言語

スーパーコンピュータの富岳に2021年4月までにJuliaがインストールされる予定らしい。

www.hpci-office.jp/pages/hardware... pic.twitter.com/RvHSO9v6My

タグ： Julia言語

posted at 10:23:24

タグ： Julia言語

posted at 10:26:31

NASAでは普通に #Julia言語を使っているようですが、日本ではどうなっているんですかね？

Juliaの良い点をNASAがどのように活かしているか分かる動画
↓
youtu.be/tQpqsmwlfY0
Modeling Spacecraft Separation Dynamics in Julia
Jonathan Diegelman twitter.com/choonggeun/sta...

タグ： Julia言語

posted at 10:44:10

#Julia言語あらためてNASA JuliaLangでググったら、こんな素敵な動画も見つかりました。

youtu.be/QdGrR0inYj8
The Julia Language 1.0 Ephemeris and Physical Constants Reader for Solar System Bodies

github.com/JuliaAstro/JPL...

タグ： Julia言語

posted at 10:53:07

#Julia言語 NASA関連

Pluto.jl
DifferentialEquations.jl
Optim.jl
GlobalSensitivity.jl
MonteCarloMeaurements.jl
ComponentArrays.jl
Unitful.jl

などを使っているみたい。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 10:57:32

#Julia言語パラメータの値をMonteCarloMeasurements.jlを使って「幅を持つ数値」で与えると、DifferentialEquations.jlはその幅を持ったパラメータ値について微分方程式を数値的に解いてくれて、Plots.jlはそれを視覚化してくれる。

幅を持つパラメータ値を扱うためのコードを自分で書く必要がない！ twitter.com/genkuroki/stat... pic.twitter.com/Ha2TkLCkfp

タグ： Julia言語

posted at 11:07:51

#Julia言語例えば微分方程式のソルバがCやFortranで書かれていて、単なる浮動小数点数にしか対応していない場合には、「幅を持つ数値」に対応するためのコードを自分で書く必要が生じる。

タグ： Julia言語

posted at 11:14:02

#Julia言語しかし、JuliaのDifferentialEquations.jlは「幅を持つ数値」を含む浮動小数点数よりも圧倒的に広い数値の型に対応している。

だから、入力のパラメータを浮動小数点数ではなく、「幅を持った数値」にすれば「幅を持った数値」について自動的に微分方程式の数値解を計算してくれます。

タグ： Julia言語

posted at 11:16:27

#Julia言語しかもその手の「魔法」は、特別に優れた人が開発したDifferentialEquations.jlのようなパッケージだけでがなく、私が書いた素朴な微分方程式の数値解を求める函数でも使えます！

Juliaが普及すれば「魔法」が大衆化されまくることになります。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 11:20:03

#Julia言語 Juliaでforループを回す素朴なコードを書けば、gccと同じ速さで計算してくれることを確認できます。

だから、gccの側で時代遅れのライブラリを使った途端に、gccは新しい無難なライブラリをデフォルトで使ってくれるJuliaに計算速度で負ける。

コードを書くのはJuliaの方が圧倒的に楽。

タグ： Julia言語

posted at 11:25:01

#Julia言語

Juliaなら気楽に書けて、gccと同じ速さで計算してくれて、ライブラリの選択で失敗しがちな人はgccではなくJuliaを使った方が無難に速く計算できる可能性が高い。

しかもJuliaなら、自分で書いた微分方程式のソルバに「幅を持った数値」をパラメータとして与える「魔法」も使える。

タグ： Julia言語

posted at 11:28:46

#Julia言語 NASAがJuliaを使うようになったことは技術的な合理性の観点からは当たり前のことでしかないと思います。

タグ： Julia言語

posted at 11:30:22

裳華房編集部 @shokabo_editors

【新刊のご案内】
もう書店に並んでいるかと思います。
見やすい新装版、お手に取っていただければ幸いです。

堀田良之『代数入門（新装版）－群と加群－』
www.shokabo.co.jp/mybooks/ISBN97... pic.twitter.com/rlnjb8nxfF

タグ：

posted at 11:34:54

#Julia言語例えば、gccの側でデフォルトのrand()を使うとJuliaに計算速度で負けます。

gcc側でメルセンヌツイスターMT19937を使ってもその後継のdSFMTがデフォルトのJuliaに計算速度で負けることになります。

gcc側でライブラリを自分で選ぶと容易に失敗できます。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 11:37:58

#Julia言語添付画像のように、シンプレクティック法の1つであるYoshida6()法で微分方程式の数値解を求め、さらにパラメータに幅を持たせた場合の挙動を視覚化することを、空行を除いて9行で実行できています。

速いだけではなく、非常に楽ちん。 pic.twitter.com/FklYwwTQLh

タグ： Julia言語

posted at 11:41:02

Lirimy @LirimyDh

@abap34 これでどうでしょうか？
docs.julialang.org/en/v1/manual/m...

タグ：

posted at 11:45:13

#Julia言語 JuliaでDifferentialEquations.jlとPlots.jlを使うと以下の仕事の仕方は時代遅れに感じられる：

❌gccやg++で微分方程式の数値解を求めるコードを書いて、計算結果はファイルに出力し、視覚化はそのファイル経由で行う。

パラメータを摂動した場合の挙動の視覚化は非常に面倒になる。

タグ： Julia言語

posted at 11:48:31

今の日本の大学教育では

❌gccやg++で微分方程式の数値解を求めるコードを書いて、計算結果はファイルに出力し、視覚化はそのファイル経由で行う。

という時代遅れな仕事の仕方が普通だと思ったまま学生を卒業させている場合が多いのではないか？

そういう学生は現代のNASAでは通用しなさそう。

タグ：

posted at 11:50:57

以下のリンク先の動画との比較:

youtu.be/QdGrR0inYj8
JuliaCon 2019 | The Julia Language Ephemeris and Physical Constants Reader for Solar System Bodies

タグ：

posted at 11:52:32

Tak itoshi @takitoshi

「Juliaで精度保証付き数値計算」更新しました。

非線形方程式の精度保証付き数値解法
taklab-blog.blogspot.com/2021/01/rigoro...

近年の区間解析の集大成ともいえる非線形方程式の解の精度保証。自動微分を使うと問題を記述するだけで、精度保証された解が出てきます。

タグ：

posted at 12:25:51

Tak itoshi @takitoshi

超越基底(佐藤陽花/エヤイヌニタㇰ/藍徽 @Distr_to_Yonder

Krawcyzk法、区間Newton法を複数回使い、与えられた範囲内の全ての解を得る。分枝限定法(branch and bound/branch and prune algorithm)、別名全解探索アルゴリズムが本当に強力なので、今後紹介したい。

タグ：

posted at 12:31:17

JuliaとDartの数値計算・機械学習ライブラリの整備に力を注いで行きたい。

(本命は前者。後者はアプリやWebシステムへの応用のため)

タグ：

posted at 12:37:11

Juliaで自動微分はどうなってるのかなと思って検索したら、うちの降旗さんのブログが出てきた。さすがはJuliaの伝道師である

タグ：

posted at 12:50:23

#Julia言語

using Statistics, StatsBase
f(x) = (mean(x), geomean(x), median(x))
[∘(fill(f, n)...)((1, 1, 2, 3, 5)) for n in 30:35]
↓
6-element Vector{～}:
～
(2.0890579497368598, 2.0890579497368598, 2.0890579497368598)
～

gist.github.com/genkuroki/e69a...

xkcd.com/2435/ twitter.com/genkuroki/stat... pic.twitter.com/i3SovM4VbR

タグ： Julia言語

posted at 13:37:12

#Julia言語

∘(fill(f, 4)...)((1, 1, 2, 3, 5))

は

f∘f∘f∘f(1, 1, 2, 3, 5)

と同じで、∘は函数の合成なので、

f(f(f(f(1, 1, 2, 3, 5))))

と同じ意味になる。

函数fのn個の合成は

∘(fill(f, n)...)

と書ける。Juliaを超高級電卓として使うときにはこれを知っているとちょっと便利。

タグ： Julia言語

posted at 13:37:13

#Julia言語実際の例は7行をちょっと超えている場合もあるのですが、以下のリンク先を見れば短い行でどれだけのことをできるかがよく分かります。

知らなかったパッケージも多数使われていた。

Seven Lines of Julia discourse.julialang.org/t/seven-lines-...

タグ： Julia言語

posted at 13:48:43

こなみひでお @konamih

200ページのJPEGデータを一発でPDFにまとめる処理を書いていて，いきなり「ImageMagickのキャッシュが限界を超えました」で落ちる悲劇に遭遇。端末からコマンドで処理したほうが早かったなあ。

タグ：

posted at 13:54:15

#Julia言語による解答例のスクショ

もとの問題文の数学的内容がほとんどそのままJuliaのコードに翻訳されていることがわかる。

難しいのはn個の函数fの合成が ∘(fill(f, n)...) と書けること。

φ(fill(f, 3)...)はφ(f, f, f) と同じ意味に、∘(f, g, h)(x,y)はf(g(h(x,y)))と同じ意味になる。 pic.twitter.com/LUR4w33cwC

タグ： Julia言語

posted at 14:00:40

#Julia言語 fの定義は

f(x) = (mean(x), geomean(x), median(x)) # Tuple

ではなく

f(x) = [mean(x), geomean(x), median(x)] # Array

でもよいが、一般に、Arrayではなく、Tupleで計算可能ならば、Tupleを使った方が計算が速くなる。

Tupleを配列のように扱うためのパッケージがStaticArrays.jl

タグ： Julia言語

posted at 14:05:15

#Julia言語巨大な行列やベクトルの数値線形代数はBLASに任せることが基本になるのですが、常微分方程式などでよく現われる小さなサイズのベクトルや行列の取り扱いについては、StaticArrays.jlの使用を検討すると良いです。これはNTupleを抽象配列化するパッケージ。

github.com/JuliaArrays/St...

タグ： Julia言語

posted at 14:11:49

#Julia言語タプルで計算できれば速いと知れば、

(1, 2, 3) + (4, 5, 6)
↓
(5, 7, 9)

のような計算もして欲しくなるのですが、タプルには足し算などが定義されていません。

StaticArrays.jlではSArray型でタプルを包んで、線形代数の演算をできるようにしています。

github.com/JuliaArrays/St...

タグ： Julia言語

posted at 14:15:20

l_ppp @ppp3141592ppp

すげー
本当にCC13、しかも17桁の方から

だいぶ前の段階で諦めたやつ
twitter.com/ppp3141592ppp/...

#カニンガム鎖 twitter.com/kazmuzik/statu...

タグ：カニンガム鎖

posted at 14:43:22

kashi @mkashi

お、ワイの講義資料が引用されておる。そしてワイの考えた謎の候補者集合の作り方も使われてる。候補者集合の最適な作り方はほんと分からん。 twitter.com/takitoshi/stat...

タグ：

posted at 15:03:41

yudai.jl @physics303

One Note の正しい使い方． pic.twitter.com/bN4dONtWhy

タグ：

posted at 15:19:12

非公開

タグ：

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

非公開

タグ：

posted at xx:xx:xx

@konamih 私が使っている板書の写真群(IMG_番号.jpg型のファイル名の形式を仮定)を2x4の画像ファイルにまとめてから、pdf化する #Julia言語スクリプトです。

-geometryで解像度を変更してサイズの爆発を防いでいます。

参考になれば幸いです。

gist.github.com/genkuroki/fe7a... twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 15:56:05

非公開

タグ：

posted at xx:xx:xx

@kikumaco #Julia言語さすがにJuliaでの自動微分でForwardDiff.jlのみを挙げるのは時代遅れな感じです。NASAの人はZygote.jlも例に挙げていました。

現在のJulia界では「自動微分よりもさらにずっと先」のことを考えるべきだという話にすでになっています↓

www.stochasticlifestyle.com/generalizing-a...

タグ： Julia言語

posted at 16:50:09

@genkuroki 5年前の記事ですから

タグ：

posted at 16:50:50

@kikumaco とはいうものの、ForwardDiff.jlは非常に便利です。

(f(x + 1e-10) - f(x - 1e-10))/(2*1e-10) とかするのはダサいし、色々悪いことが起きます。 #Julia言語

タグ： Julia言語

posted at 16:56:52

「同一のコードを別の解釈でコンパイルする」という一般論があり得ることは自明。

しかし、数値計算でその一般化が具体的にどう役に立つかは非自明。

以下は #Julia言語界から出て来たそういう方向の一般化の技術的有用性の解説。自動微分はその特殊な場合になる。

www.stochasticlifestyle.com/generalizing-a... twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 18:00:07

児玉龍彦せんせーは新型コロナ云々を言う前に、まず「チェルノブイリ膀胱炎デマ」の後始末をちゃんとするべきです

タグ：

posted at 18:09:25

でえもん @GreatDemon1701

帷子ノ辻しろっこー @ultimatile

ビジネス数学
というのは長いので、
B数学
BMath
あたりでハッシュタグ作りませんか？

タグ：

posted at 18:21:05

NumbaとCythonで迷ったらJulia(?)

タグ：

posted at 18:40:47

あおじるPPPP @kale_aojiru

受験や試験の点数、もっとはっきり言えば設問に正解するために勉強するわけじゃないよねえ。といいつつ、塾というのはそういうものを学ぶところと言えばそうなのかもしれないが、学校が「受験や試験の点数、もっとはっきり言えば設問に正解するための勉強」を率先して教えている現状、じゃあどこならば

タグ：

posted at 18:46:27

あおじるPPPP @kale_aojiru

「受験や試験の点数、もっとはっきり言えば設問に正解するため」ではない勉強を教われるのか、という問題が起きる、というか起きているのよねえ

タグ：

posted at 18:46:28

あおじるPPPP @kale_aojiru

学習塾に対する需要の多くがそういうものであるとは到底思えないが、まったくないわけでもないと思われるので、そういう方針で生き残ることは不可能ではないだろう

タグ：

posted at 18:47:43

こいなぎ @naginyan135

【級位者向け実戦型詰将棋part481(松田茂行先生作)】
これを詰ますことができた者だけが寝ることができる。 pic.twitter.com/lOCZfDkLnZ

タグ：

posted at 20:00:00

yudai.jl @physics303

見つけなくて良いと思います。

#Peing #質問箱 peing.net/ja/qs/986659573

タグ： Peing 質問箱

posted at 20:10:04

place where cat shou @catshouldnt

pic.twitter.com/2FSuOCLdfz

タグ：

posted at 20:52:04

#超算数 #移項る上で解説したものの次、数学教育現代化の方針をもつ指導要領に際して作られた指導書
文部省『中学校指導書: 数学編』MEJ 3236 大阪、大阪書籍、1970年。
等式の性質のうち加と減を【統合的に考えること】(乗除も同じ)が推奨される。同じものに同じ作用を加えるとはいかないが。 pic.twitter.com/Ma04ILjGdf

タグ：移項る超算数

posted at 20:54:37

#超算数 #移項る前回の指導要領に用語として記載された【移項】は今回の指導要領には登場しないerid.nier.go.jp/files/COFS/s44...。指導要領を解説する指導書にも登場しない。【さらに、方程式を形式的に能率よく解くことができるようにする。】という記述は移項を思わせるが(47)。【項】は前節の用語にあり。

タグ：移項る超算数

posted at 20:54:39

#超算数この指導書（文部省時代の学習指導要領解説のこと。教科書に付随する教師用指導書ではない）における偶数、倍数の定義。偶数とは【2で割り切れる】数で、2の倍数は明示的な定義がない。ただし例示は偶数と同じになっている。さらにこれらは集合の単元の学習事項で、全体集合は(0抜き)自然数。 pic.twitter.com/Jr5ruZWpOf

タグ：超算数

posted at 20:54:41

#超算数ちなみにある数の約数にはその数自体と1が含まれている(以上127ページ)。この指導書の中では、倍数はそれ自体が学習の対象になっているというよりは、集合の学習のための題材となっている。集合を論じる範囲が自然数(この当時も高校までは0を含まないようだ)なので、0は最初から考慮の外。

タグ：超算数

posted at 20:54:43

#超算数換言すると、0がある数の倍数であるかどうかを考える必要がないのだ。この状況が意図的に作り出されたかどうかはわからないが、結果として全指導要領下の指導書との継続性が保たれることになった。倍数は、59年版指導書では分数計算の道具で、70年版では集合学習の題材であったといえる。

タグ：超算数

posted at 20:54:44

#超算数 1969年指導要領とそれを解説する70年指導書は、数学教育の現代化を目指したとされる(教科調査官は前回から引き続き大野清四郎)。現代化の実態は、集合に力点が置かれたものであった。様々な経緯から数学教育に残る不具合を包括的に修正するものではなかったと思われる。

タグ：超算数

posted at 20:54:45

#超算数 1951年指導要領までの占領下指導要領では、分数の乗除が中学校まで持ち越されるなど、戦前および現代とくらべて算数・数学の学習進度が遅いという特徴があった。58年指導要領、59年指導書(中)、60年指導書(小)は、進度復旧を目したものだが、分数計算と倍数を絡めて学習すると定めてしまった。

タグ：超算数

posted at 20:54:46

#超算数次の指導要領、指導書では倍数が集合の学習と絡められたが、先回に課せられた制約を撤去するには至らなかった。分数計算のための倍数概念ならば、倍数に0がはいるかどうかは、問う必要がないという意味。戦前の中学校で0を含む倍数が扱われた形跡はあるtwitter.com/temmusu_n/stat...。

タグ：超算数

posted at 20:54:48

#超算数しかし高等小学校をふくめ、小学校で倍数がどう指導されたか、そもそも指導されたかは未調査。いまのところ、算数では倍数に0を含まない指導の遠因は、51年指導要領までの遅進度を是正しようとした58年指導要領にあると考えている。

タグ：超算数

posted at 20:54:49

#Julia言語 Julia言語で10行

nbviewer.jupyter.org/gist/genkuroki...

(1/5) 二項分布モデル内の標本分布で測った95%信頼区間にモデルのパラメーター値が含まれる確率のプロット。95%以上になる。ほとんどのパラメーターで95%より真に大きい。 twitter.com/genkuroki/stat... pic.twitter.com/c7S3KQKZ8f

タグ： Julia言語

posted at 20:59:00

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...

(2/5) 二項分布モデルの95%信頼区間 [CI_min, CI_max] を平面上の座標 (CI_min, CI_max) にプロット。丸の大きさはモデル内でその信頼区間が生じる確率の大きさに比例。赤の十字の左上側の領域ではパラメータの真の値が信頼区間に含まれている。 pic.twitter.com/kXWPTdALRC

タグ： Julia言語

posted at 21:01:39

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...

(3/5) 3次元のランダムウォーク pic.twitter.com/EedPFAjjsZ

タグ： Julia言語

posted at 21:02:10

#Julia言語続き。動画。

こういう動画を10行で作れる。 pic.twitter.com/1Fc5Gix77z

タグ： Julia言語

posted at 21:03:49

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...

(4/5) 2次元のIsingモデル。これは自作のパッケージ(公式登録済み)を使ったチート。 pic.twitter.com/mcUmqebaIg

タグ： Julia言語

posted at 21:04:41

#Julia言語動画

github.com/genkuroki/Isin...
Ising2D.jl pic.twitter.com/bbnXNIS9wH

タグ： Julia言語

posted at 21:05:53

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...

(5/5) Kuramoto model! 初期条件も周期も異なる振動子が互いの相互作用によって「ほぼ同期」するようになる。繰り返しになるが、各枡目に配置された振動子の周期は互いに異なる。それでも同期する方向の相互作用が十分に強ければ同期するようになる。 pic.twitter.com/LJuQFFRgh1

タグ： Julia言語

posted at 21:08:26

#Julia言語続き。Kuramoto model の動画。

こういう動画も10行で作れる。

訂正：1つ前のツイートの添付画像の Printf を削除。それは余計。

using Distributions, DifferentialEquations, Plots の3つのパッケージを使った自然な10行のコードでKuramoto modelの動画を作れる！ pic.twitter.com/62MS6rjaXh

タグ： Julia言語

posted at 21:11:14

#Julia言語あ、動画中の3-dimensional の綴りもおかしい！もとのソースコードどその結果の方

nbviewer.jupyter.org/gist/genkuroki...

は今から直しておきます。

タグ： Julia言語

posted at 21:13:01

#Julia言語視力の問題でこういうことになっている。こういうミスは今後増える一方になると思う。

タグ： Julia言語

posted at 21:14:13

ちょっと必要に迫られて、モンテカルロ法のチュートリアルを書いています。基礎編(9ページしかない)をざっくりと書いたので、とりあえず公開しておきます。タイトルにある「レアイベントサンプリング」は今後追加される予定。
PDFなので注意してください
www.cp.cmc.osaka-u.ac.jp/~kikuchi/texts...

タグ：

posted at 21:16:31

#Julia言語以下のリンク先のプロットは初めて作った。他の人で同様のグラフを作った人がいるかどうかは知らない。

信頼区間を平面上の点(丸)としてプロットしている。丸の大きさはモデル内標本分布で生じる確率の大きさに比例。

個人的にこれを見て少し理解が深まった。

twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 21:16:59

#Julia言語

 nbviewer.jupyter.org/gist/genkuroki...

(3/5) 3次元ランダムウォークの動画を作り直した。

ランダムウォーク自体は1行で作れます。

それの動画を作るコードの方がずっと長くなる。

それでも全部でぴったり10行におさまった。 pic.twitter.com/dRMzL6cTzr

タグ： Julia言語

posted at 21:19:32

#Julia言語 Kuramoto model の動画を再掲。

マス目の個数は16×32=512個ある。

個数を4倍にした32×64=2048個の場合も楽勝で計算できたが、動画のサイズも増えたので、この程度にした。

ソースコード↓ (10行しかありません)
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/SAQdiZDsnp

タグ： Julia言語

posted at 21:22:40

Dr. Chris Rackauckas @ChrisRackauckas

#超算数この本は【青年】に向けて書かれている。高等小学校または中学校の生徒を指し、尋常小学校(今の小学校に対応)は埒外かも。
山内太一『詳解算術: 考へ方・解き方』改訂東京、淡海堂出版部、1933年。
一応、戦前も、分数計算とともに約数や倍数の概念を教えることになっていたようだ。 pic.twitter.com/avRZUifFd3

タグ：超算数

posted at 21:26:07

100x accelerated uncertainty quantification of dynamical systems with Koopman Expectations? Come see how to use it with Pumas at the first Pumas 2.0 webinar at 12pm EST. Inviting the whole #pharma, #julialang, and #sciml community.

form.jotform.com/210605786318054

タグ： julialang pharma sciml

posted at 21:27:18

Massimo @Rainmaker1973

Watch people play ‘Pong’ on a rock wall with Augmented Climbing [source: ow.ly/E6Ug50Am51h] [more: valomotion.com] pic.twitter.com/EXNEAaybhC

タグ：

posted at 21:30:25

#Julia言語以下のリンク先は空行を除いて9行

微分方程式のソルバはシンプレクティック法の1つであるYoshida6()で、微分方程式のパラメータに幅を持たせた場合の影響を視覚化しています。

こういうのも10行必要ない。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 21:32:17

JuliaHub @JuliaHub_Inc

Julia 1.6: What has changed since Julia 1.0? In this article Lyndon White (@oxinabox_frames) lists some of the new features and changes that will be seen in the latest release. To read, click here www.oxinabox.net/2021/02/13/Jul...
#julialang pic.twitter.com/RgwIcJcoaQ

タグ： julialang

posted at 21:34:33

福島の甲状腺検査は「安心のため」って言うんだけど、がんが見つかってしまった250人は安心どころではないだろうし、B判定で安心できるかどうかは微妙だし、人によってはA2だって不安かもしれない。

よしんば30万人が安心できたとしても、その代償が200人の手術って、倫理的に許されないと思いますよ

タグ：

posted at 21:39:28

積分定数 @sekibunnteisuu

需要は少なくても、供給も少ないなら隙間産業としてなりたつかも twitter.com/kale_aojiru/st...

タグ：

posted at 21:45:57

ロム ⭐️⚡️❄️ @ROM31415

@shimeji_mattu 広い画ですが、スポンジで麻雀をした猛者もいました pic.twitter.com/XICMRcuGKW

タグ：

posted at 22:04:58

うどん伯爵 @curryking2214