JAT〜Mountaineering Year 2023

JAT〜Mountaineering Year 2023

以前道に迷って敗退した空木岳はテント泊が余儀なくされるので断念していたが諦めきれない。遂に、体力の低下が著しい昨今、早いうちに再チャレンジすることを決意。手っ取り早く駒ヶ岳ロープウェイを利用して木曽駒ヶ岳から空木岳へ南に縦走する体力の消耗出来るだけ軽減するルートとする。総距離約２５km。やれるだけやってみるつもりで望む。Mt. Utsugi, which I had to give up because I lost my way, is difficult to climb because the tent is needed for sleeping outdoors. I cannot give up on Mt. Utsugi and decide to try again as soon as possible, as my physical condition is deteriorating considerably these days. I take the Komagatake ropeway as the quickest and easiest way, and then traverse south from Mt. Kisokomagatake to Mt. Utsugi to reduce physical exhaustion as much as possible. The…

Read More Read More

SMARTSOL〜The CPR Preconditioner

SMARTSOL〜The CPR Preconditioner

前処理法(Preconditioner)としてIncomplete LU FactorizationとNested Factorizationを紹介してきた。自分の現役時代のシミュレーション環境ではこれにて充分であったが、最近より複雑な地層を対象とする大規模な油層シミュレーションの環境下ではより効果的な解法としてCPR前処理(Constrained Pressure Residual Preconditioner)なる手法が用いられているようだ。まだin-house貯留層シミュレーターSMARTのMatrix Solution Package (SMARTSOL)には実装していないが、その概要及び実装の必要性等につき考察する。

iSMART〜The Hyper-Dual Numbers in Fortran

iSMART〜The Hyper-Dual Numbers in Fortran

４ヶ月に亘る冬籠りでPython，Juliaの自動微分の実装を終え、Fortran Simulator “iSMART”をJulia言語に変換中にある。スタンフォード大の博士論文を渉猟している折、Hyper-dual numbersなる超実数を用いて厳密に一階偏導関数、加えて二階偏導関数も導出できる論文に出会う。急遽変換を中断し、FortranによるHyper-dual numbersの実装を試みることにする。あの冬籠りは一体何だったんだろうか？After four months of dedication during the winter, I finished the implementation in Python and Julia for the automatic differentiation. I am currently working on converting the in-house Fortran program “iSMART” to Julia. While browsing papers at Stanford University, I came across a Ph.D. thesis on hyper-dual numbers, which is an extension of real numbers that can derive the first-order partial derivatives as well as the second-order partial derivatives in a strict manner. Suddenly I stopped…

Read More Read More

iSMART〜Automatic Differentiation in Julia Final

iSMART〜Automatic Differentiation in Julia Final

九谷広場でJulia版ADモジュールの開発に悪戦苦闘していた日曜の午後、騒がしい排気音に集中力が途切れがちになる。当初は気に止めずにいたが女性陣が騒ぎ出したのを機に外に目を向けるとロゴを付けた旧車が何台も駆け抜けている。何事だと外に出て交通整理の係員に問いかけると、初めて北陸で開催される車イベントで堺正章さんも参画。今晩ここに泊まって明朝ゴールの京都へ向かうと聞く。目を凝らしてロゴをみると”La Festa Primavera 2023″とある。あの本邦一格式の高いイベントかと驚く。明朝出発セレモニーで往年の名車に出会えると思うと、ハイな気分になり一気にモジュールの開発が捗り完成する。さて、Julia版ADモジュールのご紹介を兼ねたJulia言語概要解説編をお届けします。On a Sunday afternoon, while I am struggling to develop Julia’s version of the AD module at Kutani Square, the loud exhaust noises cause me to lose my concentration. I don’t pay attention at first, but then I look outside and see a bunch of old cars with logo stickers driving by as the women start to make a fuss. I go outside to ask an event staffer at the traffic control what’s going on. It turns…

Read More Read More

iSMART〜Automatic Differentiation Part 2

iSMART〜Automatic Differentiation Part 2

Pythonによる自動微分クラスの構築、実装の検証編をお届けします。単なるプログラミングな知見を駆使するだけで自動微分の肝はPart 1に集約されます。Part 2 would hopefully address an implementation and verification of the Automatic Differentiation Initiative class in Python. A mere programming skill is required in this part. The core of AD could be covered in Part 1.

iSMART〜Automatic Differentiation Part 1

iSMART〜Automatic Differentiation Part 1

冬籠りの三ヶ月間どっぷり自動微分(Automatic Differentiation)に浸かりました。俺の一世代前の貯留層シミュレーターiSMARTは導入の容易な差分商で近似する数値微分法を採用しているが、演算に大幅な手間が掛かる上に最適な微小値\(\epsilon\)を見出すのが難しく誤った解に収束することも起こりうる状況にある。近年研究対象のシミュレーターでは自動微分を用いてヤコビアン行列の組込みの手間の短縮化を図っていることを知り自動微分の組込みを夢見ていた。野外活動が鈍る北陸の冬の時期、iSMARTに自動微分の実装に向けてチャレンジを開始。Three months this winter are dedicated to the automatic differentiation (AD). My primitive reservoir simulator iSMART adopts the numerical differentiation method which is easy to implement but may end up converging to a wrong solution due to the loss of the search for an optimal \(\epsilon\). I wish to incorporate an AD function into iSMART, once I learn that the latest research simulators tend to use the AD function instead of spending huge time on constructing analytical Jacobian matrices….

Read More Read More

Studio Life～2022 New Year’s Eve

Studio Life～2022 New Year’s Eve

今年は３年振りに鎌倉の本宅でファミリー全員集合し総勢９人で年越しです。愚妻は食事の用意に忙殺され、愚娘息は悪ガキの面倒に掛りきりの一方、俺は何をしていいか分からず手持ち無沙汰でボーッとして回りからヒンシュクをかっています。大晦日恒例の今年の景勝地は観光客が戻ってきたスラムダンク踏切です。This year, for the first time in three years, the whole family get together at our home in Kamakura for the New Year’s Eve with nine people totally. The wife is busy preparing meals, the daughter and son are busy taking care of their kids. I am frowned by the family when I am bored out of my mind, not knowing what to do. The scenic spot for the New Year’s Eve is the Slam Dunk railroad crossing which is…

Read More Read More

Studio Life～酒を切れない

Studio Life～酒を切れない

毎晩の飲酒はアルコール認知症・依存症のリスクを促すとの通説にとらわれて、一時白内障の術後二ヶ月間は禁酒したものの、飲酒は週４日(水、金、土、日)に控えているが依然として酒を切れない。しかし、飲酒の翌日は活力が出て体調が頗る良くなる経験をするので、体調や気分のすぐれない時などは一週間飲み通して体調の回復をはかることもある。通常焼酎乃至ウイスキー２〜３合を専ら水割りにして飲む。一番好きな樽酒日本酒やShiraz赤ワインは大事な機会に飲むことにしている。

Godzilla R32〜The Spare ECU

Godzilla R32〜The Spare ECU

以前パワーテストの測定を試みたショップでECUが突如狂ってエンジンが停止する恐れがあり、確かR33用のECUを流用する方法があると聞き及んだ。その類のショップとは縁を切ったのでR32用の純正ECUと付き合っていくしかない。ECUの高騰ぶりも顕著で俺の手が出るレベルのものは稀にしかオークションに出てこない。気長に探している矢先。