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),
112 void finalize(
int newWidth) { levelSetWidth = newWidth; }
120 grid = passedDomain->grid;
121 domain.deepCopy(grid, passedDomain->domain);
122 levelSetWidth = passedDomain->levelSetWidth;
123 pointData = passedDomain->pointData;
131 viennahrle::FillDomainWithSignedDistance(domain, pointData,
T(
NEG_VALUE),
169 return voidPointMarkers;
173 void print(std::ostream &out = std::cout) {
174 out <<
"Grid pointer: " << &grid << std::endl;
175 out <<
"lsDomain: " << &domain << std::endl;
176 out <<
"DomainSegments: " << std::endl;
178 out << &(domain.getDomainSegment(i)) << std::endl;
186 stream <<
"lsDomain";
190 stream.write(&formatVersion, 1);
193 grid.serialize(stream);
196 domain.serialize(stream);
200 const uint32_t width = levelSetWidth;
201 stream.write(
reinterpret_cast<const char *
>(&width),
sizeof(uint32_t));
205 char hasPointData = (pointData.empty()) ? 0 : 1;
206 stream.write(&hasPointData, 1);
207 if (hasPointData == 1) {
208 pointData.serialize(stream);
218 stream.read(identifier, 8);
219 if (std::string(identifier).compare(0, 8,
"lsDomain")) {
220 Logger::getInstance()
222 "Reading Domain from stream failed. Header could not be found.")
229 stream.read(&formatVersion, 1);
231 Logger::getInstance()
233 "Reading Domain of version " + std::to_string(formatVersion) +
234 " with reader of version " +
241 grid.deserialize(stream);
247 domain.deserialize(stream);
251 stream.read(
reinterpret_cast<char *
>(&width),
sizeof(uint32_t));
252 levelSetWidth = width;
256 stream.read(&hasPointData, 1);
257 if (hasPointData == 1) {
259 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:158
void deepCopy(const SmartPointer< Domain< T, D > > passedDomain)
copy all values of "passedDomain" to this Domain
Definition lsDomain.hpp:119
unsigned getNumberOfSegments() const
returns the number of segments, the levelset is split into. This is useful for algorithm parallelisat...
Definition lsDomain.hpp:148
void finalize()
this function finalizes the levelset, so it is ready for use by other algorithms
Definition lsDomain.hpp:116
void print(std::ostream &out=std::cout)
prints basic information and all memebers of the levelset structure
Definition lsDomain.hpp:173
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:112
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:144
VoidPointMarkersType & getVoidPointMarkers()
get reference to the voidPoints markers for all points
Definition lsDomain.hpp:166
const VoidPointMarkersType & getVoidPointMarkers() const
Definition lsDomain.hpp:168
const PointDataType & getPointData() const
Definition lsDomain.hpp:163
PointDataType & getPointData()
get reference to point data saved in the level set
Definition lsDomain.hpp:161
DomainType & getDomain()
get const reference to the underlying hrleDomain data structure
Definition lsDomain.hpp:142
Domain(GridType passedGrid)
Definition lsDomain.hpp:104
int getLevelSetWidth() const
Definition lsDomain.hpp:153
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:215
GridType & getGrid()
get mutable reference to the grid on which the level set is defined
Definition lsDomain.hpp:139
const GridType & getGrid() const
get reference to the grid on which the levelset is defined
Definition lsDomain.hpp:136
unsigned getNumberOfPoints() const
returns the number of defined points
Definition lsDomain.hpp:151
std::ostream & serialize(std::ostream &stream)
Serializes the Domain into a binary stream.
Definition lsDomain.hpp:184
void setLevelSetWidth(int width)
Definition lsDomain.hpp:155
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:130
static constexpr int dimensions
Definition lsDomain.hpp:48
This class holds data associated with points in space.
Definition lsPointData.hpp:20
#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:70
double T
Definition pyWrap.cpp:68