Add lunar gravity field disturbance

[200km]

Overview

Add lunar gravity field disturbance

Issue

• Expand the geopotential algorithm to Moon and/or Mars gravity potential #6

Details

I added the lunar gravity field calculation with GRGM1200A.
To correctly use this feature, we need to add the moon rotation calculation. #465

Validation results

I confirmed that the calculated disturbance is compatible with a reference code.

• n = m = 10, including the C00 coefficient.
• Moon fixed frame.

<style> </style>

	R  km	S2E m/s2	Ref m/s2	diff m/s
case 1	2000	-0.2358674	-0.2358674	-1.679E-10
	2000	-0.2358886	-0.2358886	-4.212E-10
	2000	-0.2359072	-0.2359072	1.3332E-09
case 2	-2000	0.36318726	0.36318726	-2.454E-10
	2000	-0.3632152	-0.3632152	-3.234E-10
	-1000	0.18164362	0.18164362	-1.219E-09
case 3	1000	-0.5595179	-5.60E-01	2.0208E-10
	-1500	0.83934119	0.83934119	3.7713E-09
	-1000	0.55972176	0.55972175	-5.092E-09

Scope of influence

Minor update

Supplement

NA

Note

NA

[seki-hiro]

実装ありがとうございます。approveします。
以下単純な興味です。

• 時間を計測しているようですが、referenceのツールと比較して早いとかは見えているでしょうか？

[200km]

時間を計測しているようですが、referenceのツールと比較して早いとかは見えているでしょうか？

今回の月に関しては計測していませんが、地球に関してはs2e-documentsに書いている通り、7倍程度早くなっています。ただ、アルゴリズムはほぼ同じだと思っていて、そもそもMatlabがC++に比べて遅いというのが主な原因だと想像しています。

https://github.com/ut-issl/s2e-documents/blob/develop/Specifications/Disturbance/Spec_GeoPotential.md

[seki-hiro]

matlabとc++の差が気になっていました。ありがとうございます。

[200km]

documentの比較はすごく真面目にやった訳ではなく、「実装したアルゴリズム遅すぎないよね」という確認程度なのでその点はご注意ください。
Matlabは色々最適化を試みるともっと早くできる可能性はあるとは思います。