amroc/amr/InitialCondition.h

Go to the documentation of this file.

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus
//
// Copyright (C) 2003-2007 California Institute of Technology
// Ralf Deiterding, ralf@cacr.caltech.edu

#ifndef AMROC_INITIALCONDITION_H
#define AMROC_INITIALCONDITION_H

#include <string>
#include "AMRBase.h"

template <class VectorType, int dim>
class InitialCondition : public AMRBase<VectorType,dim> {
  typedef AMRBase<VectorType,dim> base;
  typedef typename VectorType::InternalDataType DataType;
public:
  typedef typename base::vec_grid_fct_type vec_grid_fct_type;  
  typedef typename base::vec_grid_data_type vec_grid_data_type;
  typedef GridFunction<DataType,dim> grid_fct_type;
  typedef GridData<DataType,dim> grid_data_type;

  InitialCondition() : base() {}

  virtual ~InitialCondition() {}

  //******************************************************************************
  // Abstract class interface
  //******************************************************************************
  virtual void SetGrid(vec_grid_data_type& gd, grid_data_type& gdw, const int& level) = 0;   
  //******************************************************************************

  virtual void register_at(ControlDevice& Ctrl, const std::string& prefix) {
    base::LocCtrl = Ctrl.getSubDevice(prefix+"InitialCondition");
  }
  virtual void register_at(ControlDevice& Ctrl) {
    register_at(Ctrl, "");
  }

  virtual void Set(vec_grid_fct_type& u, grid_fct_type& work, const int Level) { 
    SetPhysicalTime(u,0,Level,0.0);
    forall (u,0,Level,c)
      SetGrid(u(0,Level,c),work(0,Level,c),Level);    
    end_forall
  }
};


#endif

---