equationoftime 4.1.5 → 4.1.6

Sign up to get free protection for your applications and to get access to all the features.
@@ -1,11 +1,11 @@
1
1
  # class Eot file = trigometric.rb
2
+ # methods using trig functions
2
3
  class Eot
3
- include Math
4
4
  # From trigometric.rb:
5
5
  # cosine apparent longitude
6
6
  # could be useful when dividing
7
7
  def cosine_al_sun
8
- cos_al_sun(al(@ma, @ta, Celes.faom03(@ta)))
8
+ cos_al_sun(al_sun)
9
9
  end
10
10
  alias_method :cosine_apparent_longitude, :cosine_al_sun
11
11
  alias_method :cosalsun, :cosine_al_sun
@@ -14,7 +14,7 @@ class Eot
14
14
  # cosine true longitude
15
15
  # used in solar right ascension
16
16
  def cosine_tl_sun
17
- cos_tl_sun(tl(@ma, @ta))
17
+ cos_tl_sun(tl_sun)
18
18
  end
19
19
  alias_method :cosine_true_longitude, :cosine_tl_sun
20
20
 
@@ -22,7 +22,7 @@ class Eot
22
22
  # cosine true obliquity
23
23
  # used in solar right ascension and equation of equinox
24
24
  def cosine_to_earth
25
- cos_to_earth(Celes.nut06a(@ajd, 0)[1] + Celes.obl06(@ajd, 0))
25
+ cos_to_earth(to_earth)
26
26
  end
27
27
  alias_method :cosine_true_obliquity, :cosine_to_earth
28
28
 
@@ -30,7 +30,7 @@ class Eot
30
30
  # sine apparent longitude
31
31
  # used in solar declination
32
32
  def sine_al_sun
33
- sin_al_sun(al(@ma, @ta, Celes.faom03(@ta)))
33
+ sin_al_sun(al_sun)
34
34
  end
35
35
  alias_method :sine_apparent_longitude, :sine_al_sun
36
36
 
@@ -38,7 +38,7 @@ class Eot
38
38
  # sine true longitude
39
39
  # used in solar right ascension
40
40
  def sine_tl_sun
41
- sin_tl_sun(tl(@ma, @ta))
41
+ sin_tl_sun(tl_sun)
42
42
  end
43
43
  alias_method :sine_true_longitude, :sine_tl_sun
44
44
 
@@ -46,6 +46,6 @@ class Eot
46
46
  # sine true obliquity angle of Earth
47
47
  # used in solar declination
48
48
  def sine_to_earth
49
- sin_to_earth(Celes.nut06a(@ajd, 0)[1] + Celes.obl06(@ajd, 0))
49
+ sin_to_earth(to_earth)
50
50
  end
51
51
  end
data/lib/eot/utilities.rb CHANGED
@@ -1,5 +1,5 @@
1
- # utilities.rb
2
- # class Eot
1
+ # class Eot file = utilities.rb
2
+ # methods for general use
3
3
  class Eot
4
4
  # From utilities.rb:
5
5
  # A check for default J2000
@@ -29,6 +29,19 @@ class Eot
29
29
  dt == 0 ? DT2000 : check_t_nil(dt)
30
30
  end
31
31
 
32
+ # From utilities.rb:
33
+ # formats degree string
34
+ def f_string(sgn, u, m, s, ds)
35
+ sgn +
36
+ format('%03d', u) +
37
+ ':' +
38
+ format('%02d', m) +
39
+ ':' +
40
+ format('%02d', s) +
41
+ '.' +
42
+ format('%3.3d', ds)
43
+ end
44
+
32
45
  # From utilities.rb:
33
46
  # Keeps large angles in range of 360.0
34
47
  # aliased by truncate
data/lib/eot/version.rb CHANGED
@@ -1,4 +1,6 @@
1
1
  # class Eot file = version.rb
2
+ # contains gem version
2
3
  class Eot
3
- VERSION = '4.1.5'
4
+ # current version
5
+ VERSION = '4.1.6'
4
6
  end
@@ -0,0 +1 @@
1
+ system 'bundle exec rake'
@@ -8,25 +8,26 @@ require 'eot'
8
8
 
9
9
  aliased_angles = Eot.new
10
10
 
11
- describe 'tests ajd of 2456885.0' do
11
+ describe 'tests ajd of 2456885.0 ' do
12
12
 
13
13
  before(:each) do
14
14
  aliased_angles.ajd = 2_456_885.0
15
15
  ajd = aliased_angles.ajd
16
+ aliased_angles.ma_ta_set
16
17
  # somtimes need date to check values somewhere else
17
18
  aliased_angles.date = aliased_angles.ajd_to_datetime(ajd)
18
19
  end
19
20
 
20
- it 'expected 2456885.0 for aliased_angles.ajd'do
21
+ it 'expected 2_456_885.0 for aliased_angles.ajd 'do
21
22
  assert_equal(2_456_885.0, aliased_angles.ajd)
22
23
  end
23
24
 
24
- it 'expected "2014-08-15T12:00:00+00:00" from \
25
- aliased_angles.date'.to_s do
25
+ it 'expected "2014-08-15T12:00:00+00:00" from \
26
+ aliased_angles.date.to_s ' do
26
27
  assert_equal('2014-08-15T12:00:00+00:00', aliased_angles.date.to_s)
27
28
  end
28
29
 
29
- it 'expected 3.8508003966038915 from aliased_angles.ma'do
30
+ it 'expected 3.8508003966038915 from aliased_angles.ma 'do
30
31
  assert_equal(3.8508003966038915, aliased_angles.ma)
31
32
  end
32
33
 
@@ -41,12 +42,12 @@ describe 'tests ajd of 2456885.0' do
41
42
  aliased_angles.cosine_apparent_longitude)
42
43
  end
43
44
 
44
- it 'expected -0.7943772759574919 from \
45
+ it 'expected -0.7943772759574919 from \
45
46
  aliased_angles.cosine_true_longitude()? ' do
46
47
  assert_equal(-0.7943772759574919, aliased_angles.cosine_true_longitude)
47
48
  end
48
49
 
49
- it 'expected 0.9175115346811911 from \
50
+ it 'expected 0.9175115346811911 from \
50
51
  aliased_angles.cosine_true_obliquity()? ' do
51
52
  assert_equal(0.9175115346811911, aliased_angles.cosine_true_obliquity)
52
53
  end
@@ -72,7 +73,7 @@ describe 'tests ajd of 2456885.0' do
72
73
  aliased_angles.eccentricity_earth_orbit)
73
74
  end
74
75
 
75
- it 'expected -0.021413249720702462 from \
76
+ it 'expected -0.021413249720702462 from \
76
77
  aliased_angles.equation_of_center()? ' do
77
78
  assert_equal(-0.021413249720702462, aliased_angles.equation_of_center)
78
79
  end
@@ -83,9 +84,9 @@ describe 'tests ajd of 2456885.0' do
83
84
  aliased_angles.geometric_mean_longitude)
84
85
  end
85
86
 
86
- it 'expected 1.9143229567859146 from \
87
- aliased_angles.horizon_angle()? ' do
88
- assert_equal(1.9143229567859146, aliased_angles.horizon_angle)
87
+ it 'expected 1.5857841877939605 from \
88
+ aliased_angles.horizon_angle(1)? ' do
89
+ assert_equal(1.5857841877939605, aliased_angles.horizon_angle(1))
89
90
  end
90
91
 
91
92
  it 'expected 3.8508003966038915 from \
@@ -103,13 +104,13 @@ describe 'tests ajd of 2456885.0' do
103
104
  assert_equal(0.40905940254265843, aliased_angles.mean_obliquity)
104
105
  end
105
106
 
106
- it 'expected 0.40905940254265843 from \
107
+ it 'expected 0.40905940254265843 from \
107
108
  aliased_angles.mean_obliquity_of_ecliptic()? ' do
108
109
  assert_equal(0.40905940254265843, \
109
110
  aliased_angles.mean_obliquity_of_ecliptic)
110
111
  end
111
112
 
112
- it 'expected 0.40901870461547685 from \
113
+ it 'expected 0.40901870461547685 from \
113
114
  aliased_angles.obliquity_correction()? ' do
114
115
  assert_equal(0.40901870461547685, aliased_angles.obliquity_correction)
115
116
  end
@@ -119,12 +120,12 @@ describe 'tests ajd of 2456885.0' do
119
120
  assert_equal(2.5297411654316497, aliased_angles.right_ascension)
120
121
  end
121
122
 
122
- it 'expected 0.6074784519729512 from \
123
+ it 'expected 0.6074784519729512 from \
123
124
  aliased_angles.sine_apparent_longitude()? ' do
124
125
  assert_equal(0.6074784519729512, aliased_angles.sine_apparent_longitude)
125
126
  end
126
127
 
127
- it 'expected 0.6074246812917259 from \
128
+ it 'expected 0.6074246812917259 from \
128
129
  aliased_angles.sine_true_longitude()? ' do
129
130
  assert_equal(0.6074246812917259, aliased_angles.sine_true_longitude)
130
131
  end
@@ -139,7 +140,7 @@ describe 'tests ajd of 2456885.0' do
139
140
  assert_equal(2.48877803069344, aliased_angles.true_longitude)
140
141
  end
141
142
 
142
- it 'expected 2.5101242776531474 from \
143
+ it 'expected 2.5101242776531474 from \
143
144
  aliased_angles.true_longitude_aries()? ' do
144
145
  assert_equal(2.5101242776531474, aliased_angles.true_longitude_aries)
145
146
  end
@@ -150,30 +151,31 @@ describe 'tests ajd of 2456885.0' do
150
151
  end
151
152
  end
152
153
 
153
- describe 'tests ajd of 2455055.5' do
154
+ describe 'tests ajd of 2455055.5 ' do
154
155
 
155
156
  before(:each) do
156
157
  aliased_angles.ajd = 2_455_055.0
157
158
  ajd = aliased_angles.ajd
159
+ aliased_angles.ma_ta_set
158
160
  # check date for this ajd when needed.
159
161
  aliased_angles.date = aliased_angles.ajd_to_datetime(ajd)
160
162
  end
161
163
 
162
- it 'expected 455055.0, from aliased_angles.' do
164
+ it 'expected 2_455_055.0, from aliased_angles.' do
163
165
  assert_equal(2_455_055.0, aliased_angles.ajd)
164
166
  end
165
167
 
166
- it 'expected "2009-08-11T12:00:00+00:00" from \
167
- aliased_angles.date.to_s' do
168
+ it 'expected "2009-08-11T12:00:00+00:00" from \
169
+ aliased_angles.date.to_s ' do
168
170
  assert_equal('2009-08-11T12:00:00+00:00', aliased_angles.date.to_s)
169
171
  end
170
172
 
171
- it 'expected 3.7871218188949207, from aliased_angles.ma' do
173
+ it 'expected 3.7871218188949207, from aliased_angles.ma ' do
172
174
  assert_equal(3.7871218188949207, aliased_angles.ma)
173
175
  end
174
176
 
175
177
  it 'expected 3.7871218188949207 from \
176
- aliased_angles.ma from Eot_angles.mean_anomaly()' do
178
+ aliased_angles.ma from Eot_angles.mean_anomaly() ' do
177
179
  assert_equal(3.7871218188949207, aliased_angles.mean_anomaly)
178
180
  end
179
181
 
@@ -188,22 +190,22 @@ describe 'tests ajd of 2455055.5' do
188
190
  aliased_angles.cosine_apparent_longitude)
189
191
  end
190
192
 
191
- it 'expected -0.7542060769936684 from \
193
+ it 'expected -0.7542060769936684 from \
192
194
  aliased_angles.cosine_true_longitude()? ' do
193
195
  assert_equal(-0.7542060769936684, aliased_angles.cosine_true_longitude)
194
196
  end
195
197
 
196
- it 'expected 0.9174818088112336 from \
198
+ it 'expected 0.9174818088112336 from \
197
199
  aliased_angles.cosine_true_obliquity()? ' do
198
200
  assert_equal(0.9174818088112336, aliased_angles.cosine_true_obliquity)
199
201
  end
200
202
 
201
- it 'expected 0.26426912722944046 from \
203
+ it 'expected 0.2642691272294404 from \
202
204
  aliased_angles.declination()? ' do
203
- assert_equal(0.26426912722944046, aliased_angles.declination)
205
+ assert_equal(0.2642691272294404, aliased_angles.declination)
204
206
  end
205
207
 
206
- it 'expected -0.04234904897476355 from \
208
+ it 'expected -0.04234904897476355 from \
207
209
  aliased_angles.delta_t_ecliptic()? ' do
208
210
  assert_equal(-0.04234904897476355, aliased_angles.delta_t_ecliptic)
209
211
  end
@@ -213,77 +215,77 @@ describe 'tests ajd of 2455055.5' do
213
215
  assert_equal(0.019768413456709915, aliased_angles.delta_t_elliptic)
214
216
  end
215
217
 
216
- it 'expected 0.016704576164208475 from \
218
+ it 'expected 0.016704576164208475 from \
217
219
  aliased_angles.eccentricity_earth_orbit()? ' do
218
220
  assert_equal(0.016704576164208475, \
219
221
  aliased_angles.eccentricity_earth_orbit)
220
222
  end
221
223
 
222
- it 'expected -0.019768413456709915 from \
224
+ it 'expected -0.019768413456709915 from \
223
225
  aliased_angles.equation_of_center()? ' do
224
226
  assert_equal(-0.019768413456709915, aliased_angles.equation_of_center)
225
227
  end
226
228
 
227
- it 'expected 2.445008945789877 from \
229
+ it 'expected 2.445008945789877 from \
228
230
  aliased_angles.geometric_mean_longitude()? ' do
229
231
  assert_equal(2.445008945789877, aliased_angles.geometric_mean_longitude)
230
232
  end
231
233
 
232
- it 'expected 1.9434600427973594 from \
234
+ it 'expected 1.585863261753274 from \
233
235
  aliased_angles.horizon_angle()? ' do
234
- assert_equal(1.9434600427973594, aliased_angles.horizon_angle)
236
+ assert_equal(1.585863261753274, aliased_angles.horizon_angle(1))
235
237
  end
236
238
 
237
- it 'expected 3.7871218188949207 from \
239
+ it 'expected 3.7871218188949207 from \
238
240
  aliased_angles.mean_anomaly()? ' do
239
241
  assert_equal(3.7871218188949207, aliased_angles.mean_anomaly)
240
242
  end
241
243
 
242
- it 'expected 2.444907382260759 from \
244
+ it 'expected 2.444907382260759 from \
243
245
  aliased_angles.mean_longitude_aries()? ' do
244
246
  assert_equal(2.444907382260759, aliased_angles.mean_longitude_aries)
245
247
  end
246
248
 
247
- it 'expected 0.4090707793981491 from \
249
+ it 'expected 0.4090707793981491 from \
248
250
  aliased_angles.mean_obliquity()? ' do
249
251
  assert_equal(0.4090707793981491, aliased_angles.mean_obliquity)
250
252
  end
251
253
 
252
- it 'expected 0.4090934409048494 from \
254
+ it 'expected 0.4090934409048494 from \
253
255
  aliased_angles.obliquity_correction()? ' do
254
256
  assert_equal(0.4090934409048494, aliased_angles.obliquity_correction)
255
257
  end
256
258
 
257
- it 'expected 2.467563113547267 from \
259
+ it 'expected 2.467563113547267 from \
258
260
  aliased_angles.right_ascension()? ' do
259
261
  assert_equal(2.467563113547267, aliased_angles.right_ascension)
260
262
  end
261
263
 
262
- it 'expected 0.6566577566139093 from \
264
+ it 'expected 0.6566577566139093 from \
263
265
  aliased_angles.sine_apparent_longitude()? ' do
264
266
  assert_equal(0.6566577566139093, aliased_angles.sine_apparent_longitude)
265
267
  end
266
268
 
267
- it 'expected 0.6566377946979757 from \
269
+ it 'expected 0.6566377946979757 from \
268
270
  aliased_angles.sine_true_longitude()? ' do
269
271
  assert_equal(0.6566377946979757, aliased_angles.sine_true_longitude)
270
272
  end
271
273
 
272
- it 'expected 3.767353405438211 from aliased_angles.true_anomaly()? ' do
274
+ it 'expected 3.767353405438211 from aliased_angles.true_anomaly()? ' do
273
275
  assert_equal(3.767353405438211, aliased_angles.true_anomaly)
274
276
  end
275
277
 
276
- it 'expected 2.4252405323331674 from \
278
+ it 'expected 2.4252405323331674 from \
277
279
  aliased_angles.true_longitude()? ' do
278
280
  assert_equal(2.4252405323331674, aliased_angles.true_longitude)
279
281
  end
280
282
 
281
- it 'expected 2.4449774607872907 from \
283
+ it 'expected 2.4449774607872907 from \
282
284
  aliased_angles.true_longitude_aries()? ' do
283
285
  assert_equal(2.4449774607872907, aliased_angles.true_longitude_aries)
284
286
  end
285
287
 
286
- it 'expected 0.4090934409048494 from \
288
+ it 'expected 0.4090934409048494 from \
287
289
  aliased_angles.true_obliquity()? ' do
288
290
  assert_equal(0.4090934409048494, aliased_angles.true_obliquity)
289
291
  end
@@ -13,6 +13,7 @@ describe 'aliased_displays using ajd of 2456885.0' do
13
13
  before(:each) do
14
14
  aliased_displays.ajd = 2_456_885.0
15
15
  ajd = aliased_displays.ajd
16
+ aliased_displays.ma_ta_set
16
17
  # check date for this ajd when needed.
17
18
  aliased_displays.date = aliased_displays.ajd_to_datetime(ajd)
18
19
 
@@ -36,9 +37,9 @@ describe 'aliased_displays using ajd of 2456885.0' do
36
37
  assert_equal '+142:35:33.356', aliased_displays.apparent_longitude_string
37
38
  end
38
39
 
39
- it 'expected "+013:58:51.521" \
40
+ it 'expected "+013:58:51.522" \
40
41
  from aliased_displays.declination_string()? ' do
41
- assert_equal '+013:58:51.521', aliased_displays.declination_string
42
+ assert_equal '+013:58:51.522', aliased_displays.declination_string
42
43
  end
43
44
 
44
45
  it 'expected "-04m, 29.2s" from \
@@ -81,9 +82,9 @@ describe 'aliased_displays using ajd of 2456885.0' do
81
82
  aliased_displays.julian_period_day_fraction_to_time(0)
82
83
  end
83
84
 
84
- it 'expected "+220:38:04.597" from \
85
+ it 'expected "+220:38:04.598" from \
85
86
  aliased_displays.mean_anomaly_string()? ' do
86
- assert_equal '+220:38:04.597', aliased_displays.mean_anomaly_string
87
+ assert_equal '+220:38:04.598', aliased_displays.mean_anomaly_string
87
88
  end
88
89
 
89
90
  it 'expected "+144:56:36.571" from \
@@ -91,19 +92,19 @@ describe 'aliased_displays using ajd of 2456885.0' do
91
92
  assert_equal '+144:56:36.571', aliased_displays.right_ascension_string
92
93
  end
93
94
 
94
- it 'expected "+219:24:27.797" from \
95
+ it 'expected "+219:24:27.798" from \
95
96
  aliased_displays.true_anomaly_string()? ' do
96
- assert_equal '+219:24:27.797', aliased_displays.true_anomaly_string
97
+ assert_equal '+219:24:27.798', aliased_displays.true_anomaly_string
97
98
  end
98
99
 
99
- it 'expected "+142:35:47.3184" from \
100
+ it 'expected "+142:35:47.318" from \
100
101
  aliased_displays.true_longitude_string()? ' do
101
102
  assert_equal '+142:35:47.318', aliased_displays.true_longitude_string
102
103
  end
103
104
 
104
- it 'expected "+023:26:06.163" from \
105
+ it 'expected "+023:26:06.164" from \
105
106
  aliased_displays.true_obliquity_string()? ' do
106
- assert_equal '+023:26:06.163', aliased_displays.true_obliquity_string
107
+ assert_equal '+023:26:06.164', aliased_displays.true_obliquity_string
107
108
  end
108
109
 
109
110
  end
@@ -1,32 +1,34 @@
1
1
  # aliased_utilities_spec.rb
2
-
2
+ gem 'minitest'
3
+ require 'minitest/autorun'
3
4
  lib = File.expand_path('../../../lib', __FILE__)
4
5
  $LOAD_PATH.unshift(lib) unless $LOAD_PATH.include?(lib)
5
6
  require 'eot'
6
7
 
7
- Eot_aliased_utilities = Eot.new
8
+ aliased_utilities = Eot.new
8
9
 
9
10
  describe 'tests ajd of 2456885.0' do
10
11
 
11
12
  before(:each) do
12
- Eot_aliased_utilities.ajd = 2_456_885.0
13
- ajd = Eot_aliased_utilities.ajd
13
+ aliased_utilities.ajd = 2_456_885.0
14
+ ajd = aliased_utilities.ajd
15
+ aliased_utilities.ma_ta_set
14
16
  # check date for this ajd when needed.
15
- Eot_aliased_utilities.date = Eot_aliased_utilities.ajd_to_datetime(ajd)
17
+ aliased_utilities.date = aliased_utilities.ajd_to_datetime(ajd)
16
18
  end
17
19
 
18
- it 'expected 2456885.0 for Eot_aliased_utilities.ajd'do
19
- assert_equal 2_456_885.0, Eot_aliased_utilities.ajd
20
+ it 'expected 2456885.0 for aliased_utilities.ajd'do
21
+ assert_equal 2_456_885.0, aliased_utilities.ajd
20
22
  end
21
23
 
22
- it 'expected 3.8508003966038915 for Eot_aliased_utilities.ma'do
23
- assert_equal 3.8508003966038915, Eot_aliased_utilities.ma
24
+ it 'expected 3.8508003966038915 for aliased_utilities.ma'do
25
+ assert_equal 3.8508003966038915, aliased_utilities.ma
24
26
  end
25
27
 
26
- it 'expected 0.0 returned by Eot_aliased_utilities.truncate() ' do
27
- assert_equal 0.0, Eot_aliased_utilities.truncate
28
- assert_equal 0.0, Eot_aliased_utilities.truncate(nil)
29
- assert_equal 0.0, Eot_aliased_utilities.truncate(0)
28
+ it 'expected 0.0 returned by aliased_utilities.truncate() ' do
29
+ assert_equal 0.0, aliased_utilities.truncate
30
+ assert_equal 0.0, aliased_utilities.truncate(nil)
31
+ assert_equal 0.0, aliased_utilities.truncate(0)
30
32
  end
31
33
 
32
34
  end