EMM_TensorsTest.h

Go to the documentation of this file.

#ifndef EMM_TensorsTest_H
#define EMM_TensorsTest_H

#include "CORE_Test.h"
#include "EMM_4Tensors.h"
#include "EMM_4SymmetricTensors.h"

class CORE_Run;

DEFINE_SPTR(EMM_TensorsTest);
class EMM_TensorsTest : public CORE_Test {
  
    SP_OBJECT(EMM_TensorsTest);
    // ATTRIBUTES


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

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

   

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

   
   
    
private:
    

    //tensor tests
    //=============

    tBoolean initializationTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    /* \brief test D=(A:B)
     */
    tBoolean doubleDotTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    /* \brief test D=(A:B):C
     */
    tBoolean squaredDoubleDotTests(const CORE_Run& runner,const map<tString,tString>& options) const;

    //vector tests
    //===========
    tBoolean doubleDotCrossProductTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    tBoolean doubleDotProductTests(const CORE_Run& runner,const map<tString,tString>& options) const;


    //scalar test
    //============
    /* \brief test E=(T:phi):psi
     */
    tBoolean squaredDoubleDotScalarTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    /* \brief test E=(T:UxV):phi
     */
    tBoolean doubleDotCrossDoubleDotScalarTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    /* \brief test E=(T:phi):(W:Z)
     */
    tBoolean squaredDoubleDotCrossScalarTests(const CORE_Run& runner,const map<tString,tString>& options) const;

    /* \brief test E=(T:UxV):(W:Z)
     */
    tBoolean doubleDotCrossSquaredScalarTests(const CORE_Run& runner,const map<tString,tString>& options) const;

    //2 Packed Symmetric tensors test
    tBoolean packedSymmetricTensorsTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    
    tBoolean verifySymmetry(const EMM_4Tensors& A) const;
    tBoolean verifyLeftSymmetry(const EMM_4Tensors& A) const;
    tBoolean verifyRightSymmetry(const EMM_4Tensors& A) const;
    tBoolean verifyCenterSymmetry(const EMM_4Tensors& A) const;

    
    /* \brief build A_{pqrs}=p+q+r+s+2
     */
    SP::EMM_4SymmetricTensors buildTensorS(const tString& type) const;
    /* \brief build A_{pqrs}=(p+q+1).(r+s+1)
     */
    SP::EMM_4SymmetricTensors buildTensorP(const tString& type) const;
    
};

#endif