| [28c025] | 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 |  | 
|---|
| [955051] | 19 | #include "Fragmentation/Summation/ZeroInstance.hpp" | 
|---|
|  | 20 |  | 
|---|
| [28c025] | 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 | 
|---|
| [3d9a8d] | 29 | * \param _end edge length per axis | 
|---|
|  | 30 | * \param _level number of gridpoints as \f$2^{\mathrm{level}}\f$ per unit(!) length | 
|---|
| [28c025] | 31 | */ | 
|---|
|  | 32 | SamplingGridProperties( | 
|---|
|  | 33 | const double _begin[3], | 
|---|
| [3d9a8d] | 34 | const double _end[3], | 
|---|
| [28c025] | 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 |  | 
|---|
| [c889b7] | 49 | /** Checks whether another instance is compatible with this one. | 
|---|
| [c0e8fb] | 50 | * | 
|---|
|  | 51 | *  \note Compatibility implies only that both grids have same grid spacing. | 
|---|
| [c889b7] | 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 | { | 
|---|
| [1a00bb] | 58 | return level == _props.level; | 
|---|
| [c889b7] | 59 | } | 
|---|
|  | 60 |  | 
|---|
| [3d9a8d] | 61 | /** Assignment operator. | 
|---|
|  | 62 | * | 
|---|
|  | 63 | * \param other other instance to assign ourselves to | 
|---|
|  | 64 | */ | 
|---|
|  | 65 | SamplingGridProperties& operator=(const SamplingGridProperties& other); | 
|---|
|  | 66 |  | 
|---|
| [c889b7] | 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 |  | 
|---|
| [b5f548] | 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 |  | 
|---|
| [28c025] | 115 | public: | 
|---|
|  | 116 | //!> offset of grid | 
|---|
|  | 117 | double begin[3]; | 
|---|
| [3d9a8d] | 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 | 
|---|
| [28c025] | 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]; | 
|---|
| [3d9a8d] | 133 | for (i=0; i<3; ++i) | 
|---|
|  | 134 | ar & end[i]; | 
|---|
| [28c025] | 135 | ar & level; | 
|---|
|  | 136 | } | 
|---|
|  | 137 |  | 
|---|
|  | 138 | }; | 
|---|
|  | 139 |  | 
|---|
| [955051] | 140 | template<> SamplingGridProperties ZeroInstance<SamplingGridProperties>(); | 
|---|
|  | 141 |  | 
|---|
| [28c025] | 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 |  | 
|---|
| [b5f548] | 146 | // define inline functions | 
|---|
|  | 147 | #include "SamplingGridProperties_inline.hpp" | 
|---|
|  | 148 |  | 
|---|
| [28c025] | 149 | #endif /* SAMPLINGGRIDPROPERTIES_HPP_ */ | 
|---|