clawpack/applications/euler/2d/Shockbubble/src/Problem.h

// -*- C++ -*-

// Copyright (C) 2002 Ralf Deiterding
// Brandenburgische Universitaet Cottbus

#ifndef AMROC_PROBLEM_H
#define AMROC_PROBLEM_H

#include "euler2.h"

#undef  NAUX
#define NAUX       1

#include "ClpProblem.h"

#define OWN_FLAGGING

class ScaledGradientSpecific : public ScaledGradient<VectorType,FlagType,DIM> {
  typedef ScaledGradient<VectorType,FlagType,DIM> base;
  typedef VectorType::InternalDataType DataType;
public:
  typedef base::grid_fct_type grid_fct_type;
  typedef base::flag_fct_type flag_fct_type;

  virtual bool FlagByScaledGradient(grid_fct_type& work, flag_fct_type& flags,
                                    const int& Time, const int& Level, 
                                    const DataType& TolVal, const FlagType& FlagValue) { 
    if (TolVal <= 0.0) return false;
    int TStep = TimeStep(work,Level);
    forall(work,Time,Level,c)
      work(Time+TStep,Level,c) = 0.0;
    end_forall 
    int East[2]  = { 1, 0 }; int West[2]  = { -1, 0 };
    int North[2] = { 0, 1 }; int South[2] = { 0, -1 };
    base::Difference(work, Time, Level, East, West);
    base::Difference(work, Time, Level, North, South);      
    base::FlagByValue(work, flags, Time, Level, TolVal, FlagValue);
    return true;
  }
};

class FlaggingSpecific : 
  public AMRFlagging<VectorType,FixupType,FlagType,DIM> {
  typedef AMRFlagging<VectorType,FixupType,FlagType,DIM> base;
  typedef base::solver_type solver_type;
public:
  FlaggingSpecific(solver_type& solver) : base(solver) {
      base::AddCriterion(new ScaledGradientSpecific());
      base::AddCriterion(new AbsoluteError<VectorType,FixupType,FlagType,DIM>(solver));
    }
  ~FlaggingSpecific() { DeleteAllCriterions(); }
};