ClpStdELCGFMProblem.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef AMROC_CLP_STD_ELC_GFM_PROBLEM_H
#define AMROC_CLP_STD_ELC_GFM_PROBLEM_H

#define OWN_GFMAMRSOLVER
#include "ClpStdGFMProblem.h"
#include "AMRELCGFMSolver.h"

#ifndef OWN_ELCGFMAMRSOLVER
class FluidSolverSpecific : 
  public AMRELCGFMSolver<VectorType,FixupType,FlagType,DIM> {
  typedef VectorType::InternalDataType DataType;
  typedef AMRELCGFMSolver<VectorType,FixupType,FlagType,DIM> base;
  typedef F77GFMFileOutput<VectorType,FixupType,FlagType,DIM> output_type;
public:
  FluidSolverSpecific() : base(_IntegratorSpecific, _InitialConditionSpecific, 
                               _BoundaryConditionsSpecific) {
    SetLevelTransfer(new F77LevelTransfer<VectorType,DIM>(f_prolong, f_restrict));
    SetFileOutput(new output_type(*this,f_flgout)); 
    SetFixup(new FixupSpecific(_IntegratorSpecific));
    SetFlagging(new FlaggingSpecific(*this)); 
    AddGFM(new GhostFluidMethod<VectorType,DIM>(
              new F77ELCGFMBoundary<VectorType,DIM>(f_ibndrfl,f_itrans,*this),
              new CPTLevelSet<DataType,DIM>()));
    SetCoupleGFM(0);
  }  
 
  ~FluidSolverSpecific() {
    DeleteGFM(_GFM[0]);
    delete _Flagging;
    delete _Fixup;
  }

  virtual void SendBoundaryData() {
    START_WATCH
      for (register int l=0; l<=FineLevel(base::GH()); l++) 
        ((output_type*) _FileOutput)->Transform(base::U(), base::Work(), 
                                                CurrentTime(base::GH(),l), l, 
                                                base::Dim()+4, base::t[l]);
    END_WATCH(FLUID_CPL_PRESSURE_CALCULATE)
    base::SendBoundaryData();
  }

protected:
  IntegratorSpecific _IntegratorSpecific;
  InitialConditionSpecific _InitialConditionSpecific;
  BoundaryConditionsSpecific _BoundaryConditionsSpecific;
}; 

#endif

#endif