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;
}
Example #2
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;
}