Studio Life～2023 New Year’s Eve

Studio Life～2023 New Year’s Eve

年末年始は鎌倉の本宅で年越しです。小松駅から夜行バスに乗り翌朝新宿で私鉄に乗り換え最寄り駅片瀬江ノ島駅で降車。改札口に出て太陽の光に輝く海の景色を目にする瞬間、仄かな懐かしさを覚え、関東の天候の良さを思い出して鎌倉へ戻ってきた感激に浸る。今年の年末は晴天に恵まれて寒さ知らず。生憎大晦日は天気は崩れたが正月三ヶ日も好天に恵まれる予報。I am spending the New Year’s holiday at my main residence in Kamakura. I take an overnight bus from Komatsu Station, change to a local train at Shinjuku Station next morning, and get off at Katase Enoshima Station. As soon as I step out of the exit gate and see the sea shining in the sunlight, I feel a faint sense of nostalgia and recall the excellent weather in the Kanto area. I am filled with the emotion of returning…

Read More Read More

SMARTSOL〜The AIM Implementation

SMARTSOL〜The AIM Implementation

第二回目は順応陰解法(AIM)の実装編をお届けします。AIM法の陽･陰解法への切り替えについて閾値法と局所行列安定法が提唱されている。This second thread presents an implementation of the Adaptive Implicit Method (AIM), in which the Threshold Method as well as the Local Matrix Stability Analysis Method are proposed as the switching criterion between the explicit and implicit AIM method.

SMARTSOL〜The AIM Matrix Solver

SMARTSOL〜The AIM Matrix Solver

“SMARTSOL”の重要論点の一つである順応陰解法(AIM)に関してシリーズでご紹介します。第一回目は行列解法に於けるAIM法の取り扱い方法について説明します。I would like to address about the Adaptive Implicit Method (AIM) in series, which is one of the main functions incorporated in SMARTSOL. In this first thread, we will discuss how to implement the AIM method in the matrix solver.

JAT〜深田久弥ゆかりの山々

JAT〜深田久弥ゆかりの山々

深田久弥の原点を探る久弥ゆかりの山行便りをお届けします。日本一周旅行で取り残しとなっていた岐阜は戦国時代の群雄割拠の地であったことからも何れは訪れたいと思っていた。その古戦場城巡りの帰路、能郷白山に登頂した際、余りの眺望の良さに驚き、何故深田久弥は能郷白山を落として荒島岳を百名山一山に選んだろうかとの素朴な疑問に駆られたことにある。I wish to report a mountaineering to exploring the origin of Kyuya Fukada. Since Gifu is left out of the Japan trip, I hope to visit Gifu in the future because It is the wellknown land of warlord rivalry in the Sengoku period. On the way back from visiting the old battlefield castles in Gifu, I am surprised by the panoramic view from the top of Mt. Nogo Hakusan and wonder why Kyuya Fukada choses Mt. Arashima-dake as one of…

Read More Read More

Godzilla R32〜Ultimately Impulse Buying on New Steering

Godzilla R32〜Ultimately Impulse Buying on New Steering

Godzillaに装着されているステアリングはMOMO製のDeep Cone 360mmである。本革は手に馴染んでいるがブルー色のスポークが気にいらず、純正の程度の良い中古の入手は難しかろうと何れパゴダ同様Nardi製クラシックウッドに変えたいと考えていた。日課の早朝散歩の足湯で寛いでいるとAmazonで異常な低価でNardiが販売されていることを見つけ、ついポチってしまう。The steering wheel installed in Godzilla is a MOMO Deep Cone type. The real leather is comfortable to my hand, but I do not like the blue colored spokes. Since I think that it would be difficult to get a real genuine one, I want to change the wheel to a Nardi classic wood one day, just like Pagoda. While relaxing in the foot bath during my daily early morning walk, I find Nardi on sale at Amazon…

Read More Read More

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.