| 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_ */
|
|---|