---

AMRParticleProjector.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 _AMRPARTICLEPROJECTOR_H_
#define _AMRPARTICLEPROJECTOR_H_



#include "DisjointBoxLayout.H"
#include "LevelData.H"
#include "BinFab.H"
#include "DragParticle.H"

class AMRParticleProjector
{

public:
  AMRParticleProjector();


  AMRParticleProjector(const Vector<DisjointBoxLayout>& a_vectGrids,
                       const Vector<DisjointBoxLayout>& a_vectParticleGrids,
                       const Vector<ProblemDomain>& a_vectDomain,
                       const Vector<int>& a_vectRefRatio,
                       const Vector<Real>& a_vectDx,
                       const Vector<LevelData<BinFab<DragParticle> >* >& a_vectParticles,
                       Real a_spreadingRadius,
                       Real a_correctionRadius = 2,
                       int a_gridsGrow=0);

  ~AMRParticleProjector();

  void define(const Vector<DisjointBoxLayout>& a_vectGrids,
              const Vector<DisjointBoxLayout>& a_vectParticleGrids,
              const Vector<ProblemDomain>& a_vectDomain,
              const Vector<int>& a_vectRefRatio,
              const Vector<Real>& a_vectDx, 
              const Vector<LevelData<BinFab<DragParticle> >* >& a_vectParticles,
              Real a_spreadingRadius,
              Real a_correctionRadius = 2,
              int a_gridsGrow = 0);

  void clear();
    


  void projectForce(Vector<LevelData<FArrayBox>* >& a_vectForce,
                    const Vector<LevelData<BinFab<DragParticle> >* >& a_particles);
  

  void setSpreadingRadius(const Real a_rad);
  
  void setCorrectionRadius(const Real a_rad);
  
  
  Real getSpreadingRadius() const;

  Real getCorrectionRadius() const;

  void doImages(bool m_doImages);

  bool doImages() const;


  void setGridsGrow(int a_gridsGrow);
  
  void infDomainCoarsenFactor(int a_factor);

  bool isDefined() const {return (m_isDefined && m_grownGridsDefined);};

  bool grownGridsDefined() const {return m_grownGridsDefined;};
  
  
  void setVerbosity(int a_verbosity);
  
  int verbosity() const;

  void computeD(Vector<LevelData<FArrayBox>* >& a_vectD,
                const Vector<LevelData<BinFab<DragParticle> >* >& a_particles);


  
  void solveForProjForce(Vector<LevelData<FArrayBox> *>& a_projectedForce,
                         const Vector<LevelData<FArrayBox>* >& a_D,
                         int a_lbase);



  void defineImages(List<DragParticle>& a_imageParticles,
                    const Vector<LevelData<BinFab<DragParticle> >* >& a_particles);
    



  void defineGrownGrids(const Vector<DisjointBoxLayout>& a_grids,
                        const Vector<LevelData<BinFab<DragParticle> >* >& a_vectParticles);

 

protected:
  
  void createImageParticle(DragParticle& a_image,
                           const DragParticle& a_particle,
                           Real bndryLoc, int dir, 
                           Side::LoHiSide a_side);
  


  void setDefaultValues();

  // data members -- note that we maintain two distinct grid hierachies
  // the original one, and a second one which is grown to account for
  // image particles.

  Vector<DisjointBoxLayout> m_vectGrids;

  Vector<DisjointBoxLayout> m_vectParticleGrids;

  Vector<DisjointBoxLayout> m_vectGrownGrids;

  Vector<DisjointBoxLayout> m_vectGrownParticleGrids;
  
  Vector<ProblemDomain> m_vectDomain;

  Vector<ProblemDomain> m_vectGrownDomain;

  Vector<int> m_vectRefRatio;

  Vector<int> m_vectGrownRefRatio;


  Vector<Real> m_vectDx;

  Vector<Real> m_vectGrownDx;
  
  Real m_spreadingRadius;
  
  Real m_correctionRadius;



  int m_gridsGrow;

  int m_infDomainCoarsenFactor;
  
  List<DragParticle> m_imageParticles;

  bool m_doImages;

  bool m_returnInfDomBC;

  bool m_isDefined;


  bool m_grownGridsDefined;

  int m_verbosity;

private:

  // copy constructor illegal, etc
  AMRParticleProjector (const AMRParticleProjector& );
  AMRParticleProjector operator= (const AMRParticleProjector&);

};

#endif

---