黒木玄 Gen Kuroki
- いいね数 389,756/311,170
- フォロー 995 フォロワー 14,556 ツイート 293,980
- 現在地 (^-^)/
- Web https://genkuroki.github.io/documents/
- 自己紹介 私については https://twilog.org/genkuroki と https://genkuroki.github.io と https://github.com/genkuroki と https://github.com/genkuroki/public を見て下さい。
2020年09月16日(水)
@bhaveshshrima11 @rdettogni @JuliaLanguage Juno is going to be deprecated soon and the developers are joining the VScode team, so you won't miss any feature :)
タグ：
posted at 23:48:12
#Julia言語 のリポジトリでsetroundingを検索してみた。
github.com/JuliaLang/juli...
の添付画像1の部分の情報によれば、
github.com/JuliaLang/juli...
Float{32,64}のsetrounding挙動
が怪しいので削除されたらしい。
fesetroundを直接呼ぶsetrounding_rawも発見。添付画像2
github.com/JuliaLang/juli... twitter.com/mkashi/status/... pic.twitter.com/zwxB4mViV6
タグ： Julia言語
posted at 23:45:07
fesetroundに渡す定数FE_DOWNWARDとかは、/usr/include/bits/fenv.hから取ってきたけど、この値はOSやコンパイラによって違う可能性があるし、なんならコンパイラのversionで違うかもしれない。つまり、Juliaをmakeしたときに使われたlibmとversionの揃ったheader fileから取ってくる必要がある。
タグ：
posted at 23:41:46
libmの中のfesetroundを無理やり呼び出してみたら丸めモード変わってるようだが、いいのかこれ？ (多分ダメ) pic.twitter.com/msUfZ2Rdt2
タグ：
posted at 23:26:38
It takes time to build the #JuliaLang user base, but experiences like this will surely accelerate its adoption. #programming twitter.com/bhaveshshrima1...
タグ： JuliaLang programming
posted at 23:18:59
@sekibunnteisuu ある程度の量の練習を課さないとできるようにはならないことと、あまりにも量が多いと過学習が起こり他のことに応用できなくなること、両立は機械的な指導進行では無理なんでしょうね。
タグ：
posted at 23:11:33
Juliaはやっぱり丸めモード変えられないの何とかならんのかなあ。mpfrは変えられるようだが肝心のfloatとdoubleが変えられない。version 0.5くらいまでは出来てたのが消されてしまった模様。外部の変更命令をccallで無理やり呼び出しちゃいかんのかな。 pic.twitter.com/FxG097k0GD
タグ：
posted at 23:09:08
@HW_a_pythonista @Hishinuma_t 現時点では、
習得が易しい高速な言語で
高速なライブラリ群を貼り合わせて
高速性と利便性を同居させる
という選択肢が実質Juliaしかない点は問題かもしれませんが、
技術者や科学者達にCやC++やFortranの使用を強制する方向
に進むことは
絶対にやってはいけないこと
だと思います。 twitter.com/hishinuma_t/st...
タグ：
posted at 23:07:04
@temmusu_n #超算数 #数教協 は、単元ごとにシェーマを提示して、徹底的に教え込むスタイルですね。
そういう意味では、100マスやTOSSと本質的には変わらないと思う。
どれも、短期的表面的効果はあるでしょうね。
posted at 23:03:07
@sekibunnteisuu #超算数 ああ、数教協も鍛錬主義ですね。主流派の山本喜治なんかは各教科でパターン化したようです。しかも一般教員から研究を願い出た形をとるという精神主義もおまけでついてきます。
タグ： 超算数
posted at 22:50:12
つい台所で氷を作りながら「はぁーあ……」とため息をついたところ、妻から「どうしたん、はーとかため息ついて」って言われたので、なんとなく
「あーぁ、そのパレードはどこまでも続いていく……」
と言ったら、サンホラ祭りになった。
タグ：
posted at 22:47:17
@HW_a_pythonista @Hishinuma_t 関連。以下のどちらも不十分な考え方。
❌習得が易しいが遅い言語で高速なライブラリ群を貼り合わせて使えば良い。
❌真に速度を求めるなら習得が易しい言語を使うべきではない。例えばC++を使え。
個人的な意見では特に後者は時代遅れの考え方。
前者は不十分なだけで実用的な領域があると思う。 twitter.com/genkuroki/stat...
タグ：
posted at 22:43:31
@temmusu_n 数教協の例だけど、こんだけしつこく教えれば、そりゃ特定のパターンの問題は解けるようになるでしょうね。
ameblo.jp/metameta7/entr...
過剰適応してしまって、応用は効かないでしょうが。
タグ：
posted at 22:42:10
@rdettogni @JuliaLanguage I just installed the Julia extension for VS Code and went on from there. It is pretty intuitive :) I know Juno is a powerful IDE but I like VS Code a lot
discourse.julialang.org/t/julia-vs-cod...
This is dated but still pretty helpful
タグ：
posted at 22:38:09
Re: RTs 徹底的に最適化されたライブラリは速い。
だが、そういうライブラリが例えばC++でしか使えないようだと使う人は減る。
高速なライブラリを習得し易い言語から使えるようにすると、多くの末端ユーザーに使ってもらえるようになる。
ここで末端ユーザーとは優れた技術者や科学者達のことだ。
タグ：
posted at 22:30:58
支配方程式：Euler方程式
空間離散化：有限体積法
時間積分：tvd-RK3
格子点数：半径方向100点x周方向80点
計算時間：約30分
可視化：Paraview（密度コンター）
タグ：
posted at 22:19:05
2次元軸対称FVMコードで曲線格子を用いた点爆発 (Sedov explosion test) の計算をしました．対称軸付近でも擾乱は入らず，綺麗に同心円状に衝撃波が伝搬していく様子が確認できます pic.twitter.com/oOmN8PoAE1
タグ：
posted at 22:11:09
#Julia言語 メモ
github.com/JuliaGPU/CUDA....
CUDA.jl/lib/cusparse
github.com/JuliaSparse/MK...
MKLSparse.jl
以上はどちらも最近更新されている。
github.com/gridap/SparseM...
SparseMatricesCSR.jl
もある(9ヶ月間更新されていない)。
タグ： Julia言語
posted at 22:09:32
#Python #SciPy (#Julia言語)
専門家が多くのユーザー数を持つscipyなどの疎行列のライブラリに手を出してくれるのはありがたいこと。
できればJuliaの疎行列のパッケージ(検索すると色々ある)にも手を出してくれると助かる。 twitter.com/atsushi_twi/st...
posted at 22:00:49
🆕電子版もあります コロナ社
1から始める Juliaプログラミング
www.coronasha.co.jp/np/isbn/978433...
🔗電子版 販売元リンク
kw.maruzen.co.jp/ims/itemDetail...
タグ：
posted at 22:00:27
非公開
タグ：
posted at xx:xx:xx
@sekibunnteisuu ここ数年掛算順（にかぎらずだけど）の、（いつまでたっても）ぶくぶく 湧いてくる屁理屈こねの有象無象を見てると見てると#超算数っていうより、「蛸壺算数」のネーミングの方が合ってるんじゃないかと思う。
タグ：
posted at 21:13:16
@hgn_no_otaku @genkuroki いえいえ。
私は下心という言葉をよく使います。
「なぜ対数をとるかというと、これこれこういう下心があるから」。目論見というより妖しい感じがやる気を誘うかも。うひひ。
タグ：
posted at 19:56:01
@GreatDemon1701 @genkuroki ありがとうございます。
参考になります。
6^100=10^□の形を目指したり、logをとったりする動機の説明に使えそうです。
タグ：
posted at 19:47:23
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
@HShinaoka おおっと知らなかった
いまCSRの機能だけ一通りばーっとみてきましたが，ここのコードしかないですね．．．
対角とってくる１機能だけ．．．
github.com/scipy/scipy/bl...
タグ：
posted at 18:33:00
@HShinaoka すみませんｗ一点の曇りもないのは嘘ですｗ
固有値なんかはARPACKコールしてますね
ただsparse blasとutil系は全部.pyのはず
タグ：
posted at 18:28:17
I decided to work on Javis transposing a matrix again today instead of Project Euler.
Starting in less than 6 hours!
Follow me on @Twitch to get notified:
twitch.tv/opensourcesblog
#Julialang #animation #live
Ideas on how to compare shapes are welcome :) pic.twitter.com/bssrJPzUqx
posted at 18:08:46
@hgn_no_otaku @genkuroki 既に説明されたあとのトピックだと思いますが、
私は
log10=log10¹=1
log100=log10²=2
log1000=log10³=3
・・・
「1の後に続く0の個数と同じだよね」とか、「log200は2てんいくつかになるはずだよね」とか、「6¹⁰⁰が10のナントカ乗なのかを調べるんだよ」とか言いますね。でも対数で進めるかな。
タグ：
posted at 18:04:15
ちなみに何度か言っているけど最も害悪なのはscipy.sparseです
この一点の曇りもないpythonコードの山を見て皆さんと是非とも絶望を共有したい
github.com/scipy/scipy/tr...
タグ：
posted at 18:00:02
HW_a_pythonista @HW_a_pythonista
@Hishinuma_t 速さ（実行速度）だけを求めるなら、Pythonは投げ捨ててください。応用の問題は、便利なインターフェースになっているか、どういうSolverを使うか、などを試行錯誤する必要があります。その際、素のPythonで多重ループを回したりしなければ、速度低下は、一定に抑えられます。C++やFortranとは協調。
タグ：
posted at 17:55:11
私の主張は「速さを気にするならPythonなんか投げ捨てろ」ということ．
ライブラリを呼べば速いというけどnumpy, scipyは純粋なBLAS / LAPACKのラッパーではなく独自でpythonで作っているものもある．
「ライブラリを使えばPythonは速い」という表現は誤解を生むのでやめたほうが良いと思う
タグ：
posted at 17:46:53
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
知らん記法いっぱい出てきたので
Juliaのpunctuationのページを眺めている
docs.julialang.org/en/v1/base/pun...
タグ：
posted at 16:51:01
@takasan_san_san #Julia言語 失礼しました。FFTW.jl はすでに発見済みでしたね。
私によるFFTW.jlの使用例が
nbviewer.jupyter.org/gist/genkuroki...
nbviewer.jupyter.org/github/genkuro...
にあります。 twitter.com/takasan_san_sa...
タグ： Julia言語
posted at 16:48:15
@takasan_san_san #Julia言語 でFFT情報
twilog.org/genkuroki/sear...
nbviewer.jupyter.org/gist/genkuroki...
github.com/JuliaMath/FFTW...
juliamath.github.io/FFTW.jl/stable...
(cos, sin 変換)
pkg> add FFTW
julia> using FFTW
すれば、fftなどを使えるようになります。
タグ： Julia言語
posted at 16:45:32
高三 和晃 / Kazuaki Taka @takasan_san_san
JuliaでFFTW.jlというのは見つけたけど、どうなんだろう？github.com/JuliaMath/FFTW...
タグ：
posted at 16:38:19
だーすー:||（Hiroshi Suda @suda_hiroshi
@shigerufujita バージョン（1.6.0-DEV.917 (2020-09-14)）を見て驚きました．攻めますね．
タグ：
posted at 16:35:45
高三 和晃 / Kazuaki Taka @takasan_san_san
f(t_n) みたいな離散的な時間上のデータをFourier変換したいんだけど、Juliaでそういうのにあったパッケージとかあるかな？
タグ：
posted at 16:35:14
【速報】宮城で新たに6人感染確認 3人は仙台市 累計337人に（KHB東日本放送）
#Yahooニュース
headlines.yahoo.co.jp/hl?a=20200916-...
タグ： Yahooニュース
posted at 16:21:55
Furniture lifters that use rollers can generally help move furniture and household appliances up to 200 kg and can be used on on hard wood floors, laminates and carpets buff.ly/3hBUNtW pic.twitter.com/VebUEtsFHM
タグ：
posted at 16:09:44
#Julia言語 ソースコード
using AbstractTrees
AbstractTrees.children(x::DataType) = subtypes(x)
print_tree(Number)
はよく使う。
nbviewer.jupyter.org/gist/genkuroki...
タグ： Julia言語
posted at 15:32:26
起きても興味を失っていなかったので、とりあえず出版物の総額表示義務化に反対しておきますね。出版物のみならず、ソフトウェアをパッケージ化した商品で、長い期間をかけて売る物にこれやっちゃいかんと思います。
タグ：
posted at 15:15:18
The full playlist for "Julia for dplyr useRs" is here
www.youtube.com/playlist?list=...
タグ：
posted at 15:10:10
There are many data manipulation packages/ecosystems in Julia? Don't know what to choose?
I have a video for you where I show you my favourite #julialang data manipulation packages.
www.youtube.com/watch?v=q_P2H_... (~6mins)
DataFrames.jl + DataFramesMeta.jl + DataConvenience.jl
タグ： julialang
posted at 15:10:10
JuliaでどのNNのライブラリを使えばいいのかよくわからなくなったのでスクラッチでBPTTするコードを書いた（元はPython）
compneuro-julia.github.io/chap7/7-9_bptt...
タグ：
posted at 15:04:58
@Dsuke_KATO "pure Julia" にこだわることは不合理なイデオロギーに過ぎないので、#Julia言語 から楽に使えるものはありがたく使わせてもらうのが、合理的なJuliaユーザーだと思います。
Juliaは糊言語。
ImageMagick.jlで使えるImageMagickが古過ぎという問題は報告した方がいいかも。
github.com/JuliaIO/ImageM...
タグ： Julia言語
posted at 15:03:21
The Nikon Small World in Motion competition brings together talented microscopists from all over the world. What they have been able to capture will likely astound you, even though sometimes what you're seeing is pond scum [full video: ] pic.twitter.com/mIePEMkLOd
タグ：
posted at 14:55:54
Mark Kittisopikul ht @markkitti
@genkuroki @mat_der_D julia> X = rand(0:9,5)
5-element Array{Int64,1}:
3
5
2
5
1
julia> repeat(X,3)
15-element Array{Int64,1}:
3
5
2
5
1
3
5
2
5
1
3
5
2
5
1
タグ：
posted at 14:42:04
#Julia Is Faster and Better than #R, #Matlab and #Python for #Economic #Research: Alvaro Aguirre and Jon Danielsson of the London School of Economics Systemic Risk Centre benchmarked #Julialang against R, Matlab and Python for economic research. Read more
voxeu.org/article/which-...
タグ： Economic Julia Julialang Matlab Python R Research
posted at 14:38:02
せんせーJuliaでJPEG-2000を読み込みたいです。PyCall->opencv-pythonというルートもあるのですが、Pythonアレルギーが出そうで・・・。
タグ：
posted at 14:04:37
@mat_der_D #Julia言語
① vcat(fill(X, n)...)
② Xをn個連結したものを改変せずに使ってかつ、メモリ効率の問題が気になるならば、中身がXだけのstructを作って～以下略
類似例が
nbviewer.jupyter.org/gist/genkuroki...
のMy. Catにあります(nが大きいと遅い)。Iterators.flattenも選択肢になり得る。
タグ： Julia言語
posted at 13:50:17
- Union型：複数の型を一つにまとめるもの docs.julialang.org/en/v1/manual/t...
- NTuple{N, T} は Tuple{Vararg{T,N}} のaliasで，Tの型をちょうどN要素含むタプルを定義するときに便利 docs.julialang.org/en/v1/manual/t...
- whereを使った型制約の記述（例：Array{T} where T<:Integer） docs.julialang.org/en/v1/manual/t...
タグ：
posted at 13:47:06
Julia で array X が与えられているときに、例えば縦に3つつなげたもの作るなら [X; X; X] って書けるけど、縦に n 個つなげたものを作りたい場合ってどう書けばいいんだ？
タグ：
posted at 13:38:33
#Julia言語 と #R言語 の比較についてはまず
juliadata.github.io/DataFrames.jl/...
を見て、R dplyr と Julia DataFrames.jl のコマンド対応表を眺めてみると良いかもしれません。
Juliaだけでも大抵のことができて、ユーザーが書いたコードが速く、必要があればRやPythonをJuliaで利用できます。
posted at 13:35:26
ぷりんちゃん先輩のスペクトル法ライブラリ（in #Julia言語 ）のソースを読んでいくぞ
github.com/mat-der-D/Spec...
タグ： Julia言語
posted at 13:35:19
One week of using @JuliaLanguage and I am a BIG fan. Hacked together a FE code using Gridap.jl. Extremely clean, expressive, performant and concise all at the same time. Boy those math symbols are so cool😉.
mybinder.org/v2/gh/bhaveshs... (if you want to play with it on binder) pic.twitter.com/tmDtGXAwHq
タグ：
posted at 13:26:06
#Julia言語 例えば、Juliaでの微分方程式に関する超優れもののパッケージであるDifferentialEquations.jlはRでも使用できます(Pythonでも)。
CやC++よりもJuliaの方が書き易いのでこういうのが広まるとよいと思う。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 13:23:22
#Julia言語
すでにRからJuliaを使うこともできて、改良されて使い易くなって行くと予想されます。
Rで使うライブラリで重い計算が必要なものについては、「Juliaで書いてRで使えるようにする」という選択肢が有力になって来ているかもしれません。
cran.r-project.org/web/packages/J... twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 13:15:42
言語のお勧めはよく訊かれる。
いつもjuliaを勧めてたんだけど誰も使わない。
みんな「自分のやりたい言語を勧めてくれる」人を探してるだけだなって気づいたので、何を訊かれても「なでしこでいいよ」って答えてる。
「～はどうですか」って言われたら「本質的にcと変わらんからなでしこで」と。
タグ：
posted at 13:09:21
自分では使わないのだが，#情報工学科 の学生は知っておいた方が良いし，4つ位の研究室では，既存言語やツールよりも #julia #juliaLang の方が良いのではないか？と思う．
そんな訳で，1,2年生用の環境を更新した． pic.twitter.com/sFeLdQQOoq
posted at 13:03:20
#Julia言語 翻訳【Juliaは速い～しかし、速いことの真に強力な帰結は、ユーザーがコントリビューターになるために、[CやC++のような相対的に習得が面倒な]他言語を学ぶ必要がないことである～Juliaを知っていればJuliaを使うためには十分である[.hや.cに触れる必要はない]】
mdneuzerling.com/post/first-imp... pic.twitter.com/k9JTi9FDZU
タグ： Julia言語
posted at 13:01:39
#Julia言語 遊び
Diniの曲面
3D曲面plotは、格子網を捩じったヤツの交点xyz座標を
それぞれ2次元配列で格納させるかんじ
GRバックエンドは対応してないっぽいな😢
gist.github.com/pppp314592/df7... pic.twitter.com/2rq65Fe4KK
タグ： Julia言語
posted at 13:00:37
#Julia言語 大反響の以下のリンク先のブログ
mdneuzerling.com/post/first-imp...
【Julia is fast, and the compiled code is often very similar to C. But a really powerful consequence of this is that a user doesn’t need to learn another language to become a contributor.】
本質を突いている！ twitter.com/mdneuzerling/s...
タグ： Julia言語
posted at 12:49:25
#Julia言語 Juliaの側からRを使うこともできます。
スクショのコードはよく見ないとRを直接に使っていると誤解してしまうかも。Julia内でRを使っています。
こういう感じのことたあちこちでできるようになるとみんなハッピーになれると思う。
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/m1VUbKg22Z
タグ： Julia言語
posted at 12:21:55
#Julia言語 最近、ツボにはまったJuliaネタは、#R言語 (#RStats)からJuliaを使うためのRライブラリのリポジトリで、JuliaをRのライブラリとして自動的にインストールするようにしてしまうことが提案されていること。
github.com/Non-Contradict...
Automatic Installation of Julia? #116
posted at 12:02:01
@TaleOfDuckedon @antimon2 #Julia言語 こういうのはJuliaのソースコードを見て、自分で小さなコードを書いて確認するのが速いと思います。
github.com/JuliaLang/juli...
Juliaのリポジトリでpromoteを検索
docs.julialang.org/en/v1/manual/c...
公式ドキュメントでのpromotionに関する解説
タグ： Julia言語
posted at 11:31:57
ごまふあざらし(GomahuAzaras @MathSorcerer
数学セミナーのよこえもんさんの語り口がもしろかったできゅ #Julia言語
タグ： Julia言語
posted at 11:28:49
Great post highlighting some of #julialang features! twitter.com/mdneuzerling/s...
タグ： julialang
posted at 10:28:00
@TaleOfDuckedon @genkuroki 動作原理については @genkuroki さんが詳解しているのでほぼ言うことないですが補足として、Juliaにはそもそも型キャストという概念がありません。型変換は暗黙には行われず、コードでルールが決められており（promote()関数等）、+-*/等はそれらを利用して型を合わせてから演算してます。
タグ：
posted at 09:32:38
David Neuzerling @mdneuzerling
@AriKatz20 I had no idea! I’ll update the post and credit you when I get a second. Thank you!
Looking at this blog post, it seems like the dot prefixes operators and suffixes functions? That @. looks incredibly helpful by the way.
タグ：
posted at 08:57:11
@mdneuzerling Welcome and great post! Just wanted to add that unlike the ! , the dot is less a convention and more a language feature that actually takes a function and vectorizes it on CPU, GPU etc. Aside from the convenience this has other advantages: julialang.org/blog/2017/01/m...
タグ：
posted at 08:54:13
"a user doesn’t need to learn another language to become a contributor."
Nicely put.
#JuliaLang twitter.com/mdneuzerling/s...
タグ： JuliaLang
posted at 08:44:48
David Neuzerling @mdneuzerling
✍️ First Impressions of Julia from an R User
Julia is a language that's both fast and a delight to use. Here are my first impressions, coming from an R perspective. Install Julia; you won't regret it. mdneuzerling.com/post/first-imp...
タグ：
posted at 07:34:51
#統計 自由度はパラメータ数の差という理解の方針は私だけが言っていることではない。Goodさんも同じことを言っている。
詳しい解説が欲しい人は私のツイログを検索
twilog.org/genkuroki/sear... twitter.com/genkuroki/stat...
タグ： 統計
posted at 06:07:24
@TaleOfDuckedon #Julia言語 f(x::Real, y::Integer) が定義されていても、RealとIntegerは抽象型に過ぎないので、Real型とInteger型の組み合わせを扱うネイティブコードにコンパイルされることはありません。
函数の実行時に判明している引数の具体的な型を扱うネイティブコードにコンパイルされるだけです。
タグ： Julia言語
posted at 05:03:26
@TaleOfDuckedon #Julia言語 函数の引数の具体的な型の組み合わせごとに、同じコードが異なるネイティブコードにコンパイルされる。
例：f(x,y)=x*y について、f(2,3) を実行するときには x,yがどちらもInt64型の場合のネイティブコードが生成され、f(2.0, 3.0)ではFloat64型を扱うネイティブコードが生成されます。
タグ： Julia言語
posted at 05:00:16
@TaleOfDuckedon #Julia言語
(1) Int64<:Real かつ Int64<:Integer
(2) Int64<:Integer かつ Int64<:Real
のどちらも真なので、(1),(2)のどちらのfを実行するかについて曖昧さが生じています。
Int64, Int64の組は決して他の型に変換されたりしません。
タグ： Julia言語
posted at 04:55:02
@TaleOfDuckedon #Julia言語 あともう1つ言っておくと、Realは抽象型であり、Frost64などの浮動小数点数の具体的な型ではないです。だから、Int64をRealに変換という見方そのものがおかしいです。
12と34の具体的な型はどちらもInt64です。
このInt64を使ってdispatchが行われます。続く
タグ： Julia言語
posted at 04:51:41
Im(Nyarn)＠linuxﾅﾆﾓﾜｶ @Imaginary_Nyarn
PyCallではなくrclを直に使ったライブラリならJuliaでロボットをPythonを使わずにまともに動かせそうなのでうれしい twitter.com/Imaginary_Nyar...
タグ：
posted at 04:25:55
Im(Nyarn)＠linuxﾅﾆﾓﾜｶ @Imaginary_Nyarn
ROS2でも使えるようになりそう
github.com/jdlangs/RobotO...
タグ：
posted at 04:24:06
Im(Nyarn)＠linuxﾅﾆﾓﾜｶ @Imaginary_Nyarn
JuliaのROS1クライアントライブラリ
github.com/jdlangs/RobotO...
タグ：
posted at 04:22:44
#Julia言語
sizeをローカル変数として使いたくなる場合が生じている第一の理由は、Plots.jlにおける作画サイズを意味するキーワード引数の名前がsizeであること。
例
using Plots
function plotsin()
size = (600, 300)
plot(sin; size)
end
plotsin() pic.twitter.com/UyJhR5TgvI
タグ： Julia言語
posted at 03:58:02
Julia Is Faster and Better than R, Matlab and Python for #Economic #Research: Alvaro Aguirre and Jon Danielsson of the London School of Economics Systemic Risk Centre benchmarked #Julialang against R, Matlab and Python for economic research. Read more voxeu.org/article/which-...
タグ： Economic Julialang Research
posted at 03:37:24
#Julia言語 メモ
github.com/JuliaSparse/MK...
MKLSparse.jl is a Julia package to seamlessly use the sparse functionality in MKL to speed up operations on sparse arrays in Julia. In order to use MKLSparse.jl you do not need to install Intel's MKL library...
タグ： Julia言語
posted at 03:37:07
@DocSparse @acidflask @tymwol I think @Viral_B_Shah (whose retweet bought this thread to my attention) would be the best person to discuss this with.
タグ：
posted at 03:20:35
@DocSparse @acidflask @tymwol It'd be interesting to see a list of libraries that do parallel sparse matrix multiplication on the CPU. I know of CombBLAS and SuiteSparse. Rust's sparse matrix library, sprs, just got one this summer. But that's all I'm aware of. github.com/vbarrielle/spr...
タグ：
posted at 02:03:20
@TaleOfDuckedon #Julia言語 これは実際にJuliaのコードを書いていると確かにうざいです。冷静に集合よりも複雑なタイプのインターセクションを考えなければいけない。
その辺については
youtu.be/TPuJsgyu87U
がおもろいです。
タグ： Julia言語
posted at 01:55:32
@genkuroki なるほど！（エラーメッセージちゃんと読めてませんでした）
他のプログラム言語だと、引数の順番にどの関数を当てはめるか評価が行われるものもある（→(2)が選ばれる）のですが、Juliaではエラーになるんですね！
タグ：
posted at 01:53:17
Dr. Chris Rackauckas @ChrisRackauckas
@rlmcelreath @JuliaLanguage And here's how you solve ODEs with rational numbers:
tutorials.sciml.ai/html/type_hand...
Note it's the standard rational number with the standard ODE solver, just slapped together. Multiple dispatch lets these tools compose.
タグ：
posted at 01:48:02
@TaleOfDuckedon #Julia言語
ℝ×ℤの上で定義された函数fと
ℤ×ℝの上で定義された同じ名前の付いた別の函数fが
与えられているときに、
ℤ×ℤの上での値をどちらのfで決めるかはきまっていないので、
エラーを出さざるを得なくなっているわけです。
タグ： Julia言語
posted at 01:46:03
Kohei.Sakamoto @KoheiSakamoto88
Python に慣れすぎてこのR のコードはムズムズする...
パッケージ名書かないと名前が衝突しそうで嫌だ.
一方でjulia ぐらいスッキリしてると読みやすい.
#python #rstatsj #julia twitter.com/tymwol/status/...
posted at 01:42:52
@TaleOfDuckedon #Julia言語
typeof(34) == Int64 <: Real はJuliaにおいて常にtrueです。
添付画像でエラーが出ているのは、34もReal型の一種になっていることをJuliaが認識していないからではなく、エラーメッセージに書いてあるように、(1)と(2)のfのどちらを実行するべきかが曖昧になっているからです。 pic.twitter.com/J0UNyHTiNq
タグ： Julia言語
posted at 01:42:14
@DocSparse @acidflask @tymwol I don't think there is one. My benchmarks (last summer) used a single thread for sparse matrix operations. There seems to have been an attempt to implement it but it was dropped ~5 years ago: github.com/madeleineudell...
タグ：
posted at 01:40:44
@acidflask @DocSparse @tymwol The SparseArrays.jl package handles sparse matrices (and allows overriding the plus/times operators) but as far as I'm aware, its matrix multiplication operations are single-threaded. MKLSparse.jl has parallel matvec according to this forum: discourse.julialang.org/t/parallel-spa...
タグ：
posted at 01:15:34
@genkuroki twitter.com/genkuroki/stat...
あ、すみません…。
演算結果ではなく、整数リテラルのキャストのことです。
上のTweetのf(12,34)の結果がエラーになるところで、Int64＜Realなので第1引数の12がInt64で処理されるとして、第2引数の34は暗黙的にRealにはならないのか、という。
タグ：
posted at 01:07:01
#Julia言語
gist.github.com/genkuroki/2dd0...
1. 整数と浮動小数点数の演算は整数を浮動小数点数に変換してから計算される。
2. 整数の配列の成分に1.5を代入しようとすると、1.5を整数に変換しようとしてエラーが出る。
3. 整数を整数で/で割ると浮動小数点数になるので、配列のインデックスで使えない。 twitter.com/taleofduckedon... pic.twitter.com/LFAjonwDRm
タグ： Julia言語
posted at 00:46:21
@rlmcelreath @JuliaLanguage I like that blue background as well! Looks great. #JuliaLang
タグ： JuliaLang
posted at 00:35:00
#Julia言語 私は、函数の中で局所的に size という名前のローカル変数を使ったことがあります。
Juliaのビルトインsizeその函数内で使いたくなっても、Base.sizeとして利用できるので困らない。 twitter.com/soyukkea/statu...
タグ： Julia言語
posted at 00:31:23
＞#Julia言語 ではコードを変更せずに疎行列も扱えます。
続き。蜜行列を扱う場合から疎行列を扱う場合に移行するために
* ユーザー側がコードを書き換える。
または
* コードを書き換えないために、ライブラリのインポートの仕方を変える。
のどちらかを強制するのはよくない。
タグ：
posted at 00:23:57
#Julia言語 最後の4の例では、初期条件が区間の場合の数値解が、オイラー法の同じ函数で扱えているだけではなく、グラフのプロットも含めて非常に簡単にできていることに注目。
他人の仕事を自分が書いたコードで非常に簡単に利用できるのは気分が良いです。 pic.twitter.com/IxlSm2DWpf
タグ： Julia言語
posted at 00:17:59
#Julia言語 別の例
nbviewer.jupyter.org/gist/genkuroki...
添付画像
1. シンプルなオイラー法による微分方程式の数値解を求める函数
2. 単振り子
3. 行列 U(t) に関する dU(t)/dt = A U(t) 型の微分方程式
4. 初期条件が区間の場合
2だけではなく、行列の場合の3、さらに4も1の函数で扱えている！ pic.twitter.com/zIj8BcWKow
タグ： Julia言語
posted at 00:15:24
Re: RT うけた(笑)
＞どれが疎行列を扱っているんだ？
＞#Julia言語 ではコードを変更せずに疎行列も扱えます。
Juliaのmultiple dispatchでは異なる型の対象を同じ名前の函数・同じ記号の演算子で扱える。そのおかげで、ユーザー側は完全に同じコードで異なる仕事をできる。
タグ：
posted at 00:00:07
