GFMLevelSet.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@amroc.net

#ifndef AMROC_GFMLEVELSET_H
#define AMROC_GFMLEVELSET_H

template <class DataType, int dim>
class GFMLevelSet : public controlable {

public:
  typedef GridFunction<DataType,dim> grid_fct_type;   
  typedef GridData<DataType,dim> grid_data_type;   

  GFMLevelSet() : _PlotPhi(0), _Stationary(0), _Ghosts(0), _Dim(dim) {    
    _Hierarchy = (GridHierarchy*) 0;
  }

  virtual ~GFMLevelSet() {}

  //******************************************************************************
  // Abstract class interface
  //******************************************************************************
  virtual void SetGrid(grid_data_type& gdphi, const int& Level, const double& t) = 0;   
  //******************************************************************************

  virtual void register_at(ControlDevice& Ctrl) {}
  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {}
  virtual void init() {}
  virtual void update() {}
  virtual void finish() {} 

  virtual void SetPhi(grid_fct_type& phi, const int Time, const int Level, double t) { 
    forall (phi,Time,Level,c)
      SetGrid(phi(Time,Level,c),Level,t);    
    end_forall
    START_WATCH
      Sync(phi,Time,Level);
    END_WATCH(LS_SYNC)
    START_WATCH
      ExternalBoundaryUpdate(phi,Time,Level);
    END_WATCH(BOUNDARIES_EXTERNAL)
  }

  inline const int& PlotPhi() const { return _PlotPhi; }
  inline const int& Stationary() const { return _Stationary; }

  virtual void SetupData(GridHierarchy* gh, const int& ghosts) {
    _Hierarchy = gh;
    _Ghosts = ghosts;
  }

  GridHierarchy& GH() { return *_Hierarchy; }
  const GridHierarchy& GH() const { return *_Hierarchy; }

  const int& NGhosts() const { return _Ghosts; }
  const int& Dim() const { return _Dim; }

protected:
  int _PlotPhi;
  int _Stationary;
  int _Ghosts;
  int _Dim;
  GridHierarchy* _Hierarchy;
  ControlDevice LocCtrl;
};


#endif

---