Main Page   Namespace List   Class Hierarchy   Alphabetical List   Compound List   File List   Compound Members   File Members  

QuadCFInterp.H

Go to the documentation of this file.
00001 /*   _______              __
00002     / ___/ /  ___  __ _  / /  ___
00003    / /__/ _ \/ _ \/  V \/ _ \/ _ \
00004    \___/_//_/\___/_/_/_/_.__/\___/
00005 */
00006 // CHOMBO Copyright (c) 2000-2004, The Regents of the University of
00007 // California, through Lawrence Berkeley National Laboratory (subject to
00008 // receipt of any required approvals from U.S. Dept. of Energy).  All
00009 // rights reserved.
00010 //
00011 // Redistribution and use in source and binary forms, with or without
00012 // modification, are permitted provided that the following conditions are met:
00013 //
00014 // (1) Redistributions of source code must retain the above copyright
00015 // notice, this list of conditions and the following disclaimer.
00016 // (2) Redistributions in binary form must reproduce the above copyright
00017 // notice, this list of conditions and the following disclaimer in the
00018 // documentation and/or other materials provided with the distribution.
00019 // (3) Neither the name of Lawrence Berkeley National Laboratory, U.S.
00020 // Dept. of Energy nor the names of its contributors may be used to endorse
00021 // or promote products derived from this software without specific prior
00022 // written permission.
00023 //
00024 // THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
00025 // "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
00026 // TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
00027 // PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER
00028 // OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
00029 // EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
00030 // PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
00031 // PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
00032 // LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
00033 // NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
00034 // SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
00035 //
00036 // You are under no obligation whatsoever to provide any bug fixes,
00037 // patches, or upgrades to the features, functionality or performance of
00038 // the source code ("Enhancements") to anyone; however, if you choose to
00039 // make your Enhancements available either publicly, or directly to
00040 // Lawrence Berkeley National Laboratory, without imposing a separate
00041 // written license agreement for such Enhancements, then you hereby grant
00042 // the following license: a non-exclusive, royalty-free perpetual license
00043 // to install, use, modify, prepare derivative works, incorporate into
00044 // other computer software, distribute, and sublicense such Enhancements or
00045 // derivative works thereof, in binary and source code form.
00046 //
00047 // TRADEMARKS. Product and company names mentioned herein may be the
00048 // trademarks of their respective owners.  Any rights not expressly granted
00049 // herein are reserved.
00050 //
00051 
00052 #ifndef _QUADCFINTERP_H_
00053 #define _QUADCFINTERP_H_
00054 
00055 #include <iostream>
00056 #include <cmath>
00057 #include <cassert>
00058 #include <cstdlib>
00059 #include "REAL.H"
00060 #include "IntVect.H"
00061 #include "Box.H"
00062 #include "BaseFab.H"
00063 #include "DisjointBoxLayout.H"
00064 #include "LevelData.H"
00065 #include "FArrayBox.H"
00066 #include "QuadCFStencil.H"
00067 #include "ProblemDomain.H"
00068 
00070 
00074 class QuadCFInterp
00075 {
00076 public:
00077 
00079 
00081 
00086   QuadCFInterp(
00087                const DisjointBoxLayout& a_fineBoxes,
00088                const DisjointBoxLayout* a_coarBoxes,
00089                Real  a_dxFine,
00090                int a_refRatio,
00091                int a_nComp,
00092                const Box& a_domf);
00093 
00095 
00100   QuadCFInterp(
00101                const DisjointBoxLayout& a_fineBoxes,
00102                const DisjointBoxLayout* a_coarBoxes,
00103                Real  a_dxFine,
00104                int a_refRatio,
00105                int a_nComp,
00106                const ProblemDomain& a_domf);
00107 
00109 
00113   void define(
00114               const DisjointBoxLayout& a_fineBoxes,
00115               const DisjointBoxLayout* a_coarBoxes,
00116               Real  a_dxFine,
00117               int a_refRatio,
00118               int a_nComp,
00119               const ProblemDomain& a_domf);
00120 
00122 
00125   QuadCFInterp();
00126 
00128   ~QuadCFInterp();
00129 
00131 
00134   void clear();
00135 
00137 
00140   void coarseFineInterp(
00141                         LevelData<FArrayBox>& a_phif,
00142                         const LevelData<FArrayBox>& a_phic
00143                         );
00144 
00146 
00149   bool isDefined() const;
00150 
00151 
00152 protected:
00153   //boxes for coarse buffers
00154   BoxLayout m_loCoarBoxes[SpaceDim];
00155   BoxLayout m_hiCoarBoxes[SpaceDim];
00156 
00157   //fake level identifier.
00158   int m_level;
00159 
00160   //number of components
00161   int m_nComp;
00162 
00163   //refinement ratio between levels
00164   int m_refRatio;
00165 
00166   //mesh spacing at fine level
00167   Real m_dxFine;
00168 
00169   //cfstencils on low-side faces
00170   //lives on fine input grid layout
00171   LayoutData<QuadCFStencil> m_loQCFS[SpaceDim];
00172 
00173   //cfstencils on high-side faces
00174   //lives on fine input grid layout
00175   LayoutData<QuadCFStencil> m_hiQCFS[SpaceDim];
00176 
00177   //buffer for coarse data at low edges of fine grids
00178   BoxLayoutData<FArrayBox> m_hiCoarBuffer[SpaceDim];
00179 
00180   // copiers to copy data from inputCoarLayout to m_hiCoarBuffer
00181   Copier m_hiCoarCopiers[SpaceDim];
00182 
00183   //buffer for coarse data at high edges of fine grids
00184   BoxLayoutData<FArrayBox> m_loCoarBuffer[SpaceDim];
00185 
00186   // copiers to copy data from inputCoarLayout to m_loCoarBuffer
00187   Copier m_loCoarCopiers[SpaceDim];
00188 
00189   DisjointBoxLayout m_inputFineLayout;
00190   DisjointBoxLayout m_inputCoarLayout;
00192   bool m_isDefined;
00193 
00194   ProblemDomain m_domainFine;
00195 protected:
00196 
00197   /* Coarse / Fine interpolation operator.
00198      using one sided differencing
00199      only at a particular face
00200   */
00201   void coarseFineInterp(BaseFab<Real> & a_phif,
00202                         const BaseFab<Real> & a_phic,
00203                         const QuadCFStencil& a_qcfs,
00204                         const Side::LoHiSide a_hiorlo,
00205                         const int a_idir,
00206                         const Interval& a_variables) const;
00207 
00208   //get extended phi (lives next to interpivs)
00209   void getPhiStar(BaseFab<Real> & a_phistar,
00210                   const BaseFab<Real> & a_phic,
00211                   const QuadCFStencil& a_qcfs,
00212                   const Side::LoHiSide a_hiorlo,
00213                   const int a_idir,
00214                   const Interval& a_variables) const;
00215 
00216   //interpolate over correct intvectset
00217   void interpOnIVS(BaseFab<Real> & a_phif,
00218                    const BaseFab<Real> & a_phiStar,
00219                    const QuadCFStencil& a_qcfs,
00220                    const Side::LoHiSide a_hiorlo,
00221                    const int a_idir,
00222                    const Interval& a_variables) const;
00223 };
00224 
00225 #endif

Generated on Wed Jan 19 17:51:27 2005 for Chombo&INSwithParticles by doxygen1.2.16