黒木玄 Gen Kuroki(@genkuroki)/2021年03月21日

私が大学1年で受けたコンピュータの教育はFORTRAN 77で、最初の講義でホーナー法の解説があった。

FORTRAN 77でも教えられることばかりだとつまらないので、数値計算におけるコードの自動生成の重要性なんかも教育に組み入れると面白いと思う。

#Julia言語なら完全なマクロがある。 twitter.com/Atsushi_twi/st...

タグ： Julia言語

posted at 00:20:28

@genkuroki 例のランダムウォークに関して突っ込まれたとき批判に真摯に対応しなかったことを理由にT氏をhankagosaさんがブロックしたこと，数か月前からやたら主観ベイズの某氏とT氏がべったりなことが効いてそうですねえ．W先生のお弟子さんもT氏にブロックされているようですし断絶が
twitter.com/-/status/13542...

タグ：

posted at 00:27:16

@genkuroki hankagosaさんの警鐘
twitter.com/-/status/12589...

タグ：

posted at 00:29:27

@genkuroki 無用な通知がhankagosaさんやnhayashiさんに飛ばないようなリンクになっています。
スマホアプリなら大丈夫ですが、Webブラウザを使ってる場合は-を_に置き換えることで見られます

タグ：

posted at 00:32:06

#数楽このスレッドを読んだ人の中には、自分でも sin(x) を近似する多項式を求めたい人がいると思う。

Float64に合わせて係数を作る方法の(試行錯誤そのもののナマの整理されていない)コードが、「低次の係数からFloat64での値に順次固定して求めてみる」の節にあります。

gist.github.com/genkuroki/362b...

タグ：数楽

posted at 00:36:46

#Julia言語 Juliaでも使われているFloat64でのsin(x), |x|≤π/4をよく近似する13次多項式と同等の13次多項式をJuliaを使って自力で求めてみる話が以下のリンク先スレッドにあります。

Float64に特化した係数の計算は繊細な注意が必要な問題だった。

興味がある人は自分でやってみるとよいと思う。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 00:42:31

#Julia言語私が求めた13次多項式の係数と伝統的に使われている13次多項式の係数は微小に違うのですが、相対誤差をほぼ同じです。グラフ的には区別ができない！

謎のレシピを謎のまま残さずに自力で求められるようにしておくことが、数学がらみの事柄では基本になります。

gist.github.com/genkuroki/362b... pic.twitter.com/oowxFdjsSV

タグ： Julia言語

posted at 00:49:10

#Julia言語 Float64でsin(x), |x|≤π/4を+と*だけで(実際にはmuladd, fmaを使う)計算する方法は(少なくとも私は)分かった。

Float128でもやってみて、Quadmath.jlのsinと比較してみるのは面白い問題だと思います。

こういう活きた問題が時間をかけて解けば多くのことが身に付くはずです。

タグ： Julia言語

posted at 00:54:29

訂正

訂正前：こういう活きた問題が時間をかけて解けば多くのことが身に付くはずです。

訂正後：こういう活きた問題が結構沢山あるので、適当に選んで時間をかけて解けば多くのことが身に付くはずです。

疑問1が解決すると派生の問題が必ず生じる。

タグ：

posted at 01:03:01

Kaz Muzik @kazmuzik

メタリックはんぺん @metalichanpen

DataFrames.jl github.com/JuliaData/Data... などはまだ使っていないけど、Python + Pandasから移行出来るかも‥ twitter.com/chor0sK/status... twitter.com/kamayan1128/st...

タグ：

posted at 01:25:10

地理院の標高データとpythonでとりあえずこういうことができた(白黒逆だったので再投稿) pic.twitter.com/YFZvjUlzvT

タグ：

posted at 02:09:03

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

「FastGaussQuadrature.jlで数値積分しましょう」を公開しました！！死ぬほど分かりやすくGauss求積の解説を書いたので死なないように気をつけてください…
zenn.dev/hyrodium/artic...

タグ：

posted at 02:10:54

@metalichanpen Juliaでできるかな...

タグ：

posted at 02:12:42

zennの記事は4本目で、文字数が単調増加してます pic.twitter.com/5Jt8SuVQcB

タグ：

posted at 02:12:51

高梨陣平 @jingbay

KikkaさんによるDMCAの虚偽申請に関わる記事。

読んでると救いの無さにうんざりする。 twitter.com/itm_nlab/statu...

タグ：

posted at 02:34:47

出目金 @TR_727

そして…

「千葉県のシンボル/は「カタツムリ」です。八木秀次氏が〜『正論』に詳しく書いてますが、これはカタツムリが雌雄同体であることを取り上げたもの〜」

ケント・ギルバート『やっと自虐史観のアホらしさに気づいた日本人』2016
twitter.com/tr_727/status/... pic.twitter.com/tHexP5ZoWZ

タグ：

posted at 02:34:53

Kaz Muzik @kazmuzik

【19桁】第1種カニンガム鎖、1.255e18までの処理が完了。Githubのリストを更新。raw.githubusercontent.com/kazkus/Cunning... CC13は55番目の 1231790682838960529 までが正式。CC14は3番目の 1092966255671664659 までが正式。 pic.twitter.com/IdJu5MhaYY

タグ：

posted at 02:35:06

#Julia言語 10行シリーズにマンデルブロ集合を追加した。

ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/JdwKqbNPZA

タグ： Julia言語

posted at 02:37:21

#Julia言語こういう動画を10行で作れます。

何百行もあるコードだと解読する気になれなくても、たったの10行ならば解読する気になれるかも。実際に解読できれば、こういうことを自分で自由にできるようになる。

ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/GIKHD73XsA

タグ： Julia言語

posted at 02:37:24

Julia Bloggers @juliabloggers

New post: Expand your DataFrames.jl toolbox: the flatten function - www.juliabloggers.com/expand-your-da... #julialang

タグ： julialang

posted at 02:45:18

Astellon @astellon_music

RustのサンプルだったものをCrystalで書いてみた。
pic.twitter.com/G0Hl2c9DwD

タグ：

posted at 02:54:25

Astellon @astellon_music

Julia日本語記事まとめ @julia_kizi

endler.dev/2018/fastcat/

タグ：

posted at 02:58:04

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

タグ：

posted at 03:02:40

Mark Kittisopikul ht @markkitti

OpticSim.jl: Optical Simulation software in #JuliaLang from @microsoft

github.com/microsoft/Opti...

タグ： JuliaLang

posted at 03:31:11

Mark Kittisopikul ht @markkitti

From the README: pic.twitter.com/R4j6JX1xX0

タグ：

posted at 03:31:11

非公開

タグ：

posted at xx:xx:xx

Hiroyasu Kamo @kamo_hiroyasu

@genkuroki @tsatie これ、世界を普通のトーラスでなくねじれトーラスにしたら世界がグライダーまみれになって楽しそうです。単にアスペクト比の変更でも可。

タグ：

posted at 06:07:40

#Julia言語 Mandelbrot set animation 作り直し pic.twitter.com/BjvxPjhnoo

タグ： Julia言語

posted at 06:40:56

#Julia言語作り直し2 (こちらが本命、拡大率を上げた)

これももちろん10行。 pic.twitter.com/w47s7ViI5i

タグ： Julia言語

posted at 06:41:48

#Julia言語 Game of Life

これも10行 pic.twitter.com/aJ1i2TzXS0

タグ： Julia言語

posted at 06:42:45

madfish @madfish19

凄い🤗 twitter.com/genkuroki/stat...

タグ：

posted at 06:46:16

#Julia言語 10行でライフゲーム

初期条件は5×5
2000ステップのアニメーション pic.twitter.com/b59pxDHfwG

タグ： Julia言語

posted at 07:19:09

#Julia言語

10行シリーズ全ソースコード↓
nbviewer.jupyter.org/gist/genkuroki...

添付画像は10行ライフゲームのコード

初期条件やプログレスメーターやJuliaの文法で必要なendも含めて10行！

可読性も悪くないはず。 pic.twitter.com/5yjacEYPYe

タグ： Julia言語

posted at 07:24:05

#Julia言語配列uに保持した状態を更新するコードを

for t in T
u = update(u)
end

の形式で行うと無駄なメモリアロケーションの嵐になる。

for t in T
update!(u)
end

のようにuの内容を更新するように書かないと損をする。

以下のライフゲームのコードもそうなっています。 pic.twitter.com/TjO2IbIeNy

タグ： Julia言語

posted at 07:27:07

#Julia言語配列に保持した状態を時間発展させるコードはありがちなので、ライフゲームのような遊びのプログラム作成もバカにできない。10行の中に基本が詰まっている。

10行の中にプログレスメーターを入れる余裕があるとちょっとうれしい。 pic.twitter.com/zOsM2JI69Q

タグ： Julia言語

posted at 07:30:26

@VC31415 I’d add #julia for programming. It’s faster/more flexible/closer to the math. github.com/azev77/QuantEc...

タグ： julia

posted at 08:21:29

@MapleTaffyCap @VC31415 Julia job postings are growing. Besides, Julia teaches better programming habits, so many employers who list “Python” will take Julia experience seriously. Stanford engineering classes use Julia from beginning, their grads have no problems finding jobs...

タグ：

posted at 08:52:22

eiko @EikoDonald

朝からガサガサしてるなぁと思ったら。Juliaさん脱皮っピ。なぜうちの子はイベント行く時とか脱皮するん？ちゃんとキレイに脱げました。
#レオパ #レオパードゲッコー #ヒョウモントカゲモドキ pic.twitter.com/SrlChwYHze

タグ：ヒョウモントカゲモドキレオパレオパードゲッコー

posted at 09:00:41

@MapleTaffyCap @VC31415 Agreed, but the fact that Stanford/MIT are increasingly using Julia in their curriculum gives an idea about trajectory. Good luck

タグ：

posted at 09:00:51

にるそん @ohtanilson

マイTips：複数言語でコードを書いてそれぞれからfigureとかtableを出すとき、wd直下データフォルダは共通にして、main_Rとmain_juliaとmain_pythonにメインのコードを入れて、julia_figuretableとR_figuretableのように言語ごとにフォルダ分けて出力を保存したほうがLatexに貼るとき混乱しない。

タグ：

posted at 09:20:39

@VC31415 @MapleTaffyCap Paul has some awesome #Julia packages, most of his recent stuff on Git is in jl: github.com/schrimpf?tab=r...

タグ： Julia

posted at 10:25:16

@VC31415 @MapleTaffyCap Mert posted R code to replicate your DML paper. I tried it on a big dataset & took forever. I have a strong prior it’d be much faster in Julia. I can try to translate when I have time. github.com/demirermert/ML...

タグ：

posted at 10:31:43

K @kar9222

@grant_mcdermott Agree

When picking two of R/Python/Julia, also take into account that Julia is a compiled language for performance

All else being equal, this key point makes Julia a great language to compliment R or Python
- e.g. as compared to R + Python + C when you need performance

タグ：

posted at 13:04:58

津田和俊／急激に進行した網膜剥離と闘って @kaztsuda

4点を選んでLangrange補間して3次多項式を作って数値積分するGIFです
www.desmos.com/calculator/c9r... pic.twitter.com/stoFKjDnqg

タグ：

posted at 13:46:32

非公開

タグ：

posted at xx:xx:xx

400年前にローマ法王に謁見するため太平洋を二度往復した船を復元したレプリカ。場所: 宮城県慶長使節船ミュージアム・サンファン館 www.instagram.com/p/CMqyWb2j0SV/...

タグ：

posted at 14:10:20

#Julia言語要するに、DataFrameで各列(each column)ごとの処理を高速に行いたいときの基本は、

for c in eachcol(df)
cを別の函数の引数に使っているコード
end

foreach(f, eachcol(df))

...

のような形式にすればよい。

dfの各列cを別の函数に丸ごと渡してしまうことが重要なポイント。

タグ： Julia言語

posted at 14:10:23

@m_ero_n seedを102に選んで自分のストーリーにとって都合の良いランダムウォークのサンプルを生成したという「不正の証拠」を自分自身で公開してしまっているのに、ああいう態度を取るのはまったく理解できませんでした。seed 101だと「ストーリー的に都合の悪いプロット」が得られることも容易に確認できる。 twitter.com/genkuroki/stat... pic.twitter.com/W68poKR1Hq

タグ：

posted at 14:37:27

@m_ero_n 人間、誰でも馬鹿なことをやってしまうことがあります。

TJOさんは「確かにこれだと都合の良いデータだけをピックアップして結論を出す『不正研究』になってしまっていますね。非常にまずいことをやってしまいました」と言えば済んでいた話。

科学研究の経験者ならこのまずさは認識できていたはず。

タグ：

posted at 14:40:37

@genkuroki 個人のブログなので～と庇う人が職業研究者にすらいましたが、個人のブログだからこそ謝って済むだけの話なのに変に徹底抗戦して拗らせてしまった。
ただ、ある意味かの御仁の馬脚が現れた事件でもありましたのでTJOさんの信者以外にとっては有益だったのかもしれません

タグ：

posted at 14:43:25

@m_ero_n 松浦健太郎著『StanとRでベイズ統計モデリング (Wonderful R) 』
www.amazon.co.jp/dp/4320112423

は細部に気を使って書かれている非常に良い本です。

ベイズ統計と無関係にパラパラめくってプロットの仕方を眺めるだけでも参考になる。

こういう本が売れるとよいと思います。#統計

タグ：統計

posted at 14:46:49

@m_ero_n その「信者」というのが私の視界には入っていないかも😊

タグ：

posted at 14:51:07

@genkuroki 事前ミュートが効果を発揮してますね😄

タグ：

posted at 14:53:37

松浦健太郎著『StanとRでベイズ統計モデリング (Wonderful R) 』
amazon.co.jp/dp/4320112423

が売れて、著者に興味を持った人がブログ

statmodeling.hatenablog.com/entry/bayesian...

を発見する。こういう流れができれば、みんなハッピーになれる。

この本の考え方とセンスは普遍的でStanと無関係に価値があります。

タグ：

posted at 14:58:29

#Julia言語 2次元Isingモデルのアニメーションを10行で！

10行シリーズ全ソースコード↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/UVFEihNjDU

タグ： Julia言語

posted at 17:47:50

#Julia言語 10行で作った2次元Ising模型の動画その1

200×200の盤面全体を200回更新している。 pic.twitter.com/1uPE4pIznX

タグ： Julia言語

posted at 17:50:22

#Julia言語 10行で作った2次元Ising模型の動画その2

200×200の盤面全体を200回更新するごとにプロットすることを200回繰り返している。 pic.twitter.com/BvdCT5u0c4

タグ： Julia言語

posted at 17:51:30

#Julia言語逆温度βはサイズ→∞のときの理論的な相転移点直上の値にしてあります。

タグ： Julia言語

posted at 17:52:47

Haruhiko Okumura @h_okumura

生協で買ったDell Latitude 3301（Windows 10 Pro，メモリ8G，SSD 256G）は税込99,900でとてもよかった

タグ：

posted at 17:57:20

#Julia言語 10行シリーズ

1. 信頼区間が二項分布モデルのパラメータ値を含む確率
2. 同ケースでの信頼区間の分布
3. 以上のFisher検定版
4. 3次元ランダムウォーク
5. Ising2D.jl
6. 蔵本モデル
7. ジュリア集合
8. マンデルブロ集合
9. ライフゲーム
10. 2D Ising model
↓
nbviewer.jupyter.org/gist/genkuroki...

タグ： Julia言語

posted at 17:58:11

ごまふあざらし(GomahuAzaras @MathSorcerer

@NotionHQ さん、Juliaに対応してくれると嬉しいです pic.twitter.com/z7Q2sXNmwC

タグ：

posted at 18:13:15

虚無ドットゆのっち ོ⃟⿻☀️ @ys_uec

AWS の CloudShell で wget/tar で Julia インストール・展開できたでキュ．Python入ってるんで PyPlot.plot で関数をプロットして描画結果を Actions -> Donwload file で落とすことができた.

Plots.jl はアーティファクトが重いので１GB以内に多分入らない． pic.twitter.com/IXAvyUkt69

タグ：

posted at 18:54:09

ごまふあざらし(GomahuAzaras @MathSorcerer

@h_okumura 生協の人が保証が手厚いことを売りにしてたけど、
それに対して半額の同じスペックのパソコンを買って、
壊れたらもう半額でパソコンをもう一台買う方がお得だってぼやいてた人がいて確かにその通りだなあと思った。
もう一台買い直すときにはスペックも格段に向上しているはずだし。

タグ：

posted at 18:59:10

２Core, メモリ４GB
これぐらいのとこなら遊べそう.

julia> versioninfo(verbose=true)

でもっと色々わかる．

#Julia言語 pic.twitter.com/WQmPCoLtJp

タグ： Julia言語

posted at 19:11:12

baniko @825jp

juliaのDataFramesとQueryでデータを操作する ift.tt/3s9gqrA

タグ：

posted at 19:20:15

非公開

タグ：

posted at xx:xx:xx

非公開

タグ：

posted at xx:xx:xx

k@migaki.be @k_migaki

塾も学校も課程を消して整然とした式と答えだけを残させる。真面目に有害。
矯正のために消しゴム禁止で解かせることもある。スペースの使い方やちょったした計算や図って凄い大事なのにね。 twitter.com/mimizuku_tutor...

タグ：

posted at 19:58:48

非公開

タグ：

posted at xx:xx:xx

Mattias Villani @matvil

The interactive and reactive Julia notebook system Pluto is having its own mini-conference. plutojl.org/plutocon2021
Can one of the tech giants please pour money into the development of Pluto? It is such a great tool for teaching and research communication.

#julialang pic.twitter.com/K8wx14KBiM

タグ： julialang

posted at 21:01:17

Iwao KIMURA @iwaokimura

NotionをChromeでweb越しに編集すれば，GhostTextと組み合わせて好きなエディタで編集できる．
chrome.google.com/webstore/detai... pic.twitter.com/ZbOYFj2H4y

タグ：

posted at 21:04:31

Iwao KIMURA @iwaokimura

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

Emacsの場合は以前Qiitaに書いたのとほぼ同じ．
qiita.com/iwaokimura/ite...

タグ：

posted at 21:05:56

MicrosoftにもJulianがいるのか（そりゃいるだろうが）

タグ：

posted at 21:10:37

アラミス @YfU5McAI0fpqGOg

落ち着くかな？こんなんで… pic.twitter.com/qlkXLV101s

タグ：

posted at 22:28:06

Tomoya TENDO @end_tt

吉本Cubeの新しいやつ！！
Yoshimoto CUBE GREAT No.1 youtu.be/Isg7YsnWy9g @YouTubeより

タグ：

posted at 22:28:53

o_kn @o_kn

🐶🐶🐶🐶🐶🐶
🐶周辺尤度🐶
🐶🐶🐶🐶🐶🐶

タグ：

posted at 22:30:32

Massimo @Rainmaker1973

That time a team of researchers recreated the famous Draupner freak wave measured in the North Sea in 1995, and found it bore an uncanny resemblance to "The Great Wave off Kanagawa' by the Japanese artist Katsushika Hokusai [read more: ow.ly/9iCt30np5oT] pic.twitter.com/17SvgqL90B

タグ：

posted at 22:32:34

Zenn公式 @zenn_dev

📈 ピックアップ ✨

FastGaussQuadrature.jlで数値積分しましょう
by @hyrodium

zenn.dev/hyrodium/artic...

タグ：

posted at 22:38:00

#Julia言語 10行では無理なやつ

正方形上の自由境界条件の波動方程式の離散化をDifferentialEquations.jlで数値的に解いて動画を作ってみた↓
nbviewer.jupyter.org/gist/genkuroki...

これを作るために必要な行数は40行程度でした。 twitter.com/genkuroki/stat... pic.twitter.com/0QZAtKVwpd

タグ： Julia言語

posted at 22:49:30