Defines the basic properties that a interactable object can have. More...

#include <BaseInteractable.h>

Inheritance diagram for BaseInteractable:

Public Member Functions
	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...

void	read (std::istream &is) override
	Reads a BaseInteractable from an input stream. More...

void	write (std::ostream &os) const override
	Write a BaseInteractable to an output stream. More...

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

virtual void	setIndSpecies (unsigned int indSpecies)
	Sets the index of the Species of this BaseInteractable. 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...

virtual void	move (const Vec3D &move)
	Moves this BaseInteractable by adding an amount to the position. More...

virtual void	rotate (const Vec3D &angularVelocityDt)
	Rotates 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 BaseInteraction *	getInteractionWith (BaseParticle P, unsigned timeStamp, InteractionHandler interactionHandler)=0
	Returns the interaction between this object and a given BaseParticle. 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 bool	isFixed () const =0

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 std::string	getName () const =0
	A purely virtual function. 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 Attributes
std::function< Vec3D(double)>	prescribedPosition_

std::function< Vec3D(double)>	prescribedVelocity_

std::function< Quaternion(double)>	prescribedOrientation_

std::function< Vec3D(double)>	prescribedAngularVelocity_

Vec3D	position_

Quaternion	orientation_

Vec3D	angularVelocity_

Vec3D	force_

Vec3D	torque_

std::vector< Vec3D >	forceOMP_

std::vector< Vec3D >	torqueOMP_

const ParticleSpecies *	species_

unsigned int	indSpecies_

Vec3D	velocity_

std::vector< BaseInteraction * >	interactions_

Detailed Description

Defines the basic properties that a interactable object can have.

Inherits from class BaseObject (public) Also it includes a lot of code to deal with interactable objects that have a prescibed motion. Most of the code in here is MercuryDPM internal. The only place an user will interface with this code is for setting the lambda functions that prescribe the motion of infinite mass particles.

Todo:: Check prescribed objects have infinite mass.

Definition at line 54 of file BaseInteractable.h.

Constructor & Destructor Documentation

BaseInteractable::BaseInteractable ( )

Default BaseInteractable constructor.

Simply creates an empty BaseInteractable, with all vectors relating to the positions and motions of the current object initialised to zero and all pointers to null.

Note that the function also sets the species index (indSpecies_) to zero by default, so any objects created will, by default, possess the properties associated with species 0.

Definition at line 42 of file BaseInteractable.cc.

References angularVelocity_, DEBUG, force_, indSpecies_, logger, orientation_, position_, prescribedAngularVelocity_, prescribedOrientation_, prescribedPosition_, prescribedVelocity_, Quaternion::setUnity(), Vec3D::setZero(), species_, torque_, and velocity_.

                                    :
         BaseObject()
 {
     //setting all vectors to zero
     position_.setZero();
     orientation_.setUnity();
     velocity_.setZero();
     angularVelocity_.setZero();
     force_.setZero();
     torque_.setZero();
     //setting the default species to species 0
     indSpecies_ = 0;
     //setting all relevant pointers to nullptr
     species_ = nullptr;
     prescribedPosition_ = nullptr;
     prescribedVelocity_ = nullptr;
     prescribedOrientation_ = nullptr;
     prescribedAngularVelocity_ = nullptr;
     logger(DEBUG, "BaseInteractable::BaseInteractable() finished");
 }

BaseInteractable::BaseInteractable ( const BaseInteractable & p )

Copy constructor.

Copies an existing BaseInteractable, p, and (almost) all objects it contains.

Note, that the interactions are not copied as these often require extra work. Rather, the interactions list is simply emptied using the standard C++ clear function, destroying all contents and leaving the interactions_ vector with a size of zero.

All the other properties are copied normally.

Please use this copy with care.

Definition at line 75 of file BaseInteractable.cc.

References angularVelocity_, DEBUG, force_, indSpecies_, interactions_, logger, orientation_, position_, prescribedAngularVelocity_, prescribedOrientation_, prescribedPosition_, prescribedVelocity_, species_, torque_, and velocity_.

         : BaseObject(p)
 {
     interactions_.clear();
     position_ = p.position_;
     orientation_ = p.orientation_;
     velocity_ = p.velocity_;
     angularVelocity_ = p.angularVelocity_;
     force_ = p.force_;
     torque_ = p.torque_;
     species_ = p.species_;
     indSpecies_ = p.indSpecies_;
     prescribedPosition_ = p.prescribedPosition_;
     prescribedVelocity_ = p.prescribedVelocity_;
     prescribedOrientation_ = p.prescribedOrientation_;
     prescribedAngularVelocity_ = p.prescribedAngularVelocity_;
     logger(DEBUG, "BaseInteractable::BaseInteractable(const BaseInteractable &p finished");
 }

BaseInteractable::~BaseInteractable ( )

override

Destructor, it simply destructs the BaseInteractable and all the objects it contains.

Removes all the interactions from the interactable.

Definition at line 97 of file BaseInteractable.cc.

References DEBUG, BaseObject::getId(), BaseObject::getIndex(), interactions_, logger, and VERBOSE.

 {
     logger(VERBOSE, "Deleting BaseInteractable with index= % and id = %, size = %", getIndex(), getId(),
            interactions_.size());
     while (!interactions_.empty())
     {
         interactions_.front()->removeFromHandler();
     }
     logger(DEBUG, "BaseInteractable::~BaseInteractable() finished");
 }

Member Function Documentation

void BaseInteractable::addAngularVelocity ( const Vec3D & angularVelocity )

add an increment to the angular velocity.

Parameters

[in] angularVelocity increment which to increase the angularVelocity by.

Definition at line 370 of file BaseInteractable.cc.

References angularVelocity_.

Referenced by BaseParticle::angularAccelerate().

 {
     angularVelocity_ += angularVelocity;
 }

void BaseInteractable::addForce ( const Vec3D & addForce )

Adds an amount to the force on this BaseInteractable.

Incremental version of BaseInteractable::setForce. Also see BaseInteraction::setForce for were this is used.

Parameters

[in] addForce Vec3D incremental force which is added to the total force of the interactable.

Definition at line 116 of file BaseInteractable.cc.

References force_, forceOMP_, and OMP_THREAD_NUM.

Referenced by MeshTriangle::checkInteractions(), DPMBase::computeExternalForces(), DPMBase::computeForcesDueToWalls(), DPMBase::computeInternalForce(), and BaseWall::getInteractionWith().

 {
     if (OMP_THREAD_NUM==0) {
         force_ += addForce;
     } else {
         forceOMP_[OMP_THREAD_NUM-1] += addForce;
     }
 }

void BaseInteractable::addInteraction ( BaseInteraction * I )

Adds an interaction to this BaseInteractable.

Added a new interactions to the current interactable.

Parameters

[in] I Pointer to the new interaction which is to be added to the list of interactions of this interactable.

Definition at line 292 of file BaseInteractable.cc.

References interactions_.

Referenced by InteractionHandler::addObject(), BaseInteraction::importI(), BaseInteraction::importP(), BaseInteraction::setI(), and BaseInteraction::setP().

 {
     interactions_.push_back(I);
 }

void BaseInteractable::addTorque ( const Vec3D & addTorque )

Adds an amount to the torque on this BaseInteractable.

Incremental version of BaseInteractable::setTorque. Also see BaseInteraction::setTorque for were this is used.

Parameters

[in] addTorque Vec3D incremental force which is added to the total torque of the interactable.

Definition at line 132 of file BaseInteractable.cc.

References OMP_THREAD_NUM, torque_, and torqueOMP_.

Referenced by MeshTriangle::checkInteractions(), DPMBase::computeForcesDueToWalls(), DPMBase::computeInternalForce(), and BaseWall::getInteractionWith().

 {
     if (OMP_THREAD_NUM==0) {
         torque_ += addTorque;
     } else {
         torqueOMP_[OMP_THREAD_NUM-1] += addTorque;
     }
 }

void BaseInteractable::addVelocity ( const Vec3D & velocity )

inline

adds an increment to the velocity.

See also BaseInteractable::setVelocity

Parameters

[in] velocity Vec3D containing the velocity increment which to increase the velocity by.

Definition at line 312 of file BaseInteractable.h.

References velocity_.

Referenced by BaseParticle::accelerate(), TimeDependentPeriodicBoundary::shiftAndBoostParticle(), ShearBoxBoundary::shiftHorizontalPosition(), and LeesEdwardsBoundary::shiftVerticalPosition().

313 { velocity_ += velocity; }

BaseInteractable::velocity_

Vec3D velocity_

Definition: BaseInteractable.h:500

void BaseInteractable::applyPrescribedAngularVelocity ( double time )

Computes the angular velocity from the user defined prescribed angular velocity.

This calls the prescribedAngularVelocity function if one has been defined. See also BaseInteractable::setPrescribedAngularVelocity

Parameters

[in] time double which is the current time of the simulation.

Definition at line 512 of file BaseInteractable.cc.

References prescribedAngularVelocity_, and setAngularVelocity().

Referenced by integrateAfterForceComputation(), and integrateBeforeForceComputation().

 {
     if (prescribedAngularVelocity_)
     {
         setAngularVelocity(prescribedAngularVelocity_(time));
     }
 }

void BaseInteractable::applyPrescribedOrientation ( double time )

Computes the orientation from the user defined prescribed orientation function.

This calls the prescribedOrientation function if one has been defined. See also BaseInteractable::setPrescribedOrientation

Parameters

[in] time double which is the current time of the simulation.

Definition at line 485 of file BaseInteractable.cc.

References prescribedOrientation_, and setOrientation().

Referenced by integrateBeforeForceComputation().

 {
     if (prescribedOrientation_)
     {
         setOrientation(prescribedOrientation_(time));
     }
 }

void BaseInteractable::applyPrescribedPosition ( double time )

Computes the position from the user defined prescribed position function.

This calls the prescribedPosition function if one has been defined. See also BaseInteractable::setPrescribedPosition

Parameters

[in] time double which is the current time of the simulation.

Definition at line 423 of file BaseInteractable.cc.

References prescribedPosition_, and setPosition().

Referenced by integrateBeforeForceComputation().

 {
     if (prescribedPosition_)
     {
         setPosition(prescribedPosition_(time));
     }
 }

void BaseInteractable::applyPrescribedVelocity ( double time )

Computes the velocity from the user defined prescribed velocity function.

This calls the prescribedVeclocity function if one has been defined. See also BaseInteractable::setPrescribedVelocity

Parameters

[in] time double which is the current time of the simulation.

Definition at line 454 of file BaseInteractable.cc.

References prescribedVelocity_, and setVelocity().

Referenced by integrateAfterForceComputation(), and integrateBeforeForceComputation().

 {
     if (prescribedVelocity_)
     {
         setVelocity(prescribedVelocity_(time));
     }
 }

void BaseInteractable::copyInteractionsForPeriodicParticles ( const BaseInteractable & pOriginal )

Copies interactions to this BaseInteractable whenever a periodic copy made.

This loops over all interactions of periodic (particle) and calls copySwitchPointer, which copies the interactions.

Parameters

[in] pOriginal Reference to the BaseInteractable which is to be copied to create the ghost particles.

Definition at line 386 of file BaseInteractable.cc.

References BaseInteraction::copySwitchPointer(), and interactions_.

Referenced by SubcriticalMaserBoundary::createGhostCopy(), ConstantMassFlowMaserBoundary::createGhostCopy(), PeriodicBoundary::createGhostParticle(), TimeDependentPeriodicBoundary::createGhostParticle(), ShearBoxBoundary::createHorizontalPeriodicParticles(), LeesEdwardsBoundary::createHorizontalPeriodicParticles(), AngledPeriodicBoundary::createPeriodicParticle(), ShearBoxBoundary::createVerticalPeriodicParticles(), and LeesEdwardsBoundary::createVerticalPeriodicParticles().

 {
     for (BaseInteraction* interaction : pOriginal.interactions_)
     {
         //So here this is the ghost and it is the interaction of the ghost/
         interaction->copySwitchPointer(&pOriginal, this);
     }
 }

const Vec3D & BaseInteractable::getAngularVelocity ( ) const

virtual

Returns the angular velocity of this interactable.

Returns the angular velocity of the BaseInteractbale. Note, this is the same for all BaseInteractables; it is the direction the object is moving in.

Returns: Vec3D reference which contains the current angular velocity of the current BaseInteractable.

Definition at line 341 of file BaseInteractable.cc.

References angularVelocity_.

Referenced by DeletionBoundary::checkBoundaryAfterParticleMoved(), MPISphericalParticle::copyDataFromParticleToMPIParticle(), copyVelocityFrom(), BaseParticle::getAngularMomentum(), HorizontalBaseScrew::getDistanceAndNormal(), SimpleDrumSuperquadrics::getDistanceAndNormal(), BaseParticle::getRotationalEnergy(), getVelocityAtContact(), integrateBeforeForceComputation(), BaseParticle::integrateBeforeForceComputation(), helpers::objectivenessTest(), AngledPeriodicBoundary::shiftPosition(), PeriodicBoundaryHandler::updateParticles(), and HorizontalBaseScrew::writeVTK().

 {
     return angularVelocity_;
 }

virtual Mdouble BaseInteractable::getCurvature ( const Vec3D & labFixedCoordinates ) const

inlinevirtual

returns the inverse radius, or curvature, of the surface. This value is zero for walls and gets overridden for particles that have finite radius

Todo:: should be wall-type dependent

Reimplemented in BaseParticle, and SuperQuadricParticle.

Definition at line 414 of file BaseInteractable.h.

Referenced by BaseInteraction::getEffectiveRadius().

415 { return 0.0; }

const Vec3D& BaseInteractable::getForce ( ) const

inline

Returns the force on this BaseInteractable.

Return the current force being to the BaseInteractable. Note, the code works by first computing the forces of each interaction and then it loops over all BaseInteracables applying forces to them from the interactions they are involved in.

Returns: const Vec3D reference that is the total force applied to this interactable.

Definition at line 126 of file BaseInteractable.h.

References force_.

Referenced by BaseParticle::integrateAfterForceComputation(), BaseParticle::integrateBeforeForceComputation(), BaseWall::setForceControl(), and BaseWall::setVelocityControl().

127 { return force_; }

BaseInteractable::force_

Vec3D force_

Definition: BaseInteractable.h:462

unsigned int BaseInteractable::getIndSpecies ( ) const

inline

Returns the index of the species associated with the interactable object.

Returns: Unsigned int which is the unique index of the species

Definition at line 88 of file BaseInteractable.h.

References indSpecies_.

Referenced by InteractionHandler::addInteraction(), DeletionBoundary::checkBoundaryAfterParticleMoved(), MPISphericalParticle::copyDataFromParticleToMPIParticle(), InteractionHandler::getInteraction(), ParticleHandler::readAndCreateObject(), BaseWall::setHandler(), BaseParticle::setHandler(), and BaseWall::setIndSpecies().

89 { return indSpecies_; }

BaseInteractable::indSpecies_

unsigned int indSpecies_

Definition: BaseInteractable.h:495

const std::vector<BaseInteraction*>& BaseInteractable::getInteractions ( ) const

inline

Returns a list of interactions which belong to this interactable.

Returns: An list of pointers to all the interactions which this interacable is involved in.

Definition at line 277 of file BaseInteractable.h.

References interactions_.

Referenced by MeshTriangle::checkInteractions(), PeriodicBoundaryHandler::findNewInteractions(), InteractionHandler::getExistingInteraction(), LiquidFilmParticle::getFieldVTK(), HeatFluidCoupledParticle::getFieldVTK(), InteractionHandler::getInteraction(), MeshTriangle::getInteractionWith(), BaseWall::getInteractionWith(), LiquidMigrationWilletInteraction::getNumberOfContacts(), DPMBase::removeDuplicatePeriodicParticles(), LiquidMigrationWilletInteraction::rupture(), and AngledPeriodicBoundary::shiftPosition().

278 { return interactions_; }

BaseInteractable::interactions_

std::vector< BaseInteraction * > interactions_

Definition: BaseInteractable.h:505

virtual BaseInteraction* BaseInteractable::getInteractionWith	(	BaseParticle *	P,
		unsigned	timeStamp,
		InteractionHandler *	interactionHandler
	)

pure virtual

Returns the interaction between this object and a given BaseParticle.

Todo:

TW make sure this function sets normal, distance, overlap, contact point

AT why is this a BaseParticle and not a BaseInteratable.

Implemented in BaseParticle, SuperQuadricParticle, BaseWall, MeshTriangle, TriangulatedWall, RestrictedWall, ArcWall, SineWall, Combtooth, ParabolaChute, and VChute.

Referenced by PeriodicBoundaryHandler::processReceivedInteractionData(), and Domain::processReceivedInteractionData().

virtual Mdouble BaseInteractable::getInvMass ( ) const

inlinevirtual

returns the inverse mass. This value is zero for walls and gets overridden for particles that have finite mass

Reimplemented in BaseParticle, and MeshTriangle.

Definition at line 408 of file BaseInteractable.h.

Referenced by BaseInteraction::getEffectiveMass().

409 { return 0.0; }

const Quaternion& BaseInteractable::getOrientation ( ) const

inline

Returns the orientation of this BaseInteractable.

Returns the reference to a Vec3D which contains the orientation of the interactionable. Please note the interpretation of this depends on which interactable. Please see derived objects for details.

Returns: Returns a reference to a Vec3D returns the position of the interactable.

Definition at line 230 of file BaseInteractable.h.

References orientation_.

Referenced by ScrewsymmetricIntersectionOfWalls::computeDeltaZ(), SuperQuadricParticle::computeHessianLabFixed(), SuperQuadricParticle::computeShape(), SuperQuadricParticle::computeShapeGradientLabFixed(), HorizontalBaseScrew::convertLimits(), AxisymmetricIntersectionOfWalls::convertLimits(), ScrewsymmetricIntersectionOfWalls::convertLimits(), MPISphericalParticle::copyDataFromParticleToMPIParticle(), copyPositionFrom(), InfiniteWall::createVTK(), BaseWall::getAxis(), SuperQuadricParticle::getContactPointPlanB(), InfiniteWall::getDistance(), NurbsWall::getDistanceAndNormal(), Screw::getDistanceAndNormal(), HorizontalBaseScrew::getDistanceAndNormal(), BasicIntersectionOfWalls::getDistanceAndNormal(), BasicUnionOfWalls::getDistanceAndNormal(), LevelSetWall::getDistanceAndNormal(), InfiniteWall::getDistanceAndNormal(), AxisymmetricIntersectionOfWalls::getDistanceAndNormal(), IntersectionOfWalls::getDistanceAndNormal(), InfiniteWall::getDistanceNormalOverlapSuperquadric(), BaseWall::getInteractionWith(), InfiniteWall::getNormal(), BasicIntersectionOfWalls::getVTK(), BasicUnionOfWalls::getVTK(), integrateAfterForceComputation(), BaseParticle::integrateAfterForceComputation(), integrateBeforeForceComputation(), BaseParticle::integrateBeforeForceComputation(), helpers::objectivenessTest(), CGFields::OrientationField::setFields(), AngledPeriodicBoundary::shiftPosition(), PeriodicBoundaryHandler::updateParticles(), TriangleWall::updateVertexAndNormal(), MeshTriangle::updateVertexAndNormal(), SuperQuadricParticle::writeDebugMessageMiddleOfLoop(), SuperQuadricParticle::writeDebugMessageStep1(), NurbsWall::writeVTK(), HorizontalBaseScrew::writeVTK(), LevelSetWall::writeVTK(), Screw::writeVTK(), AxisymmetricIntersectionOfWalls::writeVTK(), ScrewsymmetricIntersectionOfWalls::writeVTK(), and IntersectionOfWalls::writeVTK().

231 { return orientation_; }

BaseInteractable::orientation_

Quaternion orientation_

Definition: BaseInteractable.h:452

const Vec3D& BaseInteractable::getPosition ( ) const

inline

Returns the position of this BaseInteractable.

Returns the reference to a Vec3D which contains the position of the interactionable. Please note the interpretation of this depends on which interactable. For particles this is the centre of the particle; where for walls it is one point of the wall given \(r.n=p\)

Returns: Returns a reference to a Vec3D returns the position of the interactable.

Definition at line 218 of file BaseInteractable.h.

References position_.

219 { return position_; }

BaseInteractable::position_

Vec3D position_

Definition: BaseInteractable.h:446

const ParticleSpecies* BaseInteractable::getSpecies ( ) const

inline

Returns a pointer to the species of this BaseInteractable.

This function return a ParticleSpecies* for the current interacable. Please note, this is a ParticleSpecies; not, a BaseSpecies as interactables must have physically properties as well.

Returns: constant ParticleSpecies* pointer to the species storing the physical properties of this interactable.

Definition at line 108 of file BaseInteractable.h.

References species_.

Referenced by ThermalParticle::actionsAfterTimeStep(), ParticleHandler::addExistingObject(), ParticleHandler::addGhostObject(), WallHandler::addObject(), ParticleHandler::addObject(), IntersectionOfWalls::addObject(), SubcriticalMaserBoundary::addParticleToMaser(), ConstantMassFlowMaserBoundary::addParticleToMaser(), BaseInteraction::BaseInteraction(), SubcriticalMaserBoundary::checkBoundaryAfterParticleMoved(), ConstantMassFlowMaserBoundary::checkBoundaryAfterParticleMoved(), ChargedBondedInteraction::computeAdhesionForce(), CurvyChute::createBottom(), HeatFluidCoupledParticle::f1(), HeatFluidCoupledParticle::f2(), NurbsWall::getDistanceAndNormal(), Screw::getDistanceAndNormal(), BasicIntersectionOfWalls::getDistanceAndNormal(), BasicUnionOfWalls::getDistanceAndNormal(), IntersectionOfWalls::getDistanceAndNormal(), ChargedBondedInteraction::getElasticEnergy(), BaseParticle::getInfo(), BaseParticle::getInteractionDistance(), SuperQuadricParticle::getInteractionRadius(), SuperQuadricParticle::getInteractionWith(), BaseParticle::getMaxInteractionRadius(), BasicUnionOfWalls::getVTK(), BasicIntersectionOfWalls::getVTK(), SubcriticalMaserBoundary::isMaserParticle(), ConstantMassFlowMaserBoundary::isMaserParticle(), SubcriticalMaserBoundary::isNormalParticle(), ConstantMassFlowMaserBoundary::isNormalParticle(), DPMBase::readNextDataFile(), SubcriticalMaserBoundary::removeParticleFromMaser(), ConstantMassFlowMaserBoundary::removeParticleFromMaser(), RestrictedWall::set(), SuperQuadricParticle::setAxes(), LinearViscoelasticNormalSpecies::setCollisionTimeAndRestitutionCoefficient(), SuperQuadricParticle::setExponents(), SuperQuadricParticle::setInertia(), BaseParticle::setRadius(), BaseParticle::unfix(), and BaseInteraction::writeInteraction().

     {
         return species_;
     }

const Vec3D& BaseInteractable::getTorque ( ) const

inline

Returns the torque on this BaseInteractable.

Return the current torque being to the BaseInteractable. Note, the code works by first computing the forces of each interaction and then it loops over all BaseInteracables applying forces to them from the interactions they are involved in.

Returns: const Vec3D reference that is the total force applied to this interactable.

Definition at line 138 of file BaseInteractable.h.

References torque_.

Referenced by BaseParticle::integrateAfterForceComputation(), and BaseParticle::integrateBeforeForceComputation().

139 { return torque_; }

BaseInteractable::torque_

Vec3D torque_

Definition: BaseInteractable.h:467

const Vec3D & BaseInteractable::getVelocity ( ) const

virtual

Returns the velocity of this interactable.

Returns the velocity of the BaseInteractbale. Note, this is the same for all BaseInteractables; it is the direction the object is moving in.

Returns: Vec3D reference which contains the current velocity of the current BaseInteractable.

Definition at line 329 of file BaseInteractable.cc.

References velocity_.

Referenced by DeletionBoundary::checkBoundaryAfterParticleMoved(), HeaterBoundary::checkBoundaryAfterParticleMoved(), FrictionInteraction::computeFrictionForce(), MindlinRollingTorsionInteraction::computeFrictionForce(), MPISphericalParticle::copyDataFromParticleToMPIParticle(), copyVelocityFrom(), BaseParticle::getAngularMomentum(), BaseParticle::getKineticEnergy(), BaseParticle::getMomentum(), getVelocityAtContact(), integrateBeforeForceComputation(), BaseParticle::integrateBeforeForceComputation(), InfiniteWallWithHole::move_time(), helpers::objectivenessTest(), CircularPeriodicBoundary::rotateParticle(), CGFields::GradVelocityField::setFields(), CGFields::StandardFields::setFields(), AngledPeriodicBoundary::shiftPosition(), PeriodicBoundaryHandler::updateParticles(), and MeshTriangle::updateVerticesFromParticles().

 {
     return velocity_;
 }

const Vec3D BaseInteractable::getVelocityAtContact ( const Vec3D & contact ) const

virtual

Returns the velocity at the contact point, use by many force laws.

Reimplemented in MeshTriangle.

Definition at line 375 of file BaseInteractable.cc.

References Vec3D::cross(), getAngularVelocity(), getPosition(), and getVelocity().

 {
     return getVelocity() - Vec3D::cross(contact - getPosition(), getAngularVelocity());
 }

void BaseInteractable::integrateAfterForceComputation	(	double	time,
		double	timeStep
	)

This is part of the integration routine for objects with infinite mass.

This is the last part of time integration for interactable objects which have an infinite mass.

Parameters

[in]	time	double which is the current simulation time
[in]	timeStep	double which is the current delta time of the simulation i.e. the size of each time step.

Definition at line 611 of file BaseInteractable.cc.

References applyPrescribedAngularVelocity(), applyPrescribedVelocity(), getOrientation(), prescribedAngularVelocity_, prescribedOrientation_, prescribedPosition_, prescribedVelocity_, setAngularVelocity(), and setVelocity().

Referenced by BaseParticle::integrateAfterForceComputation(), and DPMBase::integrateAfterForceComputation().

 {
     if (prescribedPosition_)
     {
         if (prescribedVelocity_)
         {
             applyPrescribedVelocity(time + timeStep);
         }
         else
         {
             setVelocity((prescribedPosition_(time + 1.1 * timeStep) - prescribedPosition_(time + 0.9 * timeStep)) /
                         (0.2 * timeStep));
         }
     }
     else
     {
         if (prescribedVelocity_)
         {
             applyPrescribedVelocity(time + 0.5 * timeStep);
         }
     }
     if (prescribedOrientation_)
     {
         if (prescribedAngularVelocity_)
         {
             applyPrescribedAngularVelocity(time + timeStep);
         }
         else
         {
             setAngularVelocity(getOrientation().applyCInverse(
                     (prescribedOrientation_(time + 1.1 * timeStep) - prescribedOrientation_(time + 0.9 * timeStep)) /
                     (0.2 * timeStep)
             ));
         }
     }
     else
     {
         if (prescribedAngularVelocity_)
         {
             applyPrescribedAngularVelocity(time + 0.5 * timeStep);
         }
     }
 }

void BaseInteractable::integrateBeforeForceComputation	(	double	time,
		double	timeStep
	)

This is part of integrate routine for objects with infinite mass.

This does not first part of verlet integration but for objects with an infinite mass i.e. there motion is prescribed and not calculated from the applied forces. First it deals with the translation degrees of freedom and then secondly if deals with the angular degrees of freedom. Note, in both cases if the user has prescribed both positions and velocity these are used. If only only position is prescribed the velocity is computed from a finite difference. If only the velocity is prescribed the position is computed from integrating the velocity.

In the weird case they neither is set. The objects computed velocity is used to update its position.

Parameters

[in]	time	double which is the current simulation time
[in]	timeStep	double which is the current delta time of the simulation i.e. the size of each time step.

Definition at line 538 of file BaseInteractable.cc.

References applyPrescribedAngularVelocity(), applyPrescribedOrientation(), applyPrescribedPosition(), applyPrescribedVelocity(), getAngularVelocity(), getOrientation(), getVelocity(), move(), prescribedAngularVelocity_, prescribedOrientation_, prescribedPosition_, prescribedVelocity_, rotate(), setAngularVelocity(), and setVelocity().

Referenced by BaseParticle::integrateBeforeForceComputation(), and DPMBase::integrateBeforeForceComputation().

 {
     if (prescribedPosition_)
     {
         if (prescribedVelocity_)
         {
             //Both the velocity and position are defined; as we are using leap
             //frog method so the velocity is evaluated half a time later.
             applyPrescribedPosition(time + timeStep);
             applyPrescribedVelocity(time + 0.5 * timeStep);
         }
         else
         {
             //Only the position is defined.
             //Velocity is evaluated from a finite different of the Position
             //Note, we use 0.5 +- 0.1 timeStep for the velocity eval.
             applyPrescribedPosition(time + timeStep);
             setVelocity((prescribedPosition_(time + 0.6 * timeStep) - prescribedPosition_(time + 0.4 * timeStep)) /
                         (0.2 * timeStep));
         }
     }
     else
     {
         if (prescribedVelocity_)
         {
             //Only the velocity is set. The position is calculated from the
             //the integral of velocity.
             applyPrescribedVelocity(time + 0.5 * timeStep);
             move(getVelocity() * timeStep);
         }
         else
         {
             //Neither is set move based on the computed velocity of the object.
             move(getVelocity() * timeStep);
         }
     }
     if (prescribedOrientation_)
     {
         if (prescribedAngularVelocity_)
         {
             applyPrescribedOrientation(time + timeStep);
             applyPrescribedAngularVelocity(time + 0.5 * timeStep);
         }
         else
         {
             applyPrescribedOrientation(time + timeStep);
             setAngularVelocity(getOrientation().applyCInverse(
                     (prescribedOrientation_(time + 0.6 * timeStep) - prescribedOrientation_(time + 0.4 * timeStep)) /
                     (0.2 * timeStep)
             ));
         }
     }
     else
     {
         if (prescribedAngularVelocity_)
         {
             applyPrescribedAngularVelocity(time + 0.5 * timeStep);
             rotate(getAngularVelocity() * timeStep);
         }
         else
         {
             rotate(getAngularVelocity() * timeStep);
         }
     }
 }

virtual bool BaseInteractable::isFixed ( ) const

pure virtual

used to distinguish particles which belong to the flow and fixed particles/walls

Implemented in BaseWall, and BaseParticle.

void BaseInteractable::move ( const Vec3D & move )

virtual

Moves this BaseInteractable by adding an amount to the position.

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 in MeshTriangle, TriangulatedWall, and TriangleWall.

Definition at line 215 of file BaseInteractable.cc.

References position_.

Referenced by ConstantMassFlowMaserBoundary::addParticleToMaser(), ConstantMassFlowMaserBoundary::checkBoundaryAfterParticleMoved(), ConstantMassFlowMaserBoundary::createPeriodicParticle(), integrateBeforeForceComputation(), BaseParticle::integrateBeforeForceComputation(), TriangleWall::move(), TriangulatedWall::move(), MeshTriangle::move(), HopperInsertionBoundary::placeParticle(), ConstantMassFlowMaserBoundary::removeParticleFromMaser(), TimeDependentPeriodicBoundary::shiftAndBoostParticle(), ShearBoxBoundary::shiftHorizontalPosition(), LeesEdwardsBoundary::shiftHorizontalPosition(), AngledPeriodicBoundary::shiftPosition(), PeriodicBoundary::shiftPosition(), SubcriticalMaserBoundary::shiftPosition(), ConstantMassFlowMaserBoundary::shiftPosition(), ShearBoxBoundary::shiftVerticalPosition(), and LeesEdwardsBoundary::shiftVerticalPosition().

 {
     position_ += move;
 }

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

overridevirtual

Reads a BaseInteractable from an input stream.

BaseInteacatable read functions. Reads in all the information about an interacatable. Note, this can be from any istream but would normally be a file See also BaseInteractable::write

Parameters

[in] is std::istream to which the information is read from.

Implements BaseObject.

Reimplemented in BaseParticle, IntersectionOfWalls, ScrewsymmetricIntersectionOfWalls, AxisymmetricIntersectionOfWalls, InfiniteWall, TriangulatedWall, HeatFluidCoupledParticle, MeshTriangle, SphericalWall, LiquidFilmParticle, InfiniteWallWithHole, ThermalParticle, ArcWall, SuperQuadricParticle, Screw, HorizontalScrew, TriangleWall, LevelSetWall, Coil, SimpleDrumSuperquadrics, RestrictedWall, BasicIntersectionOfWalls, BasicUnionOfWalls, HorizontalBaseScrew, NurbsWall, BaseWall, CylindricalWall, SineWall, Combtooth, ParabolaChute, and VChute.

Definition at line 244 of file BaseInteractable.cc.

References angularVelocity_, force_, indSpecies_, orientation_, position_, BaseObject::read(), torque_, and velocity_.

Referenced by BaseWall::read(), and BaseParticle::read().

 {
     BaseObject::read(is);
     std::string dummy;
     is >> dummy >> indSpecies_;
     is >> dummy >> position_;
     is >> dummy;
     //this if-statement is added to read Kasper van der Vaart's data files, which contain an additional variable named positionAbsolute
     if (dummy == "positionAbsolute")
     {
         is >> dummy >> dummy >> dummy >> dummy;
     }
     is >> orientation_;
     is >> dummy >> velocity_;
     is >> dummy >> angularVelocity_;
     is >> dummy;
     if (dummy == "0")
         is >> dummy;
     is >> force_;
     is >> dummy >> torque_;
 }

bool BaseInteractable::removeInteraction ( BaseInteraction * I )

Removes an interaction from this BaseInteractable.

Removes a given interaction form the list of interactions belonging to the current interacable. This functions returns true to the interaction was found and returns false if the given interaction did not exist and the interaction was not removed. Note that this cannot be done with a range-based for, since we need the iterator to erase the interaction.

Parameters

[in] I BaseInteraction pointer which is the interaction to be removed.

Returns: bool True if the interaction was found and removed; false if the interaction did not exist for that interactable.

Definition at line 308 of file BaseInteractable.cc.

References interactions_, logger, and WARN.

Referenced by BaseInteraction::importI(), BaseInteraction::importP(), BaseInteraction::setI(), BaseInteraction::setP(), and BaseInteraction::~BaseInteraction().

 {
     for (std::vector<BaseInteraction*>::iterator it = interactions_.begin(); it != interactions_.end(); ++it)
     {
         if (I == (*it))
         {
             interactions_.erase(it);
             return true;
         }
     }
     logger(WARN, "Error in BaseInteractable::removeInteraction: Interaction could not be removed");
     return false;
 }

void BaseInteractable::resetForceTorque ( int numberOfOMPthreads )

Definition at line 141 of file BaseInteractable.cc.

References force_, forceOMP_, Vec3D::setZero(), torque_, and torqueOMP_.

Referenced by DPMBase::computeAllForces().

 {
     #ifdef MERCURY_USE_OMP
     forceOMP_.resize(numberOfOMPthreads-1);
     torqueOMP_.resize(numberOfOMPthreads-1);
     for (auto& force : forceOMP_)
     {
         force.setZero();
     }
     
     for (auto& torque : torqueOMP_)
     {
         torque.setZero();
     }
     #endif
     
     force_.setZero();
     torque_.setZero();
 }

void BaseInteractable::rotate ( const Vec3D & angularVelocityDt )

virtual

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 in MeshTriangle, TriangleWall, and Screw.

Definition at line 230 of file BaseInteractable.cc.

References orientation_, and Quaternion::updateAngularDisplacement().

Referenced by integrateBeforeForceComputation(), BaseParticle::integrateBeforeForceComputation(), Screw::rotate(), TriangleWall::rotate(), MeshTriangle::rotate(), and AngledPeriodicBoundary::shiftPosition().

 {
     //if (!angularVelocityDt.isZero()) {
     orientation_.updateAngularDisplacement(angularVelocityDt);
     //}
 }

void BaseInteractable::setAngularVelocity ( const Vec3D & angularVelocity )

set the angular velocity of the BaseInteractble.

Parameters

[in] angularVelocity Vec3D which is the angularVelocity of the interactable.

Definition at line 360 of file BaseInteractable.cc.

References angularVelocity_.

Referenced by applyPrescribedAngularVelocity(), MPISphericalParticle::copyDataFromMPIParticleToParticle(), BaseParticle::fixParticle(), integrateAfterForceComputation(), integrateBeforeForceComputation(), helpers::loadingTest(), helpers::normalAndTangentialLoadingTest(), helpers::objectivenessTest(), BaseParticle::oldRead(), DPMBase::readNextDataFile(), WallHandler::readTriangleWall(), ChuteBottom::setupInitialConditions(), and PeriodicBoundaryHandler::updateParticles().

 {
     angularVelocity_ = angularVelocity;
 }

void BaseInteractable::setForce ( const Vec3D & force )

inline

Sets the force on this BaseInteractable.

This sets the force being applied to this interactable. Note, first the code computes all forces in the interactions and then loops over all interactable objects applying the forces from the interactions to the interactables involved in the interaction.

Parameters

[in] force Vec3D which is the force to be applied.

Definition at line 149 of file BaseInteractable.h.

References force_.

150 { force_ = force; }

BaseInteractable::force_

Vec3D force_

Definition: BaseInteractable.h:462

virtual void BaseInteractable::setIndSpecies ( unsigned int indSpecies )

inlinevirtual

Sets the index of the Species of this BaseInteractable.

This set the species associated with this interactable. This function should not be used and BaseInteractable::setSpecies should be used instead. See also BaseInteractable::setSpecies

Reimplemented in BaseParticle, and BaseWall.

Definition at line 98 of file BaseInteractable.h.

References indSpecies_.

Referenced by BaseParticle::oldRead(), BaseWall::setIndSpecies(), and BaseParticle::setIndSpecies().

99 { indSpecies_ = indSpecies; }

BaseInteractable::indSpecies_

unsigned int indSpecies_

Definition: BaseInteractable.h:495

void BaseInteractable::setOrientation ( const Quaternion & orientation )

inline

Sets the orientation of this BaseInteractable.

Interpretation depends on which interactable is being considered See also BaseInteractable::getOrientation.

Parameters

[in] orientation Reference to Vec3D storing the orientation of the particle.

Definition at line 260 of file BaseInteractable.h.

References orientation_.

Referenced by IntersectionOfWalls::addObject(), applyPrescribedOrientation(), MPISphericalParticle::copyDataFromMPIParticleToParticle(), SuperQuadricParticle::getContactPointPlanB(), helpers::objectivenessTest(), InfiniteWall::oldRead(), BaseParticle::oldRead(), DPMBase::readNextDataFile(), ChuteBottom::setupInitialConditions(), TriangleWall::setVertices(), MeshTriangle::setVertices(), and PeriodicBoundaryHandler::updateParticles().

261 { orientation_ = orientation; }

BaseInteractable::orientation_

Quaternion orientation_

Definition: BaseInteractable.h:452

void BaseInteractable::setOrientationViaEuler ( Vec3D eulerAngle )

Sets the orientation of this BaseInteractable by defining the euler angles.

Definition at line 204 of file BaseInteractable.cc.

References orientation_, and Quaternion::setEuler().

Referenced by DPMBase::readNextDataFile().

 {
     orientation_.setEuler(eulerAngle);
 }

void BaseInteractable::setOrientationViaNormal ( Vec3D normal )

Sets the orientation of this BaseInteractable by defining the vector that results from the rotation of the (1,0,0) vector.

Sets the orientation of the interactable. interpretation depends on which interactable is being considered See also BaseInteractable::getOrientation.

Parameters

[in] orientation Reference to Vec3D storing the orientation of the particle.

Definition at line 199 of file BaseInteractable.cc.

References orientation_, and Quaternion::setOrientationViaNormal().

Referenced by AxisymmetricIntersectionOfWalls::AxisymmetricIntersectionOfWalls(), HorizontalBaseScrew::HorizontalBaseScrew(), Screw::Screw(), ScrewsymmetricIntersectionOfWalls::ScrewsymmetricIntersectionOfWalls(), HorizontalBaseScrew::setAxis(), SimpleDrumSuperquadrics::setAxis(), AxisymmetricIntersectionOfWalls::setAxis(), ScrewsymmetricIntersectionOfWalls::setAxis(), and InfiniteWall::setNormal().

 {
     orientation_.setOrientationViaNormal(normal);
 }

void BaseInteractable::setPosition ( const Vec3D & position )

inline

Sets the position of this BaseInteractable.

Interpretation depends on which interactable is being considered See also BaseInteractable::getPosistion.

Parameters

[in] position Reference to Vec3D storing the position of the particle.

Definition at line 239 of file BaseInteractable.h.

References position_.

Referenced by IntersectionOfWalls::addObject(), BaseWall::addParticlesAtWall(), applyPrescribedPosition(), AxisymmetricIntersectionOfWalls::AxisymmetricIntersectionOfWalls(), DropletBoundary::checkBoundaryAfterParticlesMove(), DPMBase::checkParticleForInteractionLocalPeriodic(), BaseCluster::computeInternalStructure(), MPISphericalParticle::copyDataFromMPIParticleToParticle(), Chute::createBottom(), Chute::createFlowParticle(), Membrane::createVertexParticles(), PeriodicBoundaryHandler::findTargetProcessor(), BidisperseCubeInsertionBoundary::generateParticle(), SuperQuadricParticle::getContactPointPlanB(), BasicIntersectionOfWalls::getDistanceAndNormal(), BasicUnionOfWalls::getDistanceAndNormal(), BasicIntersectionOfWalls::getVTK(), BasicUnionOfWalls::getVTK(), HorizontalBaseScrew::HorizontalBaseScrew(), InfiniteWall::InfiniteWall(), helpers::loadingTest(), helpers::normalAndTangentialLoadingTest(), helpers::objectivenessTest(), InfiniteWall::oldRead(), BaseParticle::oldRead(), BaseCluster::particleInsertionSuccessful(), ChuteInsertionBoundary::placeParticle(), FixedClusterInsertionBoundary::placeParticle(), HopperInsertionBoundary::placeParticle(), RandomClusterInsertionBoundary::placeParticle(), PolydisperseInsertionBoundary::placeParticle(), CubeInsertionBoundary::placeParticle(), FileReader::read(), DPMBase::readNextDataFile(), CircularPeriodicBoundary::rotateParticle(), ScrewsymmetricIntersectionOfWalls::ScrewsymmetricIntersectionOfWalls(), InfiniteWall::set(), ChuteBottom::setupInitialConditions(), TriangleWall::setVertices(), MeshTriangle::setVertices(), and PeriodicBoundaryHandler::updateParticles().

240 { position_ = position; }

BaseInteractable::position_

Vec3D position_

Definition: BaseInteractable.h:446

void BaseInteractable::setPrescribedAngularVelocity ( const std::function< Vec3D(double)> & prescribedAngularVelocity )

Allows the angular velocity of the infinite mass interactable to be prescribed.

Definition at line 501 of file BaseInteractable.cc.

References prescribedAngularVelocity_.

 {
     prescribedAngularVelocity_ = prescribedAngularVelocity;
 }

void BaseInteractable::setPrescribedOrientation ( const std::function< Quaternion(double)> & prescribedOrientation )

Allows the orientation of the infinite mass interactbale to be prescribed.

This is similar to the BaseInteractable::setPrescribedPosition and works the same way. See BaseInteractable::setPrescibedPosition and BaseInteractable::setPrescribedVelocity for more details how it works. Note, the rate of change of the orientation can also be set using the function BaseInteractable::setPrescribedAngularVelocity.

Parameters

[in] prescribedOrientation std::function which is the lambda function and takes a double the time and returns a Vec 3D which is the orientation of the interactable for that time.

Definition at line 474 of file BaseInteractable.cc.

References prescribedOrientation_.

 {
     prescribedOrientation_ = prescribedOrientation;
 }

void BaseInteractable::setPrescribedPosition ( const std::function< Vec3D(double)> & prescribedPosition )

Allows the position of an infinite mass interactable to be prescribed.

This functions is used to give an interactable a prescribed motion, which is defined by a std::function. This is the new moving walls interface. A demo of the use would be: setPrescribedPosition([this] (double time) { return Vec3D(getXMin(),0.0,shaker_amp * std::sin(t * 2.0 * shaker_freq * constants::pi)); } ); This example moves the wall sinusoidally with time.

Parameters

[in] prescribedPosition std::function which is the lambda function and takes a double the time and returns a Vec 3D which is the position of the interactable for that time. See also BaseInteractable::setPrescribedVelocity for more information.

Definition at line 413 of file BaseInteractable.cc.

References prescribedPosition_.

 {
     prescribedPosition_ = prescribedPosition;
 }

void BaseInteractable::setPrescribedVelocity ( const std::function< Vec3D(double)> & prescribedVelocity )

Allows the velocity of an infinite mass interactable to be prescribed.

In a similar manner to BaseInteractable::setPrescribedPosition this sets the velocity of the BaseInteactable. See also BaseInteractable::setPrescribedPosition Note, it is valid to set both the velocity and the position. No checking if these are consist is done at all. If you only set one of these function the other is automatically calculated using a by numerically differentiating or integrating the functions will is prescribed.

Parameters

[in] prescribedVelocity std::function which is the lambda function and takes a double the time and returns a Vec 3D which is the velocity of the intertable for that time.

Definition at line 444 of file BaseInteractable.cc.

References prescribedVelocity_.

Referenced by BaseWall::setForceControl(), and BaseWall::setVelocityControl().

 {
     prescribedVelocity_ = prescribedVelocity;
 }

void BaseInteractable::setSpecies ( const ParticleSpecies * species )

Sets the species of this BaseInteractable.

This function sets the species associated with this interactable object. Again this must be a ParticleSpecies. Note, it also automatically sets the index on the indSpecies_ by working up the correct index. However, index should be carefully used

Parameters

[in] species ParticleSpcies pointer which is species holding the physical properties.

Definition at line 185 of file BaseInteractable.cc.

References BaseSpecies::getHandler(), BaseObject::getIndex(), BaseHandler< T >::getObject(), indSpecies_, logger, and species_.

Referenced by BaseWall::setSpecies(), and BaseParticle::setSpecies().

 {
     logger.assert_debug(species->getHandler() != nullptr && species->getHandler()->getObject(species->getIndex())==species, "Error: Species is not part of any handler yet");
     species_ = species;
     indSpecies_ = species->getIndex();
 }

void BaseInteractable::setTorque ( const Vec3D & torque )

inline

Sets the torque on this BaseInteractable.

This sets the torque being applied to this interactable. Note, first the code computes all force/torques in the interactions and then loops over all interactable objects applying the torques from the interactions to the interactables involved in the interaction.

Parameters

[in] torque Vec3D which is the force to be applied.

Definition at line 161 of file BaseInteractable.h.

References torque_.

162 { torque_ = torque; }

BaseInteractable::torque_

Vec3D torque_

Definition: BaseInteractable.h:467

void BaseInteractable::setVelocity ( const Vec3D & velocity )

set the velocity of the BaseInteractable.