int main(){vsip_init((void*)0); { void VU_vfprintyg_d(char*,vsip_vview_d*,char*); void VU_vfreqswapIP_d(vsip_vview_d*); vsip_vview_d* Cw = vsip_vcreate_cheby_d(Nlength,ripple,0); vsip_fft_d *fft = vsip_ccfftip_create_d(Nlength,1.0,VSIP_FFT_FWD,0,0); vsip_cvview_d* FCW = vsip_cvcreate_d(Nlength,0); /*printf("CW = "); VU_vprintm_d("%6.8f ;\n",Cw); */ VU_vfprintyg_d("%6.8f\n",Cw,"Cheby_Window"); vsip_cvfill_d(vsip_cmplx_d(0,0),FCW); { vsip_vview_d *rv = vsip_vrealview_d(FCW); vsip_vcopy_d_d(Cw,rv); vsip_ccfftip_d(fft,FCW); vsip_vcmagsq_d(FCW,rv); { vsip_index ind; vsip_scalar_d max = vsip_vmaxval_d(rv,&ind); vsip_scalar_d min = max/(10e12); vsip_vclip_d(rv,min,max,min,max,rv); } vsip_vlog10_d(rv,rv); vsip_svmul_d(10,rv,rv); VU_vfreqswapIP_d(rv); VU_vfprintyg_d("%6.8f\n",rv,"Cheby_Window_Frequency_Response"); vsip_vdestroy_d(rv); } vsip_fft_destroy_d(fft); vsip_valldestroy_d(Cw); vsip_cvalldestroy_d(FCW); } vsip_finalize((void*)0); return 0; }
int main() { vsip_vview_d *dataIn, *dataClip; vsip_init((void *)0); dataIn = vsip_vcreate_d(L, VSIP_MEM_NONE); dataClip = vsip_vcreate_d(L, VSIP_MEM_NONE); /* make some data */ vsip_vramp_d(0.0, (2 * PI)/(L - 1.0), dataIn); vsip_vcos_d(dataIn,dataIn); vsip_vclip_d(dataIn,-.8,.8,-.8,.8,dataClip); printf("clip Cosine between -.8 and .8\n in => out\n "); { int i; for(i=0; i<L; i++) printf("%7.4f => %7.4f\n", vsip_vget_d(dataIn,i),vsip_vget_d(dataClip,i)); } vsip_finalize((void *)0); return 0; }