黒木玄 Gen Kuroki(@genkuroki)/2019年10月22日

#Julia言語

using PyCall
PyCall.python

でどのpythonを使う設定になっているかが分かる。私の場合には

C:\Users\genkuroki\.julia\conda\3\python.exe

これで ~/.julia 以下のpythonを使うようになっているかどうかが分かる。フルエラーメッセージの情報がないと辛い。

twitter.com/takuizum/statu...

タグ： Julia言語

posted at 23:54:05

@konamih @suokubo #Julia言語では

f(x; μ=0.0, σ=1.0) = 1/√(2π*σ^2)*exp(-(x-μ)^2/(2σ^2))

のようにして、デフォルト値を持つキーワード付き引数 μ, σ を付きの函数を定義できます。キーワード引数のリストは ; の後に書く。私は多用しています。

f(1.0)

f(3.0, σ=5.0)

のように使えます。

タグ： Julia言語

posted at 23:43:58

#Julia言語

モジュールから何がエクスポートされるのかを見るとき
names 関数が使えるできゅ！添付はセッションを起動した時の例でこの場合はMainになりますね． pic.twitter.com/nRoo1wt7LC

タグ： Julia言語

posted at 23:43:54

@konamih @suokubo #Julia言語

f(x) = 2x^2 + x +3
f(1.5) # 引数xがFloat64のf(x)がコンパイルされて実行
f(3.4) # コンパイル済みのf(x)を実行
f(4) # 引数xがInt64のf(x)がコンパイルされて実行

こんな感じ。

タグ： Julia言語

posted at 23:36:47

ロタ @Rota_JP

abcnshogiblog.com/blog-entry-143...
＞藤井猛九段
＞四間飛車はソフト的にはマイナス300でも私が指せばプラス５００。
＞なぜなら私は四間飛車の指し方を知っているから。

やばいかっこいい

タグ：

posted at 23:36:13

@konamih @suokubo #Julia言語函数にすると速くなる理由はネイティブコードにコンパイルされるから。Juliaでは

f(x) = 2x^2 + x + 3

のように定義された函数が実行されたとき、函数と引数の型の組み合わせが初出ならコンパイルしてから実行してくれます。函数にしないとコンパイルしてくれない。

タグ： Julia言語

posted at 23:34:06

#Julia言語

 docs.juliaplots.org/latest/support...

を見ると stephist というのがあります。

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

twitter.com/KB_satou/statu... pic.twitter.com/cXAwYdvJhI

タグ： Julia言語

posted at 23:15:04

あと、 ⊻ とか ≥ とか問題の定義に使われてるような unicode 文字がメソッドとしてそのまま使えるので、とりあえず brute-force するときとか問題文コピペする感覚でできるので超楽。つまり神。

タグ：

posted at 22:59:56

#Julia言語では3乗根は cbrt(x) で計算できるのですが、

\cbrt [TAB] 10 → ∛10 → 2.154434690031884

のようにも計算できる。残念ながら ∜ は定義されていない。

タグ： Julia言語

posted at 22:15:38

@tsatie #Julia言語その後に

f(1.3)

を実行すると、すでにFloat64引数のf(x)はコンパイルされているので、コンパイル結果を使って f(1.3) を評価します。

f(2)

を実行すると、Int64の引数なのでそれに合わせてf(x)が再度コンパイルされる。

函数実行時に函数と引数の組み合わせが初出ならコンパイル。

タグ： Julia言語

posted at 22:07:17

@tsatie #Julia言語は函数を実行したとき、その函数と引数の型が初出ならばその函数を引数の型に合わせてコンパイルしてから実行するという仕組みになっています。

f(x) = x + oneunit(x)

と定義して、

f(0.1)

を実行すると、引数の型Float64に合わせてf(x)がコンパイルされてから実行されます。

タグ： Julia言語

posted at 22:05:35

#Julia言語数学記号の使い方は大体において世界共通で、数学記号の入力の仕方は LaTeX 方式がスタンダードになった。そのおかげで、LaTeXで数式を入力するのと同じ感覚で、Julia対応の環境では、√, ≥, Δx, α, x₀ などを入力できる。

こういう環境はJuliaと無関係に広まって欲しいと思う。

タグ： Julia言語

posted at 22:02:07

#Julia言語 √ や ≥ などが使えることや Δx や α や x₀ などを変数や函数の名前として使うことが普及するためには、

\sqrt [TAB]
\ge [TAB]
\Delta [TAB] x
\alpha [TAB]
x\_0 [TAB]

などでそれらを気楽に入力できる環境も同時に普及させる必要がある。Juliaの開発者達は実際にそれを行った。

タグ： Julia言語

posted at 21:59:35

#Julia言語 √ や ≥ などが使えることや Δx や α や x₀ などを変数や函数の名前として使えることは、非本質的なことだと思っていたが、実際にそれが普及しているのを見ると、アルゴリズムが書いてある数学的文献と記号が揃っているのは、非常にありがたいことであることがわかった。

タグ： Julia言語

posted at 21:56:39

上原哲太郎/Tetsu. Uehara @tetsutalow

#Julia言語

Juliaでは sqrt(10) は √10 と書ける。

\sqrt [TAB] 10 で入力可能。

他にも Δx や x₀ などが変数名、函数名として使える。
(日本語も使える。絵文字も使える。)

f(x) = 2x^2 + 3x + 1

のように函数が定義できる。

twitter.com/BBFCFiro/statu...

タグ： Julia言語

posted at 21:54:38

#CSS2019 のCSSx2.0において、反PPAPに加え反PHSというのを提案しました（私ではなくテツ太郎が）。Printして、Hanko押して、Scanして送るバカバカしさ。我々はこれに代わる何かを求めているはずです。良いアイデアを皆で考えませんか。私も考えてます。 #css2019_ipsj
www.slideshare.net/tetsutalow/ppa...

タグ： CSS2019 css2019_ipsj

posted at 21:54:31

鳩 @pegion_HOLE

開発言語にJuliaが含まれている求人を教えてもらった。日本語書籍も出たし、普及するのも案外早いのかも。

タグ：

posted at 21:44:18

#Jupyter

パソコン上でJupyterLab serverを走らせておいて、

localhost:8888/tree

にアクセスすればJupyter notebookを使用可能。

ZeroTier経由でそのパソコンにアクセスすれば、出先でiPadからJupyter notebookを使用可能。Julia, Python, R, Ruby, Wolfram言語を使えるようにしてある。

タグ： Jupyter

posted at 21:41:56

#Julia言語私が行ったJulia環境の作り方の詳細は

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

で公開されている。Julia v1.1.0当時のものだが、現在でもそのまま通用するはず。

タグ： Julia言語

posted at 21:38:53

#Julia言語 Juliaは公式バイナリをダウンロードしてインストールするのが安全牌。失敗したことは皆無。

面倒なのはPythonとの連携。Python側の事情の詳細まで気にすると複雑過ぎて大変になる。

その辺の説明が面倒なので私はJulia専用のminiconda環境をごそっと新たに入れて使ってしまっている。

タグ： Julia言語

posted at 21:37:58

#Julia言語

julia> ENV["PYTHON"] = ""
julia> ]
pkg> build PyCall

すれば、Julia専用のminiconda環境をごそっと入れてくれるはず。自前のPythonを使いたいなら

julia> ENV["PYTHON"] = raw"pythonバイナリのフルパス"
julia> ]
pkg> build PyCall

twitter.com/takuizum/statu...

タグ： Julia言語

posted at 21:35:47

@konamih @suokubo #Julia言語ソースコードは

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

にあります。

タグ： Julia言語

posted at 21:25:12

@konamih @suokubo #Julia言語函数にすればコンパイルしてくれるのでもっと速いです。一瞬で計算が終わる。アスペクト比の調節もしておきました。コードは違いますが、アルゴリズムは同じ(見ればわかるはず)。異なる点は

* 函数にした
* 作画範囲をrangeで指定
* 複素数で計算
* abs2で絶対値の二乗
* heatmap pic.twitter.com/OqBnRf2Hxw

タグ： Julia言語

posted at 21:22:19

何？ #Julia言語を始める？
- まずはJuliaと検索した結果をGoogleの検索言語を日本語から英語に切り替えるのだ
- そうするとjulia って検索するとすぐにドキュメントが読めるのだ．
- GitHubのIssueなど一次情報にたどりやすいのだ．

タグ： Julia言語

posted at 21:18:56

ちなみに上の`combinations(pattern, n::Int)` は
```
combinations(pattern, n::Int) =
collect(Iterators.product((pattern for i = 1:n)...))[:]
```
のワンライナーで置き換え可能。

タグ：

posted at 21:06:24

これ参考になった。
splattingされた引数はspecializeできない(つまり高速化できない)ので、再帰とか何回も呼ばれるメソッドには使わない方がよい。

Splatting arguments causes ~30x slow down discourse.julialang.org/t/splatting-ar...

e.g.:
github.com/aviatesk/googl...

#julialang

タグ： julialang

posted at 21:05:09

juliaで競技プログラミングをやってみた感想としては, "walks like Python" よろしく非常に簡単にコードが書けるので感触はとても良かったです。
("runs like C" が生きるほど(言語面での)スピードが求められる問題はなかった)

タグ：

posted at 20:20:41

そういえばこれの julialang での解答を GitHub に push しておきました。
(The equations と Shifts の hidden testset はまだ)

github.com/aviatesk/googl... twitter.com/kdwkshh/status...

タグ：

posted at 20:14:55

KokyuHatuden @breathingpower

「EM菌はニセ科学」と社会に知られるようになったのは、この番組の功績が大きいと思います。単調な文字起こしを、のべ14万人以上の人に見ていただきました。ありがとうございます。

フジテレビEM菌報道の文字起こし - Togetter togetter.com/li/395287 @togetter_jpさんから

タグ：

posted at 20:09:43

非公開

タグ：

posted at xx:xx:xx

こなみひでお @konamih

@suokubo PC9801とN88BASICで一晩かかってやっと一枚描けたものです。今は一瞬。

タグ：

posted at 19:11:23

シキノ @sikinote

プログラムの確かめです。抗力が０になる点をy=1になるような初速度を与えた時、確かにy=1で束縛運動（青線）から解放されていることをプログラムで確かめました。 pic.twitter.com/gfSy2J126A

タグ：

posted at 19:05:53

SusumuOkubo @suokubo

kitty-poppins(=^・^=) @Tiger_Poppins

雑用事務仕事で脳みそがやる気を失っているので、何十年ぶりかMandelbrot集合など描いてみる。
学生の頃にX-windowsで描いたのに比べて、Julia楽ちんすぎる。そして計算も速すぎる。X68030 33MHzとCorei5 1.7GHzを比べてもいかんか。
問題はpyplotのアスペクト比がくるっていることか。
#Julia言語 pic.twitter.com/0Ck80c0qSM

タグ： Julia言語

posted at 19:01:00

@mihyonsong あらまあ！マークされてますね！

タグ：

posted at 18:42:24

宋美玄 @mihyonsong

@Tiger_Poppins それが、何故かブロックされてるので出来ないんです。

タグ：

posted at 18:41:13

青色1ごう @BBFCFiro

すごぃなJuliaLang。ルートを√で表記できる

タグ：

posted at 18:29:51

K.B.砂糖 @KB_satou

kitty-poppins(=^・^=) @Tiger_Poppins

中心極限定理をゾンビコーディング
(頭は動いてないけど体は動いてる) pic.twitter.com/fVlejt13hH

タグ：

posted at 18:26:25

@mihyonsong 問題を提起するならスクショで相手に見えない方法で投稿ではなく、コメRTがいいと思いますよ。

タグ：

posted at 18:22:58

#Julia言語 #Julialang twitter.com/julialangisthe...

タグ： Julialang Julia言語

posted at 17:21:33

非公開

タグ：

posted at xx:xx:xx

@julialangisthe 天才 pic.twitter.com/2FrYdn07ZA

タグ：

posted at 17:18:10

非公開

タグ：

posted at xx:xx:xx

なんか，アップデートされてないのでありまし「た」というのが正しそう． twitter.com/MathSorcerer/s...

タグ：

posted at 17:12:07

Julia 業界では \xor でタイプできるみたいですね． twitter.com/daimare_496/st... pic.twitter.com/yDso2N2DKW

タグ：

posted at 17:10:07

宋美玄 @mihyonsong

はあちゅう氏と繋がりのある方、マジで命に関わるのでアドバイスをお願いします。 twitter.com/mihyonsong/sta...

タグ：

posted at 16:52:32

宋美玄 @mihyonsong

はあちゅう氏、不適切な健康情報を拡散して炎上したら消す、というのではなくあらかじめ裏を取って欲しい。周りにヘルスリテラシーのある助言者はいないのか。医療クラスタも好きで批判してるんじゃなく止むを得ず火消ししてるんだよ。
育児関連でこういう事するのだけはマジで勘弁。ほんまに頼むわ。 pic.twitter.com/aYHlKuZL1l

タグ：

posted at 16:47:47

#Julia言語例：配列 X の各々の要素に函数 f を作用させるには f.(X) と書く。さらにその和は sum(f.(X)) と書ける。しかし、そのように書くと、メモリを余計に消費する。この点を改善するにはどうすればよいのだろうか？これについえては sum(f(x) for x in X) という逃げ道がある。

タグ： Julia言語

posted at 16:03:55

#Julia言語多くの作業をforループ抜きに可能なのだが、そうした場合に難しいのはメモリを余計に消費してしまいがちなコードになってしまうこと。

タグ： Julia言語

posted at 16:02:25

#Julia言語

似たようなことをforループ抜きに行ってみた。

1つ目のポイントは cumsum!(X, X) の部分。これを X = cumsum(X) に変えるとメモリを無駄に消費することになってはまる。
docs.julialang.org/en/v1/base/arr...!

2つ目のポイントは @. マクロを使っていること。

twitter.com/KB_satou/statu... pic.twitter.com/5HBqKkJyn0

タグ： Julia言語

posted at 15:58:59

#Julia言語 #JuliaLang

twitter.com/MathSorcerer/s...

タグ： JuliaLang Julia言語

posted at 15:51:56

www.simonwenkel.com/2018/10/25/Typ...

のコードをそのまま実行して表示できた．
dot がない場合は
brew install graphvis で pic.twitter.com/QhjjoaNMp9

タグ：

posted at 15:46:43

ラストアップデートが 5years old で絶望している

タグ：

posted at 15:23:21

Julia でも計算幾何学（Computational Geometry）
をしようぜーな動きはありますね．ただし，wrapper 多め
github.com/JuliaGeometry

タグ：

posted at 15:17:20

#Julia言語で python でいうところの if __name__=="__main__" 判定をしたい場合

は abspath(PROGRAM_FILE) == @ __FILE__
を使うと良いみたいです．
See Julia documentation

How do I check if the current file is being run as the main script?

docs.julialang.org/en/v1/manual/f... pic.twitter.com/fAuLw9oD4Y

タグ： Julia言語

posted at 15:10:34

K.B.砂糖 @KB_satou

大数の法則的な.jl
|(n 個の rand() の平均) - 0.5| の両対数のプロット pic.twitter.com/Oj9i3i0Y8b

タグ：

posted at 14:45:16

#Julia言語このスレッドの上の方でやっていることは、

* Julia言語(糊言語)からプロットのためにPythonを呼ぶ。
* Python(糊言語)はプロット時にpdflatexを呼ぶ。

という感じになっています。Juliaを使っていたのに、結果的にPythonだけではなく、pdflatexで日本語を使う方法を知る必要が生じる。

タグ： Julia言語

posted at 14:28:29

#Julia言語 matplotlib.pyplotは確かに便利で、それをPyPlot.jlやPlots.jl pyplot()経由で使う場合には、matplotlib.pyplotに関する知識が結構必要になる。

結果的にPythonにおける細かい話にもある程度詳しくなってしまう。

Juliaは糊言語の一種なのでこういう感じになりやすい。

タグ： Julia言語

posted at 14:24:38

#Julia言語注意：私の高精細表示への対応の仕方は IJulia (Jupyter)前提の対処法なので、もっとよくわかっている人が根本的な対処をしてくれると助かるのだ。

あと、Plots.jlでpyplot()を使うと思わぬところで、プロットが二重表示される。たぶん、バグ。PyPlot.clf() を適当に入れれば消せる。

タグ： Julia言語

posted at 14:22:18

#Julia言語 using Plots; pyplot(fmt=:svg) とするとPyPlotバックエンドでも高精細表示される。

このスレッドの上の方のテーマは、そのようにするとmatplotlibのpdflatex対応の機能を使ったとき不具合が生じるのでpngでの高精細表示の函数を書いたという話。

プロットで見た目の良さにこだわると沼。

タグ： Julia言語

posted at 14:18:44

#Julia言語の開発元のMITでの講義録を見ても、PyPlot.jlもPlots.jlもこだわりなく使っているようなので、私も「純粋主義」に陥らないように自分にとって便利だと思うものを自由に使うようにしている。

Juliaは「気楽にコンピューターの性能を引き出す道具」だと思っているのでそれでいい感じ。

タグ： Julia言語

posted at 14:18:44

#Julia言語

 twitter.com/astellon_music...

Plots.jl経由でGR()を使えば、legend=:outertopright などでlegendを外に追い出せます。しかし、gr()でそれをやったら縁の右端が切れたので、pyplot()を使うことにした。

添付画像1つ目がgr()版で2つ目がpyplot()版。

nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/zAmiOQwTSJ

タグ： Julia言語

posted at 14:18:44

ラーラー @sun_ozon

開発言語はC++、RUST、Juliaを足して３で割ったぐらいかな。1か月でコンパイラ周りを勉強して1か月で実装したスピードが何よりもすごい。

スーパー中学生誕生、プログラミング言語わずか数週間で開発、U-22プログラミング・コンテスト2019 - BCN＋R www.bcnretail.com/market/detail/...

タグ：

posted at 13:08:36

台風 @taifu21

@choko_midori @golgo_sardine 2017年07年10日の東京新聞が学習指導要領解説の該当箇所について文科省に問い合わせていますね。
文科省の担当者は「式が５×４でも４×５でもまったく問題ない」と言っているので、「かけ算の順序指導仕方ないよ」と言うためには、何故文科省の方針に逆らうのかの説明が求められる事になるのでは？ pic.twitter.com/lEJw5P8Gk9

タグ：

posted at 12:26:33

#Julia言語試しにJuliaで特殊函数を計算するコードの最適化の実験をやってみたい人が見るべきなのはMITでの講義の宿題の答え↓

nbviewer.jupyter.org/github/steveng...
Problem 1: Exponential integrals

Juliaによる数値計算、プロット、簡単な数式処理、マクロ、などのtutorialにもなっている。

タグ： Julia言語

posted at 12:22:58

#Julia言語

私は学生時代から「exampleは？」という質問を常に受けてきて、具体例をノータイムで出せないことは理解していない証拠だとむなされる文化圏で育ったので、例も出します。以下のリンク先。

理解度を測る最も簡単な質問は「例は？」だと思う。

関連↓

twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 11:53:27

Iwao KIMURA @iwaokimura

学園祭っぽい催しで数学科は「ポアンカレー」というカレーを出すという話はよく聞く．
一緒に「デュアリティー」というお茶を出すのはいかが．

タグ：

posted at 11:51:31

#Julia言語偏微分方程式の数値解を求めたい人にとって必須な事柄の説明は

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

にあります。具体的には = を含めた演算記号に . を付けることと、@., @ views, @ viewマクロの使い方を知ることは大事。

タグ： Julia言語

posted at 11:29:45

#Julia言語行列やベクトル達に処理を施す場合には最後の代入操作も含めて dot syntax を適用した方が良いことが多いです。用意しておいた配列内で計算するようにして、できるだけ新たな配列を確保を防ぐ方法を知っておくことは、偏微分方程式の数値解を求めたい人にとって必須の知識。続く

タグ： Julia言語

posted at 11:29:44

#Julia言語 Fortanでfor文で回していたコードをJuliaにそのまま移すことは易しい。

しかし、ベクトル化されたコードを移す場合には、dot syntaxに慣れる必要があります。Juliaでは行列A,Bについて

A * B #→行列の積

となるので、成分ごとの積にしたければ

A .* B

と . を付ける必要有り。

タグ： Julia言語

posted at 11:29:42

普段なら「原因を特定して、シンプルで合理的な解決策を探すこと」についていつも考えていそうな人が、エロい検索結果に心が影響されて「シンプルで合理的な解決策」について考えなくなる様は印象的。😅😅😅

普段モードで頭を動かせていれば、Googleの設定をいじることを自力で見つけていたと思う。

タグ：

posted at 11:17:19

日本在住者はGoogleでJuliaを検索すると困ったことになるので、Googleの設定で

地域に設定：アメリカ合衆国

にしておくべきです。Juliaの検索での「不具合」が解消します。この設定で普段の生活で困ったことはありません。 #Julia言語 pic.twitter.com/QGUwOyffg2

タグ： Julia言語

posted at 11:07:57

天むす名古屋 Temmus @temmusu_n

Julia入門には公式ドキュメント

docs.julialang.org/en/v1/

の自分に必要な部分を読むことがベストだと思うのですが、公式docは分量が多いので大変だし、英語が嫌な人もいるかもしれない。そういう場合には #Julia言語タグをつけて疑問を述べれば日本語で答えてくれる人が出て来るかも。

タグ： Julia言語

posted at 11:01:15

天むす名古屋 Temmus @temmusu_n

@metameta007 #超算数とはいえ、以降の記述も、包含除を第一、等分除を第二と称し直した上で、包含除の商に限って比という異称があることを導入したり、5÷5=1や8÷1=8を【割り算の意義】で説明せよと指示するので、完全に統合したわけではありません。この辺り、現代の算数教育界の慣行ととてもよく似ています。 pic.twitter.com/4Ogj2TE2eD

タグ：超算数

posted at 09:13:30

天むす名古屋 Temmus @temmusu_n

@metameta007 #超算数 64頁では【(甲)(乙)意義の相違を説明したる後更に(甲)(乙)を合併して次の如くに言ひ表はすことを得べし。/_割り算は積と因数の一つを知りて他の因数を見出す為めに行ふ計算なり。_】dl.ndl.go.jp/info:ndljp/pid...と甲と称される包含除と、乙と称される等分除を、統合します。 pic.twitter.com/OxLeWMUQ6F

タグ：超算数

posted at 09:01:53

@metameta007 #超算数藤沢は算術教科書では別の態度を示します。
藤沢利喜太郎編『算術教科書』上巻東京、大日本図書、1896年。
#掛算の交換法則は、35頁の6×4=4×6の例示に留めdl.ndl.go.jp/info:ndljp/pid...、62頁で割算に二つの意味を認めますdl.ndl.go.jp/info:ndljp/pid...。 pic.twitter.com/uM9DypfooX

タグ：掛算超算数

posted at 08:54:47

abap34 @abap34

Ayako Shibata@ウィメンズヘ @ayako700

Juliaは公式docみたらなんでも書いてあるからオススメだよ

タグ：

posted at 08:11:07

〈 Berger | Dillon 〉 @InertialObservr

(゜o゜;　えーー！アップルウォッチに話しかけると、そのまま電子カルテ入力されるとは。。。

未だにアップルウォッチで心電図とれない日本って、何周遅れてるんや・・・

japanese.engadget.com/2019/10/14/app...

タグ：

posted at 07:47:56

A Cycloid is a Tautochrone curve, meaning all points take the same time to reach the bottom regardless of initial height pic.twitter.com/sX94kTPWll

タグ：

posted at 02:25:19

tksm @tksmml

前処理に部分的にjulia使いたかったのでpythonから呼べるようにdockerで環境構築しました

はてなブログに投稿しました #はてなブログ
PythonからJulialangの関数を使える、Python + Julia…
tksmml.hatenablog.com/entry/2019/10/...

タグ：はてなブログ

posted at 00:15:01