EMM_OperatorsTest.h

Go to the documentation of this file.

#ifndef EMM_OperatorsTest_H
#define EMM_OperatorsTest_H

#include "EMM_Test.h"

#include "EMM_Matter.h"
#include "EMM_RealField.h"
#include "EMM_Grid3D.h"
#include "EMM_MagneticExcitationOperator.h"
#include "EMM_MatterField.h"

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

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

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

   

    virtual tBoolean test(const CORE_Run& runner,const map<tString,tString>& options) const;
    
    virtual tBoolean performanceTest(const CORE_Run& runner,const map<tString,tString>& options) const;
  
   
protected:
   
    
    
    // OPERATOR TESTS

  

    tReal testSelfAdjointness(const tUIndex& nCells,
                              const EMM_RealArray& sigma,
                              const EMM_MagneticExcitationOperator& op) const;
    
    tReal symmetricTest(const tBoolean& isSymmetric,
                        const tUIndex& nCells,
                        const EMM_RealArray& sigma,
                        const EMM_MagneticExcitationOperator& op) const;


    tReal testEDerivatives(const EMM_RealArray& sigma,const tReal& cellVolume,
                           const tReal& h,const tReal& eps_dh,
                           EMM_RealField& M,const EMM_MagneticExcitationOperator& op) const;


    tReal testHDerivatives(const EMM_RealArray& sigma,const tReal& cellVolume,
                           const tReal& h,const tReal& eps_dh,
                           EMM_RealField& M,const EMM_MagneticExcitationOperator& op) const;

    tBoolean genericOperatorTest(const SV::CORE_Object& args,
                                 const map<tString,tString>& options,
                                 const tString & opClassName,
                                 const tString& outputPath,
                                 const tReal& errMax,
                                 const tReal& dh,const tReal& eps_dh,
                                 const tReal& eps_rel,
                                 const tBoolean& isSymmetric,
                                 const tFlag& anisotropy) const;

private:
    template<class T>
    static tReal compareDiscretizedData(const tUSInt& d,const tString& fn,const EMM_MagneticExcitationOperator& op,tUIndex& indexMax);
   
    
private:
    //primary tests for operators
    tBoolean magnetostrictionPrimaryTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    tBoolean demagnetizedPrimaryTests(const CORE_Run& runner,const map<tString,tString>& options) const;
    
    //generic tests for operators
    tBoolean zeemanOperatorTest() const;
    tBoolean exchangeOperatorTest() const;
    tBoolean anisotropyOperatorTest() const;
    tBoolean demagnetizedOperatorTest() const;
    tBoolean magnetostrictionOperatorTest() const;
    tBoolean staticMagneticTensorOperatorTest() const;

 
    
  
};

#include "EMM_OperatorsTest.hpp"
#endif