[cig-commits] commit: Change some viscosities to moduli

[Mercurial]

changeset:   268:606efdc34aca
user:        Walter Landry <wlandry at caltech.edu>
date:        Thu May 03 10:40:31 2012 -0700
files:       src/StokesFACOps.h src/StokesFACOps/smooth_V_3D.C
description:
Change some viscosities to moduli


diff -r 0ba9f7511b8a -r 606efdc34aca src/StokesFACOps.h
--- a/src/StokesFACOps.h	Tue May 01 15:16:42 2012 -0700
+++ b/src/StokesFACOps.h	Thu May 03 10:40:31 2012 -0700
@@ -520,8 +520,8 @@ void smooth_V_3D(const int &ix,
   tbox::Pointer<geom::CartesianPatchGeometry> &geom,
   pdat::SideData<double> &v,
   pdat::SideData<double> &v_rhs,
-  pdat::CellData<double> &cell_viscosity,
-  pdat::EdgeData<double> &edge_viscosity,
+  pdat::CellData<double> &cell_moduli,
+  pdat::EdgeData<double> &edge_moduli,
   const pdat::CellIndex &center,
   const double Dx[3],
   const double &theta_momentum,
diff -r 0ba9f7511b8a -r 606efdc34aca src/StokesFACOps/smooth_V_3D.C
--- a/src/StokesFACOps/smooth_V_3D.C	Tue May 01 15:16:42 2012 -0700
+++ b/src/StokesFACOps/smooth_V_3D.C	Thu May 03 10:40:31 2012 -0700
@@ -13,8 +13,8 @@ void SAMRAI::solv::StokesFACOps::smooth_
  tbox::Pointer<geom::CartesianPatchGeometry> &geom,
  pdat::SideData<double> &v,
  pdat::SideData<double> &v_rhs,
- pdat::CellData<double> &cell_viscosity,
- pdat::EdgeData<double> &edge_viscosity,
+ pdat::CellData<double> &cell_moduli,
+ pdat::EdgeData<double> &edge_moduli,
  const pdat::CellIndex &center,
  const double Dx[3],
  const double &theta_momentum,
@@ -50,12 +50,12 @@ void SAMRAI::solv::StokesFACOps::smooth_
           dv_upper=v(x+offset) - v(x);
         }
 
-      double C_vx=dRm_dv_3D(cell_viscosity,edge_viscosity,center,center-pp[ix],
+      double C_vx=dRm_dv_3D(cell_moduli,edge_moduli,center,center-pp[ix],
                             edge_y+pp[iz],edge_y,edge_z+pp[iy],edge_z,
                             Dx[ix],Dx[iy],Dx[iz]);
 
       double delta_Rx=v_rhs(x)
-        - v_operator_3D(v,cell_viscosity,edge_viscosity,center,edge_y,edge_z,
+        - v_operator_3D(v,cell_moduli,edge_moduli,center,edge_y,edge_z,
                         x,y,z,pp[ix],pp[iy],pp[iz],Dx[ix],Dx[iy],Dx[iz]);
 
       /* No scaling here, though there should be. */