随分前に購入した、「計算物理I」（夏目雄平・小川建吾、朝倉書店）を参考に、4次のRunge-Kutta公式の導出を行った。卒論で一度導出はしているはずだけど、もはや記憶にない（苦笑）
上記書籍をメインに計算を追い掛けて、1階の常微分方程式の簡単な問題について、解析解との比較までは上手く行ったが、理解の程度はまだまだなので、それ以上は繰り返し読み込むしかないだろう。Simpson公式と一致するところまでは、しっかり理解したいところである。
その後は、2階の微分方程式への適用である。「数値計算による量子力学」（桜井捷海、裳華房）を参考に、もう少し考える必要がある。ブランク空きすぎだ（苦笑）
偉い人が、自社の技術者のレベルが下がっていることを嘆いていた。一つは過去の経験に固執して新たな挑戦をしないこと、もう一つは、数学力がないことである。
数学力についての話を聞いていると、本質は、数学に限らず自分の頭で考えていないことのように感じられた。過去の経験に基づいて仕事をすることも、楽をしたいという気持ちの現れだと思う。
帰宅後、速度に比例する抵抗を受けながらの落下について、紙と鉛筆で計算をした。符号の向きで考え込んだりと、初歩的なところで躓いたりしたが、一通り計算して、matplotlibで解析解をプロットするところまではできた。明日以降は複数の方法で数値計算を行って、それぞれの方法を比較しようと思う。
Pythonでプログラミングをしていると、取り敢えず本に書いてあるコードをEmacsで打ち込んでも意味はなく、真面目に理論を追い掛けたり、odeintなどの使い方を理解しないと何も残らないのだと、改めて思う。こんなことは敢えて書くまでもないかも知れないが。
小さなことからコツコツと、である。
ここ数日、平面内の荷電粒子の運動を計算をするためにPythonでコードを書いていたのだが、明かに挙動がおかしい。参考にしている書籍と比較すると、一箇所符号が逆になっているようなのだが、その理由が分からなくて困っていた。
何か符号に関する勘違いをしているのだろうと思い、上記書籍に繰り返して目を通していたら、最初のところに「質点の持っている電荷は−1であるとします」って書いてある…電荷は+1だと勘違いしてた…俺あるあるだ。
その一方で、リストやタプルを使ったfor文の扱いは参考になった。最初は愚直にインデックスをつけたりしていたのだが、そうしなくてもできることなのね。自分で考えて書いたコードよりも格段にすっきりしている。
まぁ、こうした方が知識が身に付くと考えよう。
“Don’t think. Feel!”は「燃えよドラゴン」に出てくる有名なセリフだが、”Don’t think. Move!”というのはその派生なのだろうか。いずれにしても、「ストレスフリー超大全」（樺沢紫苑）の「不安を行動で取り除く」というのはもっともな指摘である。
「Pythonによる数値計算とシミュレーション」と「なんでもPythonプログラミング」を取っ替え引っ替え、たまに「独習Python」を読みながら、放物運動のグラフを描くところまでやってみた。敢えて本に書いてあることから逸れてみることで、理解が深まったところもある。とはいえ、独立して「趣味れーしょん」に載せるレベルには程遠いので、ここで計算結果を上げるに留めておく。45度方向に投げるのが一番遠くまで飛ぶよねという、ありきたりな結果でしかない。出射角は、青が35度、赤が45度、緑が55度である。
「趣味れーしょん」の方は、10年間未更新になることを回避できれば良いと思っているが、あと2年ほどあるので、それは大丈夫だろう…多分（笑）
ついでと言っては何だが、これまでMacのメモを使っていたここの下書きを、nDiaryを使っていた頃のようにEmacsに戻した。プログラミングはEmacsですることになるので、プログラミングの時間が増えるのであれば、エディタはEmacsをメインにした方が都合が良い。
ここ数日、Pythonで文字列の処理をしようとあれこれやっている。文字列中に、全角かなから半角英数に切り替わる箇所が複数あって、そこで文字列を分割したいのだが、今のところ、切り替わる箇所を検出して、そこにセパレータになる文字を挿入するところまでは、できている。その後にsplitを使えば良いだろう。
次の問題は、今は上記の場所が1つだけある場合にしか対応していないので、それを複数箇所対応にすることだが、もう数日かかりそうだ。その後は、csvファイルを行毎に読み込んで、そうして得られるリストの特定の要素に、上記の処理を行うことが最終目標である。
10年くらい真面目にプログラミングをしていないので、この程度の処理でも結構大変だし、自宅ではEmacs、職場ではxyzzyの操作が少々ぎこちない。しかしながら、色々と心に火がついた（ような気がする）し、仕事で使うためのものなので、もう一踏ん張りである。
HDD2個とEye TV HDを、小型家電のリサイクルBOXに出してきた。片道2kmの散歩であった。
HDDの1つはTime Machineで使っていたのだが、先日動作中に異音を発したので、新しいHDDと入れ替えて、お役御免である。多分、5年以上使っていたと思う。もう一つはLinux Boxで使っていたものなので、Time Machineで使っていたものよりも古い。
Eye TV HDは、数年前にDTCP-IP対応のソフトを導入したことで不要になったが、取り敢えずSTBとMacに接続した状態で放置されていた。Time Machine用のHDDの隣に置いてあったので、丁度良い機会だと思い、処分した次第である。
また、やんごとなき事情でモニターのHDMIポートが一杯になりそうなので、Mac miniをDisplayportで接続するため、USB-C/Displayport変換ケーブルを注文した。
一本ケーブルがあったと思ったて確認したら、mini-Displayport/Displayport変換ケーブルだった。Mac miniとモニター間のHDMI接続が不調の頃に買ったのだが、変換ケーブルを買ったら不具合が収まったので放置したため、間違った買い物をしていたことに今頃気がついたという、間抜けな顛末である。
実際に使ったかどうかは記憶が定かではないのだが、多重ループを一気に抜ける場合の方法は、gotoを使うことしか知らなかった。
ところが、である。Pythonの本を読んでいたら、elseとcontinueを使う方法と、フラグ変数を使う方法の2つが出ていて、軽く衝撃を受けた。前者はelseだけを使うという意味で、後者はC言語でもgotoより一般的な方法らしいのに、今まで知らなかったことに…
年末年始で仕事関係の本を一冊読み切った。頭の整理に役立ったと思う。
Python、物理、深層学習、データ分析、ゲノム編集といった本にも手を伸ばしたが、いずれもつまみ食い程度でしかない。こういうことにもっと時間を使えたらなぁと思うのだが、現実を見なければいけないだろう（苦笑）
実家以外で年末年始を過ごすのは初めてなのだが、それゆえ初めて見る光景がいくつかある。
近所に小さな神社があるのだが、そこに参拝客が来ているのを見るのは初めてだ。「密」を心配するほどではないけれど、参拝客は途切れることなく来ているようだ。常に観測していた訳ではないけれど、三が日の間、何度か前を通り掛かった際には、そんな感じだった。
もう一つは、チンドン屋のような音楽を流しつつ歩いている集団である。その人たちの出立からは獅子舞には見えなかったが、あれは何なのだろうか。謎である。
いずれにしても、「昔ながらの」的な感じが新鮮だった。
書くだけで終わる年が多いのだが（苦笑）、明確に書いておくか否かでその後の行動にも影響が出るということなので、ここで書いておく。
１．仕事
既に少し手をつけているが、業務の中で新しいことを幾つかやろうとしている。データ分析や、それに付随するスキルを身に付けて、経営層の目に留まる結果を出したいと思っている。懸念点は、古株の人からすると、自分がやってきたことが否定されたと感じかねないことだろう。既にその傾向は見られるが、彼らの「言われたことをやってる分には責任を問われない」という発想が経営者から否定されたとも言えるので、自業自得だとも言えるんだけど、理屈ではなく感情の問題になってしまうので、そこが難しい。
もう一つの懸念点は上司だろうか。人・モノ・コトのいずれについても、その存在によるメリットだけを享受して、デメリットは私に押し付けてお仕舞いという発想なのである。この場合、古株の扱いを私に押し付けて逃げかねない。上司の上司をどう味方につけるか思案中である。
２．資格
試験に合格することは言うまでもない。例年より2ヶ月遅れになることが公表されているので、これを活用したいところではある。もっとも、やる気の問題が一番だというのは重々承知しているのだが。
３．趣味
Pythonに慣れたい。仕事のデータ分析とも関係するが、数値計算、深層学習やRaspberry Piにも手を出したい。本などは前から買ってあるんだけどねぇ（苦笑）
色々と書いたけど、どこまでできるかなぁ。