陰的時間離散化を用いたcudaおよび数値解法

（時間の離散化のために）IMPLICIT形式の有限体積法によって偏微分方程式（PDE）のセットを解決するコードを移植したいと考えています。

その結果、ADI / TDMAスキームで処理されるx、y、z方向の3重対角方程式系があります。

CUDAを使用したPDEの暗黙的な解決策については、何も見つからないようです。

ADI / TDMAスキームはCUDAで実装できますか？2D熱拡散方程式のような例はどこかにありますか？

私が見つけることができるのは、有限差分であるがEXPLICIT形式（ケンブリッジ大学）の2D熱拡散方程式のCUDAサンプルコードだけです。

ヒント/リファレンスをいただければ幸いです。

この問題は、高度にベクトル化された形式に役立ちます。お気づきのように、ADIメソッドは、3対角システムのいくつかのステップを与えます。これは線形方程式の形式であるため、CUsolverとCUblasを使用して、標準の線形代数ルーチンの並列GPUバージョンを呼び出すことができます。これらを使用して、明示的なコードを取得し、内部ループを適切なCUsolver呼び出しに変更して、コードがCPU実装とほとんど同じように見える方法で解決する必要がありますが、GPUで行列演算を実行すると、ライブラリ呼び出し。