MercuryDPM  Trunk
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
ScrewsymmetricIntersectionOfWalls.h
Go to the documentation of this file.
1 //Copyright (c) 2013-2020, The MercuryDPM Developers Team. All rights reserved.
2 //For the list of developers, see <http://www.MercuryDPM.org/Team>.
3 //
4 //Redistribution and use in source and binary forms, with or without
5 //modification, are permitted provided that the following conditions are met:
6 // * Redistributions of source code must retain the above copyright
7 // notice, this list of conditions and the following disclaimer.
8 // * Redistributions in binary form must reproduce the above copyright
9 // notice, this list of conditions and the following disclaimer in the
10 // documentation and/or other materials provided with the distribution.
11 // * Neither the name MercuryDPM nor the
12 // names of its contributors may be used to endorse or promote products
13 // derived from this software without specific prior written permission.
14 //
15 //THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
16 //ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
17 //WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18 //DISCLAIMED. IN NO EVENT SHALL THE MERCURYDPM DEVELOPERS TEAM BE LIABLE FOR ANY
19 //DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
20 //(INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
21 //LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND
22 //ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23 //(INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24 //SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25 
26 #ifndef SCREWSYMMETRICINTERSECTIONOFWALLS_H
27 #define SCREWSYMMETRICINTERSECTIONOFWALLS_H
28 
29 #include "IntersectionOfWalls.h"
30 #include "InteractionHandler.h"
31 #include "Math/Vector.h"
32 
126 {
127 public:
132 
137 
141  ScrewsymmetricIntersectionOfWalls(Vec3D position, Vec3D normal, std::vector<normalAndPosition> walls,
142  const ParticleSpecies* species);
143 
148 
153 
158 
164  bool getDistanceAndNormal(const BaseParticle& P, Mdouble& distance, Vec3D& normal_return) const final;
165 
166  void computeNormalRadialDeltaN(const Vec3D& positionLabFrame, Vec3D& normalVector, Vec3D& radialVector, Mdouble& deltaN) const;
167 
168  Mdouble computeDeltaZ(const Vec3D& positionLabFrame, Mdouble h, Mdouble pitch) const;
169 
173  void read(std::istream& is) final;
174 
178  void write(std::ostream& os) const final;
179 
183  std::string getName() const final;
184 
185  void setAxis(Vec3D a);
186 
188  void convertLimits(Vec3D& min, Vec3D& max) const;
189 
190  void writeVTK(VTKContainer& vtk) const override;
191 
192  void setRightHandedness(bool rightHandedness) {rightHandedness_=rightHandedness;}
193  bool getRightHandedness() const {return rightHandedness_;}
194 
195  void setPitch(Mdouble pitch) {pitch_=pitch;}
196  Mdouble getPitch() const {return pitch_;}
197 
198  void setThickness(Mdouble thickness) {thickness_=thickness;}
199  Mdouble getThickness() const {return thickness_;}
200 
201  void setLength(Mdouble length) {length_=length;}
202  Mdouble getLength() const {return length_;}
203 
205  bool rightHandedness_ = true;
206 
207 // void setPointsAndLines(unsigned int n) override {
208 // IntersectionOfWalls::setPointsAndLines(n);
209 // rMax_ = 0;
210 // for (auto p : C_) {
211 // rMax_ = std::max(rMax_,p.X);
212 // }
213 // logger(INFO,"rMax=%",rMax_);
214 // }
215 
216 
217 private:
218 
223 
226 
229 
230 // ///\brief The outer screw blade radius.
231 // Mdouble rMax_ = 0;
232 };
233 
234 
235 #endif
A IntersectionOfWalls is convex polygon defined as an intersection of InfiniteWall's.
double Mdouble
Definition: GeneralDefine.h:34
bool getDistanceAndNormal(const BaseParticle &P, Mdouble &distance, Vec3D &normal_return) const final
Computes the distance from the wall for a given BaseParticle and returns true if there is a collision...
Mdouble pitch_
axial length of one rotation
ScrewsymmetricIntersectionOfWalls * copy() const final
Wall copy method. It calls the copy constructor of this Wall, useful for polymorphism.
void read(std::istream &is) final
reads wall
void write(std::ostream &os) const final
outputs wall
void convertLimits(Vec3D &min, Vec3D &max) const
Use ScrewsymmetricIntersectionOfWalls to define screwsymmetric walls, such as cylinders, cones, etc.
ScrewsymmetricIntersectionOfWalls & operator=(const ScrewsymmetricIntersectionOfWalls &other)
Copy assignment operator.
Mdouble computeDeltaZ(const Vec3D &positionLabFrame, Mdouble h, Mdouble pitch) const
void computeNormalRadialDeltaN(const Vec3D &positionLabFrame, Vec3D &normalVector, Vec3D &radialVector, Mdouble &deltaN) const
Mdouble thickness_
The thickness of the screw blade.
void writeVTK(VTKContainer &vtk) const override
bool rightHandedness_
The right handedness of the screw, i.e. the direction of the screw-blade.
Definition: Vector.h:49
std::string getName() const final
Returns the name of the object.