EMM_DisplacementWaveTest.h

Go to the documentation of this file.

#ifndef EMM_DisplacementWaveTest_H
#define EMM_DisplacementWaveTest_H

#include "EMM_Test.h"


DEFINE_SPTR(EMM_DisplacementWaveTest);

class EMM_DisplacementWaveTest : public EMM_Test {
  
    SP_OBJECT(EMM_DisplacementWaveTest);
    // ATTRIBUTES


public:
   
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_DisplacementWaveTest(void);
    
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_DisplacementWaveTest(void);
    
  

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

   

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

   
private:
 
    tBoolean loadData(const tString& fn,
                      tBoolean& isTXTOutput,
                      tBoolean& isVTIOutput,
                      tUSInt& dim,
                      tBoolean periodicity[3],
                      tBoolean dirichletLC[6],
                      tFlag& neumannFaces,
                      tReal neumannLC[6],
                      tReal& C,tReal& lambda,tReal& mu,tReal& rho,
                      tReal& L100,tReal& L111,
                      tReal L[3],
                      tUInteger N[3],
                      tUIndex& nTimes,
                      tReal& cfl,
                      tString& name,
                      EMM_RealArray& U,
                      tReal M0[3],
                      tString& solverName) const;
    
    tBoolean elasticWaveTest(const tBoolean& isHyperElastic,
                             const tString& method,
                             const tBoolean& equilibrium,
                             const tString& optionsFile)  const;  
  
    
    tBoolean saveIntoFile(const tString& fileName,const tReal& Lc,const tReal L[],const tUIndex N[],const EMM_RealField& U)  const;
      
    tBoolean saveIntoFile(const tString& fileName,const tUIndex& n, const vector<const tReal*>& data)  const;
    
    tBoolean loadFromFile(const tString& fileName,const tUIndex &N,tReal* times, tReal * data)  const;
   
    tBoolean loadFromFile(const tString& fileName,const tUIndex &N,EMM_RealField& U)  const;
    
private:
    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