[cig-commits] [commit] rajesh-petsc-schur: Bug fixes: missing for loop parenthesis in Construct_arrays.c, caps_per_proc for loop variable 'm' changed to constant CPPR in General_matrix_functions.c (9be20ef)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/citcoms

On branch  : rajesh-petsc-schur
Link       : https://github.com/geodynamics/citcoms/compare/464e1b32299b15819f93efd98d969cddb84dfe51...f97ae655a50bdbd6dac1923a3471ee4dae178fbd

>---------------------------------------------------------------

commit 9be20efa6e277e8341746bf0e8d505e81a8a8d74
Author: Rajesh Kommu <rajesh.kommu at gmail.com>
Date:   Tue Sep 23 12:22:12 2014 -0700

    Bug fixes: missing for loop parenthesis in Construct_arrays.c, caps_per_proc for loop variable 'm' changed to constant CPPR in General_matrix_functions.c


>---------------------------------------------------------------

9be20efa6e277e8341746bf0e8d505e81a8a8d74
 lib/Construct_arrays.c         | 3 ++-
 lib/General_matrix_functions.c | 2 +-
 2 files changed, 3 insertions(+), 2 deletions(-)

diff --git a/lib/Construct_arrays.c b/lib/Construct_arrays.c
index d16b688..c174b65 100644
--- a/lib/Construct_arrays.c
+++ b/lib/Construct_arrays.c
@@ -563,7 +563,7 @@ void construct_sub_element(E)
      int lev,nox,noy,noz,nnn,elx,elz,ely,elzu,elxu,elt,eltu;
 
 
-  for(lev=E->mesh.levmax-1;lev>=E->mesh.levmin;lev--)
+  for(lev=E->mesh.levmax-1;lev>=E->mesh.levmin;lev--) {
           elx = E->lmesh.ELX[lev];
 	  elz = E->lmesh.ELZ[lev];
 	  ely = E->lmesh.ELY[lev];
@@ -593,6 +593,7 @@ void construct_sub_element(E)
                                  + offset[l].vector[2] * elzu * elxu;
 		      }
 		  }
+  }
 }
 
 
diff --git a/lib/General_matrix_functions.c b/lib/General_matrix_functions.c
index c8ab375..4090895 100644
--- a/lib/General_matrix_functions.c
+++ b/lib/General_matrix_functions.c
@@ -499,7 +499,7 @@ void gauss_seidel(E,d0,F,Ad,acc,cycles,level,guess)
                  Ad[eqn3] += B3[j]*UU;
                  }
 
-            if (!(E->NODE[level][m][i]&OFFSIDE))   {
+            if (!(E->NODE[level][CPPR][i]&OFFSIDE))   {
                E->temp[eqn1] = (F[eqn1] - Ad[eqn1])*E->BI[level][eqn1];
                E->temp[eqn2] = (F[eqn2] - Ad[eqn2])*E->BI[level][eqn2];
                E->temp[eqn3] = (F[eqn3] - Ad[eqn3])*E->BI[level][eqn3];