難しい

　先週から取り組んでいたGRMHDコードへのHLLD法の導入だが、半分断念した。いろいろと基礎方程式をいじくって、RMHDに適した計算方法を作ろうとしたのだが、どの計算方法も実際にコードに組むとうまくいかない。いろいろと調べたが、コードのバグではなく、もっと根本的なところに問題があるようだ。RMHDでのHLLC法からHLLD法への拡張はそう簡単ではないな。何かもっと新しい考え方が必要だ。

　私の予想だが、Mignore達も同様にHLLD法の導入を考えているだろう（少なくとも彼らの作ったコードでは非相対論MHDでHLLD法が導入されている）が、私と同じように躓いていると思う。今後、ちょっとしたときに考えてみようとは思うが、現状ではかなり厳しいな。