EMM_WaveFDMTest.h

Go to the documentation of this file.

#ifndef EMM_WaveFDMTest_H
#define EMM_WaveFDMTest_H

#include "EMM_WaveTest.h"


DEFINE_SPTR(EMM_WaveFDMTest);

class EMM_WaveFDMTest : public EMM_WaveTest {
  
    SP_OBJECT(EMM_WaveFDMTest);
    // ATTRIBUTES


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

public:
    inline static SP::EMM_WaveFDMTest New() {
        SP::EMM_WaveFDMTest p(new EMM_WaveFDMTest(),
                           EMM_WaveFDMTest::Delete());
        p->setThis(p);
        return p;
    };

   

    

   
protected:


    //signal propagation
    tBoolean computeWaveAtPreviousTime(const tBoolean& isPeriodic,
                                       const tBoolean dirichletLC[2],
                                       tReal U0(const tReal& t),//U(0,t)
                                       tReal UL(const tReal& t),//U(L,t)
                                       const EMM_RealArray& Ux_0,//U(x,0)
                                       const EMM_RealArray& dU_dt0,//dU_dt(x,0)
                                       const tReal& s2,//s=c^2*dt^2/dx^2
                                       const tReal& dt,
                                       EMM_RealArray& Ux_1) const;
    
    tBoolean computeWaveAtFirstTime(const tBoolean& isPeriodic,
                                    const tBoolean dirichletLC[2],
                                    tReal U0(const tReal& t),//U(0,t)
                                    tReal UL(const tReal& t),//U(L,t)
                                    const EMM_RealArray& Ux_0,//U(x,0)
                                    const EMM_RealArray& dU_dt0,//dU_dt(x,0)
                                    const tReal& s,//s=c^2*dt^2/dx^2
                                    const tReal& dt,
                                    EMM_RealArray& Ux_1) const;

    tBoolean computeVelocityWaveAtPreviousTime(const tBoolean& isPeriodic,
                                               const tBoolean dirichletLC[2],
                                               tReal U0(const tReal& t),//U(0,t)
                                               tReal UL(const tReal& t),//U(L,t)
                                               const EMM_RealArray& Ux_0,//U(0,x)
                                               const EMM_RealArray& Vx_0,//V(0,x)
                                               const tReal& c,
                                               const tReal& dt,
                                               const tReal& dx,
                                               EMM_RealArray& dU_dt) const;

protected:
    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;
    
    virtual tBoolean waveSystemP1Propagation(const tBoolean& isTXTOutput,
                                             const tString& prefix,
                                             const tBoolean dirichletCL[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;
    
    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;
protected:
    
   //wave examples
    virtual tBoolean  sinusoidalWave() const {
        return EMM_WaveTest::sinusoidalWave("FDM");
    }
    virtual tBoolean  triangleWave() const {
        return EMM_WaveTest::triangleWave("FDM");
    }
    virtual tBoolean  trianglePeriodicalWave() const {
        return EMM_WaveTest::trianglePeriodicalWave("FDM");
    }
    virtual tBoolean  sharkWave() const {
        
        return EMM_WaveTest::sharkWave("FDM");
    }
    
    virtual tBoolean  barWave() const {
        return EMM_WaveTest::barWave("FDM");
    }
   
};

#endif