fromamaranthimport*importnumpyasnpfromnumpyimportsin,uint32,float64,powertau=2*np.pitemp_ramp_tables=np.empty((8,1024),dtype=float64)foriinrange(0,8):base_f0=power(2,i+4)n_harmonics=np.floor(18000/base_f0)# Generate 1024 samples normalized to 1/2 waveform (hence /2*tau))# and including all harmonics up to 18kHzh=np.arange(1,n_harmonics+1)temp_ramp_tables[i]=np.sum(sin(h*base_f0/ (2*tau))/h)# compress peak-to-peak relative to the largest amplitudetemp_ramp_tables/=np.max(temp_ramp_tables)# convert to positive half of unsigned 24-bit rangeramp_tables=temp_ramp_tables*(power(2,23)-1).astype(uint32)deltemp_ramp_tables

# convert to positive half of unsigned 24-bit rangeramp_tables= [               [Const(x,unsigned(23))forxintemp_ramp_tables[y]*(power(2,23)-1).astype(uint32)]foryinrange(0,8)              ]deltemp_ramp_tables