A TriangleWall is convex polygon defined as an intersection of InfiniteWall's. More...

#include <TriangleWall.h>

Inheritance diagram for TriangleWall:

Public Member Functions
	TriangleWall ()=default
	Default constructor. More...

	TriangleWall (const TriangleWall &other)=default
	Copy constructor. More...

	~TriangleWall () override=default
	Destructor. More...

TriangleWall *	copy () const override
	Wall copy method. It calls the copy constructor of this Wall, useful for polymorphism. More...

std::string	getName () const override
	Returns the name of the object, here the string "TriangleWall". More...

void	read (std::istream &is) override
	Reads an TriangleWall from an input stream, for example a restart file. More...

void	write (std::ostream &os) const override
	Writes an TriangleWall to an output stream, for example a restart file. More...

void	setVertices (Vec3D A, Vec3D B, Vec3D C)
	Sets member variables such that the wall represents a triangle with vertices A, B, C. More...

std::array< Vec3D, 3 >	getVertices () const

void	move (const Vec3D &move) override

void	setVertices (Vec3D A, Vec3D B, Vec3D C, Vec3D position)
	Same as #setVertices(A,B,C), but sets the position explicitly. The position is important when you rotate the wall, as the wall will be rotated around this position. More...

void	writeVTK (VTKContainer &vtk) const override

const Vec3D &	getVertex (unsigned i) const
	Returns the position of a vertex. More...

bool	getDistanceAndNormal (const BaseParticle &p, Mdouble &distance, Vec3D &normal_return) const override
	Compute the distance from the wall for a given BaseParticle and return if there is a collision. If there is a collision, also return the normal vector. More...

void	rotate (const Vec3D &angularVelocity) override
	Rotates this BaseInteractable. More...

bool	isLocal (Vec3D &min, Vec3D &max) const override

bool	isInsideTriangle (const Vec3D &point) const

Public Member Functions inherited from BaseWall
	BaseWall ()
	Default constructor. More...

	BaseWall (const BaseWall &w)
	Copy constructor. More...

	~BaseWall () override
	Default destructor. More...

virtual bool	getDistanceNormalOverlap (const BaseParticle &P, Mdouble &distance, Vec3D &normal_return, Mdouble &overlap) const

virtual bool	getDistanceNormalOverlapSuperquadric (const SuperQuadricParticle &p, Mdouble &distance, Vec3D &normal_return, Mdouble &overlap) const

virtual Vec3D	getFurthestPointSuperQuadric (const Vec3D &normalBodyFixed, const Vec3D &axes, Mdouble eps1, Mdouble eps2) const

virtual void	setHandler (WallHandler *handler)
	A function which sets the WallHandler for this BaseWall. More...

WallHandler *	getHandler () const
	A function which returns the WallHandler that handles this BaseWall. More...

void	setIndSpecies (unsigned int indSpecies) override
	Define the species of this wall using the index of the species in the SpeciesHandler in this DPMBase. More...

void	setSpecies (const ParticleSpecies *species)
	Defines the species of the current wall. More...

bool	isFixed () const override

void	setForceControl (Vec3D forceGoal, Vec3D gainFactor, Vec3D baseVelocity={0, 0, 0})
	Slowly adjusts the force on a wall towards a specified goal, by adjusting (prescribing) the velocity of the wall. More...

bool	getLinePlaneIntersect (Vec3D &intersect, const Vec3D &p0, const Vec3D &p1, const Vec3D &n, const Vec3D &p)

bool	isInsideWallVTK (const Vec3D &point, const Vec3D &normal, const Vec3D &position) const

void	projectOntoWallVTK (Vec3D &point0, const Vec3D &point1, const Vec3D &normal, const Vec3D &position) const

void	intersectVTK (std::vector< Vec3D > &points, Vec3D normal, Vec3D position) const

virtual BaseInteraction *	getInteractionWithSuperQuad (SuperQuadricParticle p, unsigned timeStamp, InteractionHandler interactionHandler)

void	getVTK (std::vector< Vec3D > &points, std::vector< std::vector< double >> &triangleStrips)

const Vec3D	getAxis () const

BaseInteraction *	getInteractionWith (BaseParticle p, unsigned timeStamp, InteractionHandler interactionHandler) override
	Returns the interaction between this wall and a given particle, nullptr if there is no interaction. More...

bool	getVTKVisibility () const

void	setVTKVisibility (bool vtkVisibility)

void	addRenderedWall (BaseWall *w)

BaseWall *	getRenderedWall (size_t i) const

void	removeRenderedWalls ()

void	renderWall (VTKContainer &vtk)

void	addParticlesAtWall (unsigned numElements=50)

void	setVelocityControl (Vec3D forceGoal, Vec3D gainFactor, Vec3D baseVelocity)

Public Member Functions inherited from BaseInteractable
	BaseInteractable ()
	Default BaseInteractable constructor. More...

	BaseInteractable (const BaseInteractable &p)
	Copy constructor. More...

	~BaseInteractable () override
	Destructor, it simply destructs the BaseInteractable and all the objects it contains. More...

unsigned int	getIndSpecies () const
	Returns the index of the species associated with the interactable object. More...

const ParticleSpecies *	getSpecies () const
	Returns a pointer to the species of this BaseInteractable. More...

void	setSpecies (const ParticleSpecies *species)
	Sets the species of this BaseInteractable. More...

const Vec3D &	getForce () const
	Returns the force on this BaseInteractable. More...

const Vec3D &	getTorque () const
	Returns the torque on this BaseInteractable. More...

void	setForce (const Vec3D &force)
	Sets the force on this BaseInteractable. More...

void	setTorque (const Vec3D &torque)
	Sets the torque on this BaseInteractable. More...

void	addForce (const Vec3D &addForce)
	Adds an amount to the force on this BaseInteractable. More...

void	addTorque (const Vec3D &addTorque)
	Adds an amount to the torque on this BaseInteractable. More...

void	resetForceTorque (int numberOfOMPthreads)

void	sumForceTorqueOMP ()

const Vec3D &	getPosition () const
	Returns the position of this BaseInteractable. More...

const Quaternion &	getOrientation () const
	Returns the orientation of this BaseInteractable. More...

void	setPosition (const Vec3D &position)
	Sets the position of this BaseInteractable. More...

void	setOrientationViaNormal (Vec3D normal)
	Sets the orientation of this BaseInteractable by defining the vector that results from the rotation of the (1,0,0) vector. More...

void	setOrientationViaEuler (Vec3D eulerAngle)
	Sets the orientation of this BaseInteractable by defining the euler angles. More...

void	setOrientation (const Quaternion &orientation)
	Sets the orientation of this BaseInteractable. More...

const std::vector < BaseInteraction * > &	getInteractions () const
	Returns a list of interactions which belong to this interactable. More...

void	addInteraction (BaseInteraction *I)
	Adds an interaction to this BaseInteractable. More...

bool	removeInteraction (BaseInteraction *I)
	Removes an interaction from this BaseInteractable. More...

void	copyInteractionsForPeriodicParticles (const BaseInteractable &p)
	Copies interactions to this BaseInteractable whenever a periodic copy made. More...

void	setVelocity (const Vec3D &velocity)
	set the velocity of the BaseInteractable. More...

void	setAngularVelocity (const Vec3D &angularVelocity)
	set the angular velocity of the BaseInteractble. More...

void	addVelocity (const Vec3D &velocity)
	adds an increment to the velocity. More...

void	addAngularVelocity (const Vec3D &angularVelocity)
	add an increment to the angular velocity. More...

virtual const Vec3D &	getVelocity () const
	Returns the velocity of this interactable. More...

virtual const Vec3D &	getAngularVelocity () const
	Returns the angular velocity of this interactable. More...

void	setPrescribedPosition (const std::function< Vec3D(double)> &prescribedPosition)
	Allows the position of an infinite mass interactable to be prescribed. More...

void	applyPrescribedPosition (double time)
	Computes the position from the user defined prescribed position function. More...

void	setPrescribedVelocity (const std::function< Vec3D(double)> &prescribedVelocity)
	Allows the velocity of an infinite mass interactable to be prescribed. More...

void	applyPrescribedVelocity (double time)
	Computes the velocity from the user defined prescribed velocity function. More...

void	setPrescribedOrientation (const std::function< Quaternion(double)> &prescribedOrientation)
	Allows the orientation of the infinite mass interactbale to be prescribed. More...

void	applyPrescribedOrientation (double time)
	Computes the orientation from the user defined prescribed orientation function. More...

void	setPrescribedAngularVelocity (const std::function< Vec3D(double)> &prescribedAngularVelocity)
	Allows the angular velocity of the infinite mass interactable to be prescribed. More...

void	applyPrescribedAngularVelocity (double time)
	Computes the angular velocity from the user defined prescribed angular velocity. More...

virtual const Vec3D	getVelocityAtContact (const Vec3D &contact) const
	Returns the velocity at the contact point, use by many force laws. More...

void	integrateBeforeForceComputation (double time, double timeStep)
	This is part of integrate routine for objects with infinite mass. More...

void	integrateAfterForceComputation (double time, double timeStep)
	This is part of the integration routine for objects with infinite mass. More...

virtual Mdouble	getInvMass () const

virtual Mdouble	getCurvature (const Vec3D &labFixedCoordinates) const

Public Member Functions inherited from BaseObject
	BaseObject ()=default
	Default constructor. More...

	BaseObject (const BaseObject &p)=default
	Copy constructor, copies all the objects BaseObject contains. More...

virtual	~BaseObject ()=default
	virtual destructor More...

virtual void	moveInHandler (unsigned int index)
	Except that it is virtual, it does the same thing as setIndex() does. More...

void	setIndex (unsigned int index)
	Allows one to assign an index to an object in the handler/container. More...

void	setId (unsigned long id)
	Assigns a unique identifier to each object in the handler (container) which remains constant even after the object is deleted from the container/handler. More...

unsigned int	getIndex () const
	Returns the index of the object in the handler. More...

unsigned int	getId () const
	Returns the unique identifier of any particular object. More...

void	setGroupId (unsigned groupId)

unsigned	getGroupId () const

Private Member Functions
void	updateVertexAndNormal ()
	This function should be called after setting either position_ or vertexInLabFrame_. More...

Private Attributes
std::array< Vec3D, 3 >	vertexInLabFrame_

std::array< Vec3D, 3 >	vertex_

Vec3D	vertexMin_

Vec3D	vertexMax_

std::array< Vec3D, 3 >	edgeNormal_

std::array< Vec3D, 3 >	edge_

std::array< double, 3 >	edgeLength_

Vec3D	faceNormal_

Additional Inherited Members
Static Public Member Functions inherited from BaseWall
static void	addToVTK (const std::vector< Vec3D > &points, VTKContainer &vtk)
	Takes the points provided and adds a triangle strip connecting these points to the vtk container. More...

Detailed Description

A TriangleWall is convex polygon defined as an intersection of InfiniteWall's.

It can be defined as the intersection of a set of InfiniteWall's, defined by the normal vector into the wall and a point on the wall:

TriangleWall w;
//for each wall, specify a normal and position vector
w.addObject(Vec3D(-1.0,0.0,0.0),Vec3D(1.0,0.0,0.0));
w.addObject(Vec3D(1.0,0.0,0.0),Vec3D(0.0,0.0,0.0));
w.addObject(Vec3D(0.0,-1.0,0.0),Vec3D(0.0,1.0,0.0));
w.addObject(Vec3D(0.0,1.0,0.0),Vec3D(0.0,0.0,0.0));
wallHandler.copyAndAddObject(w);

A particle of radius r and position x touches an InfiniteWall with normal n and position p if \(p-n\cdot x\leq r\) (note 'touching particles' also includes particles that are completely enclosed inside the wall). A particle touches an TriangleWall if it touches all InfiniteWall objects (shown in the image below).

For a demonstration on how to use this class, see T8: Motion of a particle in a box with an obstacle and Flow through a 3D hourglass/silo (shown in the image below).

Definition at line 56 of file TriangleWall.h.

Constructor & Destructor Documentation

TriangleWall::TriangleWall ( )

default

Default constructor.

Referenced by copy().

TriangleWall::TriangleWall ( const TriangleWall & other )

default

Copy constructor.

TriangleWall::~TriangleWall ( )

overridedefault

Destructor.

Member Function Documentation

TriangleWall* TriangleWall::copy ( ) const

inlineoverridevirtual

Wall copy method. It calls the copy constructor of this Wall, useful for polymorphism.

Implements BaseWall.

Definition at line 78 of file TriangleWall.h.

References TriangleWall().

79 { return new TriangleWall(*this); }

TriangleWall::TriangleWall

TriangleWall()=default

Default constructor.

bool TriangleWall::getDistanceAndNormal	(	const BaseParticle &	p,
		Mdouble &	distance,
		Vec3D &	normal
	)		const

overridevirtual

Compute the distance from the wall for a given BaseParticle and return if there is a collision. If there is a collision, also return the normal vector.

Parameters

[in]	p	BaseParticle we want to calculate the distance and whether it collided of.
[out]	distance	The distance of the BaseParticle to this wall.
[out]	normal_return	If there was a collision, the normal vector to this wall will be placed here.

Returns: A boolean which says whether or not there was a collision.

This function computes whether or not there is a collision between a given BaseParticle and this TriangleWall. If there is a collision, this function also computes the distance between the BaseParticle and TriangleWall and the normal of the TriangleWall at the intersection point. It does this by calling TriangleWall::getDistanceAndNormal(const Vec3D& , Mdouble , Mdouble&, Vec3D&) const. Since this function should be called before calculating any Particle-Wall interactions, it can also be used to set the normal vector in case of curved walls.

Implements BaseWall.

Definition at line 57 of file TriangleWall.cc.

References Vec3D::dot(), edge_, edgeLength_, edgeNormal_, faceNormal_, Vec3D::getLength(), BaseInteractable::getPosition(), BaseParticle::getWallInteractionRadius(), isInsideTriangle(), and vertex_.

 {
     const Vec3D position = p.getPosition(); // note: no transfer to lab coordinates
     const Mdouble distanceMax = p.getWallInteractionRadius(this);
 
     // compute distance from face
     //get distance from particle position to the face
     const Mdouble signedDistance = Vec3D::dot(position-vertex_[0], faceNormal_);
     distance = fabs(signedDistance);
 
     // check if any contact is possible
     if (distance >= distanceMax) return false;
 
     // compute distance from edges
     const std::array<Vec3D,3> edgeBranch {position - vertex_[0], position - vertex_[1], position - vertex_[2]};
     const std::array<double,3> edgeDistance {Vec3D::dot(edgeBranch[0], edgeNormal_[0]), Vec3D::dot(edgeBranch[1], edgeNormal_[1]), Vec3D::dot(edgeBranch[2], edgeNormal_[2])};
 
     // find edge with largest distance (this will be the edge if there is a edge contact)
     const size_t id = (edgeDistance[1] > edgeDistance[2]) ?
             (edgeDistance[0] > edgeDistance[1] ? 0 : 1) : (edgeDistance[0] > edgeDistance[2] ? 0 : 2);
 
     // check if there will be contact
     if (edgeDistance[id] > distanceMax) return false;
 
     // determine if it is a face contact
     const Vec3D posProjected = position - signedDistance * faceNormal_;
     if (edgeDistance[id] <= 0 && isInsideTriangle(posProjected)){
         normal = (signedDistance >= 0) ? -faceNormal_ : faceNormal_;
         return true;
     }
 
     // determine if it is an edge or vertex contact
     const double positionAlongEdge = Vec3D::dot(edgeBranch[id], edge_[id]);
     if (positionAlongEdge < 0) {
         //possible contact with left vertex
         distance = edgeBranch[id].getLength();
         if (distance > distanceMax) return false;
         normal = edgeBranch[id] / -distance;
     } else if (positionAlongEdge > edgeLength_[id]) {
         //contact with right vertex
         const size_t idRight = (id + 1) % 3;
         distance = edgeBranch[idRight].getLength();
         if (distance > distanceMax) return false;
         normal = edgeBranch[idRight] / -distance;
     } else {
         // edge contact
         normal = edge_[id] * positionAlongEdge - edgeBranch[id];
         distance = normal.getLength();
         if (distance > distanceMax) return false;
         normal /= distance;
     }
     return true;
 }

std::string TriangleWall::getName ( ) const

inlineoverridevirtual

Returns the name of the object, here the string "TriangleWall".

Implements BaseObject.

Definition at line 84 of file TriangleWall.h.

85 { return "TriangleWall"; }

const Vec3D& TriangleWall::getVertex ( unsigned i ) const

inline

Returns the position of a vertex.

Definition at line 120 of file TriangleWall.h.

References constants::i, and vertex_.

120 {return vertex_[i];}

constants::i

const std::complex< Mdouble > i

Definition: ExtendedMath.h:50

TriangleWall::vertex_

std::array< Vec3D, 3 > vertex_

Definition: TriangleWall.h:146

std::array<Vec3D,3> TriangleWall::getVertices ( ) const

inline

Definition at line 105 of file TriangleWall.h.

References vertex_.

105 {return vertex_;}

TriangleWall::vertex_

std::array< Vec3D, 3 > vertex_

Definition: TriangleWall.h:146

bool TriangleWall::isInsideTriangle ( const Vec3D & point ) const

Definition at line 238 of file TriangleWall.cc.

References Vec3D::cross(), and vertex_.

Referenced by getDistanceAndNormal().

 {
     const Vec3D branch0 = point - vertex_[0];
     const Vec3D branch1 = point - vertex_[1];
     const Vec3D branch2 = point - vertex_[2];
 
     //compute total area
     Mdouble s = sqrt(Vec3D::cross(vertex_[1] - vertex_[0], vertex_[2] - vertex_[1]).getLengthSquared());
     //compute barycentric coordinates
     Mdouble s0 = sqrt(Vec3D::cross(branch0, branch1).getLengthSquared())/s;
     Mdouble s1 = sqrt(Vec3D::cross(branch1, branch2).getLengthSquared())/s;
     Mdouble s2 = sqrt(Vec3D::cross(branch2, branch0).getLengthSquared())/s;
     return (1 > s0 > 0 && 1 > s1 > 0 && 1 > s2 > 0);
 }

bool TriangleWall::isLocal	(	Vec3D &	min,
		Vec3D &	max
	)		const

overridevirtual

if isLocal returns true and the DPM class derived from MercuryBase, the hGrid will be used to find wall-particle contacts, using min/max.

Reimplemented from BaseWall.

Definition at line 231 of file TriangleWall.cc.

References vertexMax_, and vertexMin_.

 {
     min = vertexMin_;
     max = vertexMax_;
     return true;
 }

void TriangleWall::move ( const Vec3D & move )

overridevirtual

Moves (displaces) the interacable a given distance. Note, this just updates the position by the move.

Parameters

[in] move Reference to Vec3D which is the distance to move the interactable.

Reimplemented from BaseInteractable.

Definition at line 184 of file TriangleWall.cc.

References BaseInteractable::move(), and updateVertexAndNormal().

 {
     BaseInteractable::move(move);
     updateVertexAndNormal();
 }

void TriangleWall::read ( std::istream & is )

overridevirtual

Reads an TriangleWall from an input stream, for example a restart file.

Parameters

[in] is The input stream from which the TriangleWall is read, usually a restart file.

Reimplemented from BaseWall.

Definition at line 123 of file TriangleWall.cc.

References constants::i, BaseWall::read(), updateVertexAndNormal(), and vertexInLabFrame_.

 {
     BaseWall::read(is);
     std::string dummy;
     is >> dummy;
     for (int i = 0; i < 3; i++)
     {
         is >> vertexInLabFrame_[i];
     }
     updateVertexAndNormal();
 }

void TriangleWall::rotate ( const Vec3D & angularVelocityDt )

overridevirtual

Rotates this BaseInteractable.

Rotates the interacable a given solid angle. Note, this just updates the orientation by the angle.

This function has been declared virtual, so it can be overridden for IntersectionOfWalls.

Parameters

[in] angularVelocityDt Reference to Vec3D which is the solid angle through which the interactable is rotated.

Todo:: TW the move and rotate functions should only pass the time step, as teh velocity can be accessed directly by the object; this would simplify functions like Screw::rotate

Reimplemented from BaseInteractable.

Definition at line 111 of file TriangleWall.cc.

References Vec3D::isZero(), BaseInteractable::rotate(), and updateVertexAndNormal().

 {
     if (!angularVelocityDt.isZero())
     {
         BaseInteractable::rotate(angularVelocityDt);
         updateVertexAndNormal();
     }
 }

void TriangleWall::setVertices	(	Vec3D	A,
		Vec3D	B,
		Vec3D	C
	)

Sets member variables such that the wall represents a triangle with vertices A, B, C.

position_ is set to the center of mass of the wall
vertexInLabFrame_ is set relative to the position
updateVertexAndNormal is called to set the remaining variables

Definition at line 165 of file TriangleWall.cc.

References BaseInteractable::getPosition(), BaseInteractable::setOrientation(), BaseInteractable::setPosition(), updateVertexAndNormal(), and vertexInLabFrame_.

Referenced by WallHandler::readTriangleWall().

 {
     setPosition((A + B + C) / 3);
     setOrientation({1, 0, 0, 0});
     vertexInLabFrame_[0] = A - getPosition();
     vertexInLabFrame_[1] = B - getPosition();
     vertexInLabFrame_[2] = C - getPosition();
     //vertexInLabFrame_[0] = A;
     //vertexInLabFrame_[1] = B;
     //vertexInLabFrame_[2] = C;
     updateVertexAndNormal();
 }

void TriangleWall::setVertices	(	Vec3D	A,
		Vec3D	B,
		Vec3D	C,
		Vec3D	position
	)

Same as #setVertices(A,B,C), but sets the position explicitly. The position is important when you rotate the wall, as the wall will be rotated around this position.

Definition at line 190 of file TriangleWall.cc.

References BaseInteractable::getPosition(), BaseInteractable::setOrientation(), BaseInteractable::setPosition(), updateVertexAndNormal(), and vertexInLabFrame_.

 {
     setPosition(position);
     setOrientation({1, 0, 0, 0});
     vertexInLabFrame_[0] = A - getPosition();
     vertexInLabFrame_[1] = B - getPosition();
     vertexInLabFrame_[2] = C - getPosition();
     updateVertexAndNormal();
 }

void TriangleWall::updateVertexAndNormal ( )

private

This function should be called after setting either position_ or vertexInLabFrame_.

vertex is set to the vertex position in the real coordinate system (rotated and shifted)
vertexMin_/vertexMax_ is set to define a bounding box around the wall (for contact detection)
edge_, edgeNormal_ and faceNormal_ is updated (stored for quick computation of contact point)

Definition at line 206 of file TriangleWall.cc.

References Vec3D::cross(), edge_, edgeLength_, edgeNormal_, faceNormal_, BaseInteractable::getOrientation(), BaseInteractable::getPosition(), constants::i, Vec3D::max(), Vec3D::min(), Vec3D::normalise(), Quaternion::rotate(), vertex_, vertexInLabFrame_, vertexMax_, and vertexMin_.

Referenced by move(), read(), rotate(), and setVertices().

 {
     for (int i = 0; i < 3; i++)
     {
         vertex_[i] = vertexInLabFrame_[i];
         getOrientation().rotate(vertex_[i]);
         vertex_[i] += getPosition();
     }
     vertexMin_ = Vec3D::min(Vec3D::min(vertex_[0], vertex_[1]), vertex_[2]);
     vertexMax_ = Vec3D::max(Vec3D::max(vertex_[0], vertex_[1]), vertex_[2]);
     
     edge_ = {vertex_[1] - vertex_[0], vertex_[2] - vertex_[1], vertex_[0] - vertex_[2]};
     faceNormal_ = Vec3D::cross(edge_[0], edge_[1]);
     faceNormal_.normalise();
     
     for (int i = 0; i < 3; i++)
     {
         edgeNormal_[i] = Vec3D::cross(edge_[i], faceNormal_);
         edgeLength_[i] = edge_[i].getLength();
         edge_[i] /= edgeLength_[i];
         edgeNormal_[i].normalise();
     }
     //logger(INFO,"vertex %,%,% edge %,%,% face %",vertex_[0],vertex_[1],vertex_[2],edgeNormal_[0],edgeNormal_[1],edgeNormal_[2],faceNormal_);
 }

void TriangleWall::write ( std::ostream & os ) const

overridevirtual

Writes an TriangleWall to an output stream, for example a restart file.

Parameters

[in] os The output stream where the TriangleWall must be written to, usually a restart file.

Reimplemented from BaseWall.

Definition at line 139 of file TriangleWall.cc.

References constants::i, vertexInLabFrame_, and BaseWall::write().

 {
     BaseWall::write(os);
     os << " vertexInLabFrame ";
     for (int i = 0; i < 3; i++)
     {
         os << ' ' << vertexInLabFrame_[i];
     }
 }

void TriangleWall::writeVTK ( VTKContainer & vtk ) const

overridevirtual

adds extra information to the points and triangleStrips vectors needed to plot the wall in vtk format

Parameters

points	Coordinates of the vertices of the triangulated surfaces (in the VTK file this is called POINTS)
triangleStrips	Indices of three vertices forming one triangulated surface (in the VTK file this is called CELL)

Reimplemented from BaseWall.

Definition at line 149 of file TriangleWall.cc.

References VTKContainer::points, VTKContainer::triangleStrips, and vertex_.

 {
     const unsigned long s = vtk.points.size();
     for (auto v : vertex_)
     {
         vtk.points.push_back(v);
     }
     std::vector<double> cell;
     cell.reserve(3);
     cell.push_back(s);
     cell.push_back(s + 1);
     cell.push_back(s + 2);
     //cell.push_back(s);
     vtk.triangleStrips.push_back(cell);
 }

Member Data Documentation

std::array<Vec3D, 3> TriangleWall::edge_

private

Definition at line 158 of file TriangleWall.h.

Referenced by getDistanceAndNormal(), and updateVertexAndNormal().

std::array<double, 3> TriangleWall::edgeLength_

private

Definition at line 159 of file TriangleWall.h.

Referenced by getDistanceAndNormal(), and updateVertexAndNormal().

std::array<Vec3D, 3> TriangleWall::edgeNormal_

private

stores the wall normal n in n.x=p

Definition at line 157 of file TriangleWall.h.

Referenced by getDistanceAndNormal(), and updateVertexAndNormal().

Vec3D TriangleWall::faceNormal_

private

stores the face normal, not rotated into the lab frame; thus, if the wall rotates, this normal has to be rotated as well

Definition at line 164 of file TriangleWall.h.

Referenced by getDistanceAndNormal(), and updateVertexAndNormal().

std::array<Vec3D, 3> TriangleWall::vertex_

private

stores the position of the vertices relative to the position of the wall but not rotated into the lab frame; thus, if the wall rotates, these vertices have to be rotated as well

Definition at line 146 of file TriangleWall.h.

Referenced by getDistanceAndNormal(), getVertex(), getVertices(), isInsideTriangle(), updateVertexAndNormal(), and writeVTK().