1 # compute arctangent table for CORDIC computations in fttrigon.c
4 #units = 64*65536.0 # don't change !!
11 print "table of arctan( 1/2^n ) for PI = " + repr(units
/65536.0) + " units"
15 x
= 0.5**n
# tangent value
17 angle
= math
.atan(x
) # arctangent
18 angle2
= round(angle
*scale
) # arctangent in FT_Angle units
23 sys
.stdout
.write( comma
+ repr( int(angle2
) ) )
26 shrink
/= math
.sqrt( 1 + x
*x
)
29 print "shrink factor = " + repr( shrink
)
30 print "shrink factor 2 = " + repr( int( shrink
* (2**32) ) )
31 print "expansion factor = " + repr( 1/shrink
)