diff --git a/puremd_rc_1003/sPuReMD/src/GMRES.c b/puremd_rc_1003/sPuReMD/src/GMRES.c
index a1d8c66e7fc1b47ca7a16c1c01e79eb4d6c899dc..4aec2adea2bf820ebae9e1a02331e87748557028 100644
--- a/puremd_rc_1003/sPuReMD/src/GMRES.c
+++ b/puremd_rc_1003/sPuReMD/src/GMRES.c
@@ -46,7 +46,7 @@ static void Sparse_MatVec( const sparse_matrix * const A,
     Vector_MakeZero( b, n );
 
     #pragma omp parallel \
-        default(none) shared(n, b, b_local) private(si, ei, H, i, j, k, tid)
+        default(none) shared(n, b_local) private(si, ei, H, i, j, k, tid)
     {
 #ifdef _OPENMP
         tid = omp_get_thread_num();
@@ -131,7 +131,7 @@ static void Sparse_MatVec_Vector_Add( const sparse_matrix * const R,
     Vector_MakeZero( b, R->n );
 
     #pragma omp parallel \
-        default(none) shared(b, b_local) private(si, ei, i, k, tid)
+        default(none) shared(b_local) private(si, ei, i, k, tid)
     {
 #ifdef _OPENMP
         tid = omp_get_thread_num();
@@ -267,7 +267,7 @@ static void Jacobi_Iter( const sparse_matrix * const R, const TRIANGULARITY tri,
     static real *Dinv_b, *rp, *rp2, *rp3;
 
     #pragma omp parallel \
-        default(none) shared(b, b_local, Dinv_b, rp, rp2, rp3, iter, stderr) private(si, ei, i, k, tid)
+        default(none) shared(b_local, Dinv_b, rp, rp2, rp3, iter, stderr) private(si, ei, i, k, tid)
     {
 #ifdef _OPENMP
         tid = omp_get_thread_num();