黒木玄 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 を見て下さい。
2018年01月12日(金)
#Julia言語 よくやってしまう失敗は、=に.を付け忘れること。例えば
s = a.+b.+c.+d
と書いてしまう。すでにsを配列として確保してあれば
s .= a.+b.+c.+d
の方が速い。これは
@ . s=a+b+c+d
と書いても同じ。@の後の空白は除く。 = ではなく、.= と書かないと無駄にメモリを消費する。 pic.twitter.com/PAm6pJKfAX
タグ: Julia言語
posted at 00:38:15
#Julia言語 函数型オブジェクトの作り方は
mutable struct WeightedMean{T} w::Array{T} end
(f::WeightedMean)(x) = mean(i->f.w[i]*x[i], eachindex(f.w))
これを
mutable struct Hoge w end
(f::Hoge)(x) = mean(i->f.w[i]*x[i], eachindex(f.w))
とすると大域変数を含む函数並に遅くなる。 pic.twitter.com/2Lc4RpFv4F
タグ: Julia言語
posted at 01:49:22
#Julia言語 immutableにしても事情は同じ。函数的オブジェクトを作るときには型付きにすることが大事。
struct Poge{T} w::T end
(f::Poge)(x) = mean(i->f.w[i]*x[i], eachindex(f.w))
は速いが
struct Moge w end
(f::Moge)(x) = mean(i->f.w[i]*x[i], eachindex(f.w))
は遅い。 pic.twitter.com/HAnOL9Smyi
タグ: Julia言語
posted at 08:39:06
#Julia言語
aが大域変数のとき、函数f(x,a)をxについて数値積分するために、
quadgk(x->f(x,a), -1000, 1000)
または
g(x) = f(x,a)
quadgk(g, -1000, 1000)
と書くと遅くなる。
mutable struct Piyo{T} a::T end
(g::Piyo)(x) = f(x, g.a)
quadgk(Piyo(a), -1000, 1000)
なら速い。 pic.twitter.com/sO5r6TDABw
タグ: Julia言語
posted at 09:15:06
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
昨日のこと。まだ正負の計算がおぼつかない中1生に、①「(-2)×3」と②「2×(-3)」を答えさせた。そうしたら①は「-6」と答えられるのに、②は酷く考え込んでいる。そこでこの子が①と②のどこに違いを感じて困っているのかを探ってみると、案の定の結末・・・。
タグ:
posted at 15:09:57
#Julia言語 Julia言語での気軽な最尤法などに使うOptim.jlのoptimize函数で多次元の場合にデフォルトで使用される最適化アルゴリズムのコードは bicycle1885 さんが書いたものです。
github.com/JuliaNLSolvers...
みんな「すごい」と言っているJulia言語のチュートリアルの作者と同じ人。
タグ: Julia言語
posted at 15:11:59
#掛算 順序固定強制の被害者でも、この子みたいに小・中学生時代に呪縛を解いてもらえる子は幸せ。
呪縛が解けないまま、「でも、文章題は別です!(キリッ」と言う困った大人になってしまった人をよく見かける。 twitter.com/zatukun/status...
タグ: 掛算
posted at 15:25:55
偉くないので間違ってるかもだけど、これはループ最深部で view を呼び出した時に SubArray オブジェクトを大量に作っているのが原因みたいです。良いニュースとしては開発版ではかなり改善していて、メモリ使用量が1/5になっています。0では無いけれど。確認用コード: gist.github.com/tkf/f7bc4a78ee... twitter.com/genkuroki/stat...
タグ:
posted at 15:28:49
Juliaのドキュメント「より厳密に言えば、全ての有限の長さを持つ文字列の集合 S と文字列連結演算子 * は 自由モノイド (S, *) を形成します。」で想定する使用者が伺える
タグ:
posted at 15:30:08
こんなの自分で書けるはずないな(´・ω・`)
openlibm/e_j0.c at master · JuliaLang/openlibm github.com/JuliaLang/open...
タグ:
posted at 15:46:51
おーすごいよく調べられてる
> 早く知っておきたかったmatplotlibの基礎知識、あるいは見た目の調整が捗るArtistの話 - Qiita qiita.com/skotaro/items/...
タグ:
posted at 16:23:29
多重ディスッチの思想を知った:「…多重ディスパッチは特に数学のコードで有用です。演算がどの引数に、より「属している」かと不自然なことを考えても、意味がないからです。x + y の式の中の足し算は、 x に yよりも属してると思いますか?…」
hshindo.github.io/julia-doc-ja-v...
タグ:
posted at 16:50:00
#JuliaLang は1.0がリリースされたらボチボチ使い始めようと思ってたら到達度93%で,もう目前だった。楽しみ github.com/JuliaLang/juli...
タグ: JuliaLang
posted at 17:17:10
非公開
タグ:
posted at xx:xx:xx
藤ノ木優 アンドクター1/24発売 @kyusan_obgy
HPV池田班の報告について、調査委員会から捏造行為は無かったという結論から反ワクチン団体の人達が池田班の実験結果は正しかったなんて大誤解をしているみたいだけど、不倫会見でいう限りなく黒に近いグレーゾーンなんで誤解しないようにねwww.mhlw.go.jp/bunya/kenkou/k...池田+HPV+調査委員%27
タグ:
posted at 17:24:14
藤ノ木優 アンドクター1/24発売 @kyusan_obgy
厚生労働省の調査結果はさっき添付したリンクで全て読むことができるけど・・そもそもHPVワクチン摂取マウスの保存血清を使用した再現実験では脳への反応は認められなかった。そもそも検体数1で実験することなんて普通あり得ないんだけどね。しかもいい結果が出たら尚更再現したいはずだけど。#HPV
タグ: HPV
posted at 17:27:42
藤ノ木優 アンドクター1/24発売 @kyusan_obgy
以下厚生労働省コメント。池田教授について。NEWS23の発言については実験結果が科学的に証明されたかのような印象を与えてしまった。誤った情報が社会に広がってしまった以上再実験による科学的証明が不可欠である。科学的検証に耐えうる検体数を用意し最初から実験の再現が必要。 #HPVワクチン
タグ: HPVワクチン
posted at 17:31:22
藤ノ木優 アンドクター1/24発売 @kyusan_obgy
この内容で名誉毀損といって訴えを起こす神経の図太さには非常に感銘を受けますが・・。まあこういう経緯で宮迫のグレーゾーンや藤吉久美子の個室でマッサージと同じレベルの内容なので、反ワクチンの方々は池田班に捏造は無かった=池田班の発表は全て正しかったんだ!なんて考えはしないようにね
タグ:
posted at 17:38:21
福富さんというと、日大数学科事件を思い出す。1962年(昭和37年)、福富さんと、その引きで日大に就職していた木下素夫、銀林浩、倉田令二朗の4人がいきなりクビを宣告された。さらに福富さんは職員から暴行を受けた事件である。福富一派の教育方法が気に入らないということだったらしい。
タグ:
posted at 17:41:10
藤ノ木優 アンドクター1/24発売 @kyusan_obgy
反ワクチンにとって池田班の研究が唯一の拠り所で希望なのはわかるけど、あまりに誤解している人が多い印象だったのでツイートしておきました。ちなみに本来恨むなら村中さんやwedgeではなくて、適当な実験やデータの公表をした池田班を恨むべきなんだよ。科学的な議論ってそういうものだから
タグ:
posted at 17:41:54
非公開
タグ:
posted at xx:xx:xx
非公開
タグ:
posted at xx:xx:xx
コンフェデレーション橋はカナダ東部にある全長12.9kmの橋で、冬季に凍結する海に架かる橋として世界最長。その橋げたにより氷が規則的な形に切断されている様子の空撮。 bit.ly/2ExSk0M reddit の billerz15 氏が「フライドポテトのよう」と形容して投稿。日本人的には寒天や水羊羹が近い。 pic.twitter.com/Omz9ueuloW
タグ:
posted at 18:27:42
非公開
タグ:
posted at xx:xx:xx
ごまふあざらし(GomahuAzaras @MathSorcerer
x=cos(x) を反復法で解くコードを書いてみたよ.
#python #fortran #Julia言語 pic.twitter.com/6lMGvLgGVW
posted at 21:10:24
#JuliaLang 何故かは知らぬが今日はPackage.add(“Luxor”)が通った。これで使えるのかな。ついでにGadflyとかも入れておこう。
タグ: JuliaLang
posted at 21:58:26
ごまふあざらし(GomahuAzaras @MathSorcerer
2分法による方程式 x=cos(x) を解くプログラム.最初にPythonで書いてコピペしたコードをJuliaにそのまま記述した感じ.Juliaらしい書き方がいまいちわからぬ. pic.twitter.com/s9rspEN3mX
タグ:
posted at 22:37:33
ごまふあざらし(GomahuAzaras @MathSorcerer
@ceptree タイプ量が増えてアレな感じがします. Pythonのnotに相当するシンタックスってあります?
タグ:
posted at 22:53:05
みんなにニューロンのおもしろさが伝わったのか、ツイートが伸びてて、フォロワーも増えたので、記念にホジキン-ハクスレーモデルのアニメーションつくった pic.twitter.com/9wiXnLUMFt
タグ:
posted at 23:10:28
Jupyter はIulia, Python, Rからです。 "Second, the core programming languages supported by Jupyter are Julia, Python and R." github.com/jupyter/design...
タグ:
posted at 23:36:13
#Julia言語 Julia v0.6.2 (現行版)について次を書いて公開しました。非常に長いです。Julia言語を遅くする方法がたくさん紹介されています(解決法有り)。
nbviewer.jupyter.org/gist/genkuroki...
Julia言語で計算が遅くなった場合の解決法
v1.0が出たときにこれが時代遅れになっていれば滅茶苦茶うれしいかも。
タグ: Julia言語
posted at 23:37:58
@genkuroki 遅くなりましたが、ありがとうございます。元となったコードはあるのですが、完全に書き捨てモードで書いてたので、ちょっとそのままの形では公開できないな、、と。そのうち公開したいと思います。
タグ:
posted at 23:42:55
非公開
タグ:
posted at xx:xx:xx