LevelNodeSolver.H

Go to the documentation of this file.

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

// LevelNodeSolver.H
// adapted from LevelSolver by DTGraves, Fri, July 23, 1999
// petermc, 12 Apr 2001

#ifndef LEVELNODESOLVER_H
#define LEVELNODESOLVER_H

#include <cstdlib>
#include <iostream>
#include "SPACE.H"
#include <cmath>
#include "REAL.H"
#include "IntVectSet.H"
#include "ProblemDomain.H"
#include "NodeFArrayBox.H"
#include "DisjointBoxLayout.H"
#include "LevelData.H"
#include "NodeLevelMG.H"
#include "NodeLevelOp.H"

#include "UsingNamespace.H"

class LevelNodeSolver
{

public:



  LevelNodeSolver();


  LevelNodeSolver(const DisjointBoxLayout& a_grids,
                  const DisjointBoxLayout* a_gridsCoarsePtr,
                  const ProblemDomain&     a_domain,
                  Real                     a_dx,
                  int                      a_refToCoarse,
                  const NodeLevelOp* const a_opin,
                  int                      a_minLength = 1);


  LevelNodeSolver(const DisjointBoxLayout& a_grids,
                  const DisjointBoxLayout* a_gridsCoarsePtr,
                  const Box&               a_domain,
                  Real                     a_dx,
                  int                      a_refToCoarse,
                  const NodeLevelOp* const a_opin,
                  int                      a_minLength = 1);


  virtual ~LevelNodeSolver();


  virtual void define(
                      const DisjointBoxLayout& a_grids,
                      const DisjointBoxLayout* a_gridsCoarsePtr,
                      const ProblemDomain&     a_domain,
                      Real                     a_dx,
                      int                      a_refToCoarse,
                      const NodeLevelOp* const a_opin,
                      int                      a_minLength = 1);


  virtual void define(
                      const DisjointBoxLayout& a_grids,
                      const DisjointBoxLayout* a_gridsCoarsePtr,
                      const Box&               a_domain,
                      Real                     a_dx,
                      int                      a_refToCoarse,
                      const NodeLevelOp* const a_opin,
                      int                      a_minLength = 1);



  bool isDefined() const;




  void setnumSmoothUp(int a_numSmoothUp);


  void setnumSmoothDown(int a_numSmoothDown);


  void setBottomSmoothing(bool a_bottomSolveFlag);


  void setTolerance(Real a_tolerance);


  void setOperatorTolerance(Real a_operatorTolerance);


  void setVerbose(bool a_verbose);


  void setMaxIter(int a_maxIter);


  void setMinIter(int a_minIter);


  void clearMemory();

  void setDefaultValues();


  void levelSolveH(
                   LevelData<NodeFArrayBox>&       a_phi,
                   const LevelData<NodeFArrayBox>& a_rhs,
                   bool                            a_initializePhiToZero = true);


  void levelSolve(
                  LevelData<NodeFArrayBox>&       a_phi,
                  const LevelData<NodeFArrayBox>* a_phiCoarse,
                  const LevelData<NodeFArrayBox>& a_rhs,
                  bool                            a_initializePhiToZero = true);

protected:

  NodeLevelMG m_levelMG;

  NodeLevelOp* m_levelOpPtr;

  bool m_isDefined;

  DisjointBoxLayout m_grids;

  ProblemDomain m_domain;

  LevelData<NodeFArrayBox> m_resid;

  LevelData<NodeFArrayBox> m_scratch;

  LevelData<NodeFArrayBox> m_corr;

  int m_refToCoarse;

  Real m_dx;

  int m_maxIter;

  int m_minIter;

  Real m_tolerance;

  Real m_operatorTolerance;

  bool m_bottomSolveFlag;

  bool m_verbose;

  LayoutData< Vector<IntVectSet> > m_IVSV;

  LayoutData< Vector<IntVectSet> > m_IVSVext;

private:

  int countCoarserLevels(int a_minLength);

  void operator = (const LevelNodeSolver&){;}
  LevelNodeSolver(const LevelNodeSolver&){;}

};

#endif

---