NodeQuadCFInterp2.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

// NodeQuadCFInterp2.H
// original NodeQuadCFInterp by petermc, Mon, Apr 23, 2001
// petermc, 31 Oct 2001
#ifndef _NODEQUADCFINTERP2_H_
#define _NODEQUADCFINTERP2_H_

#include <cmath>
#include <cstdlib>
#include "BaseFab.H"
#include "NodeFArrayBox.H"
#include "LevelData.H"
#include "NodeCFIVS.H"
#include "NodeBC.H"
#include "NamespaceHeader.H"


class NodeQuadCFInterp2
{
public:



  NodeQuadCFInterp2();


  NodeQuadCFInterp2(const DisjointBoxLayout& a_grids,
                    const Box& a_domain,
                    const LayoutData<NodeCFIVS>* const a_loCFIVS,
                    const LayoutData<NodeCFIVS>* const a_hiCFIVS,
                    bool a_interfaceOnly,
                    int a_interpolationDegree,
                    int a_ncomp = 1,
                    bool a_verbose = false);


  NodeQuadCFInterp2(const DisjointBoxLayout& a_grids,
                    const ProblemDomain& a_domain,
                    const LayoutData<NodeCFIVS>* const a_loCFIVS,
                    const LayoutData<NodeCFIVS>* const a_hiCFIVS,
                    bool a_interfaceOnly,
                    int a_interpolationDegree,
                    int a_ncomp = 1,
                    bool a_verbose = false);


  ~NodeQuadCFInterp2();


  void define(
              const DisjointBoxLayout& a_grids,
              const ProblemDomain& a_domain,
              const LayoutData<NodeCFIVS>* const a_loCFIVS,
              const LayoutData<NodeCFIVS>* const a_hiCFIVS,
              bool a_interfaceOnly,
              int a_interpolationDegree,
              int a_ncomp = 1,
              bool a_verbose = false);


  void define(
              const DisjointBoxLayout& a_grids,
              const Box& a_domain,
              const LayoutData<NodeCFIVS>* const a_loCFIVS,
              const LayoutData<NodeCFIVS>* const a_hiCFIVS,
              bool a_interfaceOnly,
              int a_interpolationDegree,
              int a_ncomp = 1,
              bool a_verbose = false);




  bool isDefined() const;



  void setVerbose( bool a_verbose );



  void coarseFineInterp(
                        LevelData<NodeFArrayBox>& a_phi,
                        const LevelData<NodeFArrayBox>& a_phiCoarse);

protected:

  DisjointBoxLayout m_grids;

  ProblemDomain m_domain;

  Box m_domainCoarseNodes;

  DisjointBoxLayout m_coarsenedGrids;

  LevelData<NodeFArrayBox> m_coarseCopy;

  bool m_interfaceOnly;

  int m_interpolationDegree;

  int m_ncomp;

  bool m_isDefined;

  const LayoutData<NodeCFIVS>* m_loCFIVS;

  const LayoutData<NodeCFIVS>* m_hiCFIVS;
  // LayoutData<NodeCFIVS> m_loCFIVS[SpaceDim];
  // LayoutData<NodeCFIVS> m_hiCFIVS[SpaceDim];

  // weights for coarse/fine interpolation
  // Tuple<FArrayBox*, SpaceDim> m_wtLo, m_wtHi, m_wtC;

  bool m_verbose;

private:
  //internally useful functions

  void clearMemory();

  void setDefaultValues();

  void interpLine(
                  FArrayBox& a_fineFab,
                  const FArrayBox& a_psiFab,
                  const IntVectSet& a_psiFabNodes,
                  const IntVect& a_iv,
                  int a_idirOther);

};

#include "NamespaceFooter.H"
#endif

---