なかのにっき

電子小説を読む。

昼前に布団のクリーニングついでに大学に出てきて書類提出と準備、金曜日に残したものの回収。

午後は買い物で三鷹まで。とても暖かい。 水浴びするハトを観測 。その後夕食の仕込み、読書の続き。

調整して提出。

雨降り。昨日からうってかわって寒い。

＠ ML 学科会議:

＠ EE 会議:

＠ 教授会:

＠ 懇親会:

最終配分があって、どうやらほぼ軟着陸できそう。 Web 財務と Excel の突き合わせ確認、各種伝票が揃ったので提出。 あとは卒論の製本を出して、いまもっている領収書を出せば 2〜3 万円くらい。メモリでも買っておけばいいか。

＠ 電源修理代:

＠ ソフト代:

顧問をお引き受けしているのだけど、代替わりということで来室いただいた。よろしくお願いします。

面談日程の調整など

がんばって読んだ。心証はだいたい固まったので明日レポートを書く。

それぞれにメール。

ぼちぼち返信来て対応。

レポート送った。やれやれ。

表紙のフォントとか目次のレイアウトとかの気付いたところをいじりつつ印刷。 体調不良による遅延分 1 名を除き、ブックセンターに提出した。

decision letter の cc が届いた。ふたりめだったか。

4 月から斎藤先生の後任として着任される先生が 14:00 にいらして、斎藤先生との引き継ぎ、入る研究室の見聞、学内のご案内等々。夕方から柴田先生が合流し、 Brasserie Parisienne YG で晩ご飯。おいしかった。いろいろあって 2 名で赤を 1 本空けたので、帰宅後に頭痛が。

• 担任案件。メールのあと電話も。
• SP 部会
• 授業関連の問い合わせ
• 学会発表の打ち合わせの件
• 製本の件
• WO3 実験の件

などなど。

過去問を拾って解答つくった。

家でダラダラしていると健康によくなさそうなので歩いて買い物へ。 おにやんま、ジュンク、西友。

昨日出したメッセージに返信がきてたのでそれらに再度返信。
…こうやってメッセージが指数関数爆発しがち。抑制せねば。

解説追加。あとはシラバス見てテキスト入れたら納品、と思ったら月曜までポータルが更新作業中だったので月曜に。

そんな速いわけじゃないけど、去年導入した新 migrate が Intel(R) Xeon(R) E-2324G CPU @ 3.10GHz だったのでここで。ひとまず gfortran と libmkl-full-dev を apt install。

コピーしてボチボチコンパイルできるかテスト。 gfortran で 77 をコンパイルするときは -std=legacy をつければ警告なくコンパイルできる (参考こちら)

ifort と比較したときの gfortran のコマンドラインオプションは以下に。 いくつかのFortranコンパイラのオプションの整理

MKL をリンクしようとしたらさすがに死んだので、以下をちゃんと読もう。
Tips about Using Debian's MKL Package / FAQ

いい天気になった。買い物にいったくらいで、あとは本読み。

• 某依頼のあった書類をつくって送信
• 某問題サンプルを完成させて送信
• 新任の先生の荷物の件と python の件を確認して送信
• 花束の件振り込みして車の状況を送信
• 製本の件調整して印刷して持参

MKL を使っているのは spline 関数をつくるときの疎行列計算（だけ）なので、代替できるものはありそうなんだけど、ひとまずありものでビルドできるか試す。ChatGPT さんに聞いてみたところ、

gfortran your_code.f -o your_program \
  -I${MKLROOT}/include \
  -L${MKLROOT}/lib/intel64 \
  -lmkl_intel_lp64 -lmkl_core -lmkl_sequential \
  -lpthread -lm

でいけるでしょう、とのこと。Debian の mkl ライブラリは /usr/lib/x86_64-linux-gnu/ に入って、これは gfortran の標準的なライブラリ探索パスに入っているので (gfortran -print-search-dirs で確認可能）、

FC = gfortran
FFLAGS =  -g -Wall -std=legacy

LIBMKL = -lmkl_intel_lp64 -lmkl_core -lmkl_sequential -lpthread -lm

LFLAGS = $(LIBMKL)

STARGET2 = s2dspf
SOBJ2 = s2dspf.o
OBJ2 = spline2d.o

$(STARGET2): $(SOBJ2) $(OBJ2) $(LIBNAKANO)
        $(FC) $(FFLAGS) -o $@ $(SOBJ2) $(OBJ2) $(LFLAGS)

という感じでビルドできた。

ひとつ割れてしまった義歯があって受診。ほかにも欠けているところが3箇所あるとのことで、 しばらく通うことになりそう。とほほ。クリーニング。

青柳先生に同席いただいて 5 名の面談。

領収書をすべて提出。シュルズ先生のぶんも検品して提出まで。

ということでお願いします、というところで。

MKL と spline 周りについてはできた。デモの結果を確認するのは gnuplot だったけど、そこは ChatGPT さんの力を借りて numpy+matplotlib に変更。

関東支部に出す話、来年の電源事情・予算事情など、彼らの就職事情などなど。

もうそんなタイミングか。今年は15クラス、2行目か学生名が入っているので、 ChatGPT さんの教え の通り

=mod(row() - 2, 15) + 1

して番号を振ってコピー→数値貼り付けして納品。

ついでに名簿・名札作業もやっておく。これは 去年の記録 の通り。ただし pdftk で PW をつける作業は $f 前後の " をとらないと駄目だった。bash の仕様変わったの？

＠ ラベル屋さん:

体調不良で遅くなったお一人ぶんを刷って提出してきた。

つくる。pdf ひとまとめにするのは「pdf を結合」機能で簡単、「資料 1XX-1」みたいなラベルもコピペでいける。頁番号をふるのは以下のように。

1. ツールバー [編集] をクリック
2. 右ペイン「ヘッダーとフッター」→「追加」をクリック
3. 「ヘッダーとフッターを追加」ウィンドウの「右フッターテキスト」に "P<<1>>" と記入。
4. 同ウィンドウの「フォント名」「サイズ」を編集。SP部会資料では MS Pゴシックの 14pt。

ほんと adobe はこのへんの UI コロコロ変えよるから web の資料が探しにくくて参る。

16:00-16:20 くらいで。資料つくらないとなあ。

授業の予定を放り込み。

午前幹事会、午後研究会。5件すごく面白く聞いた。

賈先生と鈴木先生、清水先生と Soba & Co で少しだけ打ち上げ。こちらも楽しかった。

近先生が 12 日に急逝されたとの連絡をいただく。 去年退職されたばかり だったのに…

葬儀の情報 。予定ない日なので参列しよう。

10:30-11:30 on Zoom。いっぺん切れたのはなぜだ。

メールたくさんとレターパックと名刺のスキャン。

乱数周りのテストをすると謎のメモリリークが出たり SEGV が発生したり。 いろいろ調べた（迷った）結果、使ってた Mersenne Twister のバージョンが 32 bit 版だったことに気付く。駄目。

明日へ続く。

64bit MT に対応。

1. 64bit版 Mersenne Twister から mt19937-64stdint.tgz をダウンロードして展開、mt64.h と mt19937-64.c をコピー。
2. 以下のような Fortran 用のラッパー mt64-c2f.c を書いた。

   #include <stdint.h>
   
   extern void init_genrand64(unsigned long long);
   extern unsigned long long genrand64_int64(void);
   
   void sgenrand_(const int * seed){
       init_genrand64((uint64_t) (* seed));
   }
   
   /* it is an identical copy of genrand64_real1() */
   double genrand_(void){
       return (genrand64_int64() >> 11) * (1.0/9007199254740991.0);
   }
   

3. これらを

   gcc -c mt19937064.c mt64-c2f.c
   ar r libfuga.a mt19937-64.o
   ranlib libfuga.a
   

   してライブラリファイルにした
4. Fortran のソースからは

     :
         int seed = 1
         extern genrand
         call sgenrand(seed)
         write(6,*) genrand()
     :    
   

   のように書いて hoge.f とし、

   gfortran -std=legacy -o hoge hoge.f -L../lib/fuga -lfuga
   

   のようにビルド

以上で無事に呼べた。やれやれ。 ライブラリを削除してなくて 32bit のときの genrand_ sgenrand_ が残ってたりしてここでも紆余曲折。だめ。

というわけで nrnd() maxwell() のテストルーチンが確認できた。ヒストグラムと理論曲線の対応を描かせるグラフが ChatGPT + python (numy, matplotlib) で一撃でできちゃうのは本当にすごいなあ。

働こうかと思ったが、掃除・読書・買い物でおわり。公孫龍の4巻を読み終わり。あとはぶんころり先生のラノベが安かったので買ってアタマから。

?202503b#202503141 こちら の件。中央線で救護活動が入って 5 分遅刻。だめ。

喪主ご挨拶で奥様から経緯をお聞きした。慢性心不全のため、別トラブルの手術ができなくて、とのこと。 お子さんのご兄弟とも大きく独立しておられ、お孫さんもいらっしゃる…というところが多少の救いかなあ。でも早いよなあ。

吉祥寺まで富谷先生・坂本先生と電車でご一緒。

子供の前期授業料の振込。みずほに行ってみたら 1h 待ちと書いてあったので、 りそなに切り替えたら早く終わった。

春木屋で昼ごはん。

eeject を試したら、Thompson 分布と棄却法で発生させたデータのヒストグラム（正規化）が合わない。うーんうーんと唸っているうちに時間切れ。

今日から3日間。作業をしつつ聞いているとあまり入ってこないな…（すみません）。

砕けた歯のところを再建するので型取り。2週間後。

どうも python で生成させたヒストグラムがおかしいっぽい。 昔 perl で書いた、数えて作るヒストグラムだと合ったので、それでよいことにする。

16:00 過ぎから 17:30 くらいまで。当日のお仕事はなくてよさそう。

ターボ等を頂戴するなど。いったん 14-304 室に退避。

残念ながら午前は雪。卒論の製本を取りにいかないといけなかったりで、体育館での式典はオンライン中継にて。あと大学に置きネクタイがあると思っていたらなかったので、家にいっぺん取りに帰ったりもしたので(笑)

学科別の学位記授与式、部屋に戻って研究室の記念撮影、大学院の学位記授与式、大学院の授与パーティー、と参加。

3 日目。講演それぞれ大変面白かった。論文読もう。

昼は買い物くらい。床屋行こうと思っていたんだけど今日は休みだった。 夜に中学 PTA ミッション。おつかれさまでした＞誰

8:50 に登校、もろもろのミッションを無事終えた。やれやれ。

火曜に続き、齋藤先生にお声がけいただいた物品の移動。 サーモ理工の加熱装置、SEM 用の Au コーター、シュルズ先生が持っていく蒸着装置など。 今日は人手がなかったので、引き続きいったん 14-304 に保管させていただく。

整理して共有。

渡邊くんのアブストラクトを edit。

家族の送迎と自分の眼科。

調整して申し込みしてもらった。

申し送りのメモを書いて委員のみなさんに LINE で共有。

いつもの敷島館へ。家族１名が体調不良で実施が危ぶまれたがなんとか復調。

夕食の記録(笑)→ Quora: 食べたものを報告するだけのスペース。

敷島館２日目。昼も宿の部屋で蕎麦を頼み、ひたすらダラダラする予定だったが、 子供にリクエストされてちょっと外出。あとは本読み。

10 時にチェックアウト。うまい踊り子号はもうなくて、東海道線の普通列車で帰宅。

読んだ本の記録→ Quora: 読んだ？ - 温泉宿で読んだ本

• 子供の高校進学に伴う寄付金を振込して書類を送付
• PTA の引き継ぎ資料を先生方に送信
• 最後の請求書を起票
• 指導いただいた先（みっつ）への修論卒論の送付
• シュルズ先生がいらしたので鍵回収・引き継ぎ・今後の予定の打ち合わせ
• シュルズ先生に返却いただいた XPS13 の整理。しかしこりゃ無理か？
• 担任のお仕事
• SP部会賞関連の作文と調整
• 新入生・２年生ガイダンス用の資料編集と調整
• とある御縁で教えていただいたプレプリを拝読して返信

というわけで本日はおおむね時間切れ。うーむ。

ということで testing からなくなってた。ほえー。

再開。peject のチェックルーチンを整備。これは、エロージョンの動径分布の離散点 {r_i, d_i} から、{r_i, r_i d_i} をもとにスプライン関数 f(r) を作り、これを分布関数とみなして粒子放出点の r を棄却法で決め、θ は 0~2π でランダムに発生させている。

これが正しく機能しているかどうかのチェックとしては、

1. f(r) を rmin~rmax まで積分して 2π をかけた値を Ve とする。これはエロージョン領域の体積にあたる。
2. 発生させた粒子の出発位置から r を求め、rmin~rmax を bin に区切って数をカウントし、r_k~r_(k+1) の間の面積 π({r_(k+1)}^2 - {r_k}^2) と試行数 N で割り、相対密度分布 c(r) とする。
3. {r_i, d_i / Ve} と c(r) をグラフで比較する

をすれば良い、と思う。結果はこんな。たぶんだいじょうぶそう。

伊東先生にお誘いいただき、東京の寿司居酒屋で伊藤先生・相楽先生＋二次会合流村松先生、にて。 たいへん美味しゅうございました。村松先生と帰りの電車で色々話せてよかった。

村松先生以外からもいろいろヤバい話を聞いた気がするが忘れたことにしようｗ

ほんとは 2 月にしないといけなかったやつ。特に問題なく 1.5h で完了。

大学にきたら桜が咲きかけていたので 写真 を撮ってみた。

このあと ID カードを忘れたのに気付いて家まで往復することになったのであった…

床屋へ。寒い雨降りだったのでズルズルしていたら午後の遅い時間になってしまった。

買い物と家の部屋の片付けと書類整理。そろそろ回転数を上げなければ。

遅まきながら HiPIMS Today 2025 のメモから論文を pick up。昨今は Open Access なやつが多くて助かる。

謝罪のご連絡と業務プロトコル改善のご依頼。