Twitter APIの仕様変更のため、「いいね」の新規取得を終了いたしました

黒木玄 Gen Kuroki

@genkuroki

  • いいね数 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 を見て下さい。
Favolog ホーム » @genkuroki » 2019年05月12日
並び順 : 新→古 | 古→新

2019年05月12日(日)

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 その解釈のもとで x はバネ上の座標です。u_{tt} = u_{xx} - g(u) 型の偏微分方程式を「落下するバネ」の運動方程式だと解釈するとき、時刻 t におけるバネ上の点 x の高さが u=u(t,x) になっているという解釈になる.

その解釈のもとで x は空間座標ではないことに注意。

タグ: 数楽

posted at 23:50:07

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 Sine-Gordon方程式

u_{tt} = u_{xx} - sin(π u)

は高さ u について開区間 (偶数, 偶数+1) 上では重力は下向き、開区間 (奇数, 奇数+1) 上では重力は上向きのちょっと変わった状況で(自然長0の)バネを「落下」させた場合の運動方程式になっています。空間方向は上下の1次元分 u しかない。

タグ: 数楽

posted at 23:50:07

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 落下するバネで遊び始めたときには何も考えてなかった。

まさか、Sine-Gordon方程式の数値解を求めることまでやってしまうことになるとは予想していませんでした。

Sine-Gordonについて何も理解していなかったこともよくわかった。

やはりいじり倒したもの以外は全然理解していない。

タグ: 数楽

posted at 23:44:15

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語 これも

u_{tt} = u_{xx} - sin(π u)

所謂 SIne-Gordon 方程式の場合。初期状態が 2 sin(π x) です。

段々 sine カーブが歪んで行く。

x軸をもっと長くして別の解を求めることもできるのですが、まだやっていません。

ソースファイル↓
nbviewer.jupyter.org/gist/genkuroki... pic.twitter.com/tFbW4zFv6W

タグ: Julia言語 数楽

posted at 23:41:57

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語

u_{tt} = u_{xx} - sin(π u)

所謂 SIne-Gordon 方程式の場合です。

数値計算が破綻するまでプロット。 pic.twitter.com/AwxHUqkSqp

タグ: Julia言語 数楽

posted at 23:39:12

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語

u_{tt} = u_{xx} - 1/u²

これも数値計算が破綻する直前まで pic.twitter.com/ZtfEXyXwn7

タグ: Julia言語 数楽

posted at 23:39:11

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語

u_{tt} = u_{xx} - 1/u

数値計算が破綻する直前まで。 pic.twitter.com/MdZ7owHcU8

タグ: Julia言語 数楽

posted at 23:39:11

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 訂正版: u_{tt} = u_{xx} - 1 の特殊解として

u(t,x) = (f(x-t)+f(t+x))/2 - t²/2, f(x) = x²/2, -1≦x≦1

が取れる。ただし、f(x)は周期2でx軸全体に拡張しておく。

タグ: 数楽

posted at 23:36:38

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語 2つ前のツイートの動画は u_{tt} = u_{xx} - 1 の解です。これは一様重力(重力加速度1)の場合で、初期条件の x=0 での高さ(初期状態の下端)はしばらく全く動かない。

以下は非一様重力下の場合。次の動画は

u_{tt} = u_{xx} - u

の解です。数値計算が破綻するまで計算している。 pic.twitter.com/axryC0iONy

タグ: Julia言語 数楽

posted at 23:35:24

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語 1つ前のケースのexact解は

u(t,x) = (f(x-t)+f(t+x)) - t²/2, f(x) = x²/2, -1≦x≦1

です。1つ前のツイートの動画は数値解の動画なのですが、確かにこれになっていそうです。

タグ: Julia言語 数楽

posted at 23:31:04

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語

nbviewer.jupyter.org/gist/genkuroki...
非一様重力下で落下するバネの偏微分方程式のFFTを使った数値解法

u_{tt} = u_{xx} - g(u) 型の偏微分方程式(一般には非線形)を数値的に解いてみました。普通に離散化すると質点をバネで繋いだ場合と同じになりそうなので、FFTを使った別の方法を採用。 pic.twitter.com/kQAiVKXTLF

タグ: Julia言語 数楽

posted at 23:28:24

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

19年5月12日

噂のパルミジャーノなる人物のツイートをぱらぱらと眺めてみたのだけど、ようわからん。フォロワーを増やしたくて炎上を狙ってる高校生くらいにしか思えん。炎上より、プロフィールを書くほうがいいんじゃないか?

タグ:

posted at 23:25:47

らじうむ小山_PPPMP @Ra_koyama

19年5月12日

ネタかと思ったら、ほんとにこういう結果なのか…orz twitter.com/rouninsyurachi...

タグ:

posted at 22:50:48

黒木玄 Gen Kuroki @genkuroki

19年5月12日

ほぼ半々に分かれている点が何とも。

アンケート問題を読まずに無作為に答えたのと区別がつかない(笑)。

twitter.com/musicisthebest...

タグ:

posted at 22:48:14

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@kankichi57301 @sekibunnteisuu もう面倒なので私は期待の新人さんをブロックしました。

件の新人さんとやりとりする場合には返信先から私を外して下さい。

タグ:

posted at 22:45:16

k @musicisthebest_

19年5月12日

浪人界隈やばいな…。 twitter.com/rouninsyurachi...

タグ:

posted at 22:07:41

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

19年5月12日

前川の悪事は天下り斡旋に関わったことだけではなかったということです。この人物を「識者」扱いするのは愚かな行為ですよ twitter.com/genkuroki/stat...

タグ:

posted at 21:56:28

春 @sunbluesome

19年5月12日

メモリ確保を避けることと、インライン化が高速化の近道っぽい。
#Julia言語

タグ: Julia言語

posted at 21:51:03

春 @sunbluesome

19年5月12日

やっと2層NNの実装に成功。
長かった。
先人の知恵を借りた(コピペ)。
行列演算のところ、特にチューニングしていないのでnumpyの方が速くなってしまっているっぽい。
#Julia言語

タグ: Julia言語

posted at 21:50:01

Kuri TsukaΦ (くり つかお) @20_0855369231

19年5月12日

本日の演習授業にて

誰か一人くらいは「こんなん10/49に決まってますよね」と即答してほしかったのですが残念ながらそれはおらず…(条件付き確率の定義に従ってゴリゴリ計算ばかり笑)

いくつか布石は打ったつもりだったんですが難しい pic.twitter.com/vJFqeIrgxG

タグ:

posted at 21:43:38

Yossy @Yossy_K

19年5月12日

面の皮をぶ厚くして「面従腹背」でのし上がっていかないと、ああいう人は出来上がらないってことなのかも。 twitter.com/tamai1961/stat...

タグ:

posted at 21:29:41

RochejacMonmo @RochejacMonmo

19年5月12日

blogos.com/article/376481/
なんで役に立つかどうかとか数学の話をしているふうをよそおって小学校の算数みたいな例を出すんだろう?
なんで、高校数学で教えられている(or教えられていた)整数とか三角関数とか微積分とか行列とかベクトルとか、スゲー社会に役立ってるって話をしないんだろう? pic.twitter.com/4HiDWLfKL6

タグ:

posted at 21:13:09

しかのつかさ @sikano_tu

19年5月12日

権力側にいるときは天下りの斡旋もやってて、追放されると今度は左翼的に受けのいいことを言ってるわけで、要するに保身がうまい口先だけ上手な人で全く信頼できないと思うな。 twitter.com/tamai1961/stat...

タグ:

posted at 20:56:38

すむーずぷりんちゃん @mat_der_D

19年5月12日

e の定義は exp(1) 派

タグ:

posted at 20:50:27

ごまふあざらし(GomahuAzaras @MathSorcerer

19年5月12日

こういう問題N回見ていてもう疲れたので berrycondaでいいじゃないかって思ってまキュ。 twitter.com/termoshtt/stat...

タグ:

posted at 20:12:38

S (ツイートはスレッド全体をご確認く @esumii

19年5月12日

T先生が釈迦やニュートンやアインシュタインではない、は文字通りの意味では良いとして(相対的な差の比喩としては近いと思いますが)、だから自分のほうが物理を理解している、というロジックはすごい

タグ:

posted at 19:45:37

高 史明(TAKA, Fumiaki) @Fumiaki_Taka

19年5月12日

「面従できなくて1話のBパートあたりで殺されそうな人」という定評のある僕はああいう散々長い物に巻かれておいて後になって「あれは間違っていると思ってました」とか言い出す種類の人間が嫌いなんだ…。

タグ:

posted at 19:40:54

原田 実 @gishigaku

19年5月12日

前川喜平さんが「江戸しぐさ」導入の担当者だったのが意外という人がおられるようだけど『私たちの道徳』作成当時の初等中等局長が前川さんだったことはすでにわかっていたのだからむしろ当然だったわけで、この度のインタビューではそれが本人の証言で裏付けられただけ

タグ:

posted at 19:23:28

高 史明(TAKA, Fumiaki) @Fumiaki_Taka

19年5月12日

元文科の前川氏って頭はいいのかもしれんけど本当にしょーもないな…。自分が責任ある地位にいたときにバカな政治家にほいほい従っといて後になって「面従腹背だった」もクソもない。心の中で反感持つだけなら誰でもできるし「面従」したことこそが自分の遺産だろ…。

タグ:

posted at 19:20:42

非公開

タグ:

posted at xx:xx:xx

@kankichi57301 @kankichi57301

19年5月12日

@sekibunnteisuu @genkuroki 符号反転等号飛躍の儀のほうがステップが増えないってなんやそれ? #超算数

タグ: 超算数

posted at 19:13:33

黒木玄 Gen Kuroki @genkuroki

19年5月12日

後で情報をたどり易くするためのリンク

私はミュートしているのでリンク先ヘのリプライは読めない。

twitter.com/yossy_k/status...

タグ:

posted at 18:31:06

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu 「バカは黙ってろ」で済ませて問題ない案件ですね。
ミュートして大正解。

タグ:

posted at 18:18:14

みんなのメザャピョ @Mezame0318

19年5月12日

で、理想ばね想定理論を現実ばねに当てはめていいのか、に納得行く説明もらえなかったからってそこが面白い部分なのに理解してない時代遅れ大したことない無名研究者だの言いたい放題
性格に難有りだけどこの偏執的な完璧主義は学者向きな気もしますね

タグ:

posted at 18:15:36

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 別バージョンもありました pic.twitter.com/APIVDtrDFG

タグ:

posted at 18:13:04

みんなのメザャピョ @Mezame0318

19年5月12日

そのレベルの話してないやん、ぶっちゃけ
もちろん理想ばねとは違うが性質をおおまかに捉えた現象をスリンキーで見られるよ!って本文に書いてあんのに
その「おおまか」部分が異常に許せないらしい

タグ:

posted at 18:07:35

みんなのメザャピョ @Mezame0318

19年5月12日

物理クソ苦手マンだからバネの話頑張って読んだが
レッジャーノさんはスリンキーは理想ばね(ばね質量を無視する)とは違うから完全な証明とは限らないだろ 現実ばねと理想ばねのギャップを無視して理論適用できるのかよ証明しろって延々絡んでるんだなこれ

タグ:

posted at 18:06:13

黒木玄 Gen Kuroki @genkuroki

19年5月12日

そうそう、その件もあった。

【私は学長です】事件(笑)

twitter.com/sekibunnteisuu...

タグ:

posted at 18:02:06

黒木玄 Gen Kuroki @genkuroki

19年5月12日

子連れ狼の「ベビーカー」がチート性能を持っていることについては

lonewolfcub.jimdo.com/%E5%AD%90%E9%8...

などを参照せよ(笑)。

沢山の強力な武器が収納されており、拝一刀自身の手で完成させられたガトリングガンまで装備。

底には鉄板。防御にも使える。

船としても利用可能で激流も降れる。

タグ:

posted at 18:00:24

黒木玄 Gen Kuroki @genkuroki

19年5月12日

別の記録来た!

後で検索で見つけ易いようにキーワードを並べておく

パルジャミーノ
理解のレベルの差ですねー
小池一夫 子連れ狼 フィクション
移項
みんなに笑われている
大型新人

twitter.com/sekibunnteisuu...

タグ:

posted at 17:52:27

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 読みたくないでしょうが^^、資料として
twitter.com/sekibunnteisuu...
に続くツイートです pic.twitter.com/SVWPCG6SHg

タグ:

posted at 17:44:08

Akinori Ito @akinori_ito

19年5月12日

前川氏は在職中は偽歴史を教科書に載せ,天下りを斡旋して,辞職後はそれを他人事のように批判するわけね。さすがは我らの監督官庁の元次官だ

タグ:

posted at 17:38:40

黒木玄 Gen Kuroki @genkuroki

19年5月12日

とんでもニセ歴史の「江戸しぐさ」を文科省が作った道徳教材に入れる仕事をしたのは前川喜平氏(元・文部科学官僚)だった!

なんちゅう迷惑なことを!

webronza.asahi.com/national/artic... pic.twitter.com/RTRzcU4TTe

タグ:

posted at 17:35:04

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 の Plots.jl ではlegendの表示を半透明にするには、

default(:bglegend, plot_color(default(:bg), 0.5))
default(:fglegend, plot_color(ifelse(isdark(plot_color(default(:bg))), :white, :black), 0.6))

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

タグ: Julia言語

posted at 17:24:30

Kiyohiro Maeda @kmaeda

19年5月12日

面従腹背ではなく唯々諾々だったのか。 twitter.com/gishigaku/stat...

タグ:

posted at 17:21:03

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 の Plots.jl パッケージを使い倒したい人は特に

docs.juliaplots.org/latest/attribu...
Attributes

の熟読が必須。私はブラウザでブックマークしてある。

タグ: Julia言語

posted at 17:20:55

前田敦司 @maeda

19年5月12日

「「江戸しぐさ」をまだやってるんですか! あれはやめた方がいいですよてんあれはねえ、大失敗。僕が初等中等局長のとき、下村さん(下村博文・元文科相)に言われて作った。あんなインチキなもの…」他人事感が… / “親を悩ます「PTA問題」 前川喜平さんに聞いた - 前川…” htn.to/2KxrYP6P4g

タグ:

posted at 17:16:14

玉井克哉(Katsuya TAMAI) @tamai1961

19年5月12日

よくここまで無責任なことが言えるな。初等中等教育の責任者だったわけだろ。大臣が思い込みで変なこと言ったら、止めるのが仕事。それしないで辞めてから批判するとは。

大臣を止め、止まらなければ職を賭す。辞職して大臣を諌める。それをしない者が、何で局長の高給を食み、次官にまでなるのか。

タグ:

posted at 17:04:34

原田 実 @gishigaku

19年5月12日

とりあえずこの栄誉は誰に帰すべきかははっきりしたな。
togetter.com/li/853131

タグ:

posted at 16:57:34

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 ううむ、pyplot()版は線が細過ぎるのかも。切りがないのでやめます。その辺をいじって自分のために完璧な動画を作りたい人は

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

の lw=0.2 の部分を lw=0.3 とか lw=0.4 に変えて実行し直してみて下さい。変更後のファイルを私に無断で公開配布して構いません。

タグ: Julia言語

posted at 16:39:50

@kankichi57301 @kankichi57301

19年5月12日

セントアイビスに行った総数は
1+x+xy+xyz+xyzw(足し算も掛算も順不同)であり単純な話。けったいな理屈の入る余地はない。 #掛算
(なんかtypoしてるけど見逃してw)

タグ: 掛算

posted at 16:36:27

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 pyplot()バックエンドで色の設定を自前で行ったもの。plot()の中で

color_palette=(色の巡回リストの配列)

とするだけ。

2つの傾いた円盤のあいだの「相互作用」の様子が見ていて楽しい。

ソースファイル↓
nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション pic.twitter.com/VGPOKfLO5z

タグ: Julia言語

posted at 16:35:26

@kankichi57301 @kankichi57301

19年5月12日

同じくマザーグースの歌の”As I was going to St. Ives”でもおっさん一人当たりの嫁さんの数x
嫁さん一人当たりの袋の数y袋日取る当たりの猫の数z猫一匹あたりの子猫の数wとすると(続) #掛算

タグ: 掛算

posted at 16:31:36

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 gr()バックエンドで色の設定を自前で行ったもの。かなり派手になっている。確かに見易くなったような気がする。 pic.twitter.com/NmYx1rku3m

タグ: Julia言語

posted at 16:29:08

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 まず pyplot() バックエンドから gr() バックエンドに変えただけ。色の設定はデフォルトのまま。gr()バックエンドだと背景の表示が弱すぎ。しかし、動画作成時間は1/3に減る。 pic.twitter.com/g9Wp0IqKKK

タグ: Julia言語

posted at 16:29:07

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 まさにそれ。そういう細かいことをやるのが難しい。やり方を調べ上げること自体が結構面倒。だから、「デフォルトの設定を上書きしてのカスタムプロット!」がもっと公開されてくれないと不便なのだ。

やってみました。続く

twitter.com/tsatie/status/...

タグ: Julia言語

posted at 16:29:07

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

私も山川賢一氏の仲正氏批判には賛同したのだけど、山川賢一の著書は一冊も読んだことはなく、
それに対して仲正氏の本はかなり購入しているので、どちらかと言えば、仲正ブラザーズではある。
仲正ブラザーズというよりは、仲正ウォッチャーあるいは仲正ストーカーの方が正確

タグ:

posted at 16:22:26

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

注:「山川ブラザーズ」と言っても知らない人が多いと思うので、簡単に説明すると、サブカルチャー研究者の山川賢一氏が仲正氏を批判して賛同者がかなり出たので、仲正氏がまとめて「山川ブラザーズ」と命名した。

タグ:

posted at 16:18:55

非公開

タグ:

posted at xx:xx:xx

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

私は仲正氏ウォッチャーを何年もしているけど、専門よりも副業的な方が無茶苦茶忙しいように思える。
現代思想の啓蒙活動とか俳優業とか。
それはそれで人類への偉大な貢献だけど、やはり以前のように無名&匿名からの仲正批判への応答に労力と時間を割いて欲しい。

タグ:

posted at 16:14:56

らじうむ小山_PPPMP @Ra_koyama

19年5月12日

注)を付けると嫌がられる、売れないという版元が多いというツイがTLに流れているが、そんなことはないぞよ。少なくとも私は、注があって、章末・巻末に参考文献や参照URLがちゃんとあると、なんか得した気分でうれしいぞ。

タグ:

posted at 16:12:44

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

仲正氏は無名&匿名の批判にもいちいち反応する大人気なさ(褒め言葉)が良かったのに、多忙のせいかそういうことをしなくなったのは悲しい。

タグ:

posted at 16:09:30

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

仲正氏の隔月ネット連載が更新された(以前は毎月更新してたけど、多忙のため隔月に変更)。

「山川ブラザーズ」という言葉は見当たらず、「山川」だけになってしまっているのは悲しい。
雑魚には触れないという姿勢に変更してしまったんだろうか? pic.twitter.com/CTP90oMCme

タグ:

posted at 15:57:36

前田紀至子 @ki45m

19年5月12日

私の推しを!見てください! pic.twitter.com/hmKpSeY8fE

タグ:

posted at 15:51:43

前田紀至子 @ki45m

19年5月12日

ロシアの路上ライブでワガノワの生徒が飛び入りダンスするの巻。
ミュージシャン、こんなシットインされたらテンション上がりまくりだろうな〜(そしてバレエじゃなくても膝から下の動きが完璧ワガノワっ子) pic.twitter.com/mN8PHzpM6i

タグ:

posted at 15:50:03

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

@genkuroki > 仲正さんって「いいひと」ですよね。

続きがアップされました。

twitter.com/meigetsus/stat...

タグ:

posted at 15:48:08

非公開

タグ:

posted at xx:xx:xx

the_spoiler @don_jardine

19年5月12日

面従腹背とか言いながら、現役時代は上に言われた通り従って、後で自分は上の命令に従っただけと言い訳するのなら、気概も矜持も感じない。本当に面従腹背を座右の銘だというなら表向き上からの命令をハイハイと聞いておいて密かに方向を是正するとかそれくらいの芸当を見せないとただの責任逃れ。

タグ:

posted at 15:29:06

a_yam @yachts111xenon

19年5月12日

物理学会に出席してる人が
@Hal_Tasaki 先生に
あんなツイートできるんですかね? twitter.com/genkuroki/stat...

タグ:

posted at 15:28:35

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 についてググるとちょっとアレな状態になって困ると思っている人はGoogleの設定を

地域の設定:アメリカ合衆国

にするとよいです。結構おすすめ。

あと、Julia関係の検索時にキーワードとして

jl

を入れておくという技もあります。 pic.twitter.com/46gPzuIQk2

タグ: Julia言語

posted at 15:25:24

ʇɥƃıluooɯ ǝıʇɐs @tsatie

19年5月12日

元呟きの二つの動画より少し動きが早いようでより素敵感が。ただ色目にもう少し差がある方が特に別れてからの挙動が分かりやすそうに思う。 #自分でやれ #やっとこないだ環境だけ再構築しただけやし twitter.com/genkuroki/stat...

タグ: やっとこないだ環境だけ再構築しただけやし 自分でやれ

posted at 15:22:14

黒木玄 Gen Kuroki @genkuroki

19年5月12日

こっちでも記録を取っていた↓

twitter.com/ookubotact/sta...

タグ:

posted at 15:19:26

山本 卓 @sgryjp

19年5月12日

Wrote a utility (overload) function of `redirect_stdout(f, stream) which can be used for an IOBuffer:

gist.github.com/sgryjp/088f57b...

#julialang pic.twitter.com/j6LL80dtcP

タグ: julialang

posted at 15:14:24

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 初期条件を変えてもう1つ動画を作りました。N=256

nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション

こちらの方は最初から退屈せずにすむ。こういう動画を作る作業は微調節の試行錯誤の時間を結構取られる。

ぶわーっとなるところが結構快感。 pic.twitter.com/duftQo5nHe

タグ: Julia言語

posted at 15:09:59

らじうむ小山_PPPMP @Ra_koyama

19年5月12日

うーむ、「放送大学は大学じゃない→大学ですよ→もの知らずめ→学長です」の件の類似事故(笑)が、物理分野で起こっていたのかw

タグ:

posted at 15:06:56

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

辛い環境で頑張っているのに報われない日々が長く続くと、SNS世界で匿名を使って暴れたくなる人はいるからなあ

タグ:

posted at 15:06:27

OokuboTact 大久保中二病中年 @OokuboTact

19年5月12日

ツイッターの世界では、こういう無敵人がたまに現れるけど、実際に会うと普通の人だと思う。 pic.twitter.com/ydiyQM5dkm

タグ:

posted at 15:01:50

積分定数 @sekibunnteisuu

19年5月12日

地毛の「黒染め指導」は廃止を 大学生らが署名活動 www.kyobun.co.jp/news/20190508_... @kyoiku_shimbunさんから

タグ:

posted at 14:55:22

山本 卓 @sgryjp

19年5月12日

For unit testing, it is often useful to redirect stdout and/or stderr to a memory buffer so that we can examine whether the program writes expected message or not. Although I found a solution as above, there should be an easier way... Hmm... #julialang

タグ: julialang

posted at 14:53:41

山本 卓 @sgryjp

19年5月12日

Hmm, I could use another overload:

output = nothing
backup = stdout
rd, wr = redirect_stdout()
try
# Do something using stdout
finally
redirect_stdout(backup)
close(wr)
output = String(read(rd))
close(rd)
end
# Use `output` here

タグ:

posted at 14:48:01

Yossy @Yossy_K

19年5月12日

言ってることは特に異論ないけど、そう言ってる人がどんな結果を仕事で出してるのかさっぱり判らんのではどうしようもないな。 twitter.com/paruparu100/st...

タグ:

posted at 14:45:54

うみんちゅ @DEEPBLUE1219

19年5月12日

ツイッター怖い((((;゚Д゚))))))) twitter.com/genkuroki/stat...

タグ:

posted at 14:44:47

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#Julia言語 最近では Plots.jl のpyplot()バックエンドを使うことが増えている。

最近気付いたのですが、pyplot()した後には

PyPlot.plot(~)

のようにしてPyPlot.jlを普通に使えます。

matplotlibに慣れていている人にはPlots.jlのpyplot()バックエンドは結構おすすめかも。

タグ: Julia言語

posted at 14:44:24

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 初期条件は1e-8のオーダーで揺らしています。N=256です。

プロットの仕方の習得は結構重い。Exmaplesの内容がシンプル過ぎると、自分がやりたいと思った細かな作画設定の仕方が分からないので、プロットの例は色々面倒なことをやっているものが結構貴重だと思う。

タグ: 数楽

posted at 14:39:50

黒木玄 Gen Kuroki @genkuroki

19年5月12日

#数楽 #Julia言語

twitter.com/woody_egg/stat...
と同じようなことをやってみました。

nbviewer.jupyter.org/gist/genkuroki...
Lorenz方程式の数値解のアニメーション

DifferentialEquations.jl と Plots.jl を使った簡単な演習問題なのですが、数値解を求める部分よりもプロットや動画作成の方が難しい。 pic.twitter.com/HbsN5bTnwf

タグ: Julia言語 数楽

posted at 14:36:05

山本 卓 @sgryjp

19年5月12日

How to use redirect_stdout(f, [stream])? It's not easy to guess from the document...😑 #julialang

タグ: julialang

posted at 14:28:56

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu ミュートしているので読めない(笑)

タグ:

posted at 14:22:24

Ninja DAO | CryptoNi @CryptoNlnjaNFT

19年5月12日

前川さんの一連の話は、いろいろな事について「政権に問題があった」と言いながら、その後、影響力があったはずの次官退任までほとんど何もしていないので(退任はご存知天下り問題)、個人的にはどうなのかなぁ、と思っている。

タグ:

posted at 14:15:15

Ninja DAO | CryptoNi @CryptoNlnjaNFT

19年5月12日

うーん、だったら何でその後次官辞任まで何もしなかったんだ。 twitter.com/Kyoko_quetta/s...

タグ:

posted at 14:09:47

laskman @ilaskman

19年5月12日

田崎さんのバネの話、m’=m/Nとすれば色々整合しそう。連続的に質量が分布するとしたらどう扱えば良いのだろう?

タグ:

posted at 13:12:48

Akinori Ito @akinori_ito

19年5月12日

田崎先生に熱統計力学を説くとは,何者なんだ・・・?

タグ:

posted at 12:52:23

森下屋 菊助 @morishitaya

19年5月12日

@genkuroki 拝一刀のベビーカーの凄いところは、子供を守るためにガトリングガンや槍、薙刀が装備されています。

タグ:

posted at 12:49:20

梅崎直也 @unaoya

19年5月12日

昨日はベイズ統計の授業で二項分布の自然共役分布の話をして、ベータ関数の計算をしました。 pic.twitter.com/M7bYOpyEwC

タグ:

posted at 12:47:33

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

19年5月12日

本名不明、専門不明、論文の有無不明。これで学問について何か言えるのかわからん twitter.com/genkuroki/stat...

タグ:

posted at 12:39:57

黒木玄 Gen Kuroki @genkuroki

19年5月12日

きくちさんの所にはまだ来てないのか。

そのツイート1個だけじゃない点もすごいです。

小池一夫さんに【「子連れ狼」はフィクションなので、現実のベビーカー云々の根拠に使うのはやめたほうがいいと思いますよ】と言った人を見かけて以来の快挙。

伝説級の大型新人。

twitter.com/kikumaco/statu...

タグ:

posted at 12:20:20

なないお @Nanaio627

19年5月12日

特殊な環境を除いて、家族4人270万の生活で大学にいくまでの教育費が払える?満足な衣食住にも事欠く金額だよ。

タグ:

posted at 11:55:07

なないお @Nanaio627

19年5月12日

今の授業料免除を受けられる基準よりもはるかに低く設定してきて大学無償化とは。詐欺よな。

タグ:

posted at 11:52:50

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

19年5月12日

これをHal Tasakiに向かって言った強者がいるのか twitter.com/tjmlab/status/...

タグ:

posted at 11:52:03

なないお @Nanaio627

19年5月12日

大学の学費免除が家族4人で270万未満ということは、其れ位の収入の家庭の子供が大学に進学すらほぼ考えられないであろうことを知ってての設定だと思われ。つまりはほとんど国は払う気ないんだよ。
今の授業料免除を受けてる人たちの切り捨て。

タグ:

posted at 11:51:07

梶浦由記 Yuki Kajiura @Fion0806

19年5月12日

テープのり便利ですよね〜、でも案外細かいところはスティックのりの方が使いやすかったりしませんか……ヽ(´o`;。そんなこんなでまだ完全にテープのりに移行できずにいるのです……。

タグ:

posted at 11:18:08

S (ツイートはスレッド全体をご確認く @esumii

19年5月12日

33頁目や116頁目あたりなど。

タグ:

posted at 11:17:58

S (ツイートはスレッド全体をご確認く @esumii

19年5月12日

研究費あたりの論文数の話をするなら有名な www.janu.jp/report/files/2... を見れば「選択と集中」の失敗は明らかでは…

タグ:

posted at 11:10:18

某名無守 @Nanashino_kami

19年5月12日

誰も反論しないから奴らがつけあがった。誰かがガツンと反論すれば躊躇するだろう。というわけで代表的な3人をみせしめにした。つまり「舐められたら殺す」

タグ:

posted at 10:55:48

某名無守 @Nanashino_kami

19年5月12日

「史料による裏付けのない思いつきを『東大の日本史の教授』というかたがきで世間にばらまくことは、見方によっては百田氏や井沢氏よりも罪深い」
「学界で支持されていない自説を押し売りする前に、学界の主流説を世間に紹介することに努めるべき」
H郷先生は同じ立場なので容赦ない。

タグ:

posted at 10:53:51

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

19年5月12日

実験と証明は「雑談」ではないと思うんだがなあ twitter.com/curren49590380...

タグ:

posted at 10:43:49

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu 物理を分かっている人はもっとクリアにかつ地道に説明してくれます。私はすぐにミュートしました。

私は「久々の大型新人登場」だと思っていたのですが、案の定やらかしてくれました(笑)。

twitter.com/genkuroki/stat...

タグ:

posted at 10:39:23

ʇɥƃıluooɯ ǝıʇɐs @tsatie

19年5月12日

美しい。こういう事だったのか♬ とも思うた。分かってはいたけど。 twitter.com/woody_egg/stat...

タグ:

posted at 10:15:55

某名無守 @Nanashino_kami

19年5月12日

中央公論の呉座先生の記事読みました。
・百田→文筆家として知的体力が無い。口喧嘩には強いが論理的な文章は苦手
・I沢氏→あくまでも小説家
・H郷先生→学者というよりもはやタレント

三者のうち、H郷先生への批判が「ここ十年学界で言及されてない」「一般書の室も疑わしい」と一番厳しい

タグ:

posted at 06:09:26

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 物理をそれほど理解しているわけじゃないので、この人が言っていることが妥当かどうかさっぱりわからなくてスルーしていたのですが、ここから続くツイートを見て、まあ、いろいろ思いました。
twitter.com/sekibunnteisuu...

タグ:

posted at 06:07:03

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 未解決および手を付けていない部分

次数が2ずつ上がる
 ↓         
極限を取る   
 ↓         
波及速度はvを超えない

と出来ないか?

重力が一様ではない場合の、離散、連続に関しての考察

タグ:

posted at 05:58:15

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 今現在の私が理解した範囲

重力なし、一様重力で

離散モデル → 次数が2ずつ上がる
 ↓         
極限を取る   
 ↓         
連続モデル → 波及速度はvを超えない

タグ:

posted at 05:56:12

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 自然長の設定が本質的じゃないのはわかりつつも、考えるときは一応自然長を考慮しないと不安になる。何度もやるうちに「自然長は0としても構わない」と確信できるという感じです。

タグ:

posted at 05:52:23

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu #数楽 ただし、両端については慎重に考える必要があるかもしれません。

寝ます。

タグ: 数楽

posted at 05:45:25

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu #数楽 Δx=1/N, (k/m)Δx²=1, x=iΔx, L=λΔx, u(t,x)=u_i(t) とおいて、N→∞とすると

u_{tt} = u_{xx} - g(λx+u)

が得られます。y=λx+u とおくと、

y_{tt} = y_{xx} - g(y).

もとの離散版でも d²y_i/dt² の式内のLはキャンセルして消える。

結局 y_{tt} = y_{xx} - g(y) の基本形に戻ります。

タグ: 数楽

posted at 05:39:11

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu #数楽 続き。離散版は両端を除けば

d²y_i/dt² = (k/m)(y_{i+1}-y_i-L) - (k/m)(y_i-y_{i-1}-L) - g(y_i)

の形になります。Lはバネの自然長。これを

y_i= Li + u_i

で変数変換すると

d²u_i/dt² = (k/m)(u_{i+1}-u_i) - (k/m)(u_i-u_{i-1}) - g(Li+u_i).

続く

タグ: 数楽

posted at 05:34:19

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu なるほど、バネの自然長がゼロでない場合を考えているのか。確かにその場合の運動方程式は

u_{tt} = u_{xx} - g(λx+u)

の形になります。ここでλはバネ上の座標と自然長を関係付けるパラメーターです。バネの自然長が0ならλ=0になる。

そうか、この場合はまだ数値計算もしていません。続く #数楽

タグ: 数楽

posted at 05:27:02

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 離散モデルと連続モデルの関係もまだ十分に理解していません。
連続モデルでの、影響が速度vを超えて波及することない、というのを、離散モデルの次数が2ずつ上がるという結果の方から求めることは出来るのか?とかいろいろ気になります。

タグ:

posted at 05:08:48

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 物体をバネで数珠つなぎにした離散モデルの方、係数が0でない最低次数の項の次数が2ずつ増えていくというのも、ばね定数と質量が不均一の場合まではやったけど、重力が一定でない場合は、手を付けていない。

そっちを先に考えた方がいいのかな?

タグ:

posted at 05:06:15

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki バネが不均一に帯電していて、電界がある場合まで考えると

∂^2Φ/∂t^2=v^2*∂^2Φ/∂x^2-g(Φ(x,t))σ(x) まで一般化できそう。

とても手に負えそうもない。

タグ:

posted at 05:04:39

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki 結論としては、

この場合も、遠くで起こったことの影響が波及する速度は、∂^2Φ/∂t^2=v^2*∂^2Φ/∂x^2-g(Φ(x,t)) の場合、速度vを超えない

0<x、t=0で、v^2*∂^2Φ/∂x^2-g(Φ(x,0))が成り立っているなら、
vt<x、t=tで、Φ(x,t)=Φ(x,0)が成り立つ

となるものなのでしょうか?

タグ:

posted at 04:44:01

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki でも、どうせ2階微分だから、xがある位置、xからずれ、どちらをΦ(x,t)としてもよさそうですね。それなら、多少でも形が単純なg(Φ(x,t))の方がよさそう。といっても、やはり難しそう。

タグ:

posted at 04:40:05

積分定数 @sekibunnteisuu

19年5月12日

@genkuroki >∂²Φ/∂t² = ∂²Φ/∂x² - g(x)
ではなく、

そこはわかっているつもりです。

>∂²Φ/∂t² = ∂²Φ/∂x² - g(Φ)
です。

Φ(x,t)の取り方にもよると思うけど、xからのずれをΦ(x,t)としたので、
x+Φ(x,t)の位置にあることになるから、

∂^2Φ/∂t^2=v^2*∂^2Φ/∂x^2-g(x+Φ(x,t)) としました。

タグ:

posted at 04:34:26

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

19年5月12日

小池一夫【ベビーカーがしばしば問題提起されますが、江戸時代から拝一刀は大五郎をベビーカーに乗せてましたよッ!】

クソリプ【「子連れ狼」はフィクションなので、現実のベビーカー云々の根拠に使うのはやめたほうがいいと思いますよ】

この件もクッソ笑えた。

twitter.com/koikekazuo/sta... pic.twitter.com/Xx4cCT8KKQ

タグ:

posted at 02:04:29

黒木玄 Gen Kuroki @genkuroki

19年5月12日

この話を思い出した。

twitter.com/koikekazuo/sta...

タグ:

posted at 01:59:36

非公開

タグ:

posted at xx:xx:xx

黒木玄 Gen Kuroki @genkuroki

19年5月12日

あまりにも可笑しかったので記録に残しておく。

ミュート済みなので発見が遅れた。
すでに笑っている人達がいて存在に気付いた。

【理解のレベルの差ですねー】発言の返信先に注目(笑)。

また一つ伝説が生まれた!!!

この人はどこの何者? pic.twitter.com/1Ip6H8c87i

タグ:

posted at 01:42:45

闇のapj @apj

19年5月12日

@irobutsu 面白いです。中学理科,高校物理の授業を思い出しながら読んでいます。いつしか,自分がどうやってこういうややこしいのをクリアしてきたか忘れてしまっているので,振り返って考えてみることは大事かなと思いました。

タグ:

posted at 01:18:28

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu φ_{tt} = φ_{xx} - m²φ は例外的に線形でKlein-Gordon方程式。

φ_{tt} = φ_{xx} - sin φ はシャレでSine-Gordon方程式と呼ばれている。Sine-Gordon方程式はソリトン方程式なので色々特殊解の作り方が知られていますが、私は十分にいじったことがないです。今回の件に参考になる特殊解があるかも。

タグ:

posted at 01:12:48

黒木玄 Gen Kuroki @genkuroki

19年5月12日

@sekibunnteisuu 重力加速度gが高さに依存する場合のモデルは

∂²Φ/∂t² = ∂²Φ/∂x² - g(x)

ではなく、

∂²Φ/∂t² = ∂²Φ/∂x² - g(Φ)

です。xはバネ上の座標で、xに対応するバネ上の点の時刻tでの高さがΦ=Φ(t,x)なので、g=g(x)ではなく、g=g(Φ)です。

ここは結構混乱し易い所だと思う。非線形偏微分方程式。

タグ:

posted at 01:07:31

@genkurokiホーム
スポンサーリンク
▲ページの先頭に戻る
ツイート  タグ  ユーザー

User

» More...

Tag

» More...

Recent

Archive

» More...

タグの編集

掛算 統計 超算数 Julia言語 数楽 JuliaLang 十分 と教 モルグリコ 掛け算

※タグはスペースで区切ってください

送信中

送信に失敗しました

タグを編集しました