sph/docs/Domain_8h_source.html

 #pragma once


 #include "math/AffineMatrix.h"

 #include "objects/geometry/Box.h"

 #include "objects/wrappers/Optional.h"

 #include "objects/wrappers/SharedPtr.h"


 NAMESPACE_SPH_BEGIN


 enum class SubsetType {

     INSIDE,

     OUTSIDE,

 };


 struct Ghost {

     Vector position;

     Size index;

 };


 class IDomain : public Polymorphic {

 public:

     virtual Vector getCenter() const = 0;


     virtual Box getBoundingBox() const = 0;


     virtual Float getVolume() const = 0;


     virtual Float getSurfaceArea() const = 0;


     virtual bool contains(const Vector& v) const = 0;


     virtual void getSubset(ArrayView<const Vector> vs, Array<Size>& output, const SubsetType type) const = 0;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const = 0;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const = 0;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta = 2._f,

         const Float eps = 0.05_f) const = 0;


 };


 class SphericalDomain : public IDomain {

 private:

     Vector center;

     Float radius;


 public:

     SphericalDomain(const Vector& center, const Float& radius);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     INLINE bool isInsideImpl(const Vector& v) const {

         return getSqrLength(v - this->center) <= sqr(radius);

     }

 };


 class EllipsoidalDomain : public IDomain {

 private:

     Vector center;


     Vector radii;


     Float effectiveRadius;


 public:

     EllipsoidalDomain(const Vector& center, const Vector& axes);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     INLINE bool isInsideImpl(const Vector& v) const {

         return getSqrLength((v - center) / radii) <= 1._f;

     }

 };


 class BlockDomain : public IDomain {

 private:

     Box box;


 public:

     BlockDomain(const Vector& center, const Vector& edges);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;

 };


 class CylindricalDomain : public IDomain {

 private:

     Vector center;

     Float radius;

     Float height;

     bool includeBases;


 public:

     CylindricalDomain(const Vector& center, const Float radius, const Float height, const bool includeBases);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     INLINE bool isInsideImpl(const Vector& v) const {

         return getSqrLength(Vector(v[X], v[Y], center[Z]) - center) <= sqr(radius) &&

                sqr(v[Z] - center[Z]) <= sqr(0.5_f * height);

     }

 };


 class HexagonalDomain : public IDomain {

 private:

     Vector center;

     Float outerRadius; // bounding radius of the base

     Float innerRadius;

     Float height;

     bool includeBases;


 public:

     HexagonalDomain(const Vector& center, const Float radius, const Float height, const bool includeBases);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector>, Array<Float>&) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     INLINE bool isInsideImpl(const Vector& v) const {

         if (sqr(v[Z] - center[Z]) > sqr(0.5_f * height)) {

             return false;

         }

         const Float sqrLength = getSqrLength(v);

         if (sqrLength > sqr(outerRadius)) {

             return false;

         }

         if (sqrLength <= sqr(innerRadius)) {

             return true;

         }

         const Float phi = atan2(v[Y], v[X]);

         return getSqrLength(Vector(v[X], v[Y], 0._f)) <= sqr(hexagon(phi));

     }


     INLINE Float hexagon(const Float phi) const {

         return 0.5_f * SQRT_3 * 1._f / sin(phi - PI / 3._f * (floor(phi / (PI / 3._f)) - 1._f));

     }

 };


 class GaussianRandomSphere : public IDomain {

 private:

     Vector center;

     Float a;

     Float beta;


 public:

     GaussianRandomSphere(const Vector& center, const Float radius, const Float beta, const Size seed);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     Float sphericalHarmonic(const Float theta, const Float phi) const;

 };


 class HalfSpaceDomain : public IDomain {

 public:

     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;

 };

 class TransformedDomain : public IDomain {

 private:

     SharedPtr<IDomain> domain;

     AffineMatrix tm, tmInv;


 public:

     TransformedDomain(SharedPtr<IDomain> domain, const AffineMatrix& matrix);


     virtual Vector getCenter() const override;


     virtual Float getVolume() const override;


     virtual Float getSurfaceArea() const override;


     virtual Box getBoundingBox() const override;


     virtual bool contains(const Vector& v) const override;


     virtual void getSubset(ArrayView<const Vector> vs,

         Array<Size>& output,

         const SubsetType type) const override;


     virtual void getDistanceToBoundary(ArrayView<const Vector> vs, Array<Float>& distances) const override;


     virtual void project(ArrayView<Vector> vs, Optional<ArrayView<Size>> indices = NOTHING) const override;


     virtual void addGhosts(ArrayView<const Vector> vs,

         Array<Ghost>& ghosts,

         const Float eta,

         const Float eps) const override;


 private:

     Array<Vector> untransform(ArrayView<const Vector> vs) const;

 };


 NAMESPACE_SPH_END

AffineMatrix.h
Three-dimensional affine matrix.

NAMESPACE_SPH_BEGIN
NAMESPACE_SPH_BEGIN
Definition: BarnesHut.cpp:13

Box.h
Object representing a three-dimensional axis-aligned box.

radius
const float radius
Definition: CurveDialog.cpp:18

SubsetType
SubsetType
Definition: Domain.h:16

SubsetType::INSIDE
@ INSIDE
Marks all vectors inside of the domain.

SubsetType::OUTSIDE
@ OUTSIDE
Marks all vectors outside of the domain.

Size
uint32_t Size
Integral type used to index arrays (by default).
Definition: Globals.h:16

Float
double Float
Precision used withing the code. Use Float instead of float or double where precision is important.
Definition: Globals.h:13

sqr
constexpr INLINE T sqr(const T &f) noexcept
Return a squared value.
Definition: MathUtils.h:67

sin
INLINE T sin(const T f)
Definition: MathUtils.h:296

floor
INLINE auto floor(const T &f)
Definition: MathUtils.h:346

atan2
INLINE T atan2(const T y, const T x)
Definition: MathUtils.h:321

SQRT_3
constexpr Float SQRT_3
Definition: MathUtils.h:367

PI
constexpr Float PI
Mathematical constants.
Definition: MathUtils.h:361

INLINE
#define INLINE
Macros for conditional compilation based on selected compiler.
Definition: Object.h:31

NAMESPACE_SPH_END
#define NAMESPACE_SPH_END
Definition: Object.h:12

Optional.h
Wrapper of type value of which may or may not be present.

NOTHING
const NothingType NOTHING
Definition: Optional.h:16

SharedPtr.h

getSqrLength
INLINE Float getSqrLength(const Vector &v)
Definition: Vector.h:574

Vector
BasicVector< Float > Vector
Definition: Vector.h:539

Y
@ Y
Definition: Vector.h:23

X
@ X
Definition: Vector.h:22

Z
@ Z
Definition: Vector.h:24

AffineMatrix
Definition: AffineMatrix.h:14

ArrayView
Object providing safe access to continuous memory of data.
Definition: ArrayView.h:17

Array< Size >

BasicVector< Float >

BlockDomain
Block aligned with coordinate axes, defined by its center and length of each side.
Definition: Domain.h:185

BlockDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:230

BlockDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:217

BlockDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:252

BlockDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:296

BlockDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:276

BlockDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:234

BlockDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:226

BlockDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:213

BlockDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:221

BlockDomain::BlockDomain
BlockDomain(const Vector &center, const Vector &edges)
Definition: Domain.cpp:210

Box
Helper object defining three-dimensional interval (box).
Definition: Box.h:17

CylindricalDomain
Cylinder aligned with z-axis, optionally including bases (can be either open or close cylinder).
Definition: Domain.h:218

CylindricalDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:368

CylindricalDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:425

CylindricalDomain::CylindricalDomain
CylindricalDomain(const Vector &center, const Float radius, const Float height, const bool includeBases)
Definition: Domain.cpp:338

CylindricalDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:351

CylindricalDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:347

CylindricalDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:355

CylindricalDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:364

CylindricalDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:400

CylindricalDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:388

CylindricalDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:359

EllipsoidalDomain
Axis aligned ellipsoidal domain, defined by the center of sphere and lengths of three axes.
Definition: Domain.h:140

EllipsoidalDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:121

EllipsoidalDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:197

EllipsoidalDomain::EllipsoidalDomain
EllipsoidalDomain(const Vector &center, const Vector &axes)
Definition: Domain.cpp:114

EllipsoidalDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:125

EllipsoidalDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:129

EllipsoidalDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:140

EllipsoidalDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:167

EllipsoidalDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:144

EllipsoidalDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:136

EllipsoidalDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:176

GaussianRandomSphere
See Muinonen 1998.
Definition: Domain.h:320

GaussianRandomSphere::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:603

GaussianRandomSphere::GaussianRandomSphere
GaussianRandomSphere(const Vector &center, const Float radius, const Float beta, const Size seed)
Definition: Domain.cpp:588

GaussianRandomSphere::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:629

GaussianRandomSphere::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:625

GaussianRandomSphere::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:599

GaussianRandomSphere::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:612

GaussianRandomSphere::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:619

GaussianRandomSphere::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:637

GaussianRandomSphere::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:633

GaussianRandomSphere::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:607

HalfSpaceDomain
Domain representing a half-space, given by z>0.
Definition: Domain.h:360

HalfSpaceDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:687

HalfSpaceDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:666

HalfSpaceDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:683

HalfSpaceDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:671

HalfSpaceDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:675

HalfSpaceDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:679

HalfSpaceDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:714

HalfSpaceDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:730

HalfSpaceDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:707

HexagonalDomain
Similar to cylindrical domain, but bases are hexagons instead of circles.
Definition: Domain.h:263

HexagonalDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:487

HexagonalDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector >, Array< Float > &) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:516

HexagonalDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:478

HexagonalDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:496

HexagonalDomain::HexagonalDomain
HexagonalDomain(const Vector &center, const Float radius, const Float height, const bool includeBases)
Definition: Domain.cpp:464

HexagonalDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:483

HexagonalDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:474

HexagonalDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:520

HexagonalDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:546

HexagonalDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:492

IDomain
Base class for computational domains.
Definition: Domain.h:29

IDomain::getSurfaceArea
virtual Float getSurfaceArea() const =0
Returns the surface area of the domain.

IDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const =0
Returns an array of indices, marking vectors with given property by their index.

IDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const =0
Returns distances of particles lying close to the boundary.

IDomain::getBoundingBox
virtual Box getBoundingBox() const =0
Returns the bounding box of the domain.

IDomain::getVolume
virtual Float getVolume() const =0
Returns the total volume of the domain.

IDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const =0
Projects vectors outside of the domain onto its boundary.

IDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta=2._f, const Float eps=0.05_f) const =0
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...

IDomain::contains
virtual bool contains(const Vector &v) const =0
Checks if the given point lies inside the domain.

IDomain::getCenter
virtual Vector getCenter() const =0
Returns the center of the domain.

Optional
Wrapper of type value of which may or may not be present.
Definition: Optional.h:23

SharedPtr< IDomain >

SphericalDomain
Spherical domain, defined by the center of sphere and its radius.
Definition: Domain.h:102

SphericalDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:22

SphericalDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:26

SphericalDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:58

SphericalDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:86

SphericalDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:35

SphericalDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:18

SphericalDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:14

SphericalDomain::SphericalDomain
SphericalDomain(const Vector &center, const Float &radius)
Definition: Domain.cpp:10

SphericalDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:31

SphericalDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:66

TransformedDomain
Transform another domain by given transformation matrix.
Definition: Domain.h:388

TransformedDomain::getSubset
virtual void getSubset(ArrayView< const Vector > vs, Array< Size > &output, const SubsetType type) const override
Returns an array of indices, marking vectors with given property by their index.
Definition: Domain.cpp:793

TransformedDomain::getVolume
virtual Float getVolume() const override
Returns the total volume of the domain.
Definition: Domain.cpp:765

TransformedDomain::getDistanceToBoundary
virtual void getDistanceToBoundary(ArrayView< const Vector > vs, Array< Float > &distances) const override
Returns distances of particles lying close to the boundary.
Definition: Domain.cpp:799

TransformedDomain::project
virtual void project(ArrayView< Vector > vs, Optional< ArrayView< Size >> indices=NOTHING) const override
Projects vectors outside of the domain onto its boundary.
Definition: Domain.cpp:803

TransformedDomain::getSurfaceArea
virtual Float getSurfaceArea() const override
Returns the surface area of the domain.
Definition: Domain.cpp:769

TransformedDomain::addGhosts
virtual void addGhosts(ArrayView< const Vector > vs, Array< Ghost > &ghosts, const Float eta, const Float eps) const override
Duplicates positions located close to the boundary, placing copies ("ghosts") symmetrically to the ot...
Definition: Domain.cpp:807

TransformedDomain::TransformedDomain
TransformedDomain(SharedPtr< IDomain > domain, const AffineMatrix &matrix)
Definition: Domain.cpp:755

TransformedDomain::getCenter
virtual Vector getCenter() const override
Returns the center of the domain.
Definition: Domain.cpp:761

TransformedDomain::contains
virtual bool contains(const Vector &v) const override
Checks if the given point lies inside the domain.
Definition: Domain.cpp:789

TransformedDomain::getBoundingBox
virtual Box getBoundingBox() const override
Returns the bounding box of the domain.
Definition: Domain.cpp:774

Ghost
Definition: Domain.h:22

Ghost::index
Size index
Index into the original array of vectors.
Definition: Domain.h:24

Ghost::position
Vector position
Position of the ghost.
Definition: Domain.h:23

Polymorphic
Base class for all polymorphic objects.
Definition: Object.h:88