[cig-commits] commit: Remove a left over file

Mercurial hg at geodynamics.org
Fri Apr 15 16:47:04 PDT 2011


changeset:   161:799db2fc597b
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Fri Apr 15 16:46:02 2011 -0700
files:       V_Coarsen_Patch_Strategy.C
description:
Remove a left over file


diff -r 95b0f77ccfe6 -r 799db2fc597b V_Coarsen_Patch_Strategy.C
--- a/V_Coarsen_Patch_Strategy.C	Fri Apr 15 16:45:04 2011 -0700
+++ /dev/null	Thu Jan 01 00:00:00 1970 +0000
@@ -1,114 +0,0 @@
-#include "V_Coarsen_Patch_Strategy.h"
-
-void
-SAMRAI::solv::V_Coarsen_Patch_Strategy::postprocessCoarsen
-(hier::Patch& coarse,
- const hier::Patch& fine,
- const hier::Box& ,
- const hier::IntVector& )
-{
-  /* Fix up the boundary elements by iterating through the boundary
-     boxes */
-
-  /* We only care about edges, not corners, so we only iterate over
-     edge boundary boxes. */
-  const tbox::Array<hier::BoundaryBox>
-    &boundaries=coarse_fine[fine.getPatchLevelNumber()]->getEdgeBoundaries(coarse.getGlobalId());
-     
-  tbox::Pointer<pdat::SideData<double> >
-    v_fine = fine.getPatchData(v_id);
-  tbox::Pointer<pdat::SideData<double> >
-    v = coarse.getPatchData(v_id);
-
-  TBOX_ASSERT(!v.isNull());
-  TBOX_ASSERT(!v_fine.isNull());
-  TBOX_ASSERT(v_fine->getDepth() == v->getDepth());
-  TBOX_ASSERT(v->getDepth() == 1);
-
-  hier::Box gbox(v_fine->getGhostBox());
-  hier::Index ip(1,0), jp(0,1);
-  for(int mm=0; mm<boundaries.size(); ++mm)
-    {
-      hier::Box bbox=boundaries[mm].getBox();
-      int location_index=boundaries[mm].getLocationIndex();
-
-      // hier::Index unit(1,1);
-      hier::Index lower=hier::Index::coarsen(bbox.lower(),hier::Index(2,2)),
-        upper=hier::Index::coarsen(bbox.upper(),hier::Index(2,2));
-
-      // tbox::plog << "VCPS "
-      //            << mm << " "
-      //            << bbox.lower(0) << " "
-      //            << bbox.upper(0) << " "
-      //            << bbox.lower(1) << " "
-      //            << bbox.upper(1) << " "
-      //            << lower(0) << " "
-      //            << upper(0) << " "
-      //            << lower(1) << " "
-      //            << upper(1) << " "
-      //            << v_fine->getGhostBox().lower(0) << " "
-      //            << v_fine->getGhostBox().upper(0) << " "
-      //            << v_fine->getGhostBox().lower(1) << " "
-      //            << v_fine->getGhostBox().upper(1) << " "
-      //            << "\n";
-
-      for(int j=lower(1); j<=upper(1); ++j)
-        for(int i=lower(0); i<=upper(0); ++i)
-          {
-            // tbox::plog << "V Patch "
-            //            << i << " "
-            //            << j << " ";
-            /* Fix vx */
-            if(location_index==0)
-              {
-                pdat::SideIndex coarse(hier::Index(i,j),0,
-                                       pdat::SideIndex::Upper);
-                pdat::SideIndex center(coarse*2);
-                if(center[1]>=gbox.lower(1) && center[1]<gbox.upper(1))
-                  {
-                    (*v)(coarse)=((*v_fine)(center) + (*v_fine)(center+jp))/2;
-                    // tbox::plog << (*v)(coarse) << " ";
-                  }
-              }
-            else if(location_index==1)
-              {
-                pdat::SideIndex coarse(hier::Index(i,j),0,
-                                       pdat::SideIndex::Lower);
-                pdat::SideIndex center(coarse*2);
-                if(center[1]>=gbox.lower(1) && center[1]<gbox.upper(1))
-                  {
-                    (*v)(coarse)=((*v_fine)(center) + (*v_fine)(center+jp))/2;
-                    // tbox::plog << (*v)(coarse) << " ";
-                  }
-              }
-            /* Fix vy */
-            else if(location_index==2)
-              {
-                pdat::SideIndex coarse(hier::Index(i,j),1,
-                                       pdat::SideIndex::Upper);
-                pdat::SideIndex center(coarse*2);
-                if(center[0]>=gbox.lower(0) && center[0]<gbox.upper(0))
-                  {
-                    (*v)(coarse)=((*v_fine)(center) + (*v_fine)(center+ip))/2;
-                    // tbox::plog << (*v)(coarse) << " ";
-                  }
-              }
-            else if(location_index==3)
-              {
-                pdat::SideIndex coarse(hier::Index(i,j),1,
-                                       pdat::SideIndex::Lower);
-                pdat::SideIndex center(coarse*2);
-                if(center[0]>=gbox.lower(0) && center[0]<gbox.upper(0))
-                  {
-                    (*v)(coarse)=((*v_fine)(center) + (*v_fine)(center+ip))/2;
-                    // tbox::plog << (*v)(coarse) << " ";
-                  }
-              }
-            else
-              {
-                abort();
-              }
-            // tbox::plog << "\n";
-          }
-    }
-}



More information about the CIG-COMMITS mailing list