1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
|
#/*-----------------------------------------------------------------*/
#/*!
# \file mcheck.py
# \brief Check the output of the function calceph_compute with an INPOP file.
#
# \author M. Gastineau
# Astronomie et Systemes Dynamiques, IMCCE, CNRS, Observatoire de Paris.
#
# Copyright, 2016, 2017, CNRS
# email of the author : Mickael.Gastineau@obspm.fr
#*/
#/*-----------------------------------------------------------------*/
#
#/*-----------------------------------------------------------------*/
#/* License of this file :
# This file is "triple-licensed", you have to choose one of the three licenses
# below to apply on this file.
#
# CeCILL-C
# The CeCILL-C license is close to the GNU LGPL.
# ( http://www.cecill.info/licences/Licence_CeCILL-C_V1-en.html )
#
# or CeCILL-B
# The CeCILL-B license is close to the BSD.
# (http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.txt)
#
# or CeCILL v2.1
# The CeCILL license is compatible with the GNU GPL.
# ( http://www.cecill.info/licences/Licence_CeCILL_V2.1-en.html )
#
#
# This library is governed by the CeCILL-C, CeCILL-B or the CeCILL license under
# French law and abiding by the rules of distribution of free software.
# You can use, modify and/ or redistribute the software under the terms
# of the CeCILL-C,CeCILL-B or CeCILL license as circulated by CEA, CNRS and INRIA
# at the following URL "http://www.cecill.info".
#
# As a counterpart to the access to the source code and rights to copy,
# modify and redistribute granted by the license, users are provided only
# with a limited warranty and the software's author, the holder of the
# economic rights, and the successive licensors have only limited
# liability.
#
# In this respect, the user's attention is drawn to the risks associated
# with loading, using, modifying and/or developing or reproducing the
# software by the user in light of its specific status of free software,
# that may mean that it is complicated to manipulate, and that also
# therefore means that it is reserved for developers and experienced
# professionals having in-depth computer knowledge. Users are therefore
# encouraged to load and test the software's suitability as regards their
# requirements in conditions enabling the security of their systems and/or
# data to be ensured and, more generally, to use and operate it in the
# same conditions as regards security.
#
# The fact that you are presently reading this means that you have had
# knowledge of the CeCILL-C,CeCILL-B or CeCILL license and that you accept its terms.
# */
# /*-----------------------------------------------------------------*/
#/*-----------------------------------------------------------------*/
#/* main program */
#/*-----------------------------------------------------------------*/
import unittest
import openfiles
from calcephpy import CalcephBin
class TestOpen(unittest.TestCase):
def mcheck_main(self, prefetch):
peph = CalcephBin.open(openfiles.prefixsrc("../../examples/example1.dat"))
if (prefetch==0):
peph.prefetch()
f = open(openfiles.prefixsrc('../../tests/example1_tests.dat'))
d = f.readlines()
datacheck = [x.split() for x in d]
for aline in datacheck:
jd0 = float(aline[0])
target = int(aline[1])
center = int(aline[2])
PVcheck = [0,0,0,0,0,0]
for k in range(0,6):
PVcheck[k] = float(aline[k+3])
dt = jd0 - int(jd0)
jd0 = (int(jd0)) + 2.4515450000000000000E+06
PV = peph.compute(jd0, dt, target, center)
seuil = 1E-8
# /* for libration , return angles between 0 and 2*pi */
if (target == 15):
seuil = 1E-7;
j = 2
pi2 = 6.283185307179586E0
while (PV[j] >= pi2):
PV[j] = PV[j] - pi2
while (PV[j] <= 0) :
PV[j] = PV[j] + pi2;
# /* check with references coordinates */
n = 6
for j in range(0,6):
if (abs(PVcheck[j] - PV[j]) >= seuil):
print("failure for ",target, center, " at time ", jd0, dt)
print('aline=', aline)
for k in range(0,n):
print(k, "diff= ",target, center, " at time jd0, dt", PVcheck[k] - PV[k], PVcheck[k], PV[k])
print("prefetch=", prefetch)
raise RuntimeError("Internal failure")
f.close()
peph.close()
def test_mcheck(self):
self.mcheck_main(0)
self.mcheck_main(1)
if __name__ == '__main__':
unittest.main()
|
