---

NodeFArrayBox.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.
//

// NodeFArrayBox.H
// adapted from dMartin/Chombo-IAMR/util/FluxBox.H
// by Dan Martin, Fri, Jan 14, 2000
// petermc, 1 Nov 2000
// 8 Jan 2001 added 

#ifndef NODEFARRAYBOX_H
#define NODEFARRAYBOX_H

#include "ProblemDomain.H"
#include "FArrayBox.H"
#include "LevelData.H"
#include "Tuple.H"
// #include "ArrayViewData.H"
// #include "AMRIO.H"
// avoid including this because then InitWriteFAB gets multiply defined
// #include "FABView.H"


class NodeFArrayBox
{

public:



  NodeFArrayBox();


  NodeFArrayBox(const Box& a_bx, int a_nComp=1);


  ~NodeFArrayBox();


  void define(const Box& a_bx, int a_nComp=1);



  const Box& box() const;


  FArrayBox& getFab();


  const FArrayBox& getFab() const;


  Real norm(const Real a_dx,
            const Box& a_subbox, // CELL-centered
            const int a_p=2,
            const int a_startComp=0, const int a_numComp=1) const;
            

  Real norm(const Real a_dx,
            const BaseFab<int>& a_mask,
            const Box& a_subbox, // CELL-centered
            const int a_p=2,
            const int a_startComp=0, const int a_numComp=1) const;
            

  Real norm(const Real a_dx,
            const int a_p=2,
            const int a_startComp=0, const int a_numComp=1) const;


  Real norm(const Real a_dx,
            const BaseFab<int>& a_mask,
            const int a_p=2,
            const int a_startComp=0, const int a_numComp=1) const;


  Real maxnorm(const Real a_dx,
               const Box& a_subbox, // CELL-centered
               const int a_startComp=0, const int a_numComp=1) const;


  Real maxnorm(const BaseFab<int>& a_mask,
               const Box& a_subbox, // CELL-centered
               const int a_startComp=0, const int a_numComp=1) const;


  Real maxnorm(const Real a_dx,
               const int a_startComp=0, const int a_numComp=1) const;


  Real maxnorm(const BaseFab<int>& a_mask,
               const int a_startComp=0, const int a_numComp=1) const;



  void copy(const NodeFArrayBox& a_src);


  void copy(const Box& a_regionFrom,
            const Interval& a_Cdest,
            const Box& a_regionTo,
            const NodeFArrayBox& a_src,
            const Interval& a_Csrc);

  

  int size(const Box& a_R, const Interval& a_comps) const;


  void linearOut(void* a_buf, const Box& a_R, const Interval& a_comps) const;
  
  void linearIn(void* a_buf);

  void linearIn(void* a_buf, const Box& a_R, const Interval& a_comps);

  static bool preAllocatable() {return true;}
  
protected:

  // CELL-centered box. data are defined on the surrounding nodes.
  Box m_box;

  // NODE-centered data
  FArrayBox m_fab;

private:
  
  NodeFArrayBox (const NodeFArrayBox&);
  NodeFArrayBox& operator = (const NodeFArrayBox&);

};


void interiorNodes(IntVectSet& a_ivs,
                   const ProblemDomain& a_base_domain,
                   const DisjointBoxLayout& a_boxes,
                   const Box& a_box);

void interiorNodes(IntVectSet& a_ivs,
                   const Box& a_base_domain,
                   const DisjointBoxLayout& a_boxes,
                   const Box& a_box);


void interiorBoundaryNodes(LayoutData< Vector<IntVectSet> >& a_IVSV,
                           const DisjointBoxLayout& a_boxes,
                           const ProblemDomain& a_domain);

void interiorBoundaryNodes(LayoutData< Vector<IntVectSet> >& a_IVSV,
                           const DisjointBoxLayout& a_boxes,
                           const Box& a_domain);


void interiorBoundaryNodes(LayoutData< Vector<IntVectSet> >& a_IVSV,
                           const DisjointBoxLayout& a_dest,
                           const DisjointBoxLayout& a_src,
                           const ProblemDomain& a_domain);

void interiorBoundaryNodes(LayoutData< Vector<IntVectSet> >& a_IVSV,
                           const DisjointBoxLayout& a_dest,
                           const DisjointBoxLayout& a_src,
                           const Box& a_domain);


void exteriorBoundaryNodes(LayoutData< Vector<IntVectSet> >& a_exterior,
                           const LayoutData< Vector<IntVectSet> >& a_interior,
                           const DisjointBoxLayout& a_boxes);


void copyInteriorNodesOld(BoxLayoutData<NodeFArrayBox> & a_dest,
                          const LevelData<NodeFArrayBox> & a_src,
                          const Box& a_domain);


void zeroBoundaryNodes(BoxLayoutData<NodeFArrayBox>& a_dest,
                       const LayoutData< Vector<IntVectSet> >& a_IVSV);


void copyInteriorNodes(LevelData<NodeFArrayBox>& a_dest,
                       const LevelData<NodeFArrayBox>& a_src,
                       const LayoutData< Vector<IntVectSet> >& a_IVSV);


Real DotProductNodes(const BoxLayoutData<NodeFArrayBox>& a_dataOne,
                     const BoxLayoutData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const BoxLayout& a_dblIn);



Real DotProductNodes(const BoxLayoutData<NodeFArrayBox>& a_dataOne,
                     const BoxLayoutData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const BoxLayout& a_dblIn,
                     const Interval& a_comps);


Real DotProductNodes(const LevelData<NodeFArrayBox>& a_dataOne,
                     const LevelData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const ProblemDomain& a_domain,
                     const LayoutData< Vector<IntVectSet> >& a_IVSVext,
                     const Interval& a_comps);

Real DotProductNodes(const LevelData<NodeFArrayBox>& a_dataOne,
                     const LevelData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const Box& a_domain,
                     const LayoutData< Vector<IntVectSet> >& a_IVSVext,
                     const Interval& a_comps);



Real DotProductNodes(const LevelData<NodeFArrayBox>& a_dataOne,
                     const LevelData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const ProblemDomain& a_domain,
                     const Interval& a_comps);

Real DotProductNodes(const LevelData<NodeFArrayBox>& a_dataOne,
                     const LevelData<NodeFArrayBox>& a_dataTwo,
                     const BoxLayoutData< BaseFab<int> >& a_mask,
                     const Box& a_domain,
                     const Interval& a_comps);

#endif

---