黒木玄 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 を見て下さい。
2021年06月16日(水)
Juliaでフーリエ級数展開のプロット
第1,2,10,100部分和です。
授業で習ったことを応用して色んな関数で試してみました。 pic.twitter.com/Rg4Bg6NUhV
タグ：
posted at 01:08:12
#Julia言語 リンク先ではPackageCompiler.jでsysimageを作り直すことによるJuliaスクリプト起動の遅延削減について説明しました。
欠点：sysimageを作るのに3分もかかるのに、遅延は数秒～十数秒のオーダーでしか減らない。
このスレッドではパッケージのprecompileを使う別の方法を説明しましょう。 twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 01:12:11
#Julia言語 ある程度Juliaに慣れている人ならば
github.com/genkuroki/publ...
github.com/genkuroki/publ...
とスクショを見るだけで方法を理解できるはず。
パッケージのprecompileを強化することによって、heateq_calc.jl の実行時間を6～7秒から、4.7秒程度に縮めることに成功しています。 pic.twitter.com/2DMzNbLttV
タグ： Julia言語
posted at 01:17:05
#Julia言語 そのためには、Juliaスクリプトの中身の大部分をパッケージに移動して、パッケージならばスクリプトと違ってprecompile結果がキャッシュされることを使えばよいです。
熱方程式を数値的に解くコードが入っているパッケージのソースを
github.com/genkuroki/publ...
で読めます。続く
タグ： Julia言語
posted at 01:22:51
#Julia言語 その結果、Juliaスクリプト本体は
using HeatEqCalcOld
sol = HeatEqCalcOld.calc_sol()
HeatEqCalcOld. save_sol(sol)
の3行だけになっています。このスクリプトはそのスクリプトがおいてある場所で
julia --project=@. heateq_calc.jl
と実行します。
github.com/genkuroki/publ... pic.twitter.com/hBFQKe1Qfr
タグ： Julia言語
posted at 01:24:55
#Julia言語 その結果、1回目の実行ではパッケージHeatEqCalcOldのprecompileが実行され、2回目以降は0.6～1秒程度実行時間が短くなっています。 pic.twitter.com/ANiaYI9M9d
タグ： Julia言語
posted at 01:26:33
#Julia言語 さらに改良したバージョンが
github.com/genkuroki/publ...
においてあります。そこでは、heateq_calc.jlの実行のために必要なprecompile達のリストtrace.jlを
julia --project=@. --compile=all --trace-compile=src/trace.jl heateq_calc.jl
で作成して、パッケージで読み込んでいます。 pic.twitter.com/iTWb7ds6vl
タグ： Julia言語
posted at 01:29:02
#Julia言語 これによるパッケージのprecompileの強化の結果は劇的で、2回目以降のスクリプトの実行時間が4.7秒程度になっています。
改良前は6秒台でした。
この改善がうれしい人は以上と同様のことを自分でも実験して確認するとよいと思います。 pic.twitter.com/Jq8kCDZm9T
タグ： Julia言語
posted at 01:31:24
@dlfivefifty My favorite part is you can get any # of dims, with the option to lay out rows horizontally OR write everything in memory order.
タグ：
posted at 01:36:21
Charles Gaba isn't p @charles_gaba
Connecticut. The state only has 8 counties and *all* of them voted for Biden, so this doesn't really mean much, but even here there's a slight slope. pic.twitter.com/1cTv7u0Qpk
タグ：
posted at 02:12:53
#Julia言語
discourse.julialang.org/t/proper-modul...
【“goto definition” etc. in vscode is not longer working】
と言っている人がかなり攻撃的な態度を取っていると思うのですが、スクショを見ればわかるように、私の環境では VS Code での goto definition はしっかり動いています。 pic.twitter.com/RlWcAPVIPf
タグ： Julia言語
posted at 02:13:18
#Julia言語 私の環境は以下の通りです。
* Windows 10 Home
* Julia v1.6.1
* VS Code Julia Language Support (Insider) 1.2.4
現時点での私個人の不満はv1.7以上でLanguage Serverが使えないこと。安定版じゃないので仕方がない。待っていれば解決が予想される問題。
twitter.com/genkuroki/stat...
タグ： Julia言語
posted at 02:16:38
discourse デビューがいきなりバトルになるのは嫌だなあ。
どんなことであっても、「仮に○○について背に腹はかえられぬ状況に追い込まれたら、技術的にどういう選択肢があるか？」という方向の議論にすれば、白熱してもその内容は技術的な話題になって楽しいと思う。
タグ：
posted at 02:26:19
Charles Gaba isn't p @charles_gaba
Delaware, Hawaii AND Rhode Island: There's only 13 counties between all 3 states, and one of those is Kalawao, Hawaii, which only has 86 residents (none of whom appear to have been vaccinated yet), so it seems kind of pointless to look for any "trend lines" at the county level. pic.twitter.com/RgqhyeHhnp
タグ：
posted at 02:29:37
#Julia言語 #数楽 Gibbs現象が見えていて楽しい。
フーリエ級数の部分和を素直にk=1,2,…,nに関する和にせずに、例えばm≠nの場合のk=1,3,…,2m-1とk=2,4,…,2nに関する和などに変更すると、妙なところに「棘」ができて面白い場合がある。
変なことを色々やると経験値がたまりやすい。 twitter.com/1_1_2_3_5__/st...
posted at 02:40:08
Charles Gaba isn't p @charles_gaba
📣 Since this thread has gone fairly viral, I wanted to remind everyone that the vaccination levels shown are based on the TOTAL population of each county/state, not just adults. If you only include adults, every county would rank somewhat higher regardless of partisan lean.
タグ：
posted at 02:43:12
Charles Gaba isn't p @charles_gaba
📣 Also, MOST states use the % of total population *fully* vaccinated only (2 shots of Pfizer/Moderna or 1 shot of J&J). There's a few states where I'm using *total doses administered* (divided into 2x the population, w/J&J counting twice) depending on what data is available.
タグ：
posted at 02:45:10
Charles Gaba isn't p @charles_gaba
The states using total doses administered will show slightly higher rates across the board than if I'd used fully-vaxxed only, but those are still consistent across every county regardless of partisan lean (ie, the blue counties might be ~5% higher but so would the red counties).
タグ：
posted at 02:47:01
Charles Gaba isn't p @charles_gaba
I still have 18 states left to add, which I'll get to later on today if possible, including: IN, IA, KS, KY, ME, MD, MO, NE, NH, ND, OH, OK, PA, SC, SD, VA, WV, WI
Stay tuned, and remember you can support my work here: acasignups.net/support
タグ：
posted at 02:48:30
@charles_gaba @ForecasterEnten I pulled Florida this weekend. It’s not as straightforward pic.twitter.com/BiHcxNwMPL
タグ：
posted at 02:50:34
First YouTube video created with Javis www.juliabloggers.com/first-youtube-... #juliabloggers
タグ： juliabloggers
posted at 02:56:46
#JuliaLang 3×3×2 Array
a = [1;2;3;;4;5;6;;;7;8;9;;10;11;12] pic.twitter.com/R7nEOIvhRR
タグ： JuliaLang
posted at 03:11:38
#数楽 セガ社に対する好印象がうなぎのぼり！
四元数(quarternion)による回転の記述(実質的にSU(2)による記述)は、3×3の直交行列による記述よりも数学的に一段深い位置に住んでいます。(物理的にも一段深い) twitter.com/sega_official/...
タグ： 数楽
posted at 03:19:25
Julia GPU Programming with WSL2: A walkthrough on how to get started with #Julialang #GPU #programming under WSL2 on Windows. Blog by Fabian Becker (@geekproject) geekmonkey.org/julia-gpu-prog...
タグ： GPU Julialang programming
posted at 03:44:23
Finally learning #JuliaLang; it's both great and disconcerting to be away from #Python #coding
posted at 04:04:28
I have spent most of my day today trying to get code working from people’s papers where they don’t have requirements files.
The value proposition of reproducible science has never been more clear. Don’t write papers if the code is not reproducible. #AcademicChatter
タグ： AcademicChatter
posted at 06:15:31
I've never heard of #JuliaLang before. Looks pretty or is it just the font? twitter.com/Michielstock/s...
タグ： JuliaLang
posted at 06:36:15
Nathanael Wong (黄智信) @natgeo_wong
Hi everyone, I am Nathanael Wong Zhixin (黄智信, Huang Zhi Xin), Singaporean, Harvard PhD Candidate, racewalker, K-pop dancer and #JuliaLang enthusiast. I study the impact of the Maritime Continent on climate at all spatial and temporal scales. #AAPIGeoRollCall pic.twitter.com/1n6SpOQUAK
posted at 08:34:13
松田未来 「夜光雲のサリッサ」第10巻7 @macchiMC72
自転車のタイヤは幅が狭くても広くてもスピードは変わらない - GIGAZINE
gigazine.net/news/20210616-...
速く走れていると感じるのは細いタイヤが発する音による「気のせいだった」という身もふたもない結果。エネルギーロスやコーナリングの安定性を考えるとむしろ幅広の方が有利。
タグ：
posted at 08:44:49
非公開
タグ：
posted at xx:xx:xx
Ryosuke Kurihara @phys_mmc_2007
NLsolveパッケージを使えば自己無撞着方程式解けるのでは？？？ #Julia言語
タグ： Julia言語
posted at 09:22:53
#Julia言語
--projectオプション付きで起動しなくてもよいスクリプトを追加しました。
github.com/genkuroki/publ...
using Pkg; Pkg.activate(".")
using HeatEqCalc
sol = HeatEqCalc.calc_sol()
HeatEqCalc. save_sol(sol)
実行時間は4.67秒と短いまま。
twitter.com/genkuroki/stat... pic.twitter.com/v9oPw19p0c
タグ： Julia言語
posted at 11:07:47
例えば「Juliaでもx.f(y, z)をf(x,y,z)のaliasとして使えるようにしろ」とか言うから荒れる。
しかし、「仮にどうしてもx.f(y,z)のスタイルを実現しなければいけないとしたら、Juliaではどのようにするのが最良の方法なのか？」とか「PyCall.jlではどうしているのか？」なら楽しい。
タグ：
posted at 11:21:31
julia foo.jl の繰り返しが不合理だと知っていても、「julia foo.jlの実行時間を1.5秒縮めたいがどういう方法があるか？」をfoo.jlの具体例を幾つか与えてみんなで話し合えば楽しいだろう。
例↓
twitter.com/genkuroki/stat...
タグ：
posted at 11:21:31
Julia言語、グッと睨むと余計なメモリ割当がなくなって、計算速度を上げられるようになってきた。
今のところやっているのは、頻繁に更新する配列を事前に割り当てることと、要所でStaticArrayを使うことくらい。型は基本的に全部コンパイラ任せ。
タグ：
posted at 11:33:48
Pythonやってたときは、計算速度を上げようとすると、Cythonで型付けしまくる場外乱闘が始まる感じだった。一つの言語で計算の中核からインターフェースまで書けるようになってとても楽になった。タノシイヨ、Julia。
タグ：
posted at 12:08:02
@kado_judo0312 @1_1_2_3_5__ #Julia言語 Taylor展開およびFourier級数展開のコードの例
ちょっとした数値積分にquadgkは非常に便利ですよね。
TaylorSeries.jlに類似のパッケージは他にも沢山ある。
高階自動微分のコードを自分で書くと色々勉強になると思います。
github.com/genkuroki/publ... pic.twitter.com/sNSPe8H1cr
タグ： Julia言語
posted at 13:05:03
Wanted to use an array type where the columns and rows are indexed by values. Quickly found NamedArrays.jl
Easy to use and not have to worry about performance. If needed, I could've written the part of NamedArrays that I need.
The power of #julialang
タグ： julialang
posted at 13:06:38
社のインタビューを受けました。
Rが好きすぎてRでJVMを実装したnoteエンジニアに話を聞いてみた｜note株式会社 @note_corp #note #社員紹介 note.jp/n/n9e62ffd6606b
posted at 13:10:15
非公開
タグ：
posted at xx:xx:xx
#Julia言語
Pythonへの精神のoverfittingが原因の「ファイルと名前空間の無意識のうちの混同」が解消していた！よかった！
プログラミングは人類的に「精神の過剰適合」が起こり易い分野だと思う。
自分が経験して来たローカルルールを普通のことだと誤解し易い。
discourse.julialang.org/t/proper-modul...
タグ： Julia言語
posted at 13:27:08
#Julia言語 Julia discourseの閲覧によって、「精神の過剰適合」は防ぎ得ないと思うようになった。余りにも多い。
discourse.julialang.org/t/implementing...
【オブジェクトにアクセスできないと小さなプロジェクト以外のことは非常に難しくなる】
これもobject.method(x)スタイルのオブジェクト指向への過剰適合。
タグ： Julia言語
posted at 13:34:02
非公開
タグ：
posted at xx:xx:xx
これ、気色悪過ぎ。
【校則はできるだけシンプルがいい。「下着は白」4文字で表せてとても良い。】
有害な #ブラック校則 を無くせば0文字にできる。
単に0文字にするだけではなく、ブラック校則を長年放置し続けて来た人達による社会全体へのお詫び発言があればさらによい。 twitter.com/seito_sidoubu/... pic.twitter.com/Jasps8VxGe
タグ： ブラック校則
posted at 14:27:24
新谷貴司(HP制作リニューアル・SEO・ @localnavi
この理屈でいけば、下着の色についての校則をなくせば
0文字
と最高にシンプルになる。
そもそも保護者としては、教師にはこんなくだらない校則を守らせることに労力と時間を費やしてほしくない。 twitter.com/seito_sidoubu/...
タグ：
posted at 14:46:20
@KevinRxse You can copy paste it into you editor and see whether you still like it 🙂
タグ：
posted at 15:14:29
ちょろ母@‘22〜海外復活、’23はキュ @Chorouzu1017
先日、搭乗してから何気なく外を見ていたら…！
もっと最初から撮っておけば良かった。初めはハート？四つ葉？と何が始まるか分からなかったんですよ…
ものの数分で消えてしまう、でも心のこもった素晴らしいアート！
#奄美空港 pic.twitter.com/NiTvrxYw2E
タグ： 奄美空港
posted at 15:21:05
お好み焼きにはカープソースをどうぞ。 @IshidaTsuyoshi
下着の色なんか指定しないのがいちばんシンプルに決まってるだろ。
0文字だよ。 twitter.com/seito_sidoubu/...
タグ：
posted at 15:58:05
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
Rip currents are narrow currents in the surf zone that move quickly away from shore. That makes them dangerous and potentially deadly. Everything you need to know about them [source, more: buff.ly/2Ofr3Wg] pic.twitter.com/UFwA3OqYdK
タグ：
posted at 17:30:00
非公開
タグ：
posted at xx:xx:xx
#Julia言語 さらなる改良版
main.jl という名の1つのファイルの中にすべての情報を詰め込んだ。
julia main.jl PRECOMPILE を実行すると、パッケージを作成して、パッケージのprecompileが行われる。
その後は julia main.jl で小さな遅延で実行可能になる。
github.com/genkuroki/publ... pic.twitter.com/rC6CAveoUh
タグ： Julia言語
posted at 17:46:25
#Julia言語 具体的に何をやっているかは
github.com/genkuroki/publ...
を見ればわかる。Makefile, main.jl, animate.jl, screenshot.png 以外のファイルは make で自動生成されます。
「main.jl 実行時の遅延を1.5秒程度減らす」のようなことに興味がある人は理解しておくべきことだと思います。
タグ： Julia言語
posted at 17:49:30
未成年者の権利が制限されるのは「子供の利益になる」場合に限られます。
決して「大人が管理を楽にする為」ではありません。 twitter.com/seito_sidoubu/...
タグ：
posted at 17:51:53
#Julia言語 まとめ続き
さらに、main.jlというスクリプト単体で、1回目に
julia main.jl PRECOMPILE
を実行すると、パッケージを自動生成してプレコンパイルを行い、2回目以降の実行
julia main.jl
の遅延が少なくなるものを作れる。
具体例
↓
github.com/genkuroki/publ...
タグ： Julia言語
posted at 18:08:19
校則はなるべくシンプルがいいそうです。
それなら全廃をお勧めします。
生徒が授業や学習生活で困らないことは規制する権限を持ちません、あなた。 twitter.com/seito_sidoubu/...
タグ：
posted at 18:24:30
@tech_optimist @MansiiGoyal Time to stop writing papers in Python and switch to #JuliaLang :)
タグ： JuliaLang
posted at 20:13:51
改めて考えても、どう考えても頭おかしいよな。
下着の色を校則で定めて、それが守られてるかどうかどう確認するのかねえ。仮に同性の先生に限定されてチェックされるにしても、十分人権侵害なんじゃないの。 twitter.com/seito_sidoubu/...
タグ：
posted at 20:14:52
「下着の色なんてどうでもいいから」の続きが「だから下着の色は白というルールを作ったり維持したりすることは人権を害するやってはいけないことである」なら普通だと思うが、「とりあえず『下着は白』としてある」になっている。
生活指導部長として人権についてどう教えているのやら？ twitter.com/seito_sidoubu/... pic.twitter.com/eq4CMOMFuN
タグ：
posted at 20:37:15
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
非公開
タグ：
posted at xx:xx:xx
I just discovered that in @JuliaLanguage pinv() of a symmetric matrix is not symmetric. Now there is myPinv().
タグ：
posted at 21:15:56
@syutoken_sanka 私の父は白血病で余命１年でしたが臍帯血移植をし余命宣告から４年経ちましたが元気です。４年の間に私は結婚、出産をし父はおじいちゃんになりました。神戸の赤ちゃんの臍帯血と聞きました。感謝の気持ちでいっぱいです。父は移植で血液型が変わり何故か喜んでました(笑)
タグ：
posted at 22:12:49
臍帯血は捨てるものを保存して使うので、骨髄バンクのようにドナーに負担がかかったりすることがないんですよね。
ホントにもったいない twitter.com/z_zo38aika/sta...
タグ：
posted at 22:23:52
ほりたみゅ (mastodonも見てね) @Hyrodium
tiffフォーマットが64bit浮動小数点数でRGBを保存できるの知らなかった pic.twitter.com/Krw88JeEgV
タグ：
posted at 22:54:08
@kado_judo0312 #Julia言語 factorial(n)だとbigが必要になるので、
using SpecialFunctions
fact(n) = gamma(n + 1)
として、fact(n)を使うとよいです。係数cを与えて多項式を評価する函数 evalpoly(x, c) があります。cをタプルで用意しておくと、これ爆速です。 pic.twitter.com/T7jjyCFHFP
タグ： Julia言語
posted at 22:54:25
Dr. Rachel Kurchin @rachel_kurchin
Guess who has two thumbs (only one shown) and is psyched for @JuliaConOrg?!? Now to actually make my talk… pic.twitter.com/5za6Raj76a
タグ：
posted at 23:48:50
