/import/home/rettenbs/src/SWE/src/blocks/rusanov/SWE_RusanovBlock.hh

Go to the documentation of this file.

#ifndef __SWE_RUSANOVBLOCK_HH
#define __SWE_RUSANOVBLOCK_HH

#include <iostream>
#include <stdio.h>
#include <fstream>

#include "tools/help.hh"
#include "SWE_Block.hh"


using namespace std;

class SWE_RusanovBlock : public SWE_Block {

  public:
    // Constructor und Destructor
    SWE_RusanovBlock(float _offsetX = 0, float _offsetY = 0);
    virtual ~SWE_RusanovBlock();
    
  // object methods
    virtual void simulateTimestep(float dt);
    virtual float simulate(float tStart, float tEnd);
    
    virtual void computeNumericalFluxes();
    virtual void updateUnknowns(float dt);

  protected:
     
    virtual void computeBathymetrySources();

    static float computeFlux(float fLoc, float fNeigh, float xiLoc, float xiNeigh, float llf);
    float computeLocalSV(int i, int j, char dir);

    // compute the largest allowed time step for the current grid block
    virtual void computeMaxTimestep() {
       SWE_Block::computeMaxTimestep();
       // more pessimistic choice of the time step
       maxTimestep *= 0.5; 
    };

    // define additional arrays for temporary unknowns: 
    // - arrays to hold the values of the flux terms at cell edges
    Float2D Fh;
    Float2D Fhu;
    Float2D Fhv;
    Float2D Gh;
    Float2D Ghu;
    Float2D Ghv;
    // - arrays to hold the bathymetry source terms for the hu and hv equations
    Float2D Bx;
    Float2D By;
    
    // overload operator<< such that data can be written via cout <<
    // -> needs to be declared as friend to be allowed to access private data
    friend ostream& operator<< (ostream& os, const SWE_RusanovBlock& swe);
  
};

ostream& operator<< (ostream& os, const SWE_RusanovBlock& swe);



#endif