EMM_WaveFEMTest.h

Go to the documentation of this file.

#ifndef EMM_WaveFEMTest_H
#define EMM_WaveFEMTest_H

#include "EMM_WaveTest.h"
#include "MATH_FullMatrix.h"
#include "MATH_ArrayVector.h"
#include "MATH_ConjugateGradient.h"


DEFINE_SPTR(EMM_WaveFEMTest);

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


private:
    
    //insidcates if the mass matrix is diagonal
    tBoolean mIsMassMatrixCondensed;
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_WaveFEMTest(void);
    
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_WaveFEMTest(void);
    
  

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

   

public:
    inline void setIsCondensedMassMatrix(const tBoolean& isCondensed) {
        mIsMassMatrixCondensed=isCondensed;
    }
 

   
private:

    //matrix computations methods

    void computeMassMatrix(const tUIndex& nPts,
                           const tBoolean& isPeriodic,
                           const tBoolean dirichletLC[2],
                           MATH_FullMatrix& M) const;
    
    void computeCondensedMassMatrix(const tUIndex& nPts,
                                    const tBoolean& isPeriodic,
                                    const tBoolean dirichletLC[2],
                                    MATH_ArrayVector& M) const;
    
    void computeStressMatrix(const tUIndex& nPts,
                             const tBoolean& isPeriodic,
                             const tBoolean dirichletLC[2],
                             MATH_FullMatrix& S) const;
    
    //signal propagation
  
    
    tBoolean computeWaveAtFirstTime(const MATH_Matrix& M,
                                    const MATH_Matrix& X,
                                    MATH_Solver& solver,
                                    const tBoolean& isPeriodic,
                                    const tBoolean dirichletLC[2],
                                    tReal U0(const tReal& t),//U(0,t)
                                    tReal UL(const tReal& t),//U(L,t)
                                    EMM_RealArray& Ux_0,//U(x,0)
                                    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 computeWaveAtFirstTime(const MATH_ArrayVector& M,
                                    const MATH_Matrix& X,
                                    const tBoolean& isPeriodic,
                                    const tBoolean dirichletLC[2],
                                    tReal U0(const tReal& t),//U(0,t)
                                    tReal UL(const tReal& t),//U(L,t)
                                    EMM_RealArray& Ux_0,//U(x,0)
                                    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;

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 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 {
        cout << "FEM P1 System not yet implemented \n";
        return true;
    }
    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 {
        cout << "FEM P2 System not yet implemented \n";
        return true;
    };
    protected:
    
   //wave examples
    virtual tBoolean  sinusoidalWave() const {
        return mIsMassMatrixCondensed ? EMM_WaveTest::sinusoidalWave("C-FEM"):EMM_WaveTest::sinusoidalWave("FEM");
    }
    virtual tBoolean  triangleWave() const {
        return mIsMassMatrixCondensed ? EMM_WaveTest::triangleWave("C-FEM"):EMM_WaveTest::triangleWave("FEM");
    }
    virtual tBoolean  trianglePeriodicalWave() const {
        return mIsMassMatrixCondensed ? EMM_WaveTest::trianglePeriodicalWave("C-FEM"):EMM_WaveTest::trianglePeriodicalWave("FEM");
    }
    virtual tBoolean  sharkWave() const {
        
        return mIsMassMatrixCondensed ? EMM_WaveTest::sharkWave("C-FEM"):EMM_WaveTest::sharkWave("FEM");
    }
    
    virtual tBoolean  barWave() const {
        return mIsMassMatrixCondensed ? EMM_WaveTest::barWave("C-FEM"):EMM_WaveTest::barWave("FEM");
    }
    
};


#endif