---

EBArith.H

Go to the documentation of this file.

/* _______              __
  / ___/ /  ___  __ _  / /  ___
 / /__/ _ \/ _ \/  ' \/ _ \/ _ \
 \___/_//_/\___/_/_/_/_.__/\___/ 
*/
//
// This software is copyright (C) by the Lawrence Berkeley
// National Laboratory.  Permission is granted to reproduce
// this software for non-commercial purposes provided that
// this notice is left intact.
// 
// It is acknowledged that the U.S. Government has rights to
// this software under Contract DE-AC03-765F00098 between
// the U.S.  Department of Energy and the University of
// California.
//
// This software is provided as a professional and academic
// contribution for joint exchange. Thus it is experimental,
// is provided ``as is'', with no warranties of any kind
// whatsoever, no support, no promise of updates, or printed
// documentation. By using this software, you acknowledge
// that the Lawrence Berkeley National Laboratory and
// Regents of the University of California shall have no
// liability with respect to the infringement of other
// copyrights by any part of this software.
//

#ifndef _EB_ARITH_H_
#define _EB_ARITH_H_

#include "EBCellFAB.H"
#include "EBISBox.H"
#include "EBISLayout.H"
#include "EBFaceFAB.H"
#include "BoxLayout.H"
#include "BoxLayoutData.H"
#include "Stencils.H"
#include "LayoutData.H"
#include "Tuple.H"
#include "EBFluxFAB.H"
#include "EBCellFAB.H"

class EBNormType
{
public:

  enum  NormMode{OverOnlyRegular=0, OverOnlyIrregular, OverBoth};
};

class EBArith
{
public:

  static
  bool
  getAdjacentFace(FaceIndex& a_adjacentFace,
                  const FaceIndex& a_face,
                  const EBISBox& a_ebisBox, 
                  const Box& a_domain,
                  const int& a_idir,
                  const Side::LoHiSide& a_side);




  static
  void computeGradFluxStencil(VoFStencil& a_thisStencil,
                              const FaceIndex& a_thisFace,
                              const EBISBox& a_ebisBox,
                              const Box& a_domain,
                              const int& a_dir);


  static
  void computeInterpStencil(FaceStencil& a_thisStencil,
                            const FaceIndex& a_thisFace,
                            const EBISBox& a_ebisBox,
                            const Box& a_domain,
                            const int& a_dir);


  static
  void
  volWeightedSum(Real& norm, Real& volume,
                 const EBCellFAB& a_src, 
                 const Box& a_region,
                 const EBISBox& a_ebisBox,
                 const int& comp,
                 const int& p,
                 EBNormType::NormMode = EBNormType::OverBoth);



  static
  Real 
  norm(const BoxLayoutData<EBCellFAB >& a_dataOne, 
       const BoxLayout& a_layout, 
       const EBISLayout& a_ebisl, 
       const int& comp,
       const int&  p,
       EBNormType::NormMode = EBNormType::OverBoth);

  //testing interface
  static
  Real 
  norm(Real& volume, const BoxLayoutData<EBCellFAB >& a_dataOne, 
       const BoxLayout& a_layout, 
       const EBISLayout& a_ebisl, 
       const int& comp,
       const int&  p,
       EBNormType::NormMode = EBNormType::OverBoth);


  static
  Real 
  norm(const EBCellFAB& a_dataOne, 
       const Box&       a_grid, 
       const EBISBox&   a_ebisl, 
       const int&       a_comp,
       const int&       a_p,
       EBNormType::NormMode = EBNormType::OverBoth);


  static
  Real 
  norm(Real& volume, const EBCellFAB& a_dataOne, 
       const Box&       a_grid, 
       const EBISBox&   a_ebisl, 
       const int&       a_comp,
       const int&       a_p,
       EBNormType::NormMode = EBNormType::OverBoth);


  static
  void
  volWeightedSum(Real& sum, Real& volume, 
                 const BoxLayoutData<EBCellFAB >& a_dataOne, 
                 const BoxLayout& a_layout, 
                 const EBISLayout& a_ebisl, 
                 const int& comp,
                 const int& p,
                 EBNormType::NormMode = EBNormType::OverBoth);


  static
  Real 
  dotProduct(const BoxLayoutData<EBCellFAB >& a_dataOne, 
             const BoxLayoutData<EBCellFAB >& a_dataTwo,
             const BoxLayout& a_layout, 
             const EBISLayout& a_ebisl, 
             const int& a_comp);


  static
  Real 
  dotProduct(const EBCellFAB& a_dataOne, 
             const EBCellFAB& a_dataTwo,
             const Box& a_layout, 
             const EBISBox& a_ebisBox,
             const int& a_comp);



  static
  bool
  monotonePathVoFToCellVoF(VolIndex& a_vof2,
                           const VolIndex& a_vof1,
                           const IntVect& a_cell2,
                           const EBISBox& a_ebisBox);

  static
  void
  getAllVoFsInMonotonePath(Vector<VolIndex>& a_vofs,
                           const IntVect&    a_timesMoved,
                           const IntVect&    a_pathSign,
                           const VolIndex&   a_vof,
                           const EBISBox&    a_ebisBox,
                           const int&        a_redistRad);
  

private:

};

#endif

---