[cig-commits] [commit] rajesh-petsc-schur: Changed the shape of E->elt_k as part of caps_per_proc removal (509b363)

cig_noreply at geodynamics.org cig_noreply at geodynamics.org
Wed Nov 5 19:12:29 PST 2014


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

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

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

commit 509b363aa850c1cfe708999e725aa3f74ecfe445
Author: Rajesh Kommu <rajesh.kommu at gmail.com>
Date:   Tue Sep 23 18:37:49 2014 -0700

    Changed the shape of E->elt_k as part of caps_per_proc removal


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

509b363aa850c1cfe708999e725aa3f74ecfe445
 lib/Construct_arrays.c     |  6 +++---
 lib/Drive_solvers.c        |  2 +-
 lib/Element_calculations.c | 18 +++++++++---------
 lib/global_defs.h          |  2 +-
 4 files changed, 14 insertions(+), 14 deletions(-)

diff --git a/lib/Construct_arrays.c b/lib/Construct_arrays.c
index 974ac67..5eeb907 100644
--- a/lib/Construct_arrays.c
+++ b/lib/Construct_arrays.c
@@ -616,12 +616,12 @@ void construct_elt_ks(E)
 
 	for(el=1;el<=E->lmesh.NEL[lev];el++)    {
 
-	    get_elt_k(E,el,E->elt_k[lev][CPPR][el].k,lev,0);
+	    get_elt_k(E,el,E->elt_k[lev][el].k,lev,0);
 
 	    if (E->control.augmented_Lagr)
-	        get_aug_k(E,el,E->elt_k[lev][CPPR][el].k,lev);
+	        get_aug_k(E,el,E->elt_k[lev][el].k,lev);
 
-            build_diagonal_of_K(E,el,E->elt_k[lev][CPPR][el].k,lev);
+            build_diagonal_of_K(E,el,E->elt_k[lev][el].k,lev);
 
 	    }
 
diff --git a/lib/Drive_solvers.c b/lib/Drive_solvers.c
index e6dcb7c..f31c481 100644
--- a/lib/Drive_solvers.c
+++ b/lib/Drive_solvers.c
@@ -52,7 +52,7 @@ void general_stokes_solver_setup(struct All_variables *E)
     construct_node_maps(E);
   else
     for (i=E->mesh.gridmin;i<=E->mesh.gridmax;i++)
-      E->elt_k[i][CPPR]=(struct EK *)malloc((E->lmesh.NEL[i]+1)*sizeof(struct EK));
+      E->elt_k[i]=(struct EK *)malloc((E->lmesh.NEL[i]+1)*sizeof(struct EK));
 
 #ifdef USE_PETSC
   /*
diff --git a/lib/Element_calculations.c b/lib/Element_calculations.c
index 6447ea5..8ecf759 100644
--- a/lib/Element_calculations.c
+++ b/lib/Element_calculations.c
@@ -462,27 +462,27 @@ void e_assemble_del2_u(E,u,Au,level,strip_bcs)
 	  ii = (a*n+b)*dims-(dims*n+dims);
 	  /* i=1, j=1,2,3 */
 		Au[a1] +=
-		        E->elt_k[level][CPPR][e].k[ii] *
+		        E->elt_k[level][e].k[ii] *
 			u[E->ID[level][nodeb].doff[1]]
-		      + E->elt_k[level][CPPR][e].k[ii+1] *
+		      + E->elt_k[level][e].k[ii+1] *
 			u[E->ID[level][nodeb].doff[2]]
-		      + E->elt_k[level][CPPR][e].k[ii+2] *
+		      + E->elt_k[level][e].k[ii+2] *
 			u[E->ID[level][nodeb].doff[3]];
 		/* i=2, j=1,2,3 */
 		Au[a2] +=
-		        E->elt_k[level][CPPR][e].k[ii+n] *
+		        E->elt_k[level][e].k[ii+n] *
 			u[E->ID[level][nodeb].doff[1]]
-		      + E->elt_k[level][CPPR][e].k[ii+n+1] *
+		      + E->elt_k[level][e].k[ii+n+1] *
 			u[E->ID[level][nodeb].doff[2]]
-		      + E->elt_k[level][CPPR][e].k[ii+n+2] *
+		      + E->elt_k[level][e].k[ii+n+2] *
 			u[E->ID[level][nodeb].doff[3]];
 		/* i=3, j=1,2,3 */
 		Au[a3] +=
-		        E->elt_k[level][CPPR][e].k[ii+n+n] *
+		        E->elt_k[level][e].k[ii+n+n] *
 			u[E->ID[level][nodeb].doff[1]]
-		      + E->elt_k[level][CPPR][e].k[ii+n+n+1] *
+		      + E->elt_k[level][e].k[ii+n+n+1] *
 			u[E->ID[level][nodeb].doff[2]]
-		      + E->elt_k[level][CPPR][e].k[ii+n+n+2] *
+		      + E->elt_k[level][e].k[ii+n+n+2] *
 			u[E->ID[level][nodeb].doff[3]];
 
  	    }         /* end for loop b */
diff --git a/lib/global_defs.h b/lib/global_defs.h
index eb754db..d00f766 100644
--- a/lib/global_defs.h
+++ b/lib/global_defs.h
@@ -807,7 +807,7 @@ struct All_variables {
     struct SUBEL *EL[MAX_LEVELS];
     struct EG *elt_del[MAX_LEVELS];
     struct EC *elt_c[MAX_LEVELS];
-    struct EK *elt_k[MAX_LEVELS][NCS];
+    struct EK *elt_k[MAX_LEVELS];
     struct CC *cc;
     struct CCX *ccx;
     struct CC *CC[MAX_LEVELS][NCS];



More information about the CIG-COMMITS mailing list