int main (){vsip_init((void*)0);
{
int i;
vsip_cvview_d* vector = vsip_cvcreate_d(L,0);
vsip_vview_d* realIn = vsip_vrealview_d(vector);
vsip_vview_d* imagIn = vsip_vimagview_d(vector);
vsip_fft_d* fftplan = vsip_ccfftip_create_d(L,1,-1,0,0);
vsip_fft_d* fftplanI = vsip_ccfftip_create_d(L,(double)1/L,1,0,0);
char printDataOne[L][20],
printDataTwo[L][30];
/* put some data in realIn*/
vsip_vramp_d(0, 2 * VU_PI * F, realIn);
vsip_vcos_d(realIn,realIn);
/* print realIn */
for(i=0; i<L; i++){
sprintf(printDataOne[i],"%2d input -> %7.4f",i, vsip_vget_d(realIn,i));
}
/*make sure imagIn is full of zeros*/
vsip_vfill_d(0,imagIn);
/*find the fft*/
vsip_ccfftip_d(fftplan,vector);
for(i=0; i<L; i++){
sprintf(printDataTwo[i]," fft -> (%7.3f, %7.3f)",
vsip_real_d(vsip_cvget_d(vector,i)),
vsip_imag_d(vsip_cvget_d(vector,i)));
}
/*invert the fft*/
vsip_ccfftip_d(fftplanI,vector);
/*print it */
for(i=0; i<L; i++){
printf("%s %s ifft -> (%7.3f, %7.3f)\n",
printDataOne[i],
printDataTwo[i],
vsip_real_d(vsip_cvget_d(vector,i)),
vsip_imag_d(vsip_cvget_d(vector,i)));
}
{
vsip_vdestroy_d((void*)imagIn);
vsip_vdestroy_d((void*)realIn);
vsip_cblockdestroy_d(vsip_cvdestroy_d((void*)vector));
vsip_fft_destroy_d((void*)fftplan);
vsip_fft_destroy_d((void*)fftplanI);
}
} vsip_finalize((void*)0);return 0;
}
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_init((void*)0);
{
int i;
vsip_cblock_d* blockOut = vsip_cblockcreate_d(2 * L, 0);
vsip_cvview_d* vectorOut = vsip_cvbind_d(blockOut,0,2,L);
vsip_cvview_d* vectorIn = vsip_cvcreate_d(L,0);
vsip_vview_d* realIn = vsip_vrealview_d(vectorIn);
vsip_vview_d* imagIn = vsip_vimagview_d(vectorIn);
vsip_fft_d* fftplan = vsip_ccfftop_create_d(L,1,-1,0,0);
vsip_fft_d* fftplanI = vsip_ccfftop_create_d(L,(double)1/L,1,0,0);
char inputdata[L][20];
/* put some data in realIn*/
vsip_vramp_d(0, 2 * PI * F, realIn);
vsip_vcos_d(realIn,realIn);
/* print realIn */
for(i=0; i<L; i++)
sprintf(inputdata[i],"%2d input-> %7.3f",i, vsip_vget_d(realIn,i));
/*make sure imagIn is full of zeros*/
vsip_vfill_d(0,imagIn);
/*find the fft*/
vsip_ccfftop_d(fftplan,vectorIn,vectorOut);
/*invert the fft*/
vsip_ccfftop_d(fftplanI,vectorOut,vectorIn);
/*print it */
for(i=0; i<L; i++)
printf("%s fft -> (%7.3f, %7.3f) ifft -> (%7.3f, %7.3f)\n",
inputdata[i],
vsip_real_d(vsip_cvget_d(vectorOut,i)),
vsip_imag_d(vsip_cvget_d(vectorOut,i)),
vsip_real_d(vsip_cvget_d(vectorIn,i)),
vsip_imag_d(vsip_cvget_d(vectorIn,i)));
{
vsip_vdestroy_d(imagIn);
vsip_vdestroy_d(realIn);
vsip_cblockdestroy_d(vsip_cvdestroy_d(vectorIn));
vsip_cblockdestroy_d(vsip_cvdestroy_d(vectorOut));
vsip_fft_destroy_d(fftplan);
vsip_fft_destroy_d(fftplanI);
}
} vsip_finalize((void*)0);return 0;
}
int main()
{
int i;
vsip_cscalar_d z;
vsip_scalar_d data[N]; /* a public data space for I/O */
vsip_fft_d *rcfftNop;
vsip_block_d *block;
vsip_vview_d *xin;
vsip_cvview_d *yout;
vsip_init((void *)0);
rcfftNop = /* Create an out-of-place Real->Cmplx N-pt FFT */
vsip_rcfftop_create_d(N, 1.0, 1, VSIP_ALG_TIME);
/* Create a block object and bind it to the array data */
block = vsip_blockbind_d(data, N, VSIP_MEM_NONE);
xin = vsip_vbind_d(block, 0, 1, N);
/* Create another block and complex vector view for the
symmetric output */
yout = vsip_cvcreate_d((N/2)+1, VSIP_MEM_NONE);
/* Admit block to VSIPL for processing and initialize with a
linear ramp */
vsip_blockadmit_d(block, VSIP_FALSE);
vsip_vramp_d(0.0, 1.0, xin);
/* Compute an out-of-place Real->Cmplx N-pt FFT using the
rcfftNop object */
vsip_rcfftop_d(rcfftNop, xin, yout);
/* print it */
printf("Real Input Vector\n");
for(i=0; i<N; i++)
{
printf("%g\n", vsip_vget_d(xin,i));
}
printf("\nComplex Output Vector\n");
for(i=0; i<(N/2)+1; i++)
{
z = vsip_cvget_d(yout,i);
printf(SCMPLX "\n", ACMPLX(z));
}
/* Destroy the rcfftNop, blocks, and view objects */
vsip_fft_destroy_d(rcfftNop);
vsip_valldestroy_d(xin);
vsip_cvalldestroy_d(yout);
vsip_finalize((void *)0);
return(0);
}
static void destroy(fft_type *fft) { vsip_fft_destroy_d(fft);}
