[cig-commits] r12013 - in long/3D/Gale/trunk: . src/StgFEM/SLE/SystemSetup/src

walter at geodynamics.org walter at geodynamics.org
Fri May 23 14:31:48 PDT 2008


Author: walter
Date: 2008-05-23 14:31:48 -0700 (Fri, 23 May 2008)
New Revision: 12013

Modified:
   long/3D/Gale/trunk/
   long/3D/Gale/trunk/src/StgFEM/SLE/SystemSetup/src/SystemLinearEquations.c
Log:
 r2182 at earth:  boo | 2008-05-23 14:10:18 -0700
 No longer set the pressure to zero after each non-linear step



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

Modified: long/3D/Gale/trunk/src/StgFEM/SLE/SystemSetup/src/SystemLinearEquations.c
===================================================================
--- long/3D/Gale/trunk/src/StgFEM/SLE/SystemSetup/src/SystemLinearEquations.c	2008-05-23 21:31:43 UTC (rev 12012)
+++ long/3D/Gale/trunk/src/StgFEM/SLE/SystemSetup/src/SystemLinearEquations.c	2008-05-23 21:31:48 UTC (rev 12013)
@@ -675,25 +675,6 @@
         p/=nInc;
         FeVariable_SetValueAtNode(nodal_pressure,n_i,&p);
       }
-
-      /* Average the nodal pressures to get a smoothed element pressure */
-      N = FeMesh_GetElementLocalSize( pressure->feMesh);
-      for( n_i = 0; n_i < N; n_i++ ) {
-        unsigned nInc;
-        unsigned *inc;
-        unsigned nn;
-        double p,temp;
-        
-/*         Mesh_GetIncidence(pressure->feMesh, nDims, n_i, MT_VERTEX, &nInc, &inc); */
-        p=0;
-/*         for(nn=0; nn<nInc; ++nn) */
-/*           { */
-/*             FeVariable_GetValueAtNode(nodal_pressure,inc[nn],&temp); */
-/*             p+=temp; */
-/*           } */
-/*         p/=nInc; */
-        FeVariable_SetValueAtNode(pressure,n_i,&p);
-      }
     }
 }
 



More information about the cig-commits mailing list