EMM_WaveTest.h

Go to the documentation of this file.

#ifndef EMM_WaveTest_H
#define EMM_WaveTest_H

#include "EMM_Test.h"


DEFINE_SPTR(EMM_WaveTest);
class EMM_WaveTest : public EMM_Test {
  
    SP_OBJECT(EMM_WaveTest);
    // ATTRIBUTES


private:
   
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_WaveTest(void);
    
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_WaveTest(void);
    
  

public:
   
   

    
    
    virtual tBoolean test(const CORE_Run& runner,
                          const map<tString,tString>& options) const;

   
protected:

    tReal computeEnergy(const tUIndex& n,
                        const tReal& C,
                        const tReal& dx,
                        const tReal& dt,
                        const tReal* U,
                        const tReal* Ut) const;
    
    virtual tBoolean waveTE2Propagation(const tBoolean& isForward,
                                        const tBoolean& isTXTOutput,
                                        const tString& prefix,
                                        const tBoolean& isPeriodic,
                                        const tBoolean dirichletLC[2],
                                        tReal U_0(const tReal& t),//U(0,t)
                                        tReal U_L(const tReal& t),//U(L,t)
                                        const tReal& C,
                                        const tReal &L,
                                        const tReal& cfl,
                                        const tUIndex nT,
                                        EMM_RealArray& U,
                                        EMM_RealArray& V) const=0;
    
    virtual tBoolean waveSystemP1Propagation(const tBoolean& isTXTOutput,
                                             const tString& prefix,
                                             const tBoolean dirichletLC[2],
                                             tReal U_0(const tReal& t),//U(0,t)
                                             tReal U_L(const tReal& t),//U(L,t)
                                             const tReal& c,
                                             const tReal& L,
                                             const tReal& dt,
                                             const tUIndex nT,
                                             EMM_RealArray& U,
                                             EMM_RealArray& V) const=0;
    virtual tBoolean waveSystemP2Propagation(const tBoolean& isTXTOutput,
                                             const tString& prefix,
                                             const tBoolean dirichletLC[2],
                                             tReal U_0(const tReal& t),//U(0,t)
                                             tReal U_L(const tReal& t),//U(L,t)
                                             const tReal& c,
                                             const tReal& L,
                                             const tReal& cfl,
                                             const tUIndex nT,
                                             EMM_RealArray& U,
                                             EMM_RealArray& V) const=0;

    
    tBoolean save(const tString& fileName,
                  const tReal &dx,
                  const tReal *U,const tUIndex &nPts) const;
   
    tBoolean save(const tString& fileName,
                  const tReal &dx,
                  const tReal *U,const tReal * V,const tUIndex &nPts) const;

    tBoolean load(const tString& optionFile,
                  tString& path,
                  tBoolean& isTXTOutput,
                  tBoolean& isPeriodic,
                  tBoolean dirichletLC[2],
                  tReal& L,
                  tUIndex& nPoints,
                  tReal& Umax,
                  tReal& C,
                  tUIndex& nTimes,
                  tReal& cfl) const;

    //wave examples
    tBoolean  sinusoidalWave(const tString& method) const;
    tBoolean  triangleWave(const tString& method) const;
    tBoolean  trianglePeriodicalWave(const tString& method) const;
    tBoolean  sharkWave(const tString& method) const;
    tBoolean  barWave(const tString& method) const;

    //wave examples
    virtual tBoolean  sinusoidalWave() const=0;
    virtual tBoolean  triangleWave() const=0;
    virtual tBoolean  trianglePeriodicalWave() const=0;
    virtual tBoolean  sharkWave() const=0;
    virtual tBoolean  barWave() const=0;
    
    //function for x or t
    static tReal null_function(const tReal& t);
    static tReal shark_function(const tReal& t);
    static tReal triangle_function(const tReal& L,const tReal& U,const tReal& x);
    static tReal sinusoidal_function(const tReal& L,const tReal& U,const tReal& x);
    static void linearPieceWise_function(const tUInteger& n,//number of points
                                         const tReal *X,//X values
                                         const tReal* Y,//Y values
                                         const tUInteger& nEvals,//number of evaluations
                                         const tReal *Xvalues,//values of evaluations of x
                                         tReal *Yvalues);
};

#endif