7#include <hrleDomain.hpp>
8#include <hrleFillDomainWithSignedDistance.hpp>
12#include <vcLogger.hpp>
13#include <vcSmartPointer.hpp>
14#include <vcVectorType.hpp>
16#define LS_DOMAIN_SERIALIZATION_VERSION 0
20using namespace viennacore;
27template <
class T,
int D>
class Domain {
42 int levelSetWidth = 1;
51 static constexpr T POS_VALUE = std::numeric_limits<T>::max();
52 static constexpr T NEG_VALUE = std::numeric_limits<T>::lowest();
55 Domain(viennahrle::CoordType gridDelta = 1.0) {
56 viennahrle::IndexType gridMin[
D], gridMax[
D];
58 for (
unsigned i = 0; i <
D; ++i) {
61 boundaryCons[i] = BoundaryType::INFINITE_BOUNDARY;
64 grid =
GridType(gridMin, gridMax, gridDelta, boundaryCons);
69 viennahrle::CoordType gridDelta = 1.0) {
70 viennahrle::IndexType gridMin[
D], gridMax[
D];
71 for (
unsigned i = 0; i <
D; ++i) {
72 gridMin[i] = std::floor(bounds[2 * i] / gridDelta);
73 gridMax[i] = std::ceil(bounds[2 * i + 1] / gridDelta);
76 grid =
GridType(gridMin, gridMax, gridDelta, boundaryConditions);
80 Domain(std::vector<viennahrle::CoordType> bounds,
81 std::vector<unsigned> boundaryConditions,
82 viennahrle::CoordType gridDelta = 1.0) {
84 for (
unsigned i = 0; i <
D; ++i) {
85 boundaryCons[i] =
static_cast<BoundaryType>(boundaryConditions[i]);
88 SmartPointer<Domain<T, D>>
::New(bounds.data(), boundaryCons, gridDelta);
96 viennahrle::CoordType gridDelta = 1.0) {
98 SmartPointer<Domain<T, D>>
::New(bounds, boundaryConditions, gridDelta);
100 viennahrle::FillDomainWithSignedDistance(domain, pointData,
T(
NEG_VALUE),
111 template <
class... Args>
static auto New(Args &&...args) {
112 return SmartPointer<Domain>::New(std::forward<Args>(args)...);
117 void finalize(
int newWidth) { levelSetWidth = newWidth; }
125 grid = passedDomain->grid;
126 domain.deepCopy(grid, passedDomain->domain);
127 levelSetWidth = passedDomain->levelSetWidth;
128 pointData = passedDomain->pointData;
136 viennahrle::FillDomainWithSignedDistance(domain, pointData,
T(
NEG_VALUE),
174 return voidPointMarkers;
178 void print(std::ostream &out = std::cout) {
179 out <<
"Grid pointer: " << &grid << std::endl;
180 out <<
"lsDomain: " << &domain << std::endl;
181 out <<
"DomainSegments: " << std::endl;
183 out << &(domain.getDomainSegment(i)) << std::endl;
191 stream <<
"lsDomain";
195 stream.write(&formatVersion, 1);
198 grid.serialize(stream);
201 domain.serialize(stream);
205 const uint32_t width = levelSetWidth;
206 stream.write(
reinterpret_cast<const char *
>(&width),
sizeof(uint32_t));
210 char hasPointData = (pointData.empty()) ? 0 : 1;
211 stream.write(&hasPointData, 1);
212 if (hasPointData == 1) {
213 pointData.serialize(stream);
223 stream.read(identifier, 8);
224 if (std::string(identifier).compare(0, 8,
"lsDomain")) {
225 Logger::getInstance()
227 "Reading Domain from stream failed. Header could not be found.")
234 stream.read(&formatVersion, 1);
236 Logger::getInstance()
238 "Reading Domain of version " + std::to_string(formatVersion) +
239 " with reader of version " +
246 grid.deserialize(stream);
252 domain.deserialize(stream);
256 stream.read(
reinterpret_cast<char *
>(&width),
sizeof(uint32_t));
257 levelSetWidth = width;
261 stream.read(&hasPointData, 1);
262 if (hasPointData == 1) {
264 pointData.deserialize(stream);
std::vector< std::pair< viennahrle::Index< D >, T > > PointValueVectorType
Definition lsDomain.hpp:34
Domain(std::vector< viennahrle::CoordType > bounds, std::vector< unsigned > boundaryConditions, viennahrle::CoordType gridDelta=1.0)
Definition lsDomain.hpp:80
viennahrle::Domain< T, D > DomainType
Definition lsDomain.hpp:32
Domain(PointValueVectorType pointData, viennahrle::CoordType *bounds, BoundaryType *boundaryConditions, viennahrle::CoordType gridDelta=1.0)
initialise Domain with domain size "bounds", filled with point/value pairs in pointData
Definition lsDomain.hpp:94
Domain(const viennahrle::CoordType *bounds, BoundaryType *boundaryConditions, viennahrle::CoordType gridDelta=1.0)
Definition lsDomain.hpp:68
viennahrle::Grid< D > GridType
Definition lsDomain.hpp:31
static constexpr T NEG_VALUE
Definition lsDomain.hpp:52
T ValueType
Definition lsDomain.hpp:30
void clearMetaData()
Definition lsDomain.hpp:163
void deepCopy(const SmartPointer< Domain< T, D > > passedDomain)
copy all values of "passedDomain" to this Domain
Definition lsDomain.hpp:124
unsigned getNumberOfSegments() const
returns the number of segments, the levelset is split into. This is useful for algorithm parallelisat...
Definition lsDomain.hpp:153
void finalize()
this function finalizes the levelset, so it is ready for use by other algorithms
Definition lsDomain.hpp:121
void print(std::ostream &out=std::cout)
prints basic information and all memebers of the levelset structure
Definition lsDomain.hpp:178
std::vector< bool > VoidPointMarkersType
Definition lsDomain.hpp:36
PointData< T > PointDataType
Definition lsDomain.hpp:35
void finalize(int newWidth)
this function sets a new levelset width and finalizes the levelset, so it is ready for use by other a...
Definition lsDomain.hpp:117
Domain(SmartPointer< Domain > passedDomain)
Definition lsDomain.hpp:108
static constexpr T POS_VALUE
Definition lsDomain.hpp:51
BoundaryConditionEnum BoundaryType
Definition lsDomain.hpp:33
const DomainType & getDomain() const
Definition lsDomain.hpp:149
VoidPointMarkersType & getVoidPointMarkers()
get reference to the voidPoints markers for all points
Definition lsDomain.hpp:171
const VoidPointMarkersType & getVoidPointMarkers() const
Definition lsDomain.hpp:173
const PointDataType & getPointData() const
Definition lsDomain.hpp:168
PointDataType & getPointData()
get reference to point data saved in the level set
Definition lsDomain.hpp:166
DomainType & getDomain()
get const reference to the underlying hrleDomain data structure
Definition lsDomain.hpp:147
Domain(GridType passedGrid)
Definition lsDomain.hpp:104
int getLevelSetWidth() const
Definition lsDomain.hpp:158
Domain(viennahrle::CoordType gridDelta=1.0)
initalise an empty infinite Domain
Definition lsDomain.hpp:55
std::istream & deserialize(std::istream &stream)
Deserialize Domain from binary stream.
Definition lsDomain.hpp:220
GridType & getGrid()
get mutable reference to the grid on which the level set is defined
Definition lsDomain.hpp:144
static auto New(Args &&...args)
Definition lsDomain.hpp:111
const GridType & getGrid() const
get reference to the grid on which the levelset is defined
Definition lsDomain.hpp:141
unsigned getNumberOfPoints() const
returns the number of defined points
Definition lsDomain.hpp:156
std::ostream & serialize(std::ostream &stream)
Serializes the Domain into a binary stream.
Definition lsDomain.hpp:189
void setLevelSetWidth(int width)
Definition lsDomain.hpp:160
void insertPoints(PointValueVectorType pointData, bool sort=true)
re-initalise Domain with the point/value pairs in pointData This is similar to lsFromMesh with the di...
Definition lsDomain.hpp:135
static constexpr int dimensions
Definition lsDomain.hpp:48
This class holds data associated with points in space.
Definition lsPointData.hpp:21
#define LS_DOMAIN_SERIALIZATION_VERSION
Definition lsDomain.hpp:16
#define PRECOMPILE_PRECISION_DIMENSION(className)
Definition lsPreCompileMacros.hpp:24
Definition lsAdvect.hpp:36
viennahrle::BoundaryType BoundaryConditionEnum
Definition lsDomain.hpp:23
constexpr int D
Definition pyWrap.cpp:71
double T
Definition pyWrap.cpp:69