clawpack/applications/testing/2d/GFMBCTest/src/Problem.h

// -*- C++ -*-

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

#ifndef AMROC_PROBLEM_H
#define AMROC_PROBLEM_H

#include "euler2.h"
#include "ClpProblem.h"

#define OWN_INTEGRATOR
class IntegratorSpecific :  
  public ClpIntegrator<VectorType,AuxType,DIM> {
public:   
  IntegratorSpecific() : 
    ClpIntegrator<VectorType,AuxType,DIM>(NEQUSED,NWAVES,f_normalflux,f_transverseflux,
                                          f_check,f_setaux,f_source) {}
  virtual double CalculateGrid(vec_grid_data_type& NewStateVec, 
                               vec_grid_data_type& OldStateVec,
                               vec_grid_data_type* Flux[],
                               const int& level, const double& t, 
                               const double& dt, const int& mpass)
    { NewStateVec.copy(OldStateVec); return 1.0; } 
};
   
#define OWN_GFMAMRSOLVER
#include "ClpStdGFMProblem.h"
#include "StatCPTLevelSet.h"

class SolverSpecific : 
  public AMRGFMSolver<VectorType,FixupType,FlagType,DIM> {
  typedef VectorType::InternalDataType DataType;
  typedef AMRGFMSolver<VectorType,FixupType,FlagType,DIM> base;
public:
  SolverSpecific(IntegratorSpecific& integ, 
                 base::initial_condition_type& init,
                 base::boundary_conditions_type& bc) : base(integ, init, bc) {
    SetLevelTransfer(new F77LevelTransfer<VectorType,DIM>(f_prolong, f_restrict));
#ifdef f_flgout
    SetFileOutput(new F77GFMFileOutput<VectorType,FixupType,FlagType,DIM>(*this,f_flgout)); 
#else   
    SetFileOutput(new GFMFileOutput<VectorType,FixupType,FlagType,DIM>(*this)); 
#endif
    SetFixup(new FixupSpecific(integ));
    SetFlagging(new FlaggingSpecific(*this)); 
    AddGFM(new GhostFluidMethod<VectorType,DIM>(
              new F77GFMBoundary<VectorType,DIM>(f_ibndrfl,f_itrans),
              new StatCPTLevelSet<DataType,DIM>()));
  }  
 
  ~SolverSpecific() {
    DeleteGFM(_GFM[0]);
    delete _Flagging;
    delete _Fixup;
  }