44 int LDBxNRHS = LDB*NRHS;
54 if ((!A1)||(!A2)||(!B1)||(!B2)){
55 printf(
"Out of Memory \n ");
61 printf(
"-- PLASMA is initialized to run on %d cores. \n",cores);
64 LAPACKE_zlarnv_work(
IONE,
ISEED, LDAxN, A1);
65 for ( i = 0; i < N; i++)
66 for ( j = 0; j < N; j++)
67 A2[LDA*j+i] = A1[LDA*j+i];
70 LAPACKE_zlarnv_work(
IONE,
ISEED, LDBxNRHS, B1);
71 for ( i = 0; i < N; i++)
72 for ( j = 0; j < NRHS; j++)
73 B2[LDB*j+i] = B1[LDB*j+i];
82 if ((info_solution != 0)|(info != 0))
83 printf(
"-- Error in ZGESV example ! \n");
85 printf(
"-- Run of ZGESV example successful ! \n");
87 free(A1); free(A2); free(B1); free(B2); free(IPIV); free(L);
101 double Rnorm, Anorm, Xnorm, Bnorm;
103 double *work = (
double *)malloc(N*
sizeof(
double));
106 eps = LAPACKE_dlamch_work(
'e');
115 cblas_zgemm(
CblasColMajor,
CblasNoTrans,
CblasNoTrans, N, NRHS, N,
CBLAS_SADDR(alpha), A1, LDA, B2, LDB,
CBLAS_SADDR(beta), B1, LDB);
118 printf(
"============\n");
119 printf(
"Checking the Residual of the solution \n");
120 printf(
"-- ||Ax-B||_oo/((||A||_oo||x||_oo+||B||_oo).N.eps) = %e \n",Rnorm/((Anorm*Xnorm+Bnorm)*N*eps));
122 if ( isnan(Rnorm/((Anorm*Xnorm+Bnorm)*N*eps)) || (Rnorm/((Anorm*Xnorm+Bnorm)*N*eps) > 10.0) ){
123 printf(
"-- The solution is suspicious ! \n");
127 printf(
"-- The solution is CORRECT ! \n");
133 return info_solution;