黒木玄 Gen Kuroki(@genkuroki)/2021年06月

@genkuroki ありがとうございます！初歩的で恥ずかしいですが、定数をconst宣言することでかなり速くなりました。グローバル変数の引数渡しもやってなかったので試してみます。

タグ：

posted at 00:12:12

#Julia言語

ノートブックのインターフェースが日本語になってしまう症状。

D:/.julia/conda/3/Lib/site-packages/notebook/18n/ja_JP

の名前をja_JP_oldに変えた。

twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 00:25:10

kakehuton @shikihuton

チュートリアルを利用してモダンな数値計算方法へ移行するのを試みるのはありかも知れない。Juliaってgoogle colabで動かせんだっけ？

タグ：

posted at 00:28:09

@kado_judo0312 @benkyouaho #Julia言語

①プロットするだけならこれが簡単。

②失敗の再現とその理由。整数の整数乗は危ない。

③x=0.:40としてFloat64のレンジにすれば正しく計算される。

④巨大数/巨大数は危ない。全体の対数をまとめて計算してexpするのが常套手段。

物理でも n! そのものよりも log n! が重要。 pic.twitter.com/0V7XjJ9dOA

タグ： Julia言語

posted at 00:30:16

@akairamp #Julia言語のパッケージでは、

prob = (問題を記述するパラメータ達をまとめたもの)
sol = solve(prob) # 問題の解を返す函数

のスタイルをよく見ます。

solve(prob)の中で使われている函数達にもprobを引数として渡す。

続く

タグ： Julia言語

posted at 00:34:22

@akairamp #Julia言語パラメータ達のまとめ方としては、

a = 1.0
b = 2.0
c = 3.0
prob = (; a, b, c)

のような感じでよいと思います。函数の中では

function f(prob)
a, b, c = prob
a, b, c を使った計算
end

のように使う。a, b, c = prob の1行をコピペしまくる手間がちょっと増えるだけ。

タグ： Julia言語

posted at 00:37:09

@akairamp #Julia言語 v1.7 (現在のnightly build) 以上ならば

function f(prob)
(; c, a) = prob
c, a を使った計算
end

のスタイルで、probから一部のパラメータのみを抜き出せます。

私はnightly buildを使っています。nightly buildだと 2√3 が合法的なコードになります。

タグ： Julia言語

posted at 00:41:02

JuliaHub @JuliaHub_Inc

A key benefit of #Julialang is that the same language can be used for #prototyping, model building, interactive #datascience, and production deployments. In this video we will see how Julia fulfils that last promise. youtu.be/hbvC5C_4Z64

タグ： datascience Julialang prototyping

posted at 00:47:38

discourse.julialang.org/t/plots-downgr... で話題になっていた。

Jupyterなどで using Plots; plot(sin) の結果が "tiny plot" になっている人は

pkg> up GR

すればその問題を解決できます。 #Julia言語 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 00:59:21

@benao_python #Julia言語については以下のリンク先の動画と関連の沢山の動画を見ればどういう感じのものかよく分かると思います。

youtu.be/kc9HwsxE1OY

タグ： Julia言語

posted at 01:26:11

べナオ @benao_python

富谷(助教);監修シン仮面ライダー @TomiyaAkio

@genkuroki 情報提供ありがとうございます！

タグ：

posted at 01:28:41

「Python」の弱点やそのとらえ方--生みの親、グイド・ヴァンロッサム氏が語る(ZDNet Japan)
「JuliaとPythonを同じ日にコーディングしたいと考える人などいないはずだ」
してるんだよなぁ
#Yahooニュース
 news.yahoo.co.jp/articles/4fbaf...

タグ： Yahooニュース

posted at 02:25:57

Chuui @es_chuui

Python開発者(van Rossum氏)による各言語に対する雑談的な記事。

明言していないが暗黙的には、Go・Rustはアゲ、#Julia言語についてはサゲ、の模様。

>「Juliaでウェブサーバーを開発し ... など誰もいないと私は考えている」
>「1から始まり...」

japan.zdnet.com/article/351713... @zdnet_japanより

タグ： Julia言語

posted at 07:40:36

Chuui @es_chuui

Web開発にJuliaについては、私は専門外なのでなんとも。

だが、配列が1から始まる点、については、実際のところどちらが良いのだろう。
C言語から Matlab,Scilab を始めて触った時、最初は1始まりに戸惑ったものの、0始まりが絶対的に優れているとは限らないはずだ、と思ったものだが。

タグ：

posted at 07:48:17

Akinori Ito @akinori_ito

RTした記事のリプライに反ワクチンの人が結構集まっている．思ったよりも数いるのかな

タグ：

posted at 07:49:17

Akinori Ito @akinori_ito

COBOLは自然言語由来だからノーコード

タグ：

posted at 08:16:38

TaTa @Tata_oac

Ninja DAO | CryptoNi @CryptoNlnjaNFT

@ rput/rgetでR<->julia間で変数のやりとりができるので...
pythonよりRとの親和性が高い

タグ：

posted at 09:31:07

いわゆる「文系批判」の人たちの「文系イメージ」は現実より30年は遅れているので、「嫌韓」の人たちのやはり30年は遅れている「韓国イメージ」と同じで、もう放っておくしかないと思う。

タグ：

posted at 10:07:17

しんしん @zinziroge

Juliaの配列はゼロオリジンなのかぁ。
matlab寄りなのね。 twitter.com/Atsushi_twi/st...

タグ：

posted at 12:07:33

unitylist.com/p/d1j/Unity-.jl

github.com/StreetLevel/Un...

これまだ動く？ #Julia言語

タグ： Julia言語

posted at 13:06:09

ストレンジまるむし博士あるいはいかにして @fjt

@Atsushi_twi #Julia言語のコミュニティの側の反応は「こいつ、Juliaのことをろくに調べずに、印象論でディスっている」(私的要約)。

私もヴァンロッサム氏はろくに調べずにいい加減なことを言っていると思いました。

そういう態度を取らざるを得ないほど追い込まれている？(笑)

discourse.julialang.org/t/what-python-...

タグ： Julia言語

posted at 13:22:23

自作のライブラリ群をjuliaに移植することにしたよ。ぱいちょんに移植しようと考えたことは一度もないね。

タグ：

posted at 13:36:42

#Julia言語

正五角形上の波動方程式
(isin_funcを変えれば任意の領域を扱える)

ソースコードも全公開

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

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

In[5]がノイマン境界条件の場合の微分方程式を記述する函数
In[10]はディリクレ境界条件の場合 pic.twitter.com/Lvhe5FeBBT

タグ： Julia言語

posted at 13:51:01

#Julia言語

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

境界が自由に動けるノイマン境界条件の場合

山(赤)と谷(青)は境界にぶつかるとそのまま山(赤)と谷(青)のまま反射する。 pic.twitter.com/nTQacFMcJp

タグ： Julia言語

posted at 13:53:20

#Julia言語

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

境界が0に固定されているディリクレ境界条件の場合

山(赤)と谷(青)は境界にぶつかるとそれぞれ谷(青)と山(赤)に反転して反射する。

こういう類のことが1次元の波動方程式で見た方が見易い。 pic.twitter.com/vvSMpWyvQL

タグ： Julia言語

posted at 13:55:22

#Julia言語リンク先添付画像中の微分方程式を記述するg_neumanｎ!及びg_dirichlet!函数の定義が非常にシンプルであることに注意！

複雑な部分は、離散化を記述するパラメータを1つにまとめるコンストラクタParamにすべて押し付けられている。 twitter.com/genkuroki/stat... pic.twitter.com/yTSWHV1cLB

タグ： Julia言語

posted at 14:04:03

#Julia言語コンストラクタParam(isin_func, x, y)は、領域に点が入っているかどうかを判定する函数isin_funcとx軸とy軸の離散化x, yを与えると、領域内の格子点V[k]に繋がっている領域内の点の集まりE[k]などの情報を持ったオブジェクトができる。

タグ： Julia言語

posted at 14:04:04

#Julia言語本当は正方格子による近似的な領域の分割で手抜きをするのではなく、より適切な分割で離散化を記述するべきなのだが、面倒だし、サンプルコードとして複雑すぎるものになりそうなのでやめた。

誰か、より適切な分割で波動方程式を数値的に解くためのコードを見せて下さい。

タグ： Julia言語

posted at 14:07:30

#Julia言語のコミュニティを覗くと、JuliaユーザーがPythonでも楽しく仕事をしている様子も見える。例えば、最近の

discourse.julialang.org/t/what-do-you-...

では、PyCall経由だとほとんどすべてが10%速くなることを発見し、その原因の追求によって、Pythonによる計算も高速化できた、という面白い話があります。 twitter.com/genkuroki/stat...

タグ： Julia言語

posted at 14:25:50

テクニカルな問題の解決に喜びを見出すタイプの人達の話は楽しくて面白いことが多いと思います。

#Julia言語のコミュニティはそういう話題が非常に多くて良いです。

個人的な意見では #Julia言語のコミュニティは社会科学的にも非常に興味深いと思います。中心的人物達が非常にキュート。

タグ： Julia言語

posted at 14:29:22

Hidemoto Nakada @hidemotoNakada

確かにPythonのほうがずっと書きやすいんだけど、Juliaの速さは正義。Pythonのライブラリはちょっと追ってくとすぐC言語実装に潜り込むけど、Juliaだとそんなことはない。複雑なシステムプログラミングにはJuliaのほうが向いているかも。 twitter.com/Atsushi_twi/st...

タグ：

posted at 14:41:19

#Julia言語この辺も面白いです↓

discourse.julialang.org/t/julia-what-i...

なぜかJuliaのコミュニティで子供の腸内微生物の研究の話題が！(笑)

Juliaのコミュニティは「科学や技術の話題を面白がってもらえるだろう」と期待できる人達の集まりにもなっている。

タグ： Julia言語

posted at 14:41:19

USHIHARA @H_USHIHARA

以前ツイートした件( twitter.com/H_USHIHARA/sta... ), twitter.com/H_USHIHARA/sta... や黒木さんのツイート( twitter.com/genkuroki/stat... )にもあるように, Meta.parseを使えば可能だけれど, 保存の段階でDelimitedFiles.jl ( docs.julialang.org/en/v1/stdlib/D... )を使った方が良さそうだと思い始めた.
#julialang

タグ： julialang

posted at 15:04:02

#数楽

【1次元～の方が見易い】

1次元でのLaplacian (d/dx)² の離散化の話に関係がある。

離散化に境界条件も入れると、区間の両端での条件を色々変えると、古典有限型及び古典アフィン型のCartan行列が自然に出て来る。

nbviewer.jupyter.org/github/genkuro... twitter.com/genkuroki/stat...

タグ：数楽

posted at 15:07:46

#数楽要するに、古典{有限,アフィン}型のCartan行列は1次元の離散Laplacianになっている。型の違いは境界条件の違いに対応している。

さらに、それらの特性多項式は本質的にチェビシェフの多項式でもある。

数学的に基本的なオブジェクトはあらゆる場所に顔を見せる。

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

タグ：数楽

posted at 15:07:47

ケフっち＠Deterministic @__cheph__

#数楽 (d/dx)² の離散化を考えると自然に古典{有限,アフィン}型のCartan行列が出て来る。

離散化すると数学的取り扱いが一般に面倒になるのだが、このように数学者達にはよく知られている「きれいな数学」を記述する典型的なモノが自然に出て来てしまう場合もある。

タグ：数楽

posted at 15:11:05

講義でRを使っていて、大学院生はPythonを使っていて、僕自身はメインでJuliaでやってるからバグるバグる。

タグ：

posted at 17:03:35

共同通信公式 @kyodo_official

運動参加者が高須氏ら提訴－愛知、不正署名と責められ苦痛
nordot.app/77237663664968...

タグ：

posted at 17:29:00

非公開

タグ：

posted at xx:xx:xx

#Julia言語具体例

4次の古典Runge-Kuttaの時間刻み幅dt=0.05で長さ1mの単振子を数値的に解いて(添付画像①)、DifferentialEquations.jlのデフォルトのソルバのdtmax=0.05の場合(添付画像②)と比較すると、後者の方がずっと精度が高い。

全コードと計算結果↓
nbviewer.jupyter.org/github/genkuro... twitter.com/genkuroki/stat... pic.twitter.com/9b32faReXl

タグ： Julia言語

posted at 18:20:56

甘党 @amatou_blissful

60%は　3/5ですのにね twitter.com/genkuroki/stat...

タグ：

posted at 18:21:16

#Julia言語さらにDifferentialEquations.jlでは、Hamiltonianを与えただけで、自動的に(自動微分によって)正準方程式を作って数値的に解いてくれる。

単振子程度なら、ハミルトニアンと正準方程式の複雑さは同じ程度だが、一般にはハミルトニアンから正準方程式を作る計算が面倒になったりする。 pic.twitter.com/Fbsjm2Y9cM

タグ： Julia言語

posted at 18:24:26

#Julia言語手軽に使い易いYoshida6を含む大量のソルバがDifferentialEquations.jlに実装済み。

そして、「ハミルトニアンさえ与えれば十分」という機能に限らず、全体像の把握さえ不可能なほど沢山の機能がDifferentialEquations.jlにはある。

DifferentialEquations.jlは積極的に使ったほうが得。

タグ： Julia言語

posted at 18:28:14

あ〜る菊池誠(反緊縮)公式 @kikumaco

#Julia言語 DifferentialEquations.jl の使用法入門に

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

は非常に便利だと思います。

添付画像は私が書いた4次の古典Runge-Kuttaのコードですが、引数と返り値のスタイルを DifferentialEquations.jl のsolve函数とODEProblemコンストラクタの組み合わせに似せてあります。 pic.twitter.com/EFnyjSzKSe

タグ： Julia言語

posted at 18:32:11

ハミルトン系はYoshida6を使っておけばいいと思う twitter.com/genkuroki/stat...

タグ：

posted at 18:33:51

#Julia言語 DifferentialEquations.jlのような巨大パッケージをブラックボックスとして利用するのは怖い。

そういう人は、私が書いたシンプルなRunge-Kuttaの函数を理解して、それと同じようなこと(実際には圧倒的に高級なこと)をDifferentialEquations.jlでは可能だと理解すれば良いと思います。 pic.twitter.com/scvgyoIdsH

タグ： Julia言語

posted at 18:35:48

#Julia言語 DifferentialEquations.jl に限らず、他の多くのパッケージが採用している problem and solver style はJuliaと相性が非常によいと思います。

# 問題の記述
prob = FooProblem(a, b, c, d)

# アルゴリズムの指定
alg = BarAlgorithm(p, q)

# ソルバの起動
sol = solve(prob, alg)

タグ： Julia言語

posted at 18:40:00

#Julia言語函数 solve(prob, alg) は問題probの型とアルゴリズムalgの型の組み合わせごとに異なる計算法(メソッド)で解を求めてくれるようにしておく(そこでJuliaの多重ディスパッチを使う)。

このスタイルで整理できる場合は非常に多いと思います。

Juliaに限らず、広まると良いと思う。

タグ： Julia言語

posted at 18:45:17

天むす名古屋 Temmus @temmusu_n

#Julia言語既出のノート

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

では、私が実装したRunge-KuttaとDifferentialEquations.jlのadaptive=falseのRK4()アルゴリズムが同等であることの確認も行っています。

最後に Test Passed と表示されている。

ブラックボックスの利用ではこういう手続きが必須だと思う。 pic.twitter.com/MOJfuuVLbl

タグ： Julia言語

posted at 18:52:09

かわいそうだと思うけれど、言論の自由に天皇制の名誉が優越するという曲がった考えを抱いている時点で、あまり同情はできない。署名する前にこのような汚辱に塗れることはある程度予想できたのではない？

タグ：

posted at 19:09:25

💖 I'm sponsoring @JuliaLanguage because the language and its community keep bringing me immerse enthusiasm and amazing opportunities. github.com/sponsors/Julia...

タグ：

posted at 19:13:00

新卒入社以降約1年間お世話になったリクルートを退社し、今日からMITのJuliaLabにジョインしてJuliaの研究開発を稼業にしました。英語とか将来とか不安だらけですが、自分の好きな言語の開発に携わり自分が作りものにダイレクトにアプローチできることにとてもワクワクしています。 twitter.com/kdwkshh/status...

タグ：

posted at 19:15:46

2年前リクルートのインターンシップでなんとなくJuliaを使ってデータ分析をしていた時には、移り気な自分がこうも一つのプログラミング言語に夢中になっているとは思いもしませんでした。自分にとって何にも代え難い情熱と機会を与え続けてくれているこの言語への感謝を何か形として残したくて (続

タグ：

posted at 19:15:47

ごまふあざらし(GomahuAzaras @MathSorcerer

少額ですがJuliaのスポンサーになってみました。

またこれからは、開発はもちろんですが発信の量ももう少し増やしていきたいなと思っています。当面は来月末のJuliaConに向けて頑張るぞという気持ちですが、それ以外の日本語の媒体でも発信できそうなのでまた折に触れてお知らせできればと思います :)

タグ：

posted at 19:15:47

天むす名古屋 Temmus @temmusu_n

@kdwkshh すげーーーー．おめでとうございます!!!

タグ：

posted at 19:26:33

そういえば、今日初めて最寄駅でJuliaのTシャツ着た人を見ました！

タグ：

posted at 19:58:42

非公開

タグ：

posted at xx:xx:xx

積分定数 @sekibunnteisuu

天むす名古屋 Temmus @temmusu_n

@temmusu_n @genkuroki @_nagashimam @kamo_hiroyasu @IsshaKayo @teranishi #超算数　Z会は掛け算順序拘りですね。
www.youtube.com/watch?v=temWob... pic.twitter.com/cTIPae2TdF

タグ：超算数

posted at 20:22:24

@sekibunnteisuu @genkuroki @_nagashimam @kamo_hiroyasu @IsshaKayo @teranishi #超算数一応、2年生向けの動画なので、3年以降はこだわっていないの*かも*しれません。ただし「指導要領解説 by 文科省」(=/= 指導要領そのもの) が変わってから、Z会の方針を確認していないのは事実です。

タグ：超算数

posted at 20:24:50

#Julia言語デフォルトのソルバ(dtmaxを十分小さくすれば多くの場合に非常に優れている)ではなく、シンプレクティック・ソルバのYoshida6を試した場合も追加しました。

github.com/genkuroki/publ... twitter.com/kikumaco/statu... pic.twitter.com/ZclCb5pv7D

タグ： Julia言語

posted at 20:33:12

#Julia言語時間刻みを0.05にした場合(かなり粗い)

デフォルトのソルバでは全エネルギーが右肩下がりですが、シンプレクティックソルバのYoshida6では一定の値を保っています。どちらも定期的に全エネルギーの誤差がジャンプする。 pic.twitter.com/UYpcSjQbuI

タグ： Julia言語

posted at 20:33:15

#Julia言語時間刻み幅0.01の場合

デフォルトのソルバでは全エネルギーが右肩下がりだが、Yoshida6では誤差も小さく一定の値を保っている。

注意：Yoshida6が常にここまで上手く行くとは限らないので注意。しかし、Yoshida6はシンプレクティック法を使える場合の非常に有力な選択肢です。 pic.twitter.com/8XJbBRZMbN

タグ： Julia言語

posted at 20:33:17

Akinori Ito @akinori_ito

いま私が知ってる文系の先生方はバカ高いアイトラッカーで視線計測したりfMRIで脳機能計測をしたり巨大データを統計解析したりしている

タグ：

posted at 20:41:06

#Julia言語

Juliaでは、自分の所有物として公開するパッケージの中で、他人が定義した型のprobと自分が定義した型のalgを引数とするメソッドを他人が定義した函数solve(prob, alg)に自由に(無断で)追加できます。

常微分方程式の画期的な数値解法を開発した人は実際にそうすると良いと思います。

タグ： Julia言語

posted at 20:44:07

#Julia言語 #Jupyter

家庭内サーバー上での動いているJupyterLabサーバーにiPadのブラウザから接続して使っている様子のスクショ。 pic.twitter.com/IcFP121fhQ

タグ： Julia言語 Jupyter

posted at 21:10:43

#Julia言語 #Jupyter

家庭内サーバー上での動いているJupyterLabサーバーにiPadのブラウザから接続して使っている様子のスクショ。 pic.twitter.com/Si4zeM36EE

タグ： Julia言語 Jupyter

posted at 21:10:52

#Julia言語 #Jupyter

家庭内サーバー上での動いているJupyterLabサーバーにiPadのブラウザから接続して使っている様子のスクショ。

Command Paletteを左側に表示させたい場合にはAdvanced Settings EditorでCommand Paletteの項目をいじればよい。 pic.twitter.com/LQ4DC5XcUw

タグ： Julia言語 Jupyter

posted at 21:13:41

#Julia言語 #Jupyter

タッチで使う場合には、Jupytextを左に表示するようにしたcommand paletteから使用することになる。 pic.twitter.com/sXFnbhdpkl

タグ： Julia言語 Jupyter

posted at 21:31:24

おばけ @triwave33

ごまふあざらし(GomahuAzaras @MathSorcerer

JuliaをAmazon SageMaker上で使おう｜Shin＠Geo @gisp_shin #note note.com/shinkobara/n/n...

タグ： note

posted at 21:33:42

あ〜る菊池誠(反緊縮)公式 @kikumaco

最近, #Julia言語 #JuliaLang しかやってないがこれが
健全な人生とソフトウェアの付き合い方だと思っている．

タグ： JuliaLang Julia言語

posted at 22:11:16

雨宮処凛は山本太郎の本を編集しているくらいなので、太郎に肩入れというか、仲間なんですよね。

しかし、あの太郎の放射能デマを容認できるというのは弱者の味方とは言えないと思うんだなあ。そこんとこ、どう折り合いをつけてるのか、全くわからない twitter.com/qozop/status/1...

タグ：

posted at 22:12:58

athos)))))))) @athos0220

@valerauko SO上でのClojure関連のトラフィックが少なくなってるのはたしかだと思うんですけど、消えた言語の代わりにランキングに入った言語がHaskell/Julia/Objective-C/Perl/VBAとかで、Haskell, Juliaあたりはともかく他の3つが新たにランキングに入ってくるほど活発なのかは疑わしいかなと思ってます。

タグ：

posted at 22:23:55

ひらっち@ひよこエンジニア @xxhiratchxx

#Julia言語例えば九九の表を

(1:9)' .* (1:9)

で作れる。類似の使い方↓

using Plots
f(x, y) = exp(-x^2 + x*y/2 - y^2)
x = y = range(-3, 3; length=301)
z = f.(x', y)
surface(x, y, z; colorbar=false, size=(720, 540), camera=(60, 60), color=:gist_earth) twitter.com/sfitb/status/1... pic.twitter.com/SbksuCDadA

タグ： Julia言語

posted at 22:30:43

Juliaというプログラミング言語があるそうですね‼️
勉強不足で初めて知りました😅
Pythonとよく比較されているらしいのでちょっと調べてみます

タグ：

posted at 22:31:20

ろぼ / SchiRoboTech @SchiRobo

この記事のプログラムをJuliaに移植してみました。記事ではポアソン方程式のfとpが定数の場合を扱っていますが、任意の関数に拡張しようと思っています。

qiita.com/Sunset_Yuhi/it...

タグ：

posted at 22:36:09

#Julia言語 DifferentialEquations.jl のエコシステムに乗っかるときに、自分が開発したソルバのコードは自分のリポジトリにおけます。

DifferentialEquations.jl はPythonやRでも使えるので、ソルバの開発者は自分が開発したソルバをPythonやRでも使ってもらえるようにもなります。

タグ： Julia言語

posted at 22:41:49

#Julia言語

 sciml.ai/news/2018/04/3...
【DifferentialEquations.jl ～ allows for any developers to add dispatches to a common solve function. ～ it allows researchers to keep their research "their own" in their own repositories while contributing to this common cause. 】

タグ： Julia言語

posted at 22:46:23

#Julia言語ソルバの開発者は、DifferentialEquations.jl の開発者に自分が開発したソルバのコードを追加するプルリクエストを送ると必要はなく、DifferentialEquations.jl の開発者に無断で DifferentialEquations.jl にソルバを追加するパッケージを作って自分のリポジトリで配布できます。

タグ： Julia言語

posted at 22:49:58