int main(){vsip_init((void*)0); { vsip_vview_d *a = vsip_vcreate_d(N,0), *r = vsip_vcreate_d(P,0); vsip_mview_d *B = vsip_mcreate_d(P,N,VSIP_ROW,0); int i, j; vsip_vramp_d(1.0,1.0,a); for(i=0;i<P;i++) for(j=0;j<N;j++) vsip_mput_d(B,i,j,1 + i * j); printf("\n vector input \n"); vsip_mvprod_d(B,a,r); for(j=0;j<N;j++){ printf(" %6.2f ",vsip_vget_d(a,j)); } printf("\n"); printf("\n matrix input \n"); {for(i=0; i<P; i++) {for(j=0;j<N;j++){ printf(" %6.2f ",vsip_mget_d(B,i,j)); } printf(";\n"); } } printf("\n vector output \n"); for(j=0;j<P;j++){ printf("%6.2f ",vsip_vget_d(r,j)); } printf("\n"); vsip_valldestroy_d(a); vsip_valldestroy_d(r); vsip_malldestroy_d(B); }vsip_finalize((void*)0);return 0; }
void test_real_d(vsip_length N, vsip_major axis) { vsip_length rows = 5; vsip_length cols = N; vsip_length rows2 = 5; vsip_length cols2 = N/2 + 1; if (axis == VSIP_COL) { rows = N; cols = 5; rows2 = N/2 + 1; cols2 = 5; } vsip_fftm_d *f_fftm = vsip_rcfftmop_create_d(rows, cols, 1.0, axis, 0, VSIP_ALG_SPACE); vsip_fftm_d *i_fftm = vsip_crfftmop_create_d(rows, cols, 1.0/N, axis, 0, VSIP_ALG_SPACE); vsip_fftm_attr_d attr; vsip_fftm_getattr_d(f_fftm, &attr); test_assert(attr.input.r == rows && attr.input.c == cols); test_assert(attr.output.r == rows2 && attr.output.c == cols2); vsip_fftm_getattr_d(i_fftm, &attr); test_assert(attr.input.r == rows2 && attr.input.c == cols2); test_assert(attr.output.r == rows && attr.output.c == cols); vsip_mview_d *in = vsip_mcreate_d(rows, cols, VSIP_ROW, VSIP_MEM_NONE); vsip_randstate *rng = vsip_randcreate(0, 1, 1, VSIP_PRNG); vsip_mrandu_d(rng, in); vsip_randdestroy(rng); vsip_mview_d *inv = vsip_mcreate_d(rows, cols, VSIP_ROW, VSIP_MEM_NONE); vsip_cmview_d *out = vsip_cmcreate_d(rows2, cols2, VSIP_ROW, VSIP_MEM_NONE); vsip_rcfftmop_d(f_fftm, in, out); vsip_crfftmop_d(i_fftm, out, inv); test_assert(merror_db_d(inv, in) < -100.); vsip_malldestroy_d(inv); vsip_cmalldestroy_d(out); vsip_malldestroy_d(in); vsip_fftm_destroy_d(i_fftm); vsip_fftm_destroy_d(f_fftm); }
int main(){vsip_init((void*)0); { vsip_mview_d *A = vsip_mcreate_d(M,P,VSIP_ROW,0), *B = vsip_mcreate_d(P,N,VSIP_ROW,0), *R = vsip_mcreate_d(M,N,VSIP_ROW,0); int i, j; for(i=0;i<M;i++) for(j=0;j<P;j++) vsip_mput_d(A,i,j,1 + sqrt(i*j)); for(i=0;i<P;i++) for(j=0;j<N;j++) vsip_mput_d(B,i,j,1 + i * j); printf("\n A input \n"); vsip_mprod_d(A,B,R); {for(i=0;i<M;i++) {for(j=0;j<P;j++){ printf(": %6.2f ",vsip_mget_d(A,i,j)); } printf(":\n"); } } printf("\n B input \n"); {for(i=0;i<P;i++) {for(j=0;j<N;j++){ printf(": %6.2f ",vsip_mget_d(B,i,j)); } printf(":\n"); } } printf("\n R output \n"); {for(i=0;i<M;i++) {for(j=0;j<N;j++){ printf(": %6.2f ",vsip_mget_d(R,i,j)); } printf(":\n"); } } vsip_malldestroy_d(A); vsip_malldestroy_d(B); vsip_malldestroy_d(R); }vsip_finalize((void*)0);return 0; }
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; }