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