[cig-commits] r14207 - in long/3D/Gale/trunk: . src/Gale/Utils/src

[walter at geodynamics.org]

Author: walter
Date: 2009-03-03 12:26:48 -0800 (Tue, 03 Mar 2009)
New Revision: 14207

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c
Log:
 r2518 at dante:  boo | 2009-03-03 12:26:38 -0800
 Remove overcount stuff from StressBC because we are not actually overcounting



Property changes on: long/3D/Gale/trunk
___________________________________________________________________
Name: svk:merge
   - 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2516
   + 3a629746-de10-0410-b17b-fd6ecaaa963e:/cig:2518

Modified: long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c
===================================================================
--- long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c	2009-03-03 19:53:11 UTC (rev 14206)
+++ long/3D/Gale/trunk/src/Gale/Utils/src/StressBC.c	2009-03-03 20:26:48 UTC (rev 14207)
@@ -385,7 +385,6 @@
     
     if(condition)
       {
-        unsigned overcount=StressBC_get_overcount(dim,ijk,grid->sizes);
         for(entry_I=0; entry_I<self->numEntries; ++entry_I)
           {
             switch(self->_entryTbl[entry_I].type)
@@ -409,12 +408,8 @@
                                            Mesh_GetVertex(mesh,elementNodes[eNode_I]));
                 break;
               }
-            /* We have to divide by an overcount_factor, because
-               otherwise different elements will count the same node
-               more than once.  In 2D, nodes are counted twice, in 3D,
-               nodes are counted four times. */
             elForceVec[eNode_I*nodeDofCount + self->_entryTbl[entry_I].axis]+=
-              stress*area/overcount;
+              stress*area;
           }
       }
   }
@@ -532,19 +527,3 @@
     }
 }
 
-/* We have to divide by an overcount_factor, because otherwise
-   different elements will count the same node more than once.  In 2D,
-   in the interior, there are 4 elements that touch a node. On an
-   edge, 2, and on the corner 1.  Similary for 3D. */
-unsigned StressBC_get_overcount(Dimension_Index dim, IJK ijk, unsigned sizes[])
-{
-  unsigned overcount, i;
-  if(dim==2)
-    overcount=4;
-  else
-    overcount=8;
-  for(i=0;i<dim;++i)
-    if(ijk[i]==0 || ijk[i]==sizes[i]-1)
-      overcount/=2;
-  return overcount;
-}