#NEWLINE# @ jcompile ( 'float64(float64)' , nopython = True ) #NEWLINE# def heliocentric_radius_vector( jme ) : #NEWLINE# #INDENT# r0 = 0.0 #NEWLINE# r1 = 0.0 #NEWLINE# r2 = 0.0 #NEWLINE# r3 = 0.0 #NEWLINE# r4 = 0.0 #NEWLINE# for row in range ( HELIO_RADIUS_TABLE . shape [ 1 ] ) : #NEWLINE# #INDENT# r0 += ( HELIO_RADIUS_TABLE [ ( 0 , row , 0 ) ] * np . cos ( ( HELIO_RADIUS_TABLE [ ( 0 , row , 1 ) ] + ( HELIO_RADIUS_TABLE [ ( 0 , row , 2 ) ] * jme ) ) ) ) #NEWLINE# r1 += ( HELIO_RADIUS_TABLE [ ( 1 , row , 0 ) ] * np . cos ( ( HELIO_RADIUS_TABLE [ ( 1 , row , 1 ) ] + ( HELIO_RADIUS_TABLE [ ( 1 , row , 2 ) ] * jme ) ) ) ) #NEWLINE# r2 += ( HELIO_RADIUS_TABLE [ ( 2 , row , 0 ) ] * np . cos ( ( HELIO_RADIUS_TABLE [ ( 2 , row , 1 ) ] + ( HELIO_RADIUS_TABLE [ ( 2 , row , 2 ) ] * jme ) ) ) ) #NEWLINE# r3 += ( HELIO_RADIUS_TABLE [ ( 3 , row , 0 ) ] * np . cos ( ( HELIO_RADIUS_TABLE [ ( 3 , row , 1 ) ] + ( HELIO_RADIUS_TABLE [ ( 3 , row , 2 ) ] * jme ) ) ) ) #NEWLINE# r4 += ( HELIO_RADIUS_TABLE [ ( 4 , row , 0 ) ] * np . cos ( ( HELIO_RADIUS_TABLE [ ( 4 , row , 1 ) ] + ( HELIO_RADIUS_TABLE [ ( 4 , row , 2 ) ] * jme ) ) ) ) #NEWLINE# #UNINDENT# r = ( ( ( ( ( r0 + ( r1 * jme ) ) + ( r2 * ( jme ** 2 ) ) ) + ( r3 * ( jme ** 3 ) ) ) + ( r4 * ( r2 ** 4 ) ) ) / ( 10 ** 8 ) ) #NEWLINE# return r