[cig-commits] [commit] rajesh-petsc-schur: cleaned up n_assemble_del2_u (1c7d3b8)

[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 1c7d3b8a697fe031e8c25bf5600f31bbc142dcac
Author: Rajesh Kommu <rajesh.kommu at gmail.com>
Date:   Mon Sep 8 15:19:51 2014 -0700

    cleaned up n_assemble_del2_u


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

1c7d3b8a697fe031e8c25bf5600f31bbc142dcac
 lib/Element_calculations.c | 14 ++++++--------
 1 file changed, 6 insertions(+), 8 deletions(-)

diff --git a/lib/Element_calculations.c b/lib/Element_calculations.c
index 59fbc54..cbe26ea 100644
--- a/lib/Element_calculations.c
+++ b/lib/Element_calculations.c
@@ -430,6 +430,7 @@ void assemble_del2_u(E,u,Au,level,strip_bcs)
      int level;
      int strip_bcs;
 {
+
   void e_assemble_del2_u();
   void n_assemble_del2_u();
 
@@ -472,10 +473,9 @@ void e_assemble_del2_u(E,u,Au,level,strip_bcs)
 	a2 = E->ID[level][m][ii].doff[2];
 	a3 = E->ID[level][m][ii].doff[3];
 	for(b=1;b<=ends;b++) {
-	        nodeb = E->IEN[level][m][e].node[b];
-	        ii = (a*n+b)*dims-(dims*n+dims);
-			/* i=1, j=1,2 */
-	          /* i=1, j=1,2,3 */
+	  nodeb = E->IEN[level][m][e].node[b];
+	  ii = (a*n+b)*dims-(dims*n+dims);
+	  /* i=1, j=1,2,3 */
 		Au[m][a1] +=
 		        E->elt_k[level][m][e].k[ii] *
 			u[m][E->ID[level][m][nodeb].doff[1]]
@@ -541,10 +541,8 @@ void n_assemble_del2_u(E,u,Au,level,strip_bcs)
 
   for (m=1;m<=E->sphere.caps_per_proc;m++)  {
 
-     for(e=0;e<=neq;e++)
-	Au[m][e]=0.0;
-
-     u[m][neq] = 0.0;
+     for(e=0;e<neq;e++)
+      Au[m][e]=0.0;
 
      for(e=1;e<=nno;e++)     {