EMM_Output.h

Go to the documentation of this file.

#ifndef EMM_Output_H
#define EMM_Output_H

#include "EMM_IO.h"
#include "EMM_RealArray.h"
#include "EMM_VTK.h"
#include "EMM_LandauLifschitzSystem.h"

DEFINE_SPTR(EMM_Output);
class EMM_Output : public EMM_IO {
  
    SP_OBJECT(EMM_Output);
    // ATTRIBUTES
  
public:
   
  
   
private:
    //vtk readerd
    SP::EMM_VTK mVTK;
    //mModule: work array only used to save in AIM file.
    //is mutable because this modification has no incidence on class state
    mutable EMM_RealArray mModule;
    
protected: 
    // METHODS
  
    // CONSTRUCTORS 
  
    EMM_Output(void);
  
  
  
    // DESTRUCTORS 
  
    
    virtual ~EMM_Output(void);


public:
  
    inline static SP::EMM_Output New() {
        SP::EMM_Output p(new EMM_Output(),EMM_Output::Delete());
        p->setThis(p);
        return p;
    };
    
private:
    // SET methods
public:
   

  
    // GET methods 
  
public:
    // OTHERS methods

    virtual tBoolean save(const EMM_LandauLifschitzSystem& system) {
        return save(system,"");
    }
    tBoolean save(const EMM_LandauLifschitzSystem& system,const tString& suffix);

    inline tBoolean saveBackup(const EMM_LandauLifschitzSystem& system,
                               const tInteger& backupIndex) {
        tString strIndex="";
        strIndex="_back"+CORE_Integer::toString(backupIndex,getFileDigitsNumber());
        return save(system,strIndex);
    }
    static tBoolean saveVTI(const tUInteger& Nx,const tUInteger& Ny,const tUInteger& Nz,
                            const tReal& Hx,const tReal& Hy,const tReal& Hz,
                            const tString& Fname, const EMM_RealField& F,const tString& fileName) {
        tUIndex nValues;
        const tReal *Fvalues;
        F.getValues(nValues,Fvalues);
        return saveVTI(Nx,Ny,Nz,Hx,Hy,Hz,
                       Fname,F.getSize(),F.getDimension(),Fvalues,fileName);
    }
    inline static tBoolean saveVTI(const tUInteger& Nx,const tUInteger& Ny,const tUInteger& Nz,
                                   const tReal& Hx,const tReal& Hy,const tReal& Hz,
                                   const tString& Fname, const tUIndex& n, const tDimension& dim,const tReal* Fvalues,const tString& fileName) {
        //center the geometry
        tReal Ox=(Nx*(1-Hx))/2;
        tReal Oy=(Ny*(1-Hy))/2;
        tReal Oz=(Nz*(1-Hz))/2;
        
        return saveVTI(Ox,Oy,Oz,Nx,Ny,Nz,Hx,Hy,Hz,
                       Fname,n,dim,Fvalues,fileName);
        
    }
    
    static tBoolean saveVTI(const tReal& Ox,const tReal& Oy,const tReal& Oz,
                            const tUInteger& Nx,const tUInteger& Ny,const tUInteger& Nz,
                            const tReal& Hx,const tReal& Hy,const tReal& Hz,
                            const tString& Fname, const tUIndex& n, const tDimension& dim,const tReal* Fvalues,const tString& fileName);
    
    
    tBoolean saveVTI(const EMM_LandauLifschitzSystem& system,
                     const tString& fileName);

    inline tBoolean saveVTI(const EMM_LandauLifschitzSystem& system) {
        return saveVTI(system,getPath()+getPrefix()+".vti");
    }
  
    
    tBoolean saveRelaxation( const EMM_LandauLifschitzSystem& system,
                            const tString& fileName);
    
    tBoolean saveFieldToTXTFile(const EMM_RealField& M,const tString& fileName) const;

private:
    tBoolean saveFieldToAIMFile(const EMM_RealField& M,const EMM_RealField& H) const;

    
  
};

#endif