NewPoissonOp.H

Go to the documentation of this file.

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

// BVS, June 18, 2003

// We can assume that template class T has null construction.

#ifndef _NEWPOISSONOP_H_
#define _NEWPOISSONOP_H_

#include "AMRMultiGrid.H"
#include "REAL.H"
#include "Box.H"
#include "LevelDataOps.H"
#include "BCFunc.H"
#include "FArrayBox.H"
#include "NamespaceHeader.H"

class NewPoissonOp : public MGLevelOp<FArrayBox>
{
public:
  NewPoissonOp(){ ; }

  virtual ~NewPoissonOp() { ; }

  void define(const RealVect& a_dx,
              const ProblemDomain& a_domain,
              BCFunc a_bc);

  virtual void createCoarsened(    FArrayBox& a_lhs,
                                   const FArrayBox& a_rhs,
                                   const int& a_refRat);

  virtual void residual(  FArrayBox& a_lhs,
                          const FArrayBox& a_phi,
                          const FArrayBox& a_rhs,
                          bool a_homogeneous = false);
  virtual void preCond(   FArrayBox& a_correction,
                          const FArrayBox& a_residual);
  virtual void applyOp(   FArrayBox& a_lhs,
                          const FArrayBox& a_phi,
                          bool a_homogeneous = false);
  virtual void create(    FArrayBox& a_lhs,
                          const FArrayBox& a_rhs);
  virtual void assign(    FArrayBox& a_lhs,
                          const FArrayBox& a_rhs) ;
  virtual Real dotProduct(const FArrayBox& a_1,
                          const FArrayBox& a_2) ;
  virtual void incr( FArrayBox& a_lhs,
                     const FArrayBox& a_x,
                     Real a_scale) ;
  virtual void axby( FArrayBox& a_lhs, const FArrayBox& a_x,
                     const FArrayBox& a_y,
                     Real a, Real b) ;
  virtual void scale(FArrayBox& a_lhs, const Real& a_scale) ;

  virtual Real norm(const FArrayBox& a_x, int a_ord);

  virtual void setToZero(FArrayBox& a_x);


  virtual void relax(FArrayBox& a_e,
                     const FArrayBox& a_residual,
                     int iterations);

  virtual void createCoarser(FArrayBox& a_coarse,
                             const FArrayBox& a_fine,
                             bool ghost);
  virtual void restrictResidual(FArrayBox& a_resCoarse,
                                FArrayBox& a_phiFine,
                                const FArrayBox& a_rhsFine);

  virtual void prolongIncrement(FArrayBox& a_phiThisLevel,
                                const FArrayBox& a_correctCoarse);

protected:
  RealVect                m_dx;
  RealVect                m_dxCrse;
  ProblemDomain           m_domain;
  BCFunc                  m_bc;

  void levelGSRB(FArrayBox& a_e,  const FArrayBox& a_residual);

  void colorGS(FArrayBox&       a_phi,
               const FArrayBox& a_rhs,
               const IntVect&   color);

  bool nextColor(IntVect& color,
                 const IntVect& limit);

  void levelRelaxColor(FArrayBox&       a_phi,
                       const FArrayBox& a_rhs,
                       const IntVect&   a_color,
                       const Real&      a_weight,
                       const bool&      a_homogeneousPhysBC);

};

class NewPoissonOpFactory : public MGLevelOpFactory<FArrayBox>
{
public:

  NewPoissonOpFactory();
  NewPoissonOpFactory(RealVect& a_dx, BCFunc a_bc);
  void define(const RealVect& m_dx,BCFunc m_bc);
  virtual ~NewPoissonOpFactory(){;}
  virtual NewPoissonOp* MGnewOp(const ProblemDomain& a_FineindexSpace,
                                   int   a_depth,
                                   bool  a_homoOnly = true);
  virtual void MGreclaim(NewPoissonOp* a_reclaim);

  RealVect                m_dx;
  BCFunc                  m_bc;

};

#include "NamespaceFooter.H"
#endif

---