EMM_PlanarAnisotropyOperator.h

Go to the documentation of this file.

#ifndef EMM_PlanarAnisotropyOperator_H
#define EMM_PlanarAnisotropyOperator_H


#include "EMM_LinearAnisotropyOperator.h"

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

public:
   

    // ASSOCIATION
 
private:
  
    
protected: 

    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_PlanarAnisotropyOperator(void) {
    }
  
    // DESTRUCTORS 
    
    virtual ~EMM_PlanarAnisotropyOperator(void) {
    }
  
  
  
public:

    //DATA Methods
    //============
 
    /* \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}=-U_i.U_j \f$
     *
     * @return the memory size of matrix A (6)
     */
    static inline tUSInt BuildDiscretizedMatrix(const tUCInt& nU,
                                                const tReal* U,
                                                tReal* A) {
        if (nU<3) throw CORE_Exception("emicrom/core",
                                       "EMM_PlanarAnisotropyOperator::BuildDiscretizedMatrix(nU,U,A)",
                                       "nU ("+CORE_Integer::toString(nU)+") must be greater than 3");
        tDimension k,l;
        const tReal *Ul;
        const tReal *Uk=U;
        tReal *Akl=A;
        memset(A,0,6*sizeof(tReal));
        for (k=0;k<3;k++) {
            // l==k case
            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