Tweet

More from @genkuroki

黒木玄 Gen Kuroki

@genkuroki

2 Apr

#Julia言語これも面白いです。

LispとJuliaの比較のような話になっている？

以下のリンク先はJuliaでの高速多項式評価が何をやっているかの解説。係数が格納された配列を読みながらの多項式評価より速いです。

discourse.julialang.org/t/cas-benchmar…
CAS benchmarks (Symbolics.jl and Maxima)

#Julia言語【ネタ】重要情報(笑)

discourse.julialang.org/t/cas-benchmar…
【Maybe julia stands for “Jeff’s uncommon lisp is automated”?】とJeff Bezansonさんが言っていた！

#Julia言語上で紹介したスレッドは以下に紹介するスレッドから派生したっぽい。

以下のリンク先はJuliaが何をやっているかの易しい解説になっている発言を選びました。

どういう人がやって来ても、技術的な事柄を易しく解説しようとすることはやはり大事なことだな。

discourse.julialang.org/t/cas-best-pra…

Read 7 tweets

黒木玄 Gen Kuroki

@genkuroki

2 Apr

@backflip1801

ツツカナ vs. 森下卓電王戦リベンジマッチ shogidb2.com/games/07219818… #将棋 #将棋DB2 @backflip1801より

電王戦公式統一パソコン「GALLERIA電王戦」森下卓九段インタビュー - ヒューマンエラー対策森下新ルールなら技術で勝てる (1) | マイナビニュース news.mynavi.jp/kikaku/2014120…

指掛けの電王戦リベンジマッチは指継がず、森下卓九段の判定勝ちと発表。ニコ生アンケート結果は良かった85% shogi1.com/denousen-reven…

Read 5 tweets

黒木玄 Gen Kuroki

@genkuroki

2 Apr

https://twitter.com/f_sei/status/1377762238946144259

#数楽大学1年生向けの線形代数の講義では、与えられた a_{ij}とb_iが定めるx_j達に関する連立一次方程式

a_{11}x_1 + … + a_{1n}x_n = b_1
　　………
a_{m1}x_1 + … + a_{mn}x_n = b_m

をa_{ij}, b_iに関する場合分けで解く話をやったりします。

ax=bはこれのm=n=1の特別な場合。

https://twitter.com/f_sei/status/1377762238946144259

#数楽そのm=n=1の特別な場合の方程式

ax=b

は以下の3種に分けされる。

(0) a≠0の場合の ax=b
(1) a=0, b≠0の場合の 0=b
(2) a=0, b=0の場合の 0=0

(0), (2)の場合に解空間の次元はそれぞれ0, 1になり、(1)の場合に解空間は空集合になる。

これの一般化をやる。

#数楽 m=n=2の場合の方程式の中には

x + ay = p
0x + 0y = q

のようなものも含まれる、q≠0の場合は解なしになり、q=0の場合には

x = p - ay
yは任意

の形ですべての解を表せる。

こういう講義をしている人達はそれなりにいるはず。

Read 5 tweets

黒木玄 Gen Kuroki

@genkuroki

2 Apr

#数楽 #Julia言語新学期の講義のためのC^∞級函数の例の準備で、

f(x) = if x > 0 exp(-1/x) else 0 end

のような函数を扱っている人もいるはず。

現時点で最もよく使われているFloat64型の浮動小数点数の精度では

exp(-1/0.001) == 0

は true になります❗️(笑)

こういうのもネタになりそう。

#Julia言語

f(x) = if x > 0 exp(-1/x) else zero(x) end

using ForwardDiff
f′(x) = ForwardDiff.derivative(f, x)

のグラフ。最近では自動微分が普通で、数値的に微分を求めるときにも、導函数の公式を入力せずに、「正確に」かつ「効率良く」微係数を計算できます。

gist.github.com/genkuroki/ce92…

#Julia言語

f(x) = if x > 0 exp(-1/x) else zero(x) end
g(x) = if x > 0 exp(-1/x)/x^2 else zero(x) end

using ForwardDiff
f′(x) = ForwardDiff.derivative(f, x)

のgとf′の値と計算速度の比較。

導函数の公式を私が入力したgとForward.derivativeで自動微分しているf′で計算速度はほぼ同じ。

Read 7 tweets

黒木玄 Gen Kuroki

@genkuroki

2 Apr

https://twitter.com/Yosymatu/status/1377767949105127430

#Julia言語マクロは難しいです。Symbol型とExpr型がどう使われているかのイメージがないとつらい。そこでイメージを得るためのパッケージを作って以下の場所で公開しています。

github.com/genkuroki/Meta…
↓
github.com/genkuroki/Meta…

AbstractTrees.jlが神。それを使ってExprをツリー表示。

https://twitter.com/Yosymatu/status/1377767949105127430

#Julia言語イメージを得るためには、

e = :(x = 1 + 2)
dump(e)

とかしてみるとよいです。dumpの結果は見難いですが、情報がフルに入っている。見易く気軽に表示したければ、MetaUtils.jl を使えばよい。

#Julia言語自分でMetaUtils.jlを書くまで、Juliaのマクロを自分で書くことをできませんでした。

まるで「鶏と卵」。数学的な複雑さがある程度以上のものの理解ではほぼ常に「鶏と卵」状態になる。

親子丼が実は大好物。

Read 5 tweets

黒木玄 Gen Kuroki

@genkuroki

1 Apr

#Julia言語小ネタ

isless と < の違い。

それぞれの公式解説を見れば違いは明らかなのですが、違いに気付いていないユーザーは多いと思う。

< は isless にfallbackするので違いは小さいのですが、isless は sort で使用可能な全順序。

docs.julialang.org/en/v1/base/bas…
docs.julialang.org/en/v1/base/mat…

#Julia言語自分で定義した型に Base.isless を定義するときには、sort で使える全順序となるように定義するべき。

それに対して、Base.:< は全順序である必要はない。もちろん、islessにフォールバックするので、別に定義する必要が無ければ、全順序の Base.isless だけを定義しておけば足りる。

#Julia言語全順序のislessとは別に、全順序ではない < が必要なケースは、浮動小数点数の大小関係。

浮動小数点数は NaN を含んでおり、例えば、0.0 < NaN や NaN < 0.0 は false になって欲しく、実際にそうなるように定義されている。

全順序の isless で isless(Inf, NaN) は true になる。

Read 9 tweets

Share this page!

黒木玄 Gen Kuroki

Try unrolling a thread yourself!

More from @genkuroki

黒木玄 Gen Kuroki

黒木玄 Gen Kuroki

黒木玄 Gen Kuroki

黒木玄 Gen Kuroki

黒木玄 Gen Kuroki

黒木玄 Gen Kuroki

Did Thread Reader help you today?

Like this author's thread?