45 int LDBxNRHS = LDB*NRHS;
55 if ((!A1)||(!A2)||(!B1)||(!B2)){
56 printf(
"Out of Memory \n ");
62 printf(
"-- PLASMA is initialized to run on %d cores. \n",cores);
65 LAPACKE_clarnv_work(
IONE,
ISEED, LDAxN, A1);
66 for ( i = 0; i < N; i++)
67 for ( j = 0; j < N; j++)
68 A2[LDA*j+i] = A1[LDA*j+i];
71 LAPACKE_clarnv_work(
IONE,
ISEED, LDBxNRHS, B1);
72 for ( i = 0; i < N; i++)
73 for ( j = 0; j < NRHS; j++)
74 B2[LDB*j+i] = B1[LDB*j+i];
91 if ((info_solution != 0)|(info != 0))
92 printf(
"-- Error in CGETRS example ! \n");
94 printf(
"-- Run of CGETRS example successful ! \n");
96 free(A1); free(A2); free(B1); free(B2); free(IPIV); free(L);
110 float Rnorm, Anorm, Xnorm, Bnorm;
112 float *work = (
float *)malloc(N*
sizeof(
float));
115 eps = LAPACKE_slamch_work(
'e');
124 cblas_cgemm(
CblasColMajor,
CblasNoTrans,
CblasNoTrans, N, NRHS, N,
CBLAS_SADDR(alpha), A1, LDA, B2, LDB,
CBLAS_SADDR(beta), B1, LDB);
127 printf(
"============\n");
128 printf(
"Checking the Residual of the solution \n");
129 printf(
"-- ||Ax-B||_oo/((||A||_oo||x||_oo+||B||_oo).N.eps) = %e \n",Rnorm/((Anorm*Xnorm+Bnorm)*N*eps));
131 if ( isnan(Rnorm/((Anorm*Xnorm+Bnorm)*N*eps)) || (Rnorm/((Anorm*Xnorm+Bnorm)*N*eps) > 10.0) ){
132 printf(
"-- The solution is suspicious ! \n");
136 printf(
"-- The solution is CORRECT ! \n");
142 return info_solution;