AIに協力してもらって、CFDソフトを開発する試みです。 数値計算の知識は少しありますが、FORTRANでしたので、 現代的なコード（C++、Python ）あたりで開発していきたいと思ってます。 C++の基礎については以下から数回で実施しました。 今回から具体的な熱 ...

static inline IVP_DOUBLE get_closing_speed_core_i(const IVP_Impact_Solver_Long_Term *,int core_index, const IVP_U_Float_Point *rot_sped_cs, IVP_U_Float_Point *speed_ws); // based on current core speed ...

今回はPythonでの数値計算上の解析テクニックの話です。 今回の記事では、イベントの機能を使って「数値解の速度がゼロになる時刻」を求めます。 まずは結果のプロット画像です。 一連の記事では微分方程式の数値計算に solve_ivp () という関数を使っています。