@@ -461,9 +461,12 @@ def sum_mult_cos_add_mult(arr, x):
461461 for row in range (arr .shape [0 ]):
462462 s += arr [row , 0 ] * np .cos (arr [row , 1 ] + arr [row , 2 ] * x )
463463 else :
464- s = (arr [:, [0 ]] * np .cos (arr [:, [1 ]] + arr [:, [2 ]] * np .expand_dims (x , axis = 0 ))).sum (axis = 0 )
464+ s = (arr [:, [0 ]] * np .cos (
465+ arr [:, [1 ]] + arr [:, [2 ]] * np .expand_dims (x , axis = 0 )
466+ )).sum (axis = 0 )
465467 return s
466468
469+
467470@jcompile ('float64(float64)' , nopython = True )
468471def heliocentric_longitude (jme ):
469472 l0 = sum_mult_cos_add_mult (L0 , jme )
@@ -562,7 +565,6 @@ def moon_ascending_longitude(julian_ephemeris_century):
562565 nopython = True )
563566def longitude_obliquity_nutation (julian_ephemeris_century , x0 , x1 , x2 , x3 , x4 ,
564567 out ):
565-
566568 if USE_NUMBA :
567569 delta_psi_sum = 0.0
568570 delta_eps_sum = 0.0
@@ -594,9 +596,13 @@ def longitude_obliquity_nutation(julian_ephemeris_century, x0, x1, x2, x3, x4,
594596 NUTATION_YTERM_ARRAY [:, [4 ]]* np .expand_dims (x4 , axis = 0 )
595597 )
596598
597- delta_psi_sum = ((a + b * julian_ephemeris_century ) * np .sin (arg )).sum (axis = 0 )
598- delta_eps_sum = ((c + d * julian_ephemeris_century ) * np .cos (arg )).sum (axis = 0 )
599-
599+ delta_psi_sum = (
600+ (a + b * julian_ephemeris_century ) * np .sin (arg )
601+ ).sum (axis = 0 )
602+ delta_eps_sum = (
603+ (c + d * julian_ephemeris_century ) * np .cos (arg )
604+ ).sum (axis = 0 )
605+
600606 delta_psi = delta_psi_sum * 1.0 / 36000000
601607 delta_eps = delta_eps_sum * 1.0 / 36000000
602608 # seems like we ought to be able to return a tuple here instead
0 commit comments