10124913218

Committed 27 Jul 2024 04:27PM UTC coverage: 12.208%. First build

Build # 10124913218

Build Type

Pull #3794

github

Committed by

gzotti

Commit Message

Reduced to max.4 additional threads

Pull Request Pull Request #3794: Parallelize ephemeris computation with C++ parallelism.

Run Details

208 of 369 new or added lines in 24 files covered. (56.37%)

14439 of 118271 relevant lines covered (12.21%)

19314.1 hits per line

Source File
Press 'n' to go to next uncovered line, 'b' for previous

98.67

/src/core/planetsephems/marssat.cpp

1	/************************************************************************
2
3	The Ephemerides of the Martian satellites
4	(adjustement from 1877 to 2005, Version 1.0)
5	by Valery Lainey can be obtained from Valery Lainey:
6
7	V.Lainey (Lainey@oma.be)
8	ROB- 3, Avenue Circulaire, B-1180 Bruxelles (Belgium)
9	IMCCE - 77, Avenue Denfert-Rochereau 75014 Paris (France)
10
11	-----------------------------------------------------------------------
12
13	I (Johannes Gajdosik) have just taken Valery Laineys Fortran code,
14	MarsSatV1-0.f, which he kindly supplied, and rearranged it into
15	this piece of software.
16
17	I can neither allow nor forbid the usage of Valery Laineys
18	Ephemerides of the Martian satellites.
19	The copyright notice below covers not the work of Valery Lainey
20	but just my work, that is the compilation of Valery Laineys
21	Ephemerides of the Martian satellites into the software supplied in this file.
22
23
24	Copyright (c) 2006 Johannes Gajdosik
25
26	Permission is hereby granted, free of charge, to any person obtaining a
27	copy of this software and associated documentation files (the "Software"),
28	to deal in the Software without restriction, including without limitation
29	the rights to use, copy, modify, merge, publish, distribute, sublicense,
30	and/or sell copies of the Software, and to permit persons to whom the
31	Software is furnished to do so, subject to the following conditions:
32
33	The above copyright notice and this permission notice shall be included
34	in all copies or substantial portions of the Software.
35
36	THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
37	IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
38	FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
39	AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
40	LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
41	OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
42	SOFTWARE.
43
44	My modifications to the "Ephemerides of the Martian satellites" as implemented
45	in MarsSatV1-0.f are
46	1) do not calculate constant terms at runtime but beforehand
47	2) unite terms with the same frequencies
48
49	****************************************************************/
50
51	#include "marssat.hpp"
52	#include "calc_interpolated_elements.h"
53	#include "elliptic_to_rectangular.h"
54	#include "StelUtils.hpp"
55
56	#include <array>
57	#include <cmath>
58	#include <execution>
59	#include <iterator>
60	#include <numeric>
61	#include <cstring> /* memcpy */
62
63	#ifndef M_PI
64	#define M_PI 3.14159265358979323846
65	#endif
66
67
68	struct MarsSatTerm {
69	double phase,frequency,amplitude;
70	};
71
72	struct MarsSatTermList {
73	const struct MarsSatTerm *terms;
74	int size;
75	};
76
77	struct MarsSatBody {
78	const char *name;
79	double mu,l,acc;
80	double constants[6];
81	const struct MarsSatTermList lists[4];
82	//const std::vector<std::array<struct MarsSatTerm>
83	};
84
85	static const std::array<struct MarsSatTerm, 16> mars_sat_phobos_0 = {{
86	{ 5.013490350126586e+00, 2.715567382195733e+01, 4.537539999999999e-09},
87	{ 6.839910590780000e-01, 1.969446151585829e+01, 7.312000000000000e-10},
88	{ 4.514031245592586e+00, 4.073350897245015e+01, 7.785599999999993e-10},
89	{ 5.531351678889586e+00, 1.357783661756435e+01, 3.529399999999985e-10},
90	{ 5.700307292822586e+00, 4.685013393001369e+01, 2.094599999999994e-10},
91	{ 4.228062220664587e+00, 5.431134764391466e+01, 9.139999999999999e-11},
92	{ 1.167907583017500e+00, 7.461211875946264e+00, 6.204000000000000e-11},
93	{ 5.197274703601086e+00, 6.042797388545593e+01, 5.295999999999951e-11},
94	{ 5.160045430868086e+00, 3.938582310889583e+01, 1.811999999999999e-11},
95	{ 2.223277229900500e+00, 5.777452001970681e+01, 1.563999999999897e-11},
96	{ 6.921616240614999e-01, 2.103904892768837e+01, 1.535999999999999e-11},
97	{ 6.215802752555586e+00, 3.327229783044108e+01, 1.304000000000000e-11},
98	{ 1.368231494011000e+00, 3.938892179708230e+01, 1.661999999999999e-11},
99	{ 2.601539199675000e-01, 7.446011234647927e+00, 1.328000000000000e-11},
100	{ 5.903845687100000e-02, 3.941932465627423e+01, 8.559999999999978e-12},
101	{ 6.179310844209586e+00, 5.911910805492755e+01, 7.240000000000000e-12},
102	}};
103
104	static const std::array<struct MarsSatTerm, 42> mars_sat_phobos_1 = {{
105	{ 5.334169568271586e+00, 9.145914066599032e-03, 5.016848130000000e-05},
106	{ 5.334169148295587e+00, 9.145914113327495e-03, 5.016848130000000e-05},
107	{ 3.805345265960000e-01, 1.540952541904328e-03, 3.350264970000000e-05},
108	{ 3.805339484960000e-01, 1.540952613104260e-03, 3.350264903000000e-05},
109	{ 5.288363718535586e+00, 1.829243888046791e-02, 2.827839093000000e-05},
110	{ 5.288363809359586e+00, 1.829243887034939e-02, 2.827839093000000e-05},
111	{ 3.412685718051586e+00, 7.604696023359408e-03, 1.578862428000000e-05},
112	{ 3.412685998796586e+00, 7.604695992157328e-03, 1.578862427000000e-05},
113	{ 3.448999855869586e+00, 2.715567616927210e+01, 3.008844320000000e-05},
114	{ 9.194445444030001e-01, 3.101144183892481e-03, 8.178563170000001e-06},
115	{ 9.194443192780000e-01, 3.101144208596499e-03, 8.178563340000000e-06},
116	{ 5.387625479307586e+00, 1.969446275049270e+01, 2.044327843999999e-05},
117	{ 2.848996739890000e+00, 2.743829098982358e-02, 6.516951720000000e-06},
118	{ 2.848996748767000e+00, 2.743829098862537e-02, 6.516951720000000e-06},
119	{ 3.954411167839586e+00, 1.357783808463605e+01, 1.013316988000000e-05},
120	{ 2.912659283736000e+00, 4.073351425390815e+01, 7.483468759999740e-06},
121	{ 4.067858164824586e+00, 2.589727369338346e-02, 2.589359390000000e-06},
122	{ 4.067857973471586e+00, 2.589727371473152e-02, 2.589359390000000e-06},
123	{ 3.776271613850000e-01, 3.658417828853650e-02, 1.125147420000000e-06},
124	{ 3.776267985560000e-01, 3.658417832900902e-02, 1.125147420000000e-06},
125	{ 4.121936540919586e+00, 4.685014091671000e+01, 1.523532599999976e-06},
126	{ 1.231836210236000e+00, 1.674957754399965e-02, 4.976512300000000e-07},
127	{ 1.231836585135000e+00, 1.674957750247454e-02, 4.976512300000000e-07},
128	{ 2.488083095032000e+00, 1.520386207875350e-02, 6.841966800000000e-07},
129	{ 2.488083105051000e+00, 1.520386207728190e-02, 6.841966800000000e-07},
130	{ 1.576141880659500e+00, 3.504321962396548e-02, 1.124868580000000e-06},
131	{ 6.145328698142587e+00, 6.116624841770791e+00, 7.873411399999999e-07},
132	{ 2.626756105114500e+00, 5.431135233854420e+01, 1.040645720000000e-06},
133	{ 2.759579551372000e+00, 7.461213382164448e+00, 7.166789399999994e-07},
134	{ 6.005673273380586e+00, 1.314187651672942e+00, 6.022022799999936e-07},
135	{ 9.353090056600000e-02, 7.532115638546220e-04, 1.964976900000000e-07},
136	{ 9.343265452500001e-02, 7.532237916966188e-04, 1.964964500000000e-07},
137	{ 4.302928383505000e-01, 3.938582557777600e+01, 4.987307800000000e-07},
138	{ 3.587439959397586e+00, 6.042797969706294e+01, 4.804827999999994e-07},
139	{ 4.183547650369587e+00, 4.572998195120822e-02, 1.747133100000000e-07},
140	{ 4.183547815926586e+00, 4.572998193123842e-02, 1.747133100000000e-07},
141	{ 5.589444380180586e+00, 6.068411018120026e-03, 1.793076700000000e-07},
142	{ 5.589445119118587e+00, 6.068410936124942e-03, 1.793076800000000e-07},
143	{ 6.062696336024086e+00, 3.938892673561993e+01, 2.653755199999908e-07},
144	{ 4.628882088425586e+00, 3.327230135427769e+01, 2.539272800000000e-07},
145	{ 4.736462571780586e+00, 1.224558395734566e-02, 1.254715400000000e-07},
146	{ 4.736462588777586e+00, 1.224558396211886e-02, 1.254715500000000e-07},
147	}};
148
149	static const std::array<struct MarsSatTerm, 27> mars_sat_phobos_2 = {{
150	{ 1.404382124885000e+00, 7.595588511174286e-03, 1.514110912521000e-02},
151	{ 2.088385523034000e+00, 1.970205682773437e+01, 3.849620867400000e-04},
152	{ 3.895377361148586e+00, 1.069670280268190e-02, 6.903413242000000e-05},
153	{ 4.989572094750000e-01,-7.605272825959576e-03, 4.946101994000000e-05},
154	{ 8.215210221720000e-01, 4.685772969905357e+01, 3.671320788000000e-05},
155	{ 3.378076387479586e+00,-7.453616316894914e+00, 3.267983782000000e-05},
156	{ 2.772375719698000e+00, 3.939651873194720e+01, 8.750483050000000e-06},
157	{ 2.839303414769000e+00, 6.124220388371520e+00, 6.170938810000000e-06},
158	{ 3.184503170430000e-01, 6.043556479692822e+01, 6.957005940000000e-06},
159	{ 1.410329121503000e+00, 1.984322681200531e-02, 5.734088960000000e-06},
160	{ 1.336796551701000e+00, 3.327989343110183e+01, 3.431577860000000e-06},
161	{ 3.853646745987586e+00,-2.103145308554279e+01, 4.037409510000000e-06},
162	{ 2.547067366210000e-01, 1.549665594221090e-03, 2.115296220000000e-06},
163	{ 1.508224230072000e+00,-5.911151611514718e+01, 1.904827620000000e-06},
164	{ 3.462607551400000e-02,-5.165030141480156e+01, 8.696094700000001e-07},
165	{ 3.198655314207586e+00, 9.152774292423818e-03, 6.329868200000000e-07},
166	{ 5.163954670916587e+00, 1.674080796065678e-02, 6.130317800000000e-07},
167	{ 7.803937436990001e-01,-1.550610622055567e-03, 5.737217100000000e-07},
168	{ 3.264588234538586e+00, 2.716326950857758e+01, 8.181088300000000e-07},
169	{ 5.211329429707586e+00, 2.898903518276566e-02, 7.434145300000000e-07},
170	{ 2.691716763490000e+00,-2.714807769012140e+01, 5.441647800000001e-07},
171	{ 1.001733960626000e+00,-4.553367734785991e+01, 5.044515000000000e-07},
172	{ 3.212149602842586e+00,-1.968376490094933e+01, 4.284599100000000e-07},
173	{ 4.145658249560586e+00,-3.460928932187340e+01, 4.224937300000000e-07},
174	{ 5.063879137454586e+00, 2.589097180762378e-02, 3.581949800000000e-07},
175	{ 2.558639707373000e+00, 6.064928432880783e-03, 3.808881200000000e-07},
176	{ 5.292499025555586e+00, 3.075065445597794e-03, 4.055739200000000e-07},
177	}};
178
179	static const std::array<struct MarsSatTerm, 28> mars_sat_phobos_3 = {{
180	{ 2.058107128488000e+00,-7.604861328578004e-03, 9.408605183120001e-03},
181	{ 3.248856489376586e+00,-9.145863467943084e-03, 5.699538102000000e-05},
182	{ 4.557280987272586e+00, 1.829238959116374e-02, 2.474369483000000e-05},
183	{ 6.113423353213586e+00, 7.595702278066188e-03, 2.062715397000000e-05},
184	{ 2.067112903892000e+00, 2.743825478464904e-02, 5.945574160000000e-06},
185	{ 1.700689207004000e+00, 9.145729527513899e-03, 3.920190200000000e-06},
186	{ 3.577059353926586e+00,-1.829259493439447e-02, 2.167791970000000e-06},
187	{ 4.940284891770586e+00,-7.453616316894914e+00, 2.769763750000000e-06},
188	{ 2.117885747855000e+00, 3.941171892425920e+01, 1.813539170000000e-06},
189	{ 5.176055365230000e-01, 1.970205682773437e+01, 1.083429180000000e-06},
190	{ 5.864248782211586e+00, 3.658425999636452e-02, 1.066657010000000e-06},
191	{ 3.240692168340586e+00, 2.589639535977297e-02, 7.909853300000000e-07},
192	{ 7.570751994340000e-01,-1.321790869500707e+00, 8.880885099999999e-07},
193	{ 6.608910796060000e-01, 6.124220388371520e+00, 5.520758900000000e-07},
194	{ 5.519669134088586e+00, 4.685772969905357e+01, 4.397886000000000e-07},
195	{ 5.844000094873587e+00,-2.743825733444110e-02, 4.583609100000000e-07},
196	{ 4.602755395228586e+00,-1.675516082048598e-02, 3.533131100000000e-07},
197	{ 2.685649467156000e+00, 1.542518910643103e-03, 3.502988500000000e-07},
198	{ 3.353176877809586e+00, 1.225604581839110e+01, 3.828851700000000e-07},
199	{ 2.556761413073000e+00, 1.068971473943851e-02, 2.711247300000000e-07},
200	{ 4.622369718415587e+00,-2.589682634097252e-02, 1.991745400000000e-07},
201	{ 7.497690894330000e-01, 3.504362063026418e-02, 2.036897700000000e-07},
202	{ 3.965205602751586e+00, 2.714806830611846e+01, 1.726779600000000e-07},
203	{ 4.743858488241586e+00,-2.103145308554279e+01, 1.574472800000000e-07},
204	{ 3.405870070134586e+00, 4.572380100924578e-02, 1.698567700000000e-07},
205	{ 3.742380849210000e-01, 3.327989343110183e+01, 1.015004500000000e-07},
206	{ 1.374678754422000e+00,-1.970206671243377e+01, 8.259149000000001e-08},
207	{ 1.140160118874000e+00,-6.108652383295211e-03, 8.504923000000000e-08},
208	}};
209
210	static const std::array<struct MarsSatTerm, 25> mars_sat_deimos_0 = {{
211	{ 6.146803530569087e+00, 2.294413036846356e+00, 5.398360000000000e-09},
212	{ 3.432805186019586e+00, 9.935735425667586e+00, 7.095199999999997e-10},
213	{ 4.356581106602587e+00, 3.441619555269535e+00, 3.737800000000000e-10},
214	{ 5.915601483397086e+00, 9.926589650598594e+00, 2.337400000000000e-10},
215	{ 1.615243161446500e+00, 1.147206518423178e+00, 1.697000000000000e-10},
216	{ 2.114704906116000e+00, 9.917443791162608e+00, 5.314000000000000e-11},
217	{ 9.393388277920001e-01, 9.954027953247621e+00, 7.046000000000000e-11},
218	{ 3.949494792732586e+00, 9.944881541450725e+00, 4.233999999999998e-11},
219	{ 2.089322391969500e+00, 2.294097606691310e+00, 2.813999999999983e-11},
220	{ 7.795571756865000e-01, 2.294728469946876e+00, 2.549999999999971e-11},
221	{ 2.361851669664500e+00, 4.588826073692712e+00, 1.676000000000000e-11},
222	{ 2.770026883578500e+00, 9.954658773106329e+00, 9.719999999999976e-12},
223	{ 4.733522308629587e+00, 9.963173826265681e+00, 9.519999999999998e-12},
224	{ 3.258781379519086e+00, 1.472506551307767e+01, 9.739999999999990e-12},
225	{ 4.596299836113586e+00, 9.908297894321850e+00, 1.027999999999982e-11},
226	{ 3.057688443346500e+00, 2.682916339604879e+00, 6.679999999999999e-12},
227	{ 5.453932177963086e+00, 2.682285947505377e+00, 6.879999999999995e-12},
228	{ 4.328499239462586e+00, 9.936050924435195e+00, 6.199999999999996e-12},
229	{ 2.168321380562000e+00, 4.976699000854042e+00, 5.359999999999916e-12},
230	{ 2.516818588675000e+00, 9.935420104041299e+00, 3.359999999999995e-12},
231	{ 2.997827070725000e-01, 3.441304123224338e+00, 2.919999999999998e-12},
232	{ 2.143143420462500e+00, 2.682600911813873e+00, 2.979999999999999e-12},
233	{ 5.271781540988586e+00, 3.441934988580561e+00, 2.640000000000000e-12},
234	{ 5.382287205930000e-01, 9.926905039305762e+00, 2.039999999999966e-12},
235	{ 3.629187963751586e+00, 2.303558945733479e+00, 1.090000000000000e-12},
236	}};
237
238	static const std::array<struct MarsSatTerm, 21> mars_sat_deimos_1 = {{
239	{ 2.488175621789000e+00, 3.153377646687549e-04, 2.483031660190000e-03},
240	{ 2.488175634463000e+00, 3.153377641537758e-04, 2.483031659490000e-03},
241	{ 5.336889264964586e+00, 9.145915158660612e-03, 2.023426978200000e-04},
242	{ 5.336889346905586e+00, 9.145915155199243e-03, 2.023426978200000e-04},
243	{ 5.281538321363586e+00, 1.829244561727392e-02, 1.043167564200000e-04},
244	{ 5.281537958522586e+00, 1.829244563271799e-02, 1.043167564300000e-04},
245	{ 1.433490018835500e+00, 2.294413045224799e+00, 1.637133774599998e-04},
246	{ 2.852007788337500e+00, 2.743818628432718e-02, 4.820671420000000e-05},
247	{ 3.147176604035586e+00, 1.860783561096842e-02, 1.497066023000000e-05},
248	{ 3.147176307034586e+00, 1.860783562364318e-02, 1.497066034000000e-05},
249	{ 3.182516662993086e+00, 1.147206709032795e+00, 1.137028858000000e-05},
250	{ 5.917181875408586e+00, 3.441620120814550e+00, 8.350609860000001e-06},
251	{ 3.544871659744586e+00, 6.190897763118658e-04, 4.383254500000000e-06},
252	{ 3.544871649035586e+00, 6.190897769462623e-04, 4.383254500000000e-06},
253	{ 4.614193383590000e-01, 3.657933964533806e-02, 8.316176420000001e-06},
254	{ 6.754222685930000e-01, 2.775277029221846e-02, 6.447587580000001e-06},
255	{ 5.001640850015586e+00, 9.935735582783501e+00, 7.926687400000001e-06},
256	{ 4.485218623803586e+00, 8.828439970895009e-03, 3.905064210000000e-06},
257	{ 2.221013917980000e-01, 9.464721723028089e-03, 3.117801830000000e-06},
258	{ 4.485218622895586e+00, 8.828439971019056e-03, 3.905064210000000e-06},
259	{ 2.221013944660000e-01, 9.464721722885516e-03, 3.117801830000000e-06},
260	}};
261
262	static const std::array<struct MarsSatTerm, 15> mars_sat_deimos_2 = {{
263	{ 2.198649419514000e+00, 3.148401892560942e-04, 2.744131534600000e-04},
264	{ 4.366636518977586e+00, 4.977013897776327e+00, 6.015912711000000e-05},
265	{ 1.463816210370000e+00,-3.153164045300449e-04, 3.614585349000000e-05},
266	{ 1.362467588064000e+00, 2.682600831640474e+00, 2.577157903000000e-05},
267	{ 9.336792893320000e-01,-4.958721582505435e+00, 6.801997830000000e-06},
268	{ 3.152432010349586e+00, 1.535394310272536e+00, 4.452979140000000e-06},
269	{ 4.734051493918586e+00,-4.949575691510945e+00, 2.238482920000000e-06},
270	{ 5.736642853342587e+00, 2.327524649334390e-04, 1.478747570000000e-06},
271	{ 5.084945958505586e+00, 3.912335403634271e-04, 1.457720040000000e-06},
272	{ 4.232770650771586e+00, 7.271426904652134e+00, 1.347030310000000e-06},
273	{ 1.516439412109000e+00, 1.491274947692300e+01, 7.538179100000001e-07},
274	{ 3.100618466711000e+00, 1.797748798526167e-02, 8.716229700000000e-07},
275	{ 5.155432690964586e+00, 3.881877876517008e-01, 1.015648880000000e-06},
276	{ 2.250768124831000e+00,-4.940429834739583e+00, 5.090989300000000e-07},
277	{ 3.422251338868586e+00,-4.977013897776327e+00, 6.537077500000000e-07},
278	}};
279
280	static const std::array<struct MarsSatTerm, 27> mars_sat_deimos_3 = {{
281	{ 2.981506933511000e+00,-3.154811355556041e-04, 1.562693319959000e-02},
282	{ 4.557500894366586e+00, 1.829233626168517e-02, 1.321818631100000e-04},
283	{ 3.248124065112586e+00,-9.145934570587952e-03, 3.833652719000000e-05},
284	{ 1.701551338710000e+00, 9.145665986601770e-03, 2.660211842000000e-05},
285	{ 2.067662003165000e+00, 2.743821959085702e-02, 2.882438535000000e-05},
286	{ 4.805327222478586e+00, 3.158564952832469e-04, 2.713213911000000e-05},
287	{ 2.318788380092000e+00, 1.860776219134252e-02, 9.296324950000000e-06},
288	{ 3.613580709014586e+00,-1.829261290520875e-02, 4.460960290000000e-06},
289	{ 5.867810550672586e+00, 3.658405181871115e-02, 4.904156030000000e-06},
290	{ 5.496515547817586e+00,-9.461294099629817e-03, 2.366881150000000e-06},
291	{ 6.120362957845586e+00, 2.775361064625920e-02, 2.057552850000000e-06},
292	{ 3.612389208841586e+00, 8.830178627407989e-03, 2.365750140000000e-06},
293	{ 5.554163001934587e+00,-1.860800990265229e-02, 1.228814440000000e-06},
294	{ 3.561485983136586e+00, 1.797648948023003e-02, 1.215189460000000e-06},
295	{ 1.941975361869000e+00,-6.259843043363501e-04, 9.602856999999999e-07},
296	{ 5.845244776759587e+00,-2.743819348074506e-02, 1.174471770000000e-06},
297	{ 1.821048365911000e+00, 9.460829229838182e-03, 1.058441050000000e-06},
298	{ 7.587135454610000e-01, 1.682177864988152e-04, 2.383147490000000e-06},
299	{ 3.292891638604586e+00, 4.573365667404765e-02, 7.656303800000000e-07},
300	{ 3.766016108380586e+00, 3.683140215171970e-02, 4.075824300000000e-07},
301	{ 5.747244956219586e+00, 9.954342722363187e+00, 4.962340900000000e-07},
302	{ 1.724062868741000e+00,-2.774017744689255e-02, 3.003077900000000e-07},
303	{ 4.429747729642586e+00,-8.904754331211824e-03, 2.613416400000000e-07},
304	{ 8.659161909010000e-01, 2.718340176462802e-02, 2.891945200000000e-07},
305	{ 2.259603352363000e+00,-3.660152378331885e-02, 2.213788200000000e-07},
306	{ 1.108725025984000e+00, 9.935732440466172e+00, 2.449950300000000e-07},
307	{ 1.680391963791000e+00, 9.076040299352415e-03, 1.652416600000000e-07},
308	}};
309
310	static const struct MarsSatBody mars_sat_bodies[2] = {
311	{
312	"phobos", 9.549547741038312e-11, 1.970205562831390e+01, 1.657852042683113e-10,
313	{
314	0.0000626916188000,
315	2.0912973926417302,
316	-0.0000005774003141,
317	0.0000006185052790,
318	-0.0000575018919826,
319	-0.0000540676493351,
320	},
321	{
322	{mars_sat_phobos_0.data(),16},
323	{mars_sat_phobos_1.data(),42},
324	{mars_sat_phobos_2.data(),27},
325	{mars_sat_phobos_3.data(),28}
326	}
327	},
328	{
329	"deimos", 9.549547622768120e-11, 4.977013889652300e+00, -2.331793571572226e-14,
330	{
331	0.0001568134086700,
332	-1.9167537810740201,
333	-0.0000239579419518,
334	0.0000257893300229,
335	-0.0056443721379635,
336	-0.0053121806978560,
337	},
338	{
339	{mars_sat_deimos_0.data(),25},
340	{mars_sat_deimos_1.data(),21},
341	{mars_sat_deimos_2.data(),15},
342	{mars_sat_deimos_3.data(),27}
343	}
344	},
345	};
346
347	static void CalcMarsSatElem(double t,int body,double elem[6]) {	28✔
348	int j;
349	const struct MarsSatBody *bp = mars_sat_bodies + body;	28✔
350	memcpy(elem, bp->constants, 6*sizeof(double));	28✔
351
352	// for (j=0;j<2;j++) {
353	// const struct MarsSatTerm *const begin = bp->lists[j].terms;
354	// const struct MarsSatTerm *p = begin + bp->lists[j].size;
355	// while (--p >= begin) {
356	// const double d = p->phase + t*p->frequency;
357	// elem[j] += p->amplitude * cos(d);
358	// }
359	// }
360
361	// TODO: The summing should be backwards to avoid rounding issues!
362	for (j=0;j<2;j++) {	84✔
363	elem[j] += std::transform_reduce(	56✔
364	STD_EXECUTION_PAR_COMMA
365	bp->lists[j].terms, bp->lists[j].terms + bp->lists[j].size, 0.0,	56✔
366	std::plus<>(),
367	[=](const struct MarsSatTerm &trm){	1,456✔
368	const double d = trm.phase + t*trm.frequency;	1,456✔
369	return trm.amplitude * cos(d);	1,456✔
370	});
371	}
372
373
374	// for (j=2;j<4;j++) {
375	// const struct MarsSatTerm *const begin = bp->lists[j].terms;
376	// const struct MarsSatTerm *p = begin + bp->lists[j].size;
377	// while (--p >= begin) {
378	// const double d = p->phase + t*p->frequency;
379	// elem[2j-2] += p->amplitude cos(d);
380	// elem[2j-1] += p->amplitude sin(d);
381	// }
382	// }
383
384	// TODO: The summing should be backwards to avoid rounding issues!
385	for (j=2;j<4;j++) {	84✔
386	std::pair<double, double>el3456_add=
NEW 387	std::transform_reduce(	×
388	STD_EXECUTION_PAR_COMMA
389	bp->lists[j].terms, bp->lists[j].terms+bp->lists[j].size,	56✔
390	std::pair<double, double>({0.0, 0.0}),	56✔
391	[](const std::pair<double, double>&sum, const std::pair<double, double>&addon){	1,358✔
392	return std::make_pair(sum.first+addon.first, sum.second+addon.second);	1,358✔
393	},
394	[=](const struct MarsSatTerm &trm){	1,358✔
395	const double d = trm.phase + t*trm.frequency;	1,358✔
396	return std::make_pair(trm.amplitude * cos(d), trm.amplitude * sin(d));	1,358✔
397	});
398	elem[2*j-2] += el3456_add.first;	56✔
399	elem[2*j-1] += el3456_add.second;	56✔
400	}
401
402	elem[1] += (bp->l + bp->acc * t) * t;	28✔
403	}	28✔
404
405
406	static void MultMat(const double a[9],const double b[9],double c[9]) {	14✔
407	int i,j;
408	for (i=0;i<3;i++) {	56✔
409	for (j=0;j<3;j++) {	168✔
410	c[i3+j] = a[i3]b[j] + a[i3+1]b[3+j] + a[i3+2]*b[6+j];	126✔
411	}
412	}
413	}	14✔
414
415	static const double J2000_to_VSOP87[9] = {
416	9.999999999998848e-01,-4.799655442984222e-07, 0.000000000000000e+00,
417	4.403598133110236e-07, 9.174821370868568e-01, 3.977769829016507e-01,
418	-1.909192461077750e-07,-3.977769829016049e-01, 9.174821370869625e-01};
419
420	static const double ome0 = 47.68143;
421	static const double inc0 = 37.1135;
422	static const double dome = -0.1061;
423	static const double dinc = 0.0609;
424
425	static void GenerateMarsSatToVSOP87(double t,double mat_mars_sat_to_vsop87[9]) {	14✔
426	t -= 6491.5;	14✔
427	{
428	const double ome = (ome0 + dome * t / 36525.) * (M_PI/180.0);	14✔
429	const double inc = (inc0 + dinc * t / 36525.) * (M_PI/180.0);	14✔
430	const double co = cos(ome);	14✔
431	const double so = sin(ome);	14✔
432	const double ci = cos(inc);	14✔
433	const double si = sin(inc);	14✔
434	// GZ: I hope the next construct is meanwhile understood on all platforms
435	//#if defined __GNUC__ && !defined __STRICT_ANSI__
436	const double m[9] = {co,-ciso, siso,	14✔
437	so, cico,-sico,	14✔
438	0.0,si,ci};	14✔
439	//#else
440	// double m[9];
441	// m[0] = co; m[1] = -ciso; m[2] = siso;
442	// m[3] = so; m[4] = cico; m[5] = -sico;
443	// m[6] = 0.0; m[7] = si; m[8] = ci;
444	//#endif
445	MultMat(J2000_to_VSOP87,m,mat_mars_sat_to_vsop87);	14✔
446	}
447	}	14✔
448
449	static double t_0 = -1e100;
450	static double t_1 = -1e100;
451	static double t_2 = -1e100;
452	static double marssat_elem_0[2*6];
453	static double marssat_elem_1[2*6];
454	static double marssat_elem_2[2*6];
455
456	/* 1 day: */
457	#define DELTA_T 1.0
458
459	static double marssat_jd0 = -1e100;
460	static double marssat_elem[2*6];
461
462	static void CalcAllMarsSatElem(double t,double elem[12], void *user) {	14✔
463	CalcMarsSatElem(t,0,elem+(0*6));	14✔
464	CalcMarsSatElem(t,1,elem+(1*6));	14✔
465	}	14✔
466
467	static double mars_sat_to_vsop87[9];
468
469	void GetMarsSatCoor(double jd,int body,double xyz, double xyzdot) {	14✔
470	double xyz6[6];
471	GetMarsSatOsculatingCoor(jd,jd,body,xyz6);	14✔
472	xyz[0] =xyz6[0]; xyz[1] =xyz6[1]; xyz[2] =xyz6[2];	14✔
473	xyzdot[0]=xyz6[3]; xyzdot[1]=xyz6[4]; xyzdot[2]=xyz6[5];	14✔
474	}	14✔
475
476	void GetMarsSatOsculatingCoor(const double jd0,const double jd,	14✔
477	const int body,double *xyz) {
478	double x[6];
479	if (jd0 != marssat_jd0) {	14✔
480	const double t0 = jd0 - 2451545.0 + 6491.5;	14✔
481	marssat_jd0 = jd0;	14✔
482	CalcInterpolatedElements(t0,marssat_elem,12,	14✔
483	&CalcAllMarsSatElem,DELTA_T,
484	&t_0,marssat_elem_0,
485	&t_1,marssat_elem_1,
486	&t_2,marssat_elem_2,
487	NULL);
488	GenerateMarsSatToVSOP87(t0,mars_sat_to_vsop87);	14✔
489	}
490	EllipticToRectangularA(mars_sat_bodies[body].mu,marssat_elem+(body*6),jd-jd0,x);	14✔
491	xyz[0] = mars_sat_to_vsop87[0]*x[0]	14✔
492	+ mars_sat_to_vsop87[1]*x[1]	14✔
493	+ mars_sat_to_vsop87[2]*x[2];	14✔
494	xyz[1] = mars_sat_to_vsop87[3]*x[0]	14✔
495	+ mars_sat_to_vsop87[4]*x[1]	14✔
496	+ mars_sat_to_vsop87[5]*x[2];	14✔
497	xyz[2] = mars_sat_to_vsop87[6]*x[0]	14✔
498	+ mars_sat_to_vsop87[7]*x[1]	14✔
499	+ mars_sat_to_vsop87[8]*x[2];	14✔
500	// GZ This is a guess, based on the structure of other operations...
501	xyz[3] = mars_sat_to_vsop87[0]*x[3]	14✔
502	+ mars_sat_to_vsop87[1]*x[4]	14✔
503	+ mars_sat_to_vsop87[2]*x[5];	14✔
504	xyz[4] = mars_sat_to_vsop87[3]*x[3]	14✔
505	+ mars_sat_to_vsop87[4]*x[4]	14✔
506	+ mars_sat_to_vsop87[5]*x[5];	14✔
507	xyz[5] = mars_sat_to_vsop87[6]*x[3]	14✔
508	+ mars_sat_to_vsop87[7]*x[4]	14✔
509	+ mars_sat_to_vsop87[8]*x[5];	14✔
510	/*
511	printf("%d %18.9lf %15.12lf %15.12lf %15.12lf\n",
512	body,jd,xyz[0],xyz[1],xyz[2]);
513	*/
514	}	14✔

Stellarium / stellarium / 10124913218

Source File Press 'n' to go to next uncovered line, 'b' for previous

Source File
Press 'n' to go to next uncovered line, 'b' for previous