EMM_UniaxialAnisotropyOperator.h

Go to the documentation of this file.

#ifndef EMM_UniaxialAnisotropyOperator_H
#define EMM_UniaxialAnisotropyOperator_H


#include "EMM_LinearAnisotropyOperator.h"

DEFINE_SPTR(EMM_UniaxialAnisotropyOperator);
class EMM_UniaxialAnisotropyOperator : public virtual EMM_LinearAnisotropyOperator {
    
    SP_OBJECT(EMM_UniaxialAnisotropyOperator);
    // ATTRIBUTES
    
    
private:

public:
   

    // ASSOCIATION
 
private:
  
    
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_UniaxialAnisotropyOperator(void) {
    }
  
    // DESTRUCTORS 
    
    virtual ~EMM_UniaxialAnisotropyOperator(void) {
    }
  

    //DATA Methods
    //============
    
public:
    /* \brief build the symmetric matrix of size 6
     * @param nU: size of U (nU>=3)
     * @param U: anisotropy direction of size 3
     * @param A: anisotropy matrix of size 6
     *
     * It computes  the  inferior part of the  packed symmetric matrix A:
     *
     * \f$ A_{ij}=(-\delta_{ij}+U_i.U_j) \f$ where \f$ \delta_{ij}=1\f$ only and only if i=j, 0 otherwise.
     *
  
     * 
     * @return the memory size of matrix A (6)
     */
    static inline tUSInt BuildDiscretizedMatrix(const tUCInt& nU,
                                                const tReal* U,
                                                tReal* A) {
        
        tDimension k,l;
        const tReal *Ul;
        const tReal *Uk=U;
        tReal *Akl=A;
        memset(A,0,6*sizeof(tReal));
        if (nU<3) throw CORE_Exception("emicrom/core",
                                       "EMM_UniaxialAnisotropyOperator::BuildDiscretizedMatrix(nU,U,A)",
                                       "nU ("+CORE_Integer::toString(nU)+") must be greater than 3");
        for (k=0;k<3;k++) {
            // l==k case
            (*Akl)=-1;
            Ul=Uk;
            for (l=k;l<3;l++) {  
                (*Akl)+=(*Uk)*(*Ul);
                //next l
                Akl++;
                Ul++;
                
            } //end loop on l
            
            //next k
            Uk++;
        } // end loop on k
        
        return 6;
    };

};
#endif