fsi/adlib-amroc/DetonatingCanister/src/FluidProblem.h

Go to the documentation of this file.

// -*- C++ -*-

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

#ifndef AMROC_FLUIDPROBLEM_H
#define AMROC_FLUIDPROBLEM_H

#include "eulerznd3.h"

#undef  NAUX
#define NAUX      6

#include "ClpProblem.h"

#define OWN_ELCGFMAMRSOLVER
#include "ClpStdELCGFMProblem.h"

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>()));
    AddGFM(new GhostFluidMethod<VectorType,DIM>(
              new F77GFMBoundary<VectorType,DIM>(f_ibndrfl,f_itrans),
              new F77GFMLevelSet<DataType,DIM>(f_lset)));
    SetCoupleGFM(0);
  }  
 
  ~FluidSolverSpecific() {
    DeleteGFM(_GFM[1]);
    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;