RIMStdProblem.h

Go to the documentation of this file.

// -*- C++ -*-
// $Header: /home/proj/vtf3d/vtf/amroc/rim/RIMStdProblem.h,v 1.2 2004/12/01 20:51:35 ralf Exp $

#ifndef AMROC_RIM_STD_PROBLEM_H
#define AMROC_RIM_STD_PROBLEM_H

class IntegratorSpecific :  
  public RIMIntegrator<VectorType,DIM> {
public:   
  IntegratorSpecific() : 
#ifdef f_flgout
    RIMIntegrator<VectorType,DIM>(f_rimstep,f_check) {}
#else   
    RIMIntegrator<VectorType,DIM>(f_rimstep) {} 
#endif
};
  
class FixupSpecific : 
  public AMRFixup<VectorType,FixupType,DIM> {
public:
  FixupSpecific() : AMRFixup<VectorType,FixupType,DIM>() {}
};

class InitialConditionSpecific : 
  public F77InitialCondition<VectorType,DIM> {
public:    
  InitialConditionSpecific() : 
    F77InitialCondition<VectorType,DIM>(f_initial) {}
};
  
class BoundaryConditionsSpecific : 
  public F77BoundaryConditions<VectorType,DIM> {
public:    
  BoundaryConditionsSpecific() : 
    F77BoundaryConditions<VectorType,DIM>(f_boundary) {}
};

#ifndef OWN_FLAGGING
class FlaggingSpecific : 
  public AMRFlagging<VectorType,FixupType,FlagType,DIM> {
  typedef AMRFlagging<VectorType,FixupType,FlagType,DIM> base;
public:
  FlaggingSpecific(base::solver_type& solver) : base(solver) {
      base::AddCriterion(new ByValue<VectorType,FlagType,DIM>());
      base::AddCriterion(new ScaledGradient<VectorType,FlagType,DIM>());
      base::AddCriterion(new LimiterType<VectorType,FlagType,DIM>());
      base::AddCriterion(new AbsoluteError<VectorType,FixupType,FlagType,DIM>(solver));
      base::AddCriterion(new RelativeError<VectorType,FixupType,FlagType,DIM>(solver));
#ifdef f_flgout
      base::AddCriterion(new F77ByValue<VectorType,FlagType,DIM>(f_flgout));
      base::AddCriterion(new F77ScaledGradient<VectorType,FlagType,DIM>(f_flgout));
      base::AddCriterion(new F77LimiterType<VectorType,FlagType,DIM>(f_flgout));
      base::AddCriterion(new F77AbsoluteError<VectorType,FixupType,FlagType,DIM>(solver,f_flgout));
      base::AddCriterion(new F77RelativeError<VectorType,FixupType,FlagType,DIM>(solver,f_flgout));
#endif
    }

  ~FlaggingSpecific() { DeleteAllCriterions(); }
};
#endif 

#ifndef OWN_AMRSOLVER
class SolverSpecific : 
  public AMRSolver<VectorType,FixupType,FlagType,DIM> {
  typedef AMRSolver<VectorType,FixupType,FlagType,DIM> base;
public:
  SolverSpecific(IntegratorSpecific& integ, 
                 base::initial_condition_type& init,
                 base::boundary_conditions_type& bc) :
    AMRSolver<VectorType,FixupType,FlagType,DIM>(integ, init, bc) {
    SetLevelTransfer(new F77LevelTransfer<VectorType,DIM>(f_prolong, f_restrict));
#ifdef f_flgout
    SetFileOutput(new F77FileOutput<VectorType,DIM>(f_flgout)); 
#else   
    SetFileOutput(new FileOutput<VectorType,DIM>()); 
#endif
    SetFixup(new FixupSpecific());
    SetFlagging(new FlaggingSpecific(*this)); 
  }  

  ~SolverSpecific() {
    delete _LevelTransfer;
    delete _Flagging;
    delete _Fixup;
    delete _FileOutput;
  }
}; 
#endif 
 
class SolverControlSpecific : 
  public SolverControl {
public:
  SolverControlSpecific(AMRSolver<VectorType,FixupType,FlagType,DIM>& solver) : 
    SolverControl(solver) {}
}; 


#endif