ESPResSo
Extensible Simulation Package for Research on Soft Matter Systems
Loading...
Searching...
No Matches
shapes/include/shapes/Ellipsoid.hpp
Go to the documentation of this file.
1/*
2 * Copyright (C) 2010-2022 The ESPResSo project
3 * Copyright (C) 2002,2003,2004,2005,2006,2007,2008,2009,2010
4 * Max-Planck-Institute for Polymer Research, Theory Group
5 *
6 * This file is part of ESPResSo.
7 *
8 * ESPResSo is free software: you can redistribute it and/or modify
9 * it under the terms of the GNU General Public License as published by
10 * the Free Software Foundation, either version 3 of the License, or
11 * (at your option) any later version.
12 *
13 * ESPResSo is distributed in the hope that it will be useful,
14 * but WITHOUT ANY WARRANTY; without even the implied warranty of
15 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
16 * GNU General Public License for more details.
17 *
18 * You should have received a copy of the GNU General Public License
19 * along with this program. If not, see <http://www.gnu.org/licenses/>.
20 */
21
22#ifndef SHAPES_ELLIPSOID_HPP
23#define SHAPES_ELLIPSOID_HPP
24
25#include "Shape.hpp"
26
27#include <utils/Vector.hpp>
28
29namespace Shapes {
30class Ellipsoid : public Shape {
31public:
33 : m_center({0.0, 0.0, 0.0}), m_semiaxes({1.0, 1.0, 1.0}),
34 m_direction(1.0) {}
35
36 void calculate_dist(const Utils::Vector3d &pos, double &dist,
37 Utils::Vector3d &vec) const override;
38
39 void set_semiaxis_a(const double &value) { m_semiaxes[0] = value; }
40 void set_semiaxis_b(const double &value) {
41 m_semiaxes[1] = value;
42 m_semiaxes[2] = value;
43 }
44
45 // not exposed via interface; used in core unit test
46 void set_semiaxis_c(const double &value) { m_semiaxes[2] = value; }
47
48 Utils::Vector3d &center() { return m_center; }
49 double &semiaxis_a() { return m_semiaxes[0]; }
50 double &semiaxis_b() { return m_semiaxes[1]; }
51 double &semiaxis_c() { return m_semiaxes[2]; }
52 double &direction() { return m_direction; }
53
54private:
55 bool inside_ellipsoid(const Utils::Vector3d &ppos) const;
56 double newton_term(const Utils::Vector3d &ppos, const double &l) const;
57
58 Utils::Vector3d m_center;
59 Utils::Vector3d m_semiaxes;
60 double m_direction;
61};
62} // namespace Shapes
63
64#endif
Vector implementation and trait types for boost qvm interoperability.
void calculate_dist(const Utils::Vector3d &pos, double &dist, Utils::Vector3d &vec) const override
Definition Ellipsoid.cpp:30
void set_semiaxis_b(const double &value)
void set_semiaxis_a(const double &value)
void set_semiaxis_c(const double &value)