static void VI_solve_lt_d( vsip_mview_d *R, vsip_scalar_d alpha, vsip_mview_d *B) { vsip_length N = R->row_length; vsip_mview_d XX = *B; vsip_mview_d *X = &XX; vsip_vview_d xx, rr; vsip_vview_d *x = VI_mrowview_d(X,(vsip_index)0,&xx); vsip_vview_d *r_d = VI_mrowview_d(R,(vsip_index) 0,&rr); vsip_vview_d rr_m = rr; vsip_vview_d *r_m = &rr_m; vsip_stride d_s = R->row_stride + R->col_stride; r_d->length = 1; r_m->length = 0; X->col_length = 1; if(alpha == 0) return; vsip_svmul_d(alpha,x,x); vsip_svmul_d((vsip_scalar_d) 1.0/(VI_VGET_D(r_d,0)),x,x); N--; while(N-- >0){ r_d->offset += d_s; r_m->length++; r_m->offset += R->col_stride; x->offset += X->col_stride; vsip_svmul_d(alpha,x,x); VI_vsubvmprodIP_d(r_m,X,x); vsip_svmul_d((vsip_scalar_d) 1.0/(VI_VGET_D(r_d,0)),x,x); X->col_length++; } return; }
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; }
static void VI_svmprod_d( vsip_scalar_d beta, vsip_vview_d *v, vsip_mview_d *A, vsip_vview_d *w) { vsip_vmprod_d(v,A,w); vsip_svmul_d(beta,w,w); return; }
int main(){vsip_init((void*)0); { vsip_mview_d *A = vsip_mcreate_d(N,N,VSIP_COL,0); vsip_vview_d *x0 = vsip_vcreate_d(N,0); vsip_vview_d *b = vsip_vcreate_d(N,0); vsip_mview_d *X = vsip_mcreate_d(N,3,VSIP_ROW,0); vsip_mview_d *XT = vsip_mcreate_d(N,3,VSIP_COL,0); vsip_vramp_d(1,1,x0); vsip_vmul_d(x0,x0,x0); vsip_mput_d(A,0,0,-3); vsip_mput_d(A,0,1,7); vsip_mput_d(A,0,2,10); vsip_mput_d(A,0,3,12); vsip_mput_d(A,1,0,0); vsip_mput_d(A,1,1,13); vsip_mput_d(A,1,2,18); vsip_mput_d(A,1,3,6); vsip_mput_d(A,2,0,2); vsip_mput_d(A,2,1,-9); vsip_mput_d(A,2,2,6); vsip_mput_d(A,2,3,3); vsip_mput_d(A,3,0,1); vsip_mput_d(A,3,1,2); vsip_mput_d(A,3,2,3); vsip_mput_d(A,3,3,4); {int i,j; printf("A = [\n"); for(i=0; i<N; i++){ for(j=0; j<N; j++) printf("%9.2f%s",vsip_mget_d(A,i,j),(j == N-1) ? "":","); printf(";\n"); } printf("]\n"); } { int k; vsip_vview_d *x; vsip_length L = vsip_mgetrowlength_d(X); for(k=0; k<L; k++){ x = vsip_mcolview_d(X,k); vsip_mvprod_d(A,x0,b); vsip_vcopy_d_d(b,x); vsip_svmul_d(2.0,x0,x0); vsip_vdestroy_d(x); } {int i,j; printf("X = [\n");for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(X,i,j),(j == 2) ? "":","); printf(";\n"); } printf("]\n"); } { vsip_lu_d* luAop = vsip_lud_create_d(N); vsip_mcopy_d_d(X,XT); if(luAop == NULL) exit(1); vsip_lud_d(luAop,A); vsip_lusol_d(luAop,0,X); vsip_lusol_d(luAop,1,XT); vsip_lud_destroy_d(luAop); } } {int i,j; printf("A\\X\n"); for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(X,i,j),(j == 2) ? "":","); printf(";\n"); } } {int i,j; printf("A'\\X\n"); for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(XT,i,j),(j == 2) ? "":","); printf(";\n"); } } { vsip_valldestroy_d(b); vsip_valldestroy_d(x0); vsip_malldestroy_d(X); vsip_malldestroy_d(A); } }vsip_finalize((void*)0);return 1; }
int main(){vsip_init((void*)0); { vsip_mview_d *A = vsip_mcreate_d(N,N,VSIP_ROW,0); vsip_vview_d *x0 = vsip_vcreate_d(N,0); vsip_vview_d *b = vsip_vcreate_d(N,0); vsip_mview_d *X = vsip_mcreate_d(N,3,VSIP_ROW,0); vsip_mview_d *XT = vsip_mcreate_d(N,3,VSIP_ROW,0); vsip_vramp_d(1,1,x0); vsip_vmul_d(x0,x0,x0); { int i; vsip_vview_d *r = vsip_mrowview_d(A,0); vsip_offset o = vsip_vgetoffset_d(r); for(i=0; i<N; i++){ vsip_vputoffset_d(r,o); vsip_vramp_d(i,1,r); o +=N; } vsip_vdestroy_d(r); } { vsip_vview_d *r = vsip_mdiagview_d(A,0); vsip_vfill_d(0,r); vsip_vdestroy_d(r); } {int i,j; printf("A = [\n"); for(i=0; i<N; i++){ for(j=0; j<N; j++) printf("%9.2f%s",vsip_mget_d(A,i,j),(j == N-1) ? "":","); printf(";\n"); } printf("]\n"); } { int k; vsip_vview_d *x; vsip_length L = vsip_mgetrowlength_d(X); for(k=0; k<L; k++){ x = vsip_mcolview_d(X,k); vsip_mvprod_d(A,x0,b); vsip_vcopy_d_d(b,x); vsip_svmul_d(2.0,x0,x0); vsip_vdestroy_d(x); } {int i,j; printf("X = [\n"); for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(X,i,j),(j == 2) ? "":","); printf(";\n"); } printf("]\n"); } {vsip_lu_d* luAop = vsip_lud_create_d(N); vsip_mcopy_d_d(X,XT); if(luAop == NULL) exit(1); vsip_lud_d(luAop,A); vsip_lusol_d(luAop,0,X); vsip_lusol_d(luAop,1,XT); vsip_lud_destroy_d(luAop); } } {int i,j; printf("A\\X\n"); for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(X,i,j),(j == 2) ? "":","); printf(";\n"); } } {int i,j; printf("A'\\X\n"); for(i=0; i<N; i++){ for(j=0; j<3; j++) printf("%9.2f%s",vsip_mget_d(XT,i,j),(j == 2) ? "":","); printf(";\n"); } } { vsip_valldestroy_d(b); vsip_valldestroy_d(x0); vsip_malldestroy_d(X); vsip_malldestroy_d(A); } }vsip_finalize((void*)0);return 1; }