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