| 1 | /* | 
|---|
| 2 | * SamplingGridProperties.hpp | 
|---|
| 3 | * | 
|---|
| 4 | *  Created on: 25.07.2012 | 
|---|
| 5 | *      Author: heber | 
|---|
| 6 | */ | 
|---|
| 7 |  | 
|---|
| 8 | #ifndef SAMPLINGGRIDPROPERTIES_HPP_ | 
|---|
| 9 | #define SAMPLINGGRIDPROPERTIES_HPP_ | 
|---|
| 10 |  | 
|---|
| 11 | // include config.h | 
|---|
| 12 | #ifdef HAVE_CONFIG_H | 
|---|
| 13 | #include <config.h> | 
|---|
| 14 | #endif | 
|---|
| 15 |  | 
|---|
| 16 | #include "boost/serialization/export.hpp" | 
|---|
| 17 | #include "boost/serialization/array.hpp" | 
|---|
| 18 |  | 
|---|
| 19 | #include "Fragmentation/Summation/ZeroInstance.hpp" | 
|---|
| 20 |  | 
|---|
| 21 | /** This class stores a sample function on a three-dimensional grid. | 
|---|
| 22 | * | 
|---|
| 23 | */ | 
|---|
| 24 | class SamplingGridProperties { | 
|---|
| 25 | public: | 
|---|
| 26 | /** Constructor for class SamplingGridProperties. | 
|---|
| 27 | * | 
|---|
| 28 | * \param _begin offset of grid | 
|---|
| 29 | * \param _end edge length per axis | 
|---|
| 30 | * \param _level number of gridpoints as \f$2^{\mathrm{level}}\f$ per unit(!) length | 
|---|
| 31 | */ | 
|---|
| 32 | SamplingGridProperties( | 
|---|
| 33 | const double _begin[3], | 
|---|
| 34 | const double _end[3], | 
|---|
| 35 | const int _level); | 
|---|
| 36 |  | 
|---|
| 37 | /** Copy constructor for class SamplingGridProperties. | 
|---|
| 38 | * | 
|---|
| 39 | * \param _props instance to copy from | 
|---|
| 40 | */ | 
|---|
| 41 | SamplingGridProperties(const SamplingGridProperties &_props); | 
|---|
| 42 |  | 
|---|
| 43 | /** Default constructor. | 
|---|
| 44 | */ | 
|---|
| 45 | SamplingGridProperties(); | 
|---|
| 46 |  | 
|---|
| 47 | virtual ~SamplingGridProperties(); | 
|---|
| 48 |  | 
|---|
| 49 | /** Checks whether another instance is compatible with this one. | 
|---|
| 50 | * | 
|---|
| 51 | *  \note Compatibility implies only that both grids have same grid spacing. | 
|---|
| 52 | * | 
|---|
| 53 | * \param _props other properties to check against | 
|---|
| 54 | * \return true - are compatible, false - else | 
|---|
| 55 | */ | 
|---|
| 56 | bool isCompatible(const SamplingGridProperties &_props) const | 
|---|
| 57 | { | 
|---|
| 58 | return level == _props.level; | 
|---|
| 59 | } | 
|---|
| 60 |  | 
|---|
| 61 | /** Assignment operator. | 
|---|
| 62 | * | 
|---|
| 63 | * \param other other instance to assign ourselves to | 
|---|
| 64 | */ | 
|---|
| 65 | SamplingGridProperties& operator=(const SamplingGridProperties& other); | 
|---|
| 66 |  | 
|---|
| 67 | /** Equality operator for class SamplingGridProperties. | 
|---|
| 68 | * | 
|---|
| 69 | * \param _props other object to compare to | 
|---|
| 70 | */ | 
|---|
| 71 | bool operator==(const SamplingGridProperties &_props) const; | 
|---|
| 72 |  | 
|---|
| 73 | /** Inequality operator for class SamplingGridProperties. | 
|---|
| 74 | * | 
|---|
| 75 | * \param _props other object to compare to | 
|---|
| 76 | */ | 
|---|
| 77 | bool operator!=(const SamplingGridProperties &_props) const | 
|---|
| 78 | { | 
|---|
| 79 | return (!(*this == _props)); | 
|---|
| 80 | } | 
|---|
| 81 |  | 
|---|
| 82 | /** Returns the volume of the domain for this sampled function. | 
|---|
| 83 | * | 
|---|
| 84 | * @return volume | 
|---|
| 85 | */ | 
|---|
| 86 | const double getVolume() const; | 
|---|
| 87 |  | 
|---|
| 88 | /** Returns the total number of gridpoints of the discrete mesh covering the volume. | 
|---|
| 89 | * | 
|---|
| 90 | * @return number of gridpoints sampled_values should have | 
|---|
| 91 | */ | 
|---|
| 92 | const size_t getTotalGridPoints() const; | 
|---|
| 93 |  | 
|---|
| 94 | /** Returns the number of grid points per axis. | 
|---|
| 95 | * | 
|---|
| 96 | * @return number of grid points per unit length | 
|---|
| 97 | */ | 
|---|
| 98 | const size_t getGridPointsPerAxis() const; | 
|---|
| 99 |  | 
|---|
| 100 | /** Returns the length of the domain for the given \a axis. | 
|---|
| 101 | * | 
|---|
| 102 | * \param axis axis for which to get step length | 
|---|
| 103 | * \return domain length for the given axis, i.e. end - begin | 
|---|
| 104 | */ | 
|---|
| 105 | const double getTotalLengthPerAxis(const size_t axis) const; | 
|---|
| 106 |  | 
|---|
| 107 | /** Returns the real step length from one discrete grid point to the next. | 
|---|
| 108 | * | 
|---|
| 109 | * \param axis axis for which to get step length | 
|---|
| 110 | * \return step length for the given axis, as domain length over getGridPointsPerAxis() | 
|---|
| 111 | */ | 
|---|
| 112 | const double getDeltaPerAxis(const size_t axis) const; | 
|---|
| 113 |  | 
|---|
| 114 |  | 
|---|
| 115 | public: | 
|---|
| 116 | //!> offset of grid | 
|---|
| 117 | double begin[3]; | 
|---|
| 118 | //!> size of grid, i.e. edge length per axis of domain | 
|---|
| 119 | double end[3]; | 
|---|
| 120 | //!> level, i.e. \f$2^{\mathrm{level}}\f$ grid points per axis per unit(!) length | 
|---|
| 121 | int level; | 
|---|
| 122 |  | 
|---|
| 123 | private: | 
|---|
| 124 |  | 
|---|
| 125 | friend class boost::serialization::access; | 
|---|
| 126 | // serialization | 
|---|
| 127 | template <typename Archive> | 
|---|
| 128 | void serialize(Archive& ar, const unsigned int version) | 
|---|
| 129 | { | 
|---|
| 130 | int i; | 
|---|
| 131 | for (i=0; i<3; ++i) | 
|---|
| 132 | ar & begin[i]; | 
|---|
| 133 | for (i=0; i<3; ++i) | 
|---|
| 134 | ar & end[i]; | 
|---|
| 135 | ar & level; | 
|---|
| 136 | } | 
|---|
| 137 |  | 
|---|
| 138 | }; | 
|---|
| 139 |  | 
|---|
| 140 | template<> SamplingGridProperties ZeroInstance<SamplingGridProperties>(); | 
|---|
| 141 |  | 
|---|
| 142 | // we need to give this class a unique key for serialization | 
|---|
| 143 | // its is only serialized through its base class FragmentJob | 
|---|
| 144 | BOOST_CLASS_EXPORT_KEY(SamplingGridProperties) | 
|---|
| 145 |  | 
|---|
| 146 | // define inline functions | 
|---|
| 147 | #include "SamplingGridProperties_inline.hpp" | 
|---|
| 148 |  | 
|---|
| 149 | #endif /* SAMPLINGGRIDPROPERTIES_HPP_ */ | 
|---|