FM3DMagnetoStaticH5Block.h

Go to the documentation of this file.

#ifndef CLASSIC_FIELDMAP3DMAGNETOSTATICH5BLOCK_HH
#define CLASSIC_FIELDMAP3DMAGNETOSTATICH5BLOCK_HH

#include "Fields/Fieldmap.h"
#include "hdf5.h"
#include "H5hut.h"

class FM3DMagnetoStaticH5Block: public Fieldmap {

public:
    virtual bool getFieldstrength(const Vector_t &R, Vector_t &E, Vector_t &B) const;
    virtual void getFieldDimensions(double &zBegin, double &zEnd, double &rBegin, double &rEnd) const;
    virtual void getFieldDimensions(double &xIni, double &xFinal, double &yIni, double &yFinal, double &zIni, double &zFinal) const;
    virtual bool getFieldDerivative(const Vector_t &R, Vector_t &E, Vector_t &B, const DiffDirection &dir) const;
    virtual void swap();
    virtual void getInfo(Inform *msg);
    virtual double getFrequency() const;
    virtual void setFrequency(double freq);

    virtual bool isInside(const Vector_t &r) const;
private:
    FM3DMagnetoStaticH5Block(std::string aFilename);
    ~FM3DMagnetoStaticH5Block();

    virtual void readMap();
    virtual void freeMap();

    std::vector<h5_float64_t> FieldstrengthEz_m;    
    std::vector<h5_float64_t> FieldstrengthEx_m;    
    std::vector<h5_float64_t> FieldstrengthEy_m;    
    std::vector<h5_float64_t> FieldstrengthBz_m;    
    std::vector<h5_float64_t> FieldstrengthBx_m;    
    std::vector<h5_float64_t> FieldstrengthBy_m;    
    h5_float64_t frequency_m;

    h5_float64_t xbegin_m;
    h5_float64_t xend_m;
    //     int xcentral_idx_m;

    h5_float64_t ybegin_m;
    h5_float64_t yend_m;
    //     int ycentral_idx_m;

    h5_float64_t zbegin_m;
    h5_float64_t zend_m;

    h5_float64_t hx_m;                   
    h5_float64_t hy_m;                   
    h5_float64_t hz_m;                   
    h5_int64_t num_gridpx_m;              
    h5_int64_t num_gridpy_m;              
    h5_int64_t num_gridpz_m;              
    bool swap_m;
    friend class Fieldmap;
};

inline bool FM3DMagnetoStaticH5Block::isInside(const Vector_t &r) const
{
    return ((r(0) >= xbegin_m && r(0) < xend_m) &&
            (r(1) >= ybegin_m && r(1) < yend_m) &&
            (r(2) >= zbegin_m && r(2) < zend_m));
}

#endif