added missing mkl call in parcsrmv

Author: liruipeng

SRC/parcsrmv.c

@@ -33,10 +33,10 @@ void pEVSL_ParcsrMatvecCommBegin(pevsl_Parcsr *A, double *x) {
 void pEVSL_ParcsrMatvecCommEnd(pevsl_Parcsr *A) {
     int err;
     err = MPI_Waitall(A->comm_handle->num_proc_send_to, A->comm_handle->send_requests, \
-                A->comm_handle->send_status);
+                      A->comm_handle->send_status);
     PEVSL_CHKERR(err != MPI_SUCCESS);
     err = MPI_Waitall(A->comm_handle->num_proc_recv_from, A->comm_handle->recv_requests, \
-                A->comm_handle->recv_status);
+                      A->comm_handle->recv_status);
     PEVSL_CHKERR(err != MPI_SUCCESS);
 }
 

SRC/spmat.c

@@ -104,7 +104,14 @@ void dcsrgemv(char trans, int nrow, int ncol, double alp, double *a,
 
 // y = alp*A*x + bet*y
 int pEVSL_MatvecGen(double alp, pevsl_Csr *A, double *x, double bet, double *y) {
+
+#ifdef USE_MKL
+  char cN = 'N';
+  mkl_dcsrmv(&cN, &(A->nrows), &(A->ncols), &alp, "GXXCXX", A->a, A->ja, A->ia, 
+             A->ia+1, x, &bet, y);
+#else
   dcsrgemv('N', A->nrows, A->ncols, alp, A->a, A->ia, A->ja, x, bet, y);
+#endif
 
   return 0;
 }